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)
