def result_extent(self): dataset = util.gdal_open(self.scenario.lg) shape = dataset.RasterXSize, dataset.RasterYSize geotransform = dataset.GetGeoTransform() return util.google_extent_from_geotransform( shape, geotransform)
def save_raster(self, fieldname): self.save_file(fieldname) path = self.metadata[fieldname]['path'] dataset = gdal_open(path) if dataset is None: self.add_field_error(fieldname, 'Kan rasterbestand niet openen') return self.cleaned_data.get(fieldname) shape = dataset.GetRasterBand(1).ReadAsArray().shape self.metadata[fieldname] = { 'path': path, 'shape': shape, 'geotransform': dataset.GetGeoTransform() } return self.cleaned_data.get(fieldname)