print 'Resolution of auto grid:', resol(grid_auto)
# %% Interpolation on auto grid
print 'Interpolation:'
gridded_auto = xyz.togrid()
# equivalent a :
# >>> gridded_auto = xyz.togrid(xyz.grid)
# %% Interpolation on manual grid
print 'Interpolation and masking, then extraction'
# - defintion
grid_manual = create_grid((-5.3, -4.91, .01), (48.1, 48.41, .01))
# - interpolation
gridded_manual = xyz.togrid(grid_manual, mask='h')
# - extraction with margin
xyz_up = xyz.clip(zone=(None, None, None, 48.3), margin=2)
# - si None, valeurs limites internes (i.e xyz.xmin(), ...)
# - margin : marge relative en unite de resolution
# -> ici : ymax = 48.3 + xyz.resol()[1]*2
# %% Save
print 'Save'
prefix = __file__[:-2] + 'up'
xyz_up.save(prefix + '.xyz') # ascii
xyz_up.save(prefix + '.nc') # netcdf/grd
# %% Plots
print 'Plots'
# - init
P.figure(figsize=(4.5, 8))
P.rc('font', size=8)
print 'Resolution grille auto :', resol(grid_auto)
# Interpolation sur grille auto
print 'Interpolation auto'
gridded_auto = xyz.togrid()
# equivalent a :
# >>> gridded_auto = xyz.togrid(xyz.grid)
# Interpolation sur grille manuelle
print 'Interpolation et masquage manuels puis extraction'
# - defintion de la grille
grid_manual = create_grid((-5.3, -4.91, .01), (48.1, 48.41, .01))
# - interpolation
gridded_manual = xyz.togrid(grid_manual, mask='h')
# Extraction d'une sous-zone (xmin,ymin...) avec marge
xyz_up = xyz.clip(zone=(None, None, None, 48.3), margin=2)
# - si None, valeurs limites internes (i.e xyz.xmin(), ...)
# - margin : marge relative en unite de resolution
# -> ici : ymax = 48.3 + xyz.resol()[1]*2
# Sauvegarde
print 'Sauvegarde'
prefix = __file__[:-2]+'up'
xyz_up.save(prefix+'.xyz') # ascii
xyz_up.save(prefix+'.nc') # netcdf/grd
# Plots
print 'Plots'
# - init
P.figure(figsize=(4.5, 8))
P.rc('font', size=8)
