French RGF93 - Lambert 93 will be available in gvSIG Mobile

Thanks to Hervé Quinquenel (from the French Ecole Nationale des Sciences Géographiques, ENSG, http://www.ensg.eu and the Institut Géographique National, IGN, http://www.ign.fr) we now have the official formulas for the conversion between WGS84 coordinates (EPSG:4326, used by the GNSS systems such as GPS) and the new official French projection (RGF93 - Lambert 93, EPSG:2154). You can see the formulas in the source code SVN repository, searching for the tag "AJOUT HERVE" in this file. Of course this conversion will also be available in the upcoming gvSIG Mobile 1.0.



To check the correctness of this new Java code, we can use the coordinates of the Point zéro des routes de France (France road system zero point):

EPSG:2154
X = 652215.52 m
Y = 6861681.77 m

EPSG:4326
48º 51' 12.24845'' N (48.8534023472)
2º 20' 55.62563'' E (2.3487848972)

If we use those EPSG:4326 coordinates as the input for the formulas sent by Hervé, we get:

X = 652215.5199688374
Y = 6861681.770095885

The divergence is less than one milimeter, which means that we are actually using the official formulas. For those who are a bit skeptical, the coordinates of Paris are not part of the conversion formula, so in terms of testing, the coordinates we have used are a random sample.



Point zéro des routes de France (Source: Wikipedia).

Comments

Post a Comment

Popular posts from this blog

Computing and applying an affine transformation to a vector layer using gvSIG and Quantum GIS

Image
Geomatics Supervisor and Egyptian Archaeology PhD student Anna Hodgkinson told me that she is trying to apply an affine transformation to a set of coordinates. These are originally in an unknown local coordinate system. She has provided 4 reference points for which both the local coordinates and the EPSG:32636 coordinates are known: ID ; EPSG_32636_X  ; EPSG_32636_Y   ; Local_X     ; Local_Y 1  ; 300646.885431 ; 3231905.223285 ; -144.333424 ; -50.678221 4  ; 300939.099675 ; 3232334.126771 ;  100.477740 ; 407.439814 7  ; 300519.458438 ; 3232235.904115 ; -304.216081 ; 262.236955 8  ; 300717.256565 ; 3232552.638076 ; -141.107317 ; 602.476276 First we need to compute which affine transformation takes the coordinates from the local grid to EPSG:32636. To achieve this, we can use the raster georeferencing tool in gvSIG since it lets you input by hand the source and target coordinates and then creates a WLD file whose 6 parameter...
Read more

Transformación de CANOA (Venezuela) a WGS84 en gvSIG mediante una rejilla que simula el método Molodensky-Badekas

Image
Puesto que gvSIG 1.9 no permite utilizar la transformación Molodensky-Badekas de 10 parámetros , una alternativa consiste en crear una rejilla o grilla de transformación en formato NTv2 utilizando los resultados de la transformación en un gran número de puntos de control repartidos por todo el territorio (en este caso, Venezuela). Hermógenes Suárez y José Nicolás Pérez Duin han publicado en la lista de usuarios de gvSIG una hoja de cálculo con las fórmulas y los parámetros necesarios para realizar la transformación. Para realizar la rejilla de transformación, partimos del rectángulo mínimo que contiene a todo el territorio de Venezuela: Longitud: [-73.5, -59.75] Latitud: [0.5, 12.25] La distancia entre los nodos de la rejilla debe, por un lado, ser lo suficientemente pequeña como para proporcionar una precisión adecuada, pero por otro lado no debemos generar una rejilla con un número desorbitado de nodos. En el primer intento, generaremos una rejilla con una distancia entre nodos de ...
Read more