python - 卫星地面轨迹的 Python 插值

Question

我正在尝试插入以下卫星地面轨道：

问题是 (lon,lat)->(0,0) 处的不连续性，这会导致该区域的拟合曲线不佳：

我不确定参数插值是否是最好的方法，现在我刚刚应用了线性插值：

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from scipy import interpolate

df = pd.read_table('groundtracks', sep='\t')

x = df['lon']
y = df['lat']

f = interpolate.interp1d(x, y, kind='linear')

xnew = np.linspace(x.min(), x.max(), num=x.count()*2)
ynew = f(xnew)

Answer 1

您不应该在 0 附近进行插值，因为卫星在两个非常不同的时间点存在。

您应该将 360° 添加到数据的后半部分（可能对应于左上角曲线），以便获得从大约 -30° 到 360° 的连续曲线。然后插值。然后执行相反的操作：lon[lon>180] -= 360。

Answer 2

作为替代方案，您也可以尝试将轨道分成两个不同的部分，并分别对它们进行插值。