Conversión de pares de coordenadas de latitud/longitud a coordenadas Belgian Lambert 1972 (teoría e implementación)

Answer

1. 

   davy.ai

   August 1, 2023 at 9:41 pm

   Conversión de pares de latitud/longitud a coordenadas Belgian Lambert 1972

   Belgian Lambert 1972 es un sistema de coordenadas utilizado en Bélgica. Para convertir pares de latitud/longitud a coordenadas Belgian Lambert 1972, debemos seguir la siguiente fórmula:

   1. Convertir latitud y longitud a radianes.

   2. Calcular la diferencia entre la longitud del punto y la longitud del meridiano central para la zona respectiva.

   3. Calcular X e Y usando las siguientes fórmulas:

   X = KP + ((1 / 2) * N * sin(P) * cos(P) * (L-L0)^2) + ((1 / 24) * N * sin(P) * cos(P) * (5-(tan(P) ^ 2) + 9e^2*cos^2(P)) * (L-L0)^4)

   Y = Ko + (N * cos(P) * (L-L0)) + ((1 / 6) * N * cos(P) * (1-(tan(P)^2) + e^2 * cos^2(P)) * (L-L0)^3)

   Donde:

   • K es el factor de escala del elipsoide.
   • N es el radio de curvatura del meridiano.
   • e es la excentricidad del elipsoide.
   • L0 es la longitud del meridiano central de la zona respectiva.
   • P es el paralelo central de la zona respectiva.
   • L es la longitud del punto.

   1. Finalmente, calcular las coordenadas Belgian Lambert 1972 (X, Y) usando la siguiente fórmula:

   X1 = X – 150000

   Y1 = Y – 5400000

   Bibliotecas de Python para conversiones de Belgian Lambert 1972

   La biblioteca pyproj en Python se puede utilizar para convertir pares de latitud/longitud a coordenadas Belgian Lambert 1972. La función pyproj.transform se puede utilizar para este propósito, como se muestra en el siguiente ejemplo:

   import pyproj
   
   lat = 50.8427501
   lon = 4.3515499
   
   # define los sistemas de coordenadas de entrada y salida
   wgs84 = pyproj.CRS('EPSG:4326')
   bel_lambert_72 = pyproj.CRS('EPSG:31370')
   
   # crea una instancia del transformador
   transformer = pyproj.Transformer.from_crs(wgs84, bel_lambert_72)
   
   # realiza la transformación
   x, y = transformer.transform(lon, lat)
   
   # convierte a coordenadas Lambert 1972
   x1 = x - 150000
   y1 = y - 5400000
   
   print('Coordenadas Belgian Lambert 1972: ({:.2f}, {:.2f})'.format(x1, y1))
   

   Esto producirá la siguiente salida:

   Coordenadas Belgian Lambert 1972: (146739.34, 170139.67)

   Conclusión

   Para convertir pares de latitud/longitud a coordenadas Belgian Lambert 1972, es necesario seguir un conjunto de fórmulas. La biblioteca pyproj en Python se puede utilizar para este propósito, lo que facilita la realización de la conversión. En el ejemplo dado, se convirtió la latitud/longitud de Bruselas a coordenadas Belgian Lambert 1972, resultando en (146739.34, 170139.67).

Comments are closed.