def plot_grid(lon, lat, dhdt, cell=None, plot_cell=False, contourf=True): fig = plt.figure() m = viz.make_proj_stere(BBOX_REG) plon, plat = lon[GRID_CELL[1]], lat[GRID_CELL[0]] lon, lat = np.meshgrid(lon, lat) xx, yy = m(lon, lat) x, y = m(plon, plat) print 'point cell (lon lat x y):', plon, plat, x, y dhdt = np.ma.masked_invalid(dhdt) m.pcolormesh(xx, yy, dhdt, cmap=CMAP, vmin=-ABS_VAL, vmax=ABS_VAL, alpha=1) m.plot(x, y, 'mo', markersize=8) if plot_cell: lon, lat = util.box(cell) x, y = m(lon, lat) m.plot(x, y, 'k', linewidth=2) if contourf: fig2 = plt.figure() m.contourf(xx, yy, dhdt, 25, cmap=CMAP) #plt.colorbar().set_label(LABEL) viz.colorbar(fig, CMAP, (-ABS_VAL, ABS_VAL)) ''' w = 50.0 / (axis[1] - axis[0]) rect = 0.142 - w, 0.08, 2 * w, 0.02 cmap = CMAP viz.colorbar(fig, cmap, colorlim, legend, rect, ticks, ticklabels, size=7, weight='bold', color='k', edgecolor='w') leg = fig.add_axes([0, 0, 1, 1]) leg.set_axis_off() viz.text(leg, 0.70, 0.90, title, ha='center', va='top', size=10, weight='bold', color='k', edgecolor='w') ''' return fig
def plot_grid(lon, lat, dhdt, cell=None, plot_cell=False, contourf=True): fig = plt.figure() m = viz.make_proj_stere(BBOX_REG) plon, plat = lon[GRID_CELL[1]], lat[GRID_CELL[0]] lon, lat = np.meshgrid(lon, lat) xx, yy = m(lon, lat) x, y = m(plon, plat) print 'point cell (lon lat x y):', plon, plat, x, y dhdt = np.ma.masked_invalid(dhdt) m.pcolormesh(xx, yy, dhdt, cmap=CMAP, vmin=-ABS_VAL, vmax=ABS_VAL, alpha=1) m.plot(x, y, 'mo', markersize=8) if plot_cell: lon, lat = util.box(cell) x, y = m(lon, lat) m.plot(x, y, 'k', linewidth=2) if contourf: fig2 = plt.figure() m.contourf(xx, yy, dhdt, 25, cmap=CMAP) #plt.colorbar().set_label(LABEL) viz.colorbar(fig, CMAP, (-ABS_VAL,ABS_VAL)) ''' w = 50.0 / (axis[1] - axis[0]) rect = 0.142 - w, 0.08, 2 * w, 0.02 cmap = CMAP viz.colorbar(fig, cmap, colorlim, legend, rect, ticks, ticklabels, size=7, weight='bold', color='k', edgecolor='w') leg = fig.add_axes([0, 0, 1, 1]) leg.set_axis_off() viz.text(leg, 0.70, 0.90, title, ha='center', va='top', size=10, weight='bold', color='k', edgecolor='w') ''' return fig
def plot_grid(lon, lat, dhdt, cell=None, plot_cell=False, contourf=True, savefig=None): fig = plt.figure() m = viz.make_proj_stere(BBOX_REG) plon, plat = lon[GRID_CELL[1]], lat[GRID_CELL[0]] # shift grid half degree lon -= (lon[1] - lon[0]) / 2. lat -= (lat[1] - lat[0]) / 2. lon, lat = np.meshgrid(lon, lat) xx, yy = m(lon, lat) x, y = m(plon, plat) print 'point cell (lon lat x y):', plon, plat, x, y dhdt = np.ma.masked_invalid(dhdt) m.pcolormesh(xx, yy, dhdt, cmap=CMAP, vmin=-ABS_VAL, vmax=ABS_VAL, alpha=1) p_ = m.drawparallels(np.arange(-90., -60, 6), labels=[0, 0, 0, 1], color='0.8') m_ = m.drawmeridians(np.arange(-180, 180., 12), labels=[1, 0, 0, 0], color='0.8') #m.plot(x, y, 'mo', markersize=8) if plot_cell: lon, lat = util.box(cell) x, y = m(lon, lat) m.plot(x, y, 'k', linewidth=2) if contourf: fig2 = plt.figure() m.contourf(xx, yy, dhdt, 25, cmap=CMAP) ### ''' x2, y2 = np.loadtxt('/Users/fpaolo/data/masks/scripps/scripps_iceshelves_v1_geod.txt', usecols=(0,1), unpack=True, comments='%') x2, y2 = m(x2, y2) m.plot(x2, y2, 'o', markersize=0.8) ''' ### #plt.colorbar().set_label(LABEL) #viz.colorbar(fig, CMAP, (-ABS_VAL,ABS_VAL), title='m/yr') ''' w = 50.0 / (axis[1] - axis[0]) rect = 0.142 - w, 0.08, 2 * w, 0.02 cmap = CMAP viz.colorbar(fig, cmap, colorlim, legend, rect, ticks, ticklabels, size=7, weight='bold', color='k', edgecolor='w') leg = fig.add_axes([0, 0, 1, 1]) leg.set_axis_off() viz.text(leg, 0.70, 0.90, title, ha='center', va='top', size=10, weight='bold', color='k', edgecolor='w') ''' if savefig is not None: plt.savefig(savefig + '_map.png') return fig
def plot_map(lon, lat, RR, bbox, mfile, mres=1): fig1 = plt.figure() ax = fig1.add_subplot(111) m = viz.make_proj_stere(bbox) m = viz.plot_grid_proj(m, lon, lat, RR, contourf=False, shift=True) plt.colorbar(orientation='vertical', shrink=.7) x, y, mask = viz.get_gtif_subreg(m, mfile, mres) mask = np.ma.masked_values(mask, 4) m.imshow(mask, zorder=1, cmap=plt.cm.gray_r) #m = plot_boundaries(m, (lon.min(),lon.max(),lat.min()-10,lat.max()+3)) p_ = m.drawparallels(np.arange(-90.,-70.,3), labels=[0,0,1,0], color='0.6') m_ = m.drawmeridians(np.arange(-80.,-20.,5), labels=[1,0,0,0], color='0.6') return fig1
def plot_map(lon, lat, grid, bbox, mfile, mres=1, **kw): fig1 = plt.figure() ax = fig1.add_subplot(111) m = viz.make_proj_stere(bbox) m = viz.plot_grid_proj(m, lon, lat, grid, contourf=False, shift=True, **kw) #plt.colorbar(orientation='vertical', shrink=.7) #x, y, mask = viz.get_gtif_subreg(m, mfile, mres) #mask = np.ma.masked_values(mask, 4) #m.imshow(mask, zorder=0, cmap=plt.cm.gray_r) #m = plot_boundaries(m, (lon.min(),lon.max(),lat.min()-10,lat.max()+3)) #p_ = m.drawparallels(np.arange(-90.,-70.,3), labels=[0,0,0,1], color='0.6') #m_ = m.drawmeridians(np.arange(-80.,-20.,5), labels=[1,0,0,0], color='0.6') #viz.colorbar(fig1, CMAP, (-ABS_VAL,ABS_VAL)) return fig1
def plot_grid(lon, lat, dhdt, cell=None, plot_cell=False, contourf=True, savefig=None): fig = plt.figure() m = viz.make_proj_stere(BBOX_REG) plon, plat = lon[GRID_CELL[1]], lat[GRID_CELL[0]] # shift grid half degree lon -= (lon[1] - lon[0])/2. lat -= (lat[1] - lat[0])/2. lon, lat = np.meshgrid(lon, lat) xx, yy = m(lon, lat) x, y = m(plon, plat) print 'point cell (lon lat x y):', plon, plat, x, y dhdt = np.ma.masked_invalid(dhdt) m.pcolormesh(xx, yy, dhdt, cmap=CMAP, vmin=-ABS_VAL, vmax=ABS_VAL, alpha=1) p_ = m.drawparallels(np.arange(-90.,-60, 2), labels=[1,0,0,0], color='0.3') m_ = m.drawmeridians(np.arange(-180,180., 15), labels=[0,0,1,0], color='0.3') #m.plot(x, y, 'mo', markersize=8) if plot_cell: lon, lat = util.box(cell) x, y = m(lon, lat) m.plot(x, y, 'k', linewidth=2) if contourf: fig2 = plt.figure() m.contourf(xx, yy, dhdt, 25, cmap=CMAP) ### ''' x2, y2 = np.loadtxt('/Users/fpaolo/data/masks/scripps/scripps_iceshelves_v1_geod.txt', usecols=(0,1), unpack=True, comments='%') x2, y2 = m(x2, y2) m.plot(x2, y2, 'o', markersize=1.2) ''' ### #plt.colorbar().set_label(LABEL) rect = 0.25, 0.11, 0.5, 0.02 viz.colorbar(fig, CMAP, (-ABS_VAL,ABS_VAL), title='Elevation change rate (m/yr)', rect=rect) ''' w = 50.0 / (axis[1] - axis[0]) rect = 0.142 - w, 0.08, 2 * w, 0.02 cmap = CMAP viz.colorbar(fig, cmap, colorlim, legend, rect, ticks, ticklabels, size=7, weight='bold', color='k', edgecolor='w') leg = fig.add_axes([0, 0, 1, 1]) leg.set_axis_off() viz.text(leg, 0.70, 0.90, title, ha='center', va='top', size=10, weight='bold', color='k', edgecolor='w') ''' if savefig is not None: plt.savefig(savefig+'_map.png') return fig
def plot_map(lon, lat, grid, bbox, mfile, mres=1, **kw): """ **kw : keyword args contourf=[True/False] vmin=int vmax=int """ fig1 = plt.figure() ax = fig1.add_subplot(111) m = viz.make_proj_stere(bbox) m = viz.plot_grid_proj(m, lon, lat, grid, shift=True, **kw) plt.colorbar(orientation='vertical', shrink=.65) x, y, mask = viz.get_gtif_subreg(m, mfile, mres) mask = np.ma.masked_values(mask, 4) m.imshow(mask, zorder=1, cmap=plt.cm.gray_r) #m = plot_boundaries(m, (lon.min(),lon.max(),lat.min()-10,lat.max()+3)) p_ = m.drawparallels(np.arange(-90.,-60,3), labels=[0,0,0,1], color='0.2') m_ = m.drawmeridians(np.arange(-180,180.,5), labels=[1,0,0,0], color='0.2') return fig1
def plot_grid(lon, lat, dhdt, cell=None, plot_cell=False, contourf=True): fig = plt.figure() m = viz.make_proj_stere(BBOX_REG) lon, lat = np.meshgrid(lon, lat) xx, yy = m(lon, lat) dhdt = np.ma.masked_invalid(dhdt) m.pcolormesh(xx, yy, dhdt, cmap=CMAP, vmin=-ABS_VAL, vmax=ABS_VAL, alpha=1) p_ = m.drawparallels(np.arange(-90.,-60, 5), labels=[0,0,0,1], color='0.8') m_ = m.drawmeridians(np.arange(-180,180., 10), labels=[1,0,0,0], color='0.8') ### ''' xx, yy = np.loadtxt('/Users/fpaolo/data/masks/scripps/scripps_iceshelves_v1_geod.txt', usecols=(0,1), unpack=True, comments='%') xx, yy = m(xx, yy) m.plot(xx, yy, 'o', markersize=0.1) ''' ### #m.plot(x, y, 'mo', markersize=8) if plot_cell: lon, lat = util.box(cell) x, y = m(lon, lat) m.plot(x, y, 'k', linewidth=2) if contourf: fig2 = plt.figure() m.contourf(xx, yy, dhdt, 25, cmap=CMAP) #plt.colorbar().set_label(LABEL) viz.colorbar(fig, CMAP, (-ABS_VAL,ABS_VAL)) ''' w = 50.0 / (axis[1] - axis[0]) rect = 0.142 - w, 0.08, 2 * w, 0.02 cmap = CMAP viz.colorbar(fig, cmap, colorlim, legend, rect, ticks, ticklabels, size=7, weight='bold', color='k', edgecolor='w') leg = fig.add_axes([0, 0, 1, 1]) leg.set_axis_off() viz.text(leg, 0.70, 0.90, title, ha='center', va='top', size=10, weight='bold', color='k', edgecolor='w') ''' return fig
def plot_map(lon, lat, grid, bbox, mfile, mres=1, **kw): """ **kw : keyword args contourf=[True/False] vmin=int vmax=int """ fig1 = plt.figure() ax = fig1.add_subplot(111) m = viz.make_proj_stere(bbox) m = viz.plot_grid_proj(m, lon, lat, grid, shift=True, **kw) plt.colorbar(orientation='vertical', shrink=.65) x, y, mask = viz.get_gtif_subreg(m, mfile, mres) mask = np.ma.masked_values(mask, 4) m.imshow(mask, zorder=1, cmap=plt.cm.gray_r) #m = plot_boundaries(m, (lon.min(),lon.max(),lat.min()-10,lat.max()+3)) p_ = m.drawparallels(np.arange(-90., -60, 3), labels=[0, 0, 0, 1], color='0.2') m_ = m.drawmeridians(np.arange(-180, 180., 5), labels=[1, 0, 0, 0], color='0.2') return fig1