Trasformazione coordinate Python da ECI a ECEF

Question

Sto usando il pacchetto python sgp4 (sviluppato da Brandon Rhodes); Sto inserendo le informazioni TLE per ISS (Zarya). Ricevo la posizione del satellite ad una certa data, nel sistema di coordinate ECI. Tuttavia, voglio convertirlo in sistema di coordinate ECEF.

Esistono almeno due istanze di funzioni MATLAB implementate (ecitoecef, eci2ecef) e una funzione C (eci2ecef) implementata per il progetto RTKLIB. (Non posso pubblicare collegamenti perché ho una restrizione di due collegamenti per domanda dal questa è la mia prima domanda e ho una reputazione 1)

Non sono riuscito a trovare una funzione python implementata. Tuttavia, in questa domanda precedente, Python satellite tracking with spg4, pyephem- positions not matching l'utente dice "Ho trovato almeno una parte del problema. Spg84.propagate() restituisce la posizione in ECI, non ECEF. Esegui rapidamente eci2ecef e si allinea perfettamente con la risposta di previsione ".

È che eci2ecef si riferisce a un'implementazione Python? È disponibile online? C'è un documento online che spiega esattamente il passo da compiere se dovessi essere implementato da me? In questa domanda How do I convert ECI coordinates to longitude latitude and altitude to display on a Map? in un commento si dice "Si differenziano per la velocità di rotazione della Terra in modo che la trasformazione sia solo un sinonimo 2-D: -sin trasform di tipo cos". Qual è l'angolo della trasformazione? È il tempo di Greenwhich Sideral? Quando usiamo SGP4, è anche la GST che specifichiamo, o altre volte?

Grazie mille per il vostro aiuto!

Answer 1

Per prima cosa ho studiato il documento AIAA 2006-6753 che è on-line su Celestrak, e che ha costituito la base per rendere la libreria sgp4 disponibili in Python, in primo luogo:

http://celestrak.com/publications/AIAA/2006-6753/

Sulla base delle sue raccomandazioni - e in particolare sulle informazioni nella sua Appendice - Sto tentando di integrare i satelliti della Terra nella mia nuova libreria di astronomia Skyfield. Una nuova versione con questo lavoro dovrebbe essere presto fuori, ma nel frattempo è possibile studiare le trasformazioni che sto usando (in base alla carta oltre) a:

https://github.com/brandon-rhodes/python-skyfield/blob/master/skyfield/sgp4lib.py

Non sono sicuro che ho la Trasforma correttamente, ma controlleremo i risultati rispetto a tutte le fonti ufficiali di posizioni che posso trovare per determinare se ho commesso errori sistematici.

Answer 2

La trasformazione da ECI a ECEF è anche implementata in satellite-js derivata da PyEphem.

Potrebbe essere necessario verificarne l'esattezza.