2014-11-21 4 views
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

1

다른 방법으로 트랙을 두 개의 별개의 반으로 분할하고 개별적으로 보간 할 수도 있습니다.