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))