2014-11-21 54 views
0

我想要内插以下卫星地面轨迹:的Python插值跟踪

enter image description here

的问题是在(LON,LAT)的不连续性 - >(0,0),其导致该地区一个贫穷的拟合曲线:

enter image description here

我不知道,如果一个参数插值将是最好的办法,现在我刚申请的线性插值:

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) 

回答

1

您不应该在0附近插值,因为卫星在两个不同的时间点存在。

您应该在数据的后半部分添加360°(可能对应于左上角的曲线),以便您可以获得从-30°到360°的连续曲线。然后插入。然后执行相反的操作:lon[lon>180] -= 360

+0

谢谢,这是一个很好的解决方法,但是数据需要事先进行一些操作,以按照您的建议按象限分割。 – npires 2014-11-21 16:43:55

1

作为替代方案,您也可以尝试将曲目分成两个不同的半部分,然后分别对它们进行插值。