예제 #1
0
파일: gdal_raster.py 프로젝트: ujjwaln/namo
    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)
예제 #2
0
파일: gdal_raster.py 프로젝트: ujjwaln/namo
    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)
예제 #3
0
 def __init__(self, file_name):
     self.file_name = check_zip(file_name)