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)