def test_area_def2basemap(self): area_def = utils.parse_area_file(os.path.join(os.path.dirname(__file__), 'test_files', 'areas.cfg'), 'ease_sh')[0] bmap = plot.area_def2basemap(area_def) self.assertTrue(bmap.rmajor == bmap.rminor and bmap.rmajor == 6371228.0, 'Failed to create Basemap object')
def saveDebug_pyresample(weight_sum, filename): """ Save image for debugging onto map with grid and coastlines using pyresample """ import matplotlib matplotlib.use('AGG', warn=False) from pyresample import plot import matplotlib.pyplot as plt import matplotlib.cm as cm new_image = np.array(255.0 * weight_sum / np.max(weight_sum), 'uint8') bmap = plot.area_def2basemap(pc(Satellites.outwidth, Satellites.outheight), resolution='c') bmap.drawcoastlines(linewidth=0.2, color='red') bmap.drawmeridians(np.arange(-180, 180, 10), linewidth=0.2, color='red') bmap.drawparallels(np.arange(-90, 90, 10), linewidth=0.2, color='red') bmap.imshow(new_image, origin='upper', vmin=0, vmax=255, cmap=cm.Greys_r) # @UndefinedVariable # plt.show() i = datetime.datetime.utcnow() plt.text(0, -75, "%s UTC" % i.isoformat(), color='green', size=15, horizontalalignment='center') plt.savefig(filename, bbox_inches='tight', pad_inches=0, dpi=400) plt.close()
def test_area_def2basemap(self): area_def = utils.parse_area_file( os.path.join(os.path.dirname(__file__), 'test_files', 'areas.cfg'), 'ease_sh')[0] bmap = plot.area_def2basemap(area_def) self.assertTrue( bmap.rmajor == bmap.rminor and bmap.rmajor == 6371228.0, 'Failed to create Basemap object')
def test_area_def2basemap(self): """Test the area to Basemap object conversion function.""" from pyresample import plot from pyresample import parse_area_file area_def = parse_area_file(os.path.join(os.path.dirname(__file__), 'test_files', 'areas.yaml'), 'ease_sh')[0] bmap = plot.area_def2basemap(area_def) self.assertTrue(bmap.rmajor == bmap.rminor and bmap.rmajor == 6371228.0, 'Failed to create Basemap object')
def saveDebug(new_image, filename): import matplotlib matplotlib.use('AGG', warn=False) from pyresample import plot import matplotlib.pyplot as plt import matplotlib.cm as cm bmap = plot.area_def2basemap(SatelliteData.pc(), resolution='c') bmap.drawcoastlines() bmap.drawmeridians(np.arange(-180, 180, 45)) bmap.drawparallels(np.arange(-90, 90, 10)) col = bmap.imshow(new_image, origin='upper', vmin=None, vmax=None, cmap = cm.Greys_r) plt.savefig(filename, bbox_inches='tight', pad_inches=0, dpi = 200) plt.close()
def saveDebug_pyresample(weight_sum, filename): """ Save image for debugging onto map with grid and coastlines using pyresample """ import matplotlib matplotlib.use('AGG', warn=False) from pyresample import plot import matplotlib.pyplot as plt import matplotlib.cm as cm new_image = np.array(255.0 * weight_sum / np.max(weight_sum), 'uint8') bmap = plot.area_def2basemap(SatelliteData.pc(), resolution='c') bmap.drawcoastlines(linewidth=0.2, color='red') bmap.drawmeridians(np.arange(-180, 180, 45), linewidth=0.2, color='red') bmap.drawparallels(np.arange(-90, 90, 10), linewidth=0.2, color='red') bmap.imshow(new_image, origin='upper', vmin=0, vmax=255, cmap=cm.Greys_r) # @UndefinedVariable # plt.show() plt.savefig(filename, bbox_inches='tight', pad_inches=0, dpi=400) plt.close()
print(area_def) # ### plot the reprojected image # # pyresample can take its projection structure and turn it into a basemap instance using # [area_def2basemap](https://github.com/pytroll/pyresample/blob/3fce7fc832e0b86369d452d504a21fb65205435c/pyresample/plot.py#L91) # # This works because both pyresample and basemap use the same underlying projection code called # [pyproj](https://github.com/jswhit/pyproj). # In[4]: plt.close('all') fig,ax = plt.subplots(1,1, figsize=(8,8)) bmap=area_def2basemap(area_def,ax=ax,resolution='c') num_meridians=180 num_parallels = 90 vmin=None; vmax=None col = bmap.imshow(result, origin='upper', vmin=0, vmax=0.4) label='channel 1' bmap.drawmeridians(np.arange(-180, 180, num_meridians),labels=[True,False,False,True]) bmap.drawparallels(np.arange(-90, 90, num_parallels),labels=[False,True,True,False]) bmap.drawcoastlines() fig.colorbar(col, shrink=0.5, pad=0.05).set_label(label) fig.canvas.draw() plt.show() # ## Why is this image clipped on the right? #
800, \ [x_ll, x_ur, y_ll, y_ur]) swath_def = geometry.SwathDefinition(lon, lat) grid_def = geometry.GridDefinition(lons=lon, lats=lat) overlap_fraction = swath_def.overlaps(area_def) overlap_fraction = swath_def.overlap_rate(area_def) result = kd_tree.resample_nearest(swath_def, S_VV_ABS, area_def, \ radius_of_influence=200, fill_value=None) bmap = plot.area_def2basemap(area_def) m = Basemap(projection='ortho',lat_0=lat.mean(),lon_0=lon.mean(),resolution='i') # compute map projection coordinates of grid. lons, lats = m.makegrid(RasterYSize+1, RasterXSize+1) x, y = m(lons, lats)