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