Conversion coordonnées gps
Bonjour
Je dispose d'un tableau numpy de dimension 17000*3 formaté comme suivant :
longitude ; latitude ; année de saisie;
et je souhaiterais étudier l'évolution des saisies selon les années.
Pour cela j'aurai besoin de convertir mes coordonnées longitude, latitude en données de format (x,y) sur un plan (sans afficher nécessairement ce plan). Et je me fiche de connaître l'origine des abscisses et des ordonnées, puisque ce qui va m'intéresser c'est l'éloignement des données entre elles et la façon dont la surface la plus petite contenant toutes les saisies s'étend chaque année.
Pour cela j'ai essayé d'utiliser le programme suivant :
Quelqu'un aurait il une piste s'il vous plaît ?
Merci d'avance.
Je dispose d'un tableau numpy de dimension 17000*3 formaté comme suivant :
longitude ; latitude ; année de saisie;
et je souhaiterais étudier l'évolution des saisies selon les années.
Pour cela j'aurai besoin de convertir mes coordonnées longitude, latitude en données de format (x,y) sur un plan (sans afficher nécessairement ce plan). Et je me fiche de connaître l'origine des abscisses et des ordonnées, puisque ce qui va m'intéresser c'est l'éloignement des données entre elles et la façon dont la surface la plus petite contenant toutes les saisies s'étend chaque année.
Pour cela j'ai essayé d'utiliser le programme suivant :
def lat_long_to_cartesiennes(lat,longi): R=6371 #globe radius (km) x=R*longi print(np.tan(((np.pi/4)+(lat/2)))) y=R*np.log(np.tan(((np.pi/4)+(lat/2)))) return(x,y)Mais cela ne fonctionne pas car je me retrouve avec des valeurs négatives dans le log par moment, et quand je n'en ai pas, les résultats me semblent aberrants (j'ai des distances énormes dès la première année alors que ma surface de répartition est très faible).
Quelqu'un aurait il une piste s'il vous plaît ?
Merci d'avance.
Connectez-vous ou Inscrivez-vous pour répondre.
Réponses
$x = R\cos(lat).\cos(longi)$
$y = R\cos(lat).\sin(longi) $
$z= R\sin(lat) $
https://forums.futura-sciences.com/mathematiques-superieur/905320-coordonnees-gps-mille-nautique-cercle-probleme-reel-urgent.html
en particulier les posts de Chadok
je n'ai pas trop compris tes formules, Chadok en donne qui fonctionnent
peut-être vérifier les unités, la longitude et la latitude sont en degrés, en radians dans ton programme ?
tes formules tiennent elles compte du fait que la circonférence d'un parallèle se rétrécit plus on se rapproche d'un pôle