# - regrillage de la large vers une grille moins fine from vacumm.misc.grid import regridding, resol, create_grid xr, yr = resol(var_large.getGrid()) grid_large = create_grid((-7., -1, xr*4.5), (46., 49, yr*4.5)) var_large = regridding.regrid2d(var_large, grid_large) # 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)
# - regrillage de la large vers une grille moins fine from vacumm.misc.grid import regridding, resol, create_grid xr, yr = resol(var_large.getGrid()) grid_large = create_grid((-7., -1, xr * 4.5), (46., 49, yr * 4.5)) var_large = regridding.regrid2d(var_large, grid_large) # 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)
# 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')