# 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')
