这里是我生成的一些随机数据:
这给我这个图形: 然后我建立了一个等高线图:
import scipy.signal as sgn
import scipy.interpolate as intr
import numpy.ma as ma
x = np.linspace(-100,0,500)
y = sgn.sawtooth(2 * np.pi * .2 * x)
y = (sgn.sawtooth(2 * np.pi * .2 * x)+1)/2
y = (y+1)*25
plt.plot(x,y)
z = np.sin(2*np.pi*.1*x)+np.sin(2*np.pi*.1*y)
这给我这个图形: 然后我建立了一个等高线图:
xi,yi = np.meshgrid(np.linspace(x.min(),x.max(),200),np.linspace(y.min(),y.max(),200))
zi = intr.griddata((x,y), z, (xi, yi) , method='cubic')
plt.contourf( xi,yi,zi,100); plt.colorbar()
对于一个方形的遮罩,我这样做:xi,yi = np.meshgrid(np.linspace(x.min(),x.max(),200),np.linspace(y.min(),y.max(),200))
mask =(yi> 25) & (yi< 35) & (xi > -55) & (xi < -25)
zi = intr.griddata((x,y), z, (xi, yi) , method='cubic')
zi = ma.masked_array(zi, mask = mask )
plt.contourf( xi,yi,zi,100); plt.colorbar()
我的问题是,如何在一个正方形的同一基准位置上放置一个三角形遮罩,给定三角形的三个坐标,即(-50,0),(-25,0),(-37,25)。
mask =(yi> 0) & (yi< 25) & (xi > -53) & (xi < -25) & ((xi+53) < yi)
这样的东西? - Danmask =(yi> 25) & (yi< 35) & (xi > -55) & (xi < -25) & ((xi+95) > 2*yi)
,大致类似于此。 - Dan