2014-11-21 38 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

作爲替代方案,您也可以嘗試將曲目分成兩個不同的半部分,然後分別對它們進行插值。