(lat,lon)之后的成像lat,lon,alt数据 - > (x,y) - > (lat,lon)转换

时间:2016-03-30 22:35:42

标签: python scipy coordinate-transformation linear-interpolation

我有(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提供更规则的网格。这样做有什么好办法?

0 个答案:

没有答案
相关问题