def asgridarray(self):
xdata = self.dims[1].getDimValue()
ydata = self.dims[0].getDimValue()
gdata = GridArray(self.array, xdata, ydata, self.fill_value, self.proj)
return gdata
def savegrid(self, x, y, fname, format='surfer', **kwargs):
gdata = GridArray(self.array, x.array, y.array, -9999.0)
if format == 'surfer':
gdata.saveAsSurferASCIIFile(fname)
elif format == 'bil':
gdata.saveAsBILFile(fname)
elif format == 'esri_ascii':
gdata.saveAsESRIASCIIFile(fname)
elif format == 'micaps4':
desc = kwargs.pop('description', 'var')
date = kwargs.pop('date', datetime.datetime.now())
date = miutil.jdate(date)
hours = kwargs.pop('hours', 0)
level = kwargs.pop('level', 0)
smooth = kwargs.pop('smooth', 1)
boldvalue = kwargs.pop('boldvalue', 0)
proj = kwargs.pop('proj', None)
if proj is None:
gdata.saveAsMICAPS4File(fname, desc, date, hours, level,
smooth, boldvalue)
else:
if proj.isLonLat():
gdata.saveAsMICAPS4File(fname, desc, date, hours, level,
smooth, boldvalue)
else:
gdata.saveAsMICAPS4File(fname, desc, date, hours, level,
smooth, boldvalue, proj)
def savegrid(self, x, y, fname, format='surfer', **kwargs):
gdata = GridArray(self.array, x.array, y.array, -9999.0)
if format == 'surfer':
gdata.saveAsSurferASCIIFile(fname)
elif format == 'bil':
gdata.saveAsBILFile(fname)
elif format == 'esri_ascii':
gdata.saveAsESRIASCIIFile(fname)
elif format == 'micaps4':
desc = kwargs.pop('description', 'var')
date = kwargs.pop('date', datetime.datetime.now())
date = miutil.jdate(date)
hours = kwargs.pop('hours', 0)
level = kwargs.pop('level', 0)
smooth = kwargs.pop('smooth', 1)
boldvalue =kwargs.pop('boldvalue', 0)
proj = kwargs.pop('proj', None)
if proj is None:
gdata.saveAsMICAPS4File(fname, desc, date, hours, level, smooth, boldvalue)
else:
if proj.isLonLat():
gdata.saveAsMICAPS4File(fname, desc, date, hours, level, smooth, boldvalue)
else:
gdata.saveAsMICAPS4File(fname, desc, date, hours, level, smooth, boldvalue, proj)
