def __init__(self, srs=0): """Create Spatial Reference System from input parameter""" if sys.version_info.major == 2: str_types = (str, unicode) else: str_types = str # create SRS osr.SpatialReference.__init__(self) # parse input parameters status = 1 if srs is 0: # generate default WGS84 SRS status = self.ImportFromEPSG(4326) elif isinstance(srs, str_types): # parse as proj4 string status = self.ImportFromProj4(str(srs)) if status > 0: # parse as WKT string status = self.ImportFromWkt(str(srs)) if status > 0: raise NansatProjectionError('Proj4 or WKT (%s) is wrong' % srs) elif isinstance(srs, int): # parse as EPSG code status = self.ImportFromEPSG(srs) if status > 0: raise NansatProjectionError('EPSG %d is wrong' % srs) elif type(srs) in [osr.SpatialReference, NSR]: # parse from input Spatial Reference status = self.ImportFromWkt(srs.ExportToWkt()) if status > 0: raise NansatProjectionError('NSR %s is wrong' % srs)
def __init__(self, srs=0): """Create Spatial Reference System from input parameter""" # create SRS osr.SpatialReference.__init__(self) # parse input parameters status = 1 if srs is 0: # generate default WGS84 SRS status = self.ImportFromWkt(osr.SRS_WKT_WGS84) elif isinstance(srs, str): # parse as proj4 string status = self.ImportFromProj4(str(srs)) if status > 0: # parse as WKT string status = self.ImportFromWkt(str(srs)) if status > 0: raise NansatProjectionError('Proj4 or WKT (%s) is wrong' % srs) # TODO: catch long in python 3 elif isinstance(srs, int): # parse as EPSG code status = self.ImportFromEPSG(srs) if status > 0: raise NansatProjectionError('EPSG %d is wrong' % srs) elif type(srs) in [osr.SpatialReference, NSR]: # parse from input Spatial Reference status = self.ImportFromWkt(srs.ExportToWkt()) if status > 0: raise NansatProjectionError('NSR %s is wrong' % srs)
def __init__(self, srs=None, ext=None, ds=None, **kwargs): """Create Domain from GDALDataset or string options or lat/lon grids""" # If too much information is given raise error if ds is not None and srs is not None and ext is not None: raise ValueError( 'Ambiguous specification of both dataset, srs- and ext-strings.' ) # choose between input opitons: # ds # ds and srs # srs and ext # if only a dataset is given: # copy geo-reference from the dataset if ds is not None and srs is None: self.vrt = VRT.from_gdal_dataset(ds) # If dataset and srs are given (but not ext): # use AutoCreateWarpedVRT to determine bounds and resolution elif ds is not None and srs is not None: srs = NSR(srs) tmp_vrt = gdal.AutoCreateWarpedVRT(ds, None, srs.wkt) if tmp_vrt is None: raise NansatProjectionError( 'Could not warp the given dataset to the given SRS.') else: self.vrt = VRT.from_gdal_dataset(tmp_vrt) # If SpatialRef and extent string are given (but not dataset) elif srs is not None and ext is not None: srs = NSR(srs) # create full dictionary of parameters extent_dict = Domain._create_extent_dict(ext) # convert -lle to -te if 'lle' in extent_dict.keys(): extent_dict = self._convert_extentDic(srs, extent_dict) # get size/extent from the created extent dictionary geo_transform, raster_x_size, raster_y_size = self._get_geotransform( extent_dict) # create VRT object with given geo-reference parameters self.vrt = VRT.from_dataset_params(x_size=raster_x_size, y_size=raster_y_size, geo_transform=geo_transform, projection=srs.wkt, gcps=[], gcp_projection='') elif 'lat' in kwargs and 'lon' in kwargs: warnings.warn( 'Domain(lon=lon, lat=lat) will be deprectaed!' 'Use Domain.from_lonlat()', NansatFutureWarning) # create self.vrt from given lat/lon self.vrt = VRT.from_lonlat(kwargs['lon'], kwargs['lat']) else: raise ValueError('"dataset" or "srsString and extentString" ' 'or "dataset and srsString" are required')
def __init__(self, srs=None, ext=None, ds=None, lon=None, lat=None, name='', logLevel=None): """Create Domain from GDALDataset or string options or lat/lon grids""" # set default attributes self.logger = add_logger('Nansat', logLevel) self.name = name self.logger.debug('ds: %s' % str(ds)) self.logger.debug('srs: %s' % srs) self.logger.debug('ext: %s' % ext) # If too much information is given raise error if ds is not None and srs is not None and ext is not None: raise ValueError( 'Ambiguous specification of both dataset, srs- and ext-strings.' ) # choose between input opitons: # ds # ds and srs # srs and ext # lon and lat # if only a dataset is given: # copy geo-reference from the dataset if ds is not None and srs is None: self.vrt = VRT.from_gdal_dataset(ds) # If dataset and srs are given (but not ext): # use AutoCreateWarpedVRT to determine bounds and resolution elif ds is not None and srs is not None: srs = NSR(srs) tmp_vrt = gdal.AutoCreateWarpedVRT(ds, None, srs.wkt) if tmp_vrt is None: raise NansatProjectionError( 'Could not warp the given dataset to the given SRS.') else: self.vrt = VRT.from_gdal_dataset(tmp_vrt) # If SpatialRef and extent string are given (but not dataset) elif srs is not None and ext is not None: srs = NSR(srs) # create full dictionary of parameters extent_dict = Domain._create_extent_dict(ext) # convert -lle to -te if 'lle' in extent_dict.keys(): extent_dict = self._convert_extentDic(srs, extent_dict) # get size/extent from the created extent dictionary geo_transform, raster_x_size, raster_y_size = self._get_geotransform( extent_dict) # create VRT object with given geo-reference parameters self.vrt = VRT.from_dataset_params(x_size=raster_x_size, y_size=raster_y_size, geo_transform=geo_transform, projection=srs.wkt, gcps=[], gcp_projection='') elif lat is not None and lon is not None: # create self.vrt from given lat/lon self.vrt = VRT.from_lonlat(lon, lat) else: raise ValueError('"dataset" or "srsString and extentString" ' 'or "dataset and srsString" are required') self.logger.debug('vrt.dataset: %s' % str(self.vrt.dataset))