# -*- coding: utf8 -*- # Definition de la region de travail Bretagne zone = (-6.5, 47.2, -2, 49.05) # Recuperation du trait Europeen from vacumm.bathy.shorelines import GSHHS tc = GSHHS(input='h', clip=zone) # Interpolation du niveau de la mer sur le trait xyz = tc.bathy() # Plot xyz.plot(size=10, savefigs=__file__, show=False)#, savefigs_pdf=True)
# -*- coding: utf8 -*- import pylab as P # Definition de la region de travail Bretagne zone = (-6.5, 47.2, -2, 49.05) # Recuperation du trait Europeen from vacumm.bathy.shorelines import GSHHS tc = GSHHS(input='h', clip=zone) # Interpolation du niveau de la mer sur le trait xyz = tc.bathy() # Plot xyz.plot(size=10, savefigs=__file__, show=False)#, savefigs_pdf=True) P.close()
#!/usr/bin/env python
# -*- coding: utf8 -*-
"""Manipulation de traits de cote (:mod:`vacumm.bathy.shorelines`)"""
from vcmq import N, os, merc, create_grid, resol, map2, P
from vacumm.bathy.shorelines import GSHHS, Histolitt
# Lecture de divers traits sur zone reduite
zone = (-6.5, 47.2, -2, 49.05)
gf = GSHHS('h', clip=zone)
gl = GSHHS('l', clip=zone) # -> essayer Histolitt
# Trace
kwpt = dict(fill=False, points=True, s=2., alpha=.7, points_linewidth=0)
gf.plot(color='r', zorder=12, m_left=.1, label='High resolution', show=False, m_figsize=(5.5, 6), **kwpt)
gl.plot(color='g', zorder=11, label='Low resolution', **kwpt)
# Infos
print gf.resol()
print gl.resol()
print gf.xmax
xy = gl.xy # -> essayer get_point() et get_xy()
print xy.shape
# Bathy cotiere
xyz = gf.bathy()
xyz.plot(size=10)
