Пример #1
0
 def __init__(self,input_raster,data_type=np.uint16,band_selection=0,projection=None):
     self.input_raster = input_raster
     self.band_list = read_image(input_raster, data_type, band_selection)
     self.rows, self.cols, self.nbands, self.geo_transform, self.projection = read_image_parameters(input_raster)
     if projection and self._parser_projection != projection:
         self.__change_projection(projection)
     self.resolution = abs(self.geo_transform[1])
Пример #2
0
 def __band_list(self,std_fix=True):
     raster_arrays_str = " ".join([raster.input_raster for raster in self.__rasters])
     self.tmpfile = tempfile.NamedTemporaryFile().name
     cmd = RasterGdalExtent.GDAL_MERGE_CMD+str(" -init \"0 0 255\" -n 0 -o {} ".format(self.tmpfile))+raster_arrays_str
     executeGdal(cmd)
     output_band = read_image(self.tmpfile,0,0)
     return output_band