我有(lonvec,latvec,altvec)数据,是从基于equirectangular投影的数据集中提取的。我正在将坐标转换为(xvec,yvec,altvec)。我使用此数据进行插值:
X = np.linspace(xvec.min(), xvec.max(), 300)
Y = np.linspace(yvec.min(), yvec.max(), 300)
Z = np.linspace(zvec.min(), zvec.max(), 300)
X2, Y2 = np.meshgrid(X, Y)
interp = scipy.interpolate.LinearNDInterpolator(cartcoord, mesvec0, fill_value=0)
Z0 = interp(X2, Y2)
此时,我有一个(xvec,yvec)点的网格网格和插值数据集(Z0)。这些变量(X2,Y2和Z0)中的每一个都是二维的。
我想保留插值点(及其对应的x,y位置),同时在显示最终图像之前将x,y坐标转换回lat,lon坐标。最终图像应为(lat,lon,alt)。从(lat,lon)到(x,y)进行变换的原因是为LinearNDInterpolator提供更规则的网格。这样做有什么好办法?