def __init__(self, raster, geot, nodata_value=np.nan, fill_value=-<NUM_LIT>, projection=None, datatype=None):
super(GeoRaster, self).__init__()<EOL>if isinstance(raster, np.ma.core.MaskedArray):<EOL><INDENT>self.raster = raster<EOL><DEDENT>else:<EOL><INDENT>self.raster = np.ma.masked_array(raster, mask=raster == nodata_value,<EOL>fill_value=fill_value)<EOL><DEDENT>self.geot = geot<EOL>self.nodata_value = nodata_value<EOL>self.shape = raster.shape<EOL>self.x_cell_size = geot[<NUM_LIT:1>]<EOL>self.y_cell_size = geot[-<NUM_LIT:1>]<EOL>self.xmin = geot[<NUM_LIT:0>]<EOL>self.ymax = geot[<NUM_LIT:3>]<EOL>self.xmax = self.xmin + self.x_cell_size * self.shape[<NUM_LIT:1>]<EOL>self.ymin = self.ymax + self.y_cell_size * self.shape[<NUM_LIT:0>]<EOL>self.bounds = (self.xmin, self.ymin, self.xmax, self.ymax)<EOL>self.projection = projection<EOL>self.datatype = datatype<EOL>self.mcp_cost = None<EOL>self.weights = None<EOL>self.G = None<EOL>self.Gamma = None<EOL>self.Join_Counts = None<EOL>self.Moran = None<EOL>self.Geary = None<EOL>self.Moran_Local = None<EOL>
Initialize Georaster
Usage:
    geo=GeoRaster(raster, geot, nodata_value=ndv)
where
    raster: Numpy masked array with the raster data,
            which could be loaded with from_file(file) or load_tiff(file)
    geot: GDAL Geotransformation
    nodata_value: No data value in raster, optional