def __init__(self, filename=None):
# set default spatial reference
sr = osr.SpatialReference()
sr.ImportFromEPSG(4326)
self.srs = sr.ExportToWkt()
Base.__init__(self,filename)
def __init__(self, filename=None):
Base.__init__(self, filename)
# set default spatial reference
sr = osr.SpatialReference()
sr.ImportFromEPSG(6933) # EASE-2.0
#sr.ImportFromEPSG(4326) # WGS84
self.srs = sr.ExportToWkt()
def __init__(self, filename=None):
Base.__init__(self, filename)
# set default spatial reference
sr = osr.SpatialReference()
sr.ImportFromEPSG(4326)
self.srs = sr.ExportToWkt()
self.doy = 0 # day of the year to retrieve
self.ssh = None
def get_lon_lat(self, ds):
path = ds.GetDescription()
if path.endswith('1km'):
lon = self.get_data(Base.get_dataset(self, LON1K))
lat = self.get_data(Base.get_dataset(self, LAT1K))
elif path.endswith('3km'):
lon = self.get_data(Base.get_dataset(self, LON3K))
lat = self.get_data(Base.get_dataset(self, LAT3K))
else:
raise
return (lon, lat)
def get_dataset(self, name):
ds = Base.get_dataset(self, name)
if ds:
if not ds.GetProjection():
ds.SetProjection(self.wgs84)
ds.SetGeoTransform([-180, 0.1, 0, -90, 0, 0.1])
return ds
def get_dataset(self, name):
ds = Base.get_dataset(self, name)
if ds:
self.set_projection(ds)
return ds
def connect(self, host=NTSG, port=0, timeout=-999):
Base.connect(self, host, port, timeout)
def get_dataset(self, name):
ds = Base.get_dataset(self, name)
if ds:
if not ds.GetProjection():
ds.SetProjection(self.wgs84)
return ds