Le portail francophone de la géomatique

Bonjour et bienvenu

Théoriquement, vous devriez trouver votre bonheur là :
http://geodesie.ign.fr/contenu/fichiers … lg0073.pdf

D'une manière générale, la plupart des algorithmes concernant la géodésie Française se trouve sur le site geodesie.ign.fr

Cordialement

Merci François pour ta réponse, cela correspond parfaitement à ce que je cherchais !

Pour ceux que cela intéresserait, voilà ce que donne l'algorithme décrit, prenant en entrée les coordonnées ETRS-LAEA (X,Y) et donnant en sortie longitude et latitude (lambda, phi) :

conversion_degre <- function(X,Y){
 
  a <- 6378137.0
  e <- 0.081819191043
  lambda0 <- 0.174532925199
  phi1 <- 0.907571211037
  X0 <- 4321000
  Y0 <- 3210000
  epsilon <- 10^(-11)
 
  x <- X-X0
  y <- Y-Y0
  qp <- (1-e^2)*(1/(1-e^2) - 1/(2*e) * log((1-e)/(1+e)))
 
  q1 <- (1-e^2)*(sin(phi1)/(1-e^2*sin(phi1)^2) - (1/(2*e)) * log((1-e*sin(phi1))/(1+e*sin(phi1))))
  beta1 <- asin(q1/qp)
  m1 <- cos(phi1)/sqrt(1-e^2*sin(phi1)^2)
  Rq <- a*sqrt(qp/2)
  D <- (a*m1)/(Rq*cos(beta1))
  rho <- sqrt((x/D)^2 + (D*y)^2)
  ce <- 2*asin(rho/(2*Rq))
 
  q <- qp*(cos(ce)*sin(beta1) + D*y*sin(ce)*cos(beta1)/rho)
 
  lambda <- lambda0 + atan((x*sin(ce))/((D*rho*cos(beta1)*cos(ce)) - D^2*y*sin(beta1)*sin(ce)))
  phi0 <- asin(q/2)
 
  phi_iplus1 <- Inf
  phi_i <- phi0
  ecart <- Inf
  while(ecart>epsilon){
    phi_iplus1 <- phi_i + ((1-e^2*sin(phi_i)^2)^2/(2*cos(phi_i)))*( q/(1-e^2)-sin(phi1)/(1-e^2*sin(phi_i)^2)+(1/(2*e))*log((1-e*sin(phi_i)) / (1+e*sin(phi_i))))
    ecart <- abs(phi_iplus1-phi_i)
    phi_i <- phi_iplus1
  }
  phi <- phi_i
 
  return(c(lambda*(180/pi),phi*(180/pi)))
}