def __init__(self, dataset_name, srid, bottom_up_data=False): dataset_name = check_zip(dataset_name) ds = gdal.Open(dataset_name, GA_ReadOnly) if ds is None: raise Exception("Could not open dataset %s" % dataset_name) self.dsname = dataset_name self.ds = ds self.RasterCount = ds.RasterCount gt = ds.GetGeoTransform() size = (ds.RasterXSize, ds.RasterYSize) ul = (gt[0], gt[3]) scale = (gt[1], gt[5]) skew = (gt[2], gt[4]) projection = ds.GetProjectionRef() #set band 1 as default self.raster_band = self.ds.GetRasterBand(1) gdal_datatype = self.raster_band.DataType nodata_value = self.raster_band.GetNoDataValue() super(GDALRaster, self).__init__(size=size, ul=ul, scale=scale, skew=skew, srid=srid, gdal_datatype=gdal_datatype, nodata_value=nodata_value, nodata_range=None, bottom_up=bottom_up_data)
def __init__(self, file_name): self.file_name = check_zip(file_name)