def test_gmap(): g = GoogleCenterMap(center_ll=(10.762660, 46.794221), zoom=13, size_x=640, size_y=640) m = Map(g.grid, countries=False, nx=640) m.set_lonlat_countours(interval=0.025) m.set_shapefile(get_demo_file('Hintereisferner.shp'), linewidths=2, edgecolor='darkred') m.set_rgb(g.get_vardata()) m.visualize(addcbar=False) dem = salem.GeoTiff(get_demo_file('hef_srtm.tif')) dem.set_subset(margin=-100) dem = salem.grids.local_mercator_grid(center_ll=(10.76, 46.798444), extent=(10000, 7000)) i, j = dem.ij_coordinates g = GoogleVisibleMap(x=i, y=j, src=dem, size_x=500, size_y=400) img = g.get_vardata() m = Map(dem, countries=False) assert_raises(ValueError, m.set_data, img) m.set_lonlat_countours(interval=0.025) m.set_shapefile(get_demo_file('Hintereisferner.shp'), linewidths=2, edgecolor='darkred') m.set_rgb(img, g.grid) m.visualize(addcbar=False)
def test_center(self): import matplotlib as mpl gm = GoogleCenterMap(center_ll=(10.762660, 46.794221), zoom=13, size_x=500, size_y=500, use_cache=False) gm.set_roi(shape=get_demo_file('Hintereisferner.shp')) gm.set_subset(toroi=True, margin=10) img = gm.get_vardata()[..., :3] img[np.nonzero(gm.roi == 0)] /= 2. # from PIL import Image # Image.fromarray((img * 255).astype(np.uint8)).save( # get_demo_file('hef_google_roi.png')) ref = mpl.image.imread(get_demo_file('hef_google_roi.png')) rmsd = np.sqrt(np.mean((ref - img)**2)) self.assertTrue(rmsd < 0.2) # assert_allclose(ref, img, atol=2e-2) gm = GoogleCenterMap(center_ll=(10.762660, 46.794221), zoom=13, size_x=500, size_y=500) gm.set_roi(shape=get_demo_file('Hintereisferner.shp')) gm.set_subset(toroi=True, margin=10) img = gm.get_vardata()[..., :3] img[np.nonzero(gm.roi == 0)] /= 2. rmsd = np.sqrt(np.mean((ref - img)**2)) self.assertTrue(rmsd < 0.2) gm = GoogleCenterMap(center_ll=(10.762660, 46.794221), zoom=13, size_x=500, size_y=500) gm2 = GoogleCenterMap(center_ll=(10.762660, 46.794221), zoom=13, size_x=500, size_y=500, scale=2) assert (gm.grid.nx * 2) == gm2.grid.nx assert gm.grid.extent == gm2.grid.extent
def test_center(self): gm = GoogleCenterMap(center_ll=(10.762660, 46.794221), zoom=13, size_x=500, size_y=500, use_cache=False) gm.set_roi(shape=get_demo_file('Hintereisferner.shp')) gm.set_subset(toroi=True, margin=10) img = gm.get_vardata() img[np.nonzero(gm.roi == 0)] /= 2. # from scipy.misc import toimage # toimage(img).save(get_demo_file('hef_google_roi.png')) ref = mpl.image.imread(get_demo_file('hef_google_roi.png')) rmsd = np.sqrt(np.mean((ref - img)**2)) self.assertTrue(rmsd < 0.1) # assert_allclose(ref, img, atol=2e-2) gm = GoogleCenterMap(center_ll=(10.762660, 46.794221), zoom=13, size_x=500, size_y=500) gm.set_roi(shape=get_demo_file('Hintereisferner.shp')) gm.set_subset(toroi=True, margin=10) img = gm.get_vardata() img[np.nonzero(gm.roi == 0)] /= 2. rmsd = np.sqrt(np.mean((ref - img)**2)) self.assertTrue(rmsd < 0.1)
def test_center(self): gm = GoogleCenterMap(center_ll=(10.762660, 46.794221), zoom=13, size_x=500, size_y=500, use_cache=False) gm.set_roi(shape=get_demo_file("Hintereisferner.shp")) gm.set_subset(toroi=True, margin=10) img = gm.get_vardata() img[np.nonzero(gm.roi == 0)] /= 2.0 # from scipy.misc import toimage # toimage(img).save(get_demo_file('hef_google_roi.png')) ref = mpl.image.imread(get_demo_file("hef_google_roi.png")) rmsd = np.sqrt(np.mean((ref - img) ** 2)) self.assertTrue(rmsd < 0.1) # assert_allclose(ref, img, atol=2e-2) gm = GoogleCenterMap(center_ll=(10.762660, 46.794221), zoom=13, size_x=500, size_y=500) gm.set_roi(shape=get_demo_file("Hintereisferner.shp")) gm.set_subset(toroi=True, margin=10) img = gm.get_vardata() img[np.nonzero(gm.roi == 0)] /= 2.0 rmsd = np.sqrt(np.mean((ref - img) ** 2)) self.assertTrue(rmsd < 0.1)