def plotMap(self):
"""Plot return period wind speed maps"""
lon, lat, years, inputData = self.loadFile(self.inputFile, 'wspd')
if lon.min() > 180.:
lon = lon - 360.
[xgrid, ygrid] = np.meshgrid(lon, lat)
dmask = inputData.mask
inputData = metutils.convert(inputData, 'mps', self.plotUnits.units)
inputData = ma.array(inputData, mask=dmask)
map_kwargs = dict(llcrnrlon=xgrid.min(),
llcrnrlat=ygrid.min(),
urcrnrlon=xgrid.max(),
urcrnrlat=ygrid.max(),
projection='merc',
resolution='i')
for i, year in enumerate(years):
log.debug("Plotting %d-year return period hazard map", year)
title = '%d-Year ARI Cyclonic Wind Hazard' % (year)
imageFilename = '%d_yrRP_hazard_map.png' % (year)
filename = pjoin(self.plotPath, imageFilename)
cbarlab = "Wind speed (%s)"%self.plotUnits.units
levels = self.plotUnits.levels
if self.smooth:
dx = np.mean(np.diff(xgrid))
data = smooth(inputData[i, :, :], int(1/dx))
else:
data = inputData[i, :, :]
saveHazardMap(data, xgrid, ygrid, title, levels,
cbarlab, map_kwargs, filename)
self.progressbar.update((i + 1) / float(len(years)), 0.0, 0.9)
def plotMap(self):
"""Plot return period wind speed maps"""
lon, lat, years, inputData = self.loadFile(self.inputFile, 'wspd')
if lon.min() > 180.:
lon = lon - 360.
[xgrid, ygrid] = np.meshgrid(lon, lat)
inputData = metutils.convert(inputData, 'mps', self.plotUnits.units)
map_kwargs = dict(llcrnrlon=xgrid.min(),
llcrnrlat=ygrid.min(),
urcrnrlon=xgrid.max(),
urcrnrlat=ygrid.max(),
projection='merc',
resolution='i')
for i, year in enumerate(years):
log.debug("Plotting %d-year return period hazard map"%(year))
title = '%d-Year Return Period Cyclonic Wind Hazard' % (year)
imageFilename = '%d_yrRP_hazard_map.png' % (year)
filename = pjoin(self.plotPath, imageFilename)
cbarlab = "Wind speed (%s)"%self.plotUnits.units
levels = self.plotUnits.levels
saveHazardMap(inputData[i, :, :], xgrid, ygrid, title, levels,
cbarlab, map_kwargs, filename)
self.progressbar.update((i + 1) / float(len(years)), 0.0, 0.9)