# On ajoute un traît de côte pour le masquage from vacumm.bathy.bathy import GriddedBathy, GriddedBathyMerger bathy_large = GriddedBathy(var_large, shoreline='i') bathy_small = var_small # Création de la grille finale de résolution intermédiaire final_grid = create_grid((-6., .5, xr*2.5), (45, 47., yr*2.5)) # On crée maintenant le merger merger = GriddedBathyMerger(final_grid) # - ajout de la bathy basse resolution en premier (en dessous) merger += bathy_large # - puis ajout de celle haute résolution merger += bathy_small # On définit le traît de côte pour le masquage merger.set_shoreline('h') # Fusion vers la grille finale bathy = merger.merge() # Plot merger.plot(show=False) from vacumm.misc.plot import savefigs, add_grid kwgrid = dict(linewidth=.5, alpha=.5, samp=2) add_grid(grid_large, color='r', **kwgrid) add_grid(var_small.getGrid(), color='#00ff00', **kwgrid) savefigs(__file__)
# On ajoute un traît de côte pour le masquage from vacumm.bathy.bathy import GriddedBathy, GriddedBathyMerger bathy_large = GriddedBathy(var_large, shoreline='i') bathy_small = var_small # Création de la grille finale de résolution intermédiaire final_grid = create_grid((-6., .5, xr * 2.5), (45, 47., yr * 2.5)) # On crée maintenant le merger merger = GriddedBathyMerger(final_grid) # - ajout de la bathy basse resolution en premier (en dessous) merger += bathy_large # - puis ajout de celle haute résolution merger += bathy_small # On définit le traît de côte pour le masquage merger.set_shoreline('h') # Fusion vers la grille finale bathy = merger.merge() # Plot merger.plot(show=False) from vacumm.misc.plot import savefigs, add_grid kwgrid = dict(linewidth=.5, alpha=.5, samp=2) add_grid(grid_large, color='r', **kwgrid) add_grid(var_small.getGrid(), color='#00ff00', **kwgrid) savefigs(__file__) P.close()
# Plot bathy.plot(title='Direct') # Bathy standard from vacumm.bathy.bathy import NcGriddedBathy, bathy_list print bathy_list().keys() sbathy = NcGriddedBathy(lon=(-5.2, -3), lat=(48., 48.9), name='etopo2') sbathy.plot(title='Bathy par defaut') # Fusion de bathy mgrid = create_grid((-6.1, -3, 0.01), (47.8, 48.9, .006)) # -> changer la grille from vacumm.bathy.bathy import GriddedBathyMerger merger = GriddedBathyMerger(mgrid) merger += sbathy merger += bathy merger.plot(title='Fusion / plot direct') # Plot via plot_bathy mbathy = merger.merge() plot_bathy(mbathy, title='Fusion / plot via plot_bathy') # Regrillage rgrid = create_grid((-6.1, -3, 0.05), (47.8, 48.9, .02)) Mbathy = GriddedBathy(mbathy) rbathy = Mbathy.regrid(rgrid) plot_bathy(rbathy, title='Regrillee')