miércoles, 1 de junio de 2011

Diferencias entre Psad56 Mean y Psad56 Perú, Bolivia, Guyana, etc...

Quizá en nuestro trabajo con mapas, planos, información georeferenciada, tenemos que ver en que proyección y Datum se encuentran, y notamos que están en una proyección por ejemplo:

UTM, Datum PSAD 56, Zona 18S

…Pero nos podemos topar que nos digan más información como:

UTM, Datum PSAD 56 Mean, Zona 18S (que por lo general en mi país, Perú, es el más común en usarlo).

Y esta otra…

UTM, Datum PSAD 56 Perú, Zona 18S, UTM, Datum PSAD 56 Bolivia, Zona 18S, UTM, Datum PSAD 56 Guyana, Zona 18S, entre otras...

Entonces estamos en la interrogante… ¿En qué se diferencian? En este momento tratare de explicar esto.

Quizá una forma para diferenciarlos es usando como base al Datum WGS 84: ya que los “Diferenciales X Y Z” o” Shift X Y Z” de estos Datums PSAD 56 se diferencian entre ellos como vemos a continuación:


Diferenciales o Shift para WGS 84

Psad 56

X

Y

Z

Mean

-288

175

-376

Bolivia

-270

188

-388

Chile N

-270

183

-390

Chile S

-305

243

-442

Colombia

-282

169

-371

Ecuador

-278

171

-367

Guyana

-298

159

-369

Perú

-279

175

-379

Ahora veremos cómo pasar de un Datum PSAD 56 Mean a PSAD 56 Perú mediante el Datum WGS 84:

En este ejemplo usaremos el ArcGIS para hacer este cambio:

1) Tenemos un cuadrado proyectado llamado “PSAD56_MEAN” en PSAD_1956_UTM_18S (En el ArcGIS no hay un Datum que diga PSAD_1956_UTM_18S MEAN o PSAD_1956_UTM_18S PERU) y estamos asumiendo primeramente que sea PSAD 56 MEAN con las siguientes coordenadas:

a. 210000, 8910000

b. 210000, 8920000

c. 220000, 8920000

d. 220000, 8910000

2) Pasamos nuestro cuadrado de PSAD 56 MEAN a WGS 84 (recuérdese que estamos asumiendo que nuestro Datum primeramente es PSAD 56 MEAN) y acá comienza nuestro paso de PSAD56 MEAN a PSAD 56 PERU, en la parte en rojo de la imagen de abajo nos piden que opcionalmente (aunque la verdad para mi es algo esencial en este tipo de trabajos) el tipo de transformación geográfica (o Geographic Transformation) y seleccionamos el primero de la lista desglosable que dice “PSAD_1956_To_WGS_1984_1” que contiene los Diferenciales o Shifts para pasar de PSAD 56 MEAN a WGS 84 y viceversa (con el botón con el símbolo + lo añadimos en el recuadro debajo de la lista desglosable) y pulsamos el botón OK.

3) Y tenemos nuestro cuadrado en WGS 84 llamado “WGS84” desde el PSAD 56 MEAN.

4) Ahora a nuestro cuadrado “WGS84” lo transformamos a PSAD56 PERU y lo llamaremos “PSAD56_PERU” y repetimos lo que hicimos en el punto 2 con la diferencia que en la lista desglosable de Geographic Transformation en la lista desglosable seleccionamos “PSAD_1956_To_WGS_1984_8” que contiene los Diferenciales o Shifts para pasar de WGS 84 a PSAD 56 PERU y viceversa (con el botón con el símbolo + lo añadimos en el recuadro debajo de la lista desglosable) y pulsamos el botón OK.


5) Y tenemos nuestro cuadrado llamado PSAD56_PERU en el Datum PSAD 56 PERU.

6) Para verificar que todo lo dicho anteriormente se cumple, usaremos el “Geotrans V2.4” para demostrar que estamos bien con los Datums de nuestros cuadrados


7) Usaremos los vértices Superiores Izquierdos de nuestros cuadrados proyectados en PSAD 56 MEAN, WGS 84 y PSAD 56 PERU y rotulamos con sus coordenadas.

8) Solo nos enfocaremos entre los Datums PSAD 56 MEAN y PSAD 56 PERU ya que el “Geotrans V2.4” contiene directamente esta transformación. Seleccionamos en la primera parte superior del “Geotrans V2.4” las coordenadas en PSAD 56 MEAN y lo configuramos con los parámetros de dicho Datum.

9) Y en la parte inferior del “Geotrans V2.4” configuramos para que sean proyectados en PSAD 56 PERU.

10) Presionamos el botón “Convert : Upper - > Lower” y tenemos las coordenadas proyectadas de PSAD56 MEAN a PSAD56 PERU y vemos que son las mismas coordenadas que están rotuladas.

Para hacerlo en Mapinfo, tenemos que editar el archivo Mapinfo.prj, lo que haremos es introducir el siguiente código (en la imagen siguiente esta al costado de la flecha en rojo):

"PSAD56 / UTM zona 18S Peru", 8, 999, 4, -279, 175, -379, 7, -75, 0, 0.9996, 500000, 10000000

Donde:

8: Proyección Transversal Mercator
999: Formato para la adición de una nueva proyección con un Datum definido por el usuario
4: Elipsoide Internacional 1924
Enlace-279,175-379: shifts WGS 84 (ver cuadro más arriba)
7: metros
-75: Meridiano central zona 18S
0: Origen Latitud
0.9996: Factor de escala
500000: Falso este
10000000: Falso norte

Este tema se explico en la entrada ¿Cómo crear nuevas proyecciones en el Mapinfo?

Y la comprobación será así (es en algo parecido al anterior):

11) Creamos un nuevo archivo y creamos un recuadro con las siguientes coordenadas:

a. 210000, 8910000

b. 210000, 8920000

c. 220000, 8920000

d. 220000, 8910000

Luego lo proyectamos “Psad 56 / 18 s” (en el Mapinfo no indica que sea Mean, Perú u otro Datum del tipo Psad 56), así que asumiremos que sea el Psad 56 Mean y lo llamaremos “Psad56_Mean_18s”.


12) Luego lo guardamos como un archivo de Mapinfo con nombre “Psad56_Peru_18s” y lo proyectamos con la proyección que recién acabamos de crear "PSAD56 / UTM zona 18S Peru"

13) Abrimos este nuevo archivo “Psad56_Peru_18s” en una nueva ventana y comprobamos sus coordenadas para ver si la transformación es la correcta. Para ello usaremos de nuevo el "Geotrans V2.4." Y observamos que la transformación es correcta.

Esperando que sea útil esta información, si tuvieran alguna consulta o comentario, lo pueden hacer al correo recursos.web.gis@gmail.com ó a info@ingeografos.com.pe. Gracias y saludos.

ª