print xyz1.xmin # -> essayer get_xmin avec mask
print xyz1.resol() # -> essayer avec deg=...
# Plot
kwp = dict(size=40, map_res=None, masked_alpha=.1)
xyz1.plot(mode='both', **kwp) # -> essayer autre mode
# Fichier
fbathy = __file__[:-2] + 'xyz'
xyz1.save(fbathy)
xyz3 = XYZBathy(fbathy)
xyz3.plot(title='XYZ3', **kwp)
# Fusions
xyz = xyz1 + xyz2
merger = XYZBathyMerger()
merger += xyz1 # -> supprimer/ajouter
merger.append(xyz2)
for i, b in enumerate(merger):
b.long_name = 'Niveau : %i' % i # Pour les legendes
b.set_transp(False) # Opacite
merger.plot(mode='cluster', **kwp) # -> essayer autre mode
xyz = merger.get_xyz(long_name='Merged')
xyz.plot(**kwp)
# Regrillage
grid_auto = xyz.grid
grid_manual = create_grid((-5.3, -4.91, .01), (48.1, 48.41, .01))
gridded = xyz.togrid(grid_manual)
map2(gridded, res=None)
P.rcParams['font.size'] = 9
P.subplots_adjust(bottom=.03, top=.97, hspace=.25)
kwplot = dict(show=False,
colorbar=False,
map_res=None,
margin=0.,
map_autoresize=0,
xmin=xyz.xmin,
xmax=xyz.xmax,
ymin=xyz.ymin,
ymax=xyz.ymax)
xyz.plot(title='Fusion directe', **kwplot)
# Utilisation d'un fusionneur
# - init
merger = XYZBathyMerger()
# - ajout de xyz direct
merger += xyz1
merger.append(xyz2)
merger += xyz3
merger += xyz4
# - ajout a partir d'un fichier
merger += fxyz5
# - suppression d'un xyz
merger -= xyz2
# - suppression du dernier dataset
del merger[-1]
# - changement d'autres attributs
for i, xyz in enumerate(merger):
xyz.long_name = 'Niveau : %i' % i # Pour les legendes
xyz.set_transp(False) # Opacite
N.savetxt(fxyz5, gene_bathy(-5.15, 48.2, .2, .1, amp=15).transpose()) # center
# Fusion directe
xyz = xyz1 + xyz2 + xyz3 + xyz4 + fxyz5
# Plot
import pylab as P ; P.figure(figsize=(5., 8.5)) ; P.subplot(311)
P.rcParams['font.size'] = 9
P.subplots_adjust(bottom=.03, top=.97, hspace=.25)
kwplot = dict(show=False, colorbar=False,map_res=None, margin=0., map_autoresize=0,
xmin=xyz.xmin, xmax=xyz.xmax, ymin=xyz.ymin, ymax=xyz.ymax)
xyz.plot(title='Fusion directe', **kwplot)
# Utilisation d'un fusionneur
# - init
merger = XYZBathyMerger()
# - ajout de xyz direct
merger += xyz1
merger.append(xyz2)
merger += xyz3
merger += xyz4
# - ajout a partir d'un fichier
merger += fxyz5
# - suppression d'un xyz
merger -= xyz2
# - suppression du dernier dataset
del merger[-1]
# - changement d'autres attributs
for i, xyz in enumerate(merger):
xyz.long_name = 'Niveau : %i'%i # Pour les legendes
xyz.set_transp(False) # Opacite
print xyz1.resol() # -> essayer avec deg=...
# Plot
kwp = dict(size=40, map_res=None, masked_alpha=0.1)
xyz1.plot(mode="both", **kwp) # -> essayer autre mode
# Fichier
fbathy = __file__[:-2] + "xyz"
xyz1.save(fbathy)
xyz3 = XYZBathy(fbathy)
xyz3.plot(title="XYZ3", **kwp)
# Fusions
xyz = xyz1 + xyz2
merger = XYZBathyMerger()
merger += xyz1 # -> supprimer/ajouter
merger.append(xyz2)
for i, b in enumerate(merger):
b.long_name = "Niveau : %i" % i # Pour les legendes
b.set_transp(False) # Opacite
merger.plot(mode="cluster", **kwp) # -> essayer autre mode
xyz = merger.get_xyz(long_name="Merged")
xyz.plot(**kwp)
# Regrillage
grid_auto = xyz.grid
grid_manual = create_grid((-5.3, -4.91, 0.01), (48.1, 48.41, 0.01))
gridded = xyz.togrid(grid_manual)
map2(gridded, res=None)
