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 :

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.