Изменение координат вектора магнитного поля со специей ⇐ Python

[Anonymous]

У меня есть данные о магнитном поле космического корабля "Вояджер-1" в координатах IHG (инерциальная гелиографическая, IAU_SUN в спайсе). Я хочу преобразовать эти данные в координаты в системе S3 (координаты, связанные с Юпитером, IAU_JUPITER в spice).
Поэтому я хочу перейти от (Bx, By, Bz) в IHG к (Br, Btheta, Bphi) в S3.
Поэтому я написал функцию, которая использует функцию pxform для вычисления матрицы вращения между двумя интересующими меня кадрами. Затем он применяет эту функцию к моим известным данным через произведение матрицы на вектор, а затем меняет координаты с декартовых на сферические.
Я пробовал эту функцию на известных точках в обеих системах, и это не очень хорошо, но мне кажется, что это не так. чтобы найти проблему.
def transform_vector(dataframe, ET_1979, index):

# Get the rotation matrix from IHG to S3
transform_matrix = spice.pxform("IAU_SUN", "IAU_JUPITER", ET_1979[index])

# Magnetic field vector in Heliographic frame
B_IHG = np.array([dataframe['Bx'][index], dataframe['By'][index], dataframe['Bz'][index]])

# Transform to System III
B_tamp = spice.mxv(transform_matrix, B_IHG)

# Give the coordinates Br, Btheta, Bphi
B_S3 = spice.recsph(B_tamp)

# Return transformed Cartesian and spherical components as a pandas Series
return pd.Series({
'days' : VG_1979['days_1979'][index],
'et' : ET_1979[index],
'Br': B_S3[0],
'Btheta': B_S3[1],
'Bphi': B_S3[2],
})


Подробнее здесь: https://stackoverflow.com/questions/792 ... with-spice