map2(mvmask, fillcontinents=False, drawcoastlines_color='r', drawcoastlines_zorder=100,
res=resolution, proj='merc', fill='pcolor', drawcoastlines_linewidth=2,
cmap=cmap_rs(['1', '0.6']), contour=False, colorbar=False,
savefig ='courses_masking_0.png', show=True)
# Cote via mask
indices = get_coastal_indices(mask)
dist = get_dist_to_coast(grid, mask)
# Erosion de la cote
f = cdms2.open(data_sample('mars3d.xy.nc'))
temp = f('temp')
f.close()
temp2 = erode_coast(temp, maxiter=3) # changez maxiter et tracez
# Enveloppe de points
xy = P.randn(2, 500)*0.5+N.array([-5,48]).reshape(2,-1)
xe, ye = convex_hull(xy, method='delaunay') # changez la methode
# Creer des plolygones
xpoly = [-5.8, -5, -5, -5.8]
ypoly = [48, 48, 48.5, 48.5]
polys = polygons([[xpoly, ypoly]]) # meme chose avec des min/max + tracez
# -> appliquez a polygon_mask
# Selection de points
drawcoastlines_linewidth=2,
cmap=cmap_rs(['1', '0.6']),
contour=False,
colorbar=False,
savefig='courses_masking_0.png',
show=True)
# Cote via mask
indices = get_coastal_indices(mask)
dist = get_dist_to_coast(grid, mask)
# Erosion de la cote
f = cdms2.open(data_sample('mars3d.xy.nc'))
temp = f('temp')
f.close()
temp2 = erode_coast(temp, maxiter=3) # changez maxiter et tracez
# Enveloppe de points
xy = P.randn(2, 500) * 0.5 + N.array([-5, 48]).reshape(2, -1)
xe, ye = convex_hull(xy, method='delaunay') # changez la methode
# Creer des plolygones
xpoly = [-5.8, -5, -5, -5.8]
ypoly = [48, 48, 48.5, 48.5]
polys = polygons([[xpoly, ypoly]]) # meme chose avec des min/max + tracez
# -> appliquez a polygon_mask
# Selection de points
xsel, ysel = polygon_select(xy[0], xy[1], polys) # essay l'option mask
P.plot(xy[0], xy[1], 'o')
P.plot(xpoly, ypoly, 'k-', lw=2)
xx, yy = N.indices((50, 100), 'f')
x0 = y0 = 25
dxy = 30
var = MV2.exp(-((xx - x0)**2 + (yy - y0)**2) / dxy**2)
# Masques
# - reference
mask = var.filled() > .9
mask[:, 50:] = True
mask[15:35, 65:85] = False
# - variable
var[:] = MV2.masked_greater(var, .8)
var[:, 50:] = MV2.masked
# Erode
vare = erode_coast(var, mask)
# Plots
P.figure(figsize=(6, 9))
P.subplots_adjust(hspace=.2)
P.subplot(311)
P.pcolormesh(var.asma())
P.title('Original variable')
add_key(1, color='w')
P.subplot(312)
P.pcolormesh(mask)
P.title('Reference mask')
add_key(2, color='w')
P.subplot(313)
P.pcolormesh(vare.asma())
P.title('With coastal erosion')
xx, yy = N.indices((50, 100), 'f')
x0 = y0 = 25
dxy = 30
var = MV2.exp(-((xx-x0)**2+(yy-y0)**2)/dxy**2)
# Masques
# - reference
mask = var.filled()>.9
mask[:, 50:] = True
mask[15:35, 65:85] = False
# - variable
var[:] = MV2.masked_greater(var, .8)
var[:, 50:] = MV2.masked
# Erode
vare = erode_coast(var, mask)
# Plots
P.figure(figsize=(6, 9))
P.subplots_adjust(hspace=.2)
P.subplot(311)
P.pcolormesh(var.asma())
P.title('Original variable') ; add_key(1, color='w')
P.subplot(312)
P.pcolormesh(mask)
P.title('Reference mask') ; add_key(2, color='w')
P.subplot(313)
P.pcolormesh(vare.asma())
P.title('With coastal erosion') ; add_key(3, color='w')
savefigs(__file__)
