Evaluación del error asociado a la transformación de datum del
Sistema Geodésico Mundial (WGS84) a Ocotepeque-CR
utilizando el método de tres parámetros de Molodensky
Jorge Fallas
Laboratorio de Teledetección y Sistemas de Información Geográfica
PRMVS-EDECA, UNIVERSIDAD NACIONAL
2003
i
INDICE
Introducción ...................................................................................................................... 1
Metodología ...................................................................................................................... 2
Resultados y discusión....................................................................................................... 5
Error de los parámetros disponibles para realizar la transformación de datum .................... 6
Ajustes empíricos .............................................................................................................. 9
Traslación de ejes de Rody Rodríguez ............................................................................... 9
Ajuste ecuación lineal...................................................................................................... 10
Error en Osa y Manzanillo ............................................................................................... 12
La extensión CR_proyecciones_datum ............................................................................ 13
Error esperado al proyectar datos de CRTM a Lambert y viceversa.................................. 13
Conclusiones ................................................................................................................... 14
Referencias...................................................................................................................... 14
1
Introducción
El Sistema de Posicionamiento Global es un sistema de navegación y posicionamiento satelital que
ha estado en uso en el área militar desde 1978 (Bennett, 1990; [Link]/news/;
[Link]/news/); sin embargo su uso se popularizó en el sector civil hasta el año 1995; cuando el
sistema fue declarado operacional por el Departamento de Defensa de los Estados Unidos de
América (Enge y Misra, 1999). En la actualidad los receptores del SPG se utilizan en áreas tan
diversas como la navegación terrestre, área y marítima; el levantamiento de fincas; en la
actualización de cartografía temática; en la agricultura de precisión; en el estudio de deformaciones
de estructuras tales como puntes, represas y placas tectónicas y en la navegación espacial (Bolstad,
1993, Enge y Misra, 1999).
El usuario del SPG tiene acceso a 6 o más satélites en un 96% del tiempo, a 8 o más satélites en un
32% del tiempo y a 9 o más satélites un 5% del tiempo. Si el receptor tiene capacidad de recibir
señales del SPG y del sistema homólogo de la Federación Rusa denominado GLONASS (GLObal
Navigation Satellite Systems) su receptor tendrá acceso a 9 o más satélites un 99% del tiempo (Enge
y Misra, 1999). La reducción en el costo de los receptores ha popularizado su uso en el área de
medio ambiente y recursos naturales. Sin embargo no siempre el uso de una tecnología moderna
implica una mayor exactitud en los datos que se registran en el campo. La mayoría de los receptores
de bajo costo (US$100 a 300) tienen un error en posicionamiento de 10 a 15 m. Los nuevos
receptores con capacidad de recibir señales del WAAS (Wide Area Augmentation System) pueden
reducir dicho error a 5m. Sin embargo, si no se configura y utiliza el receptor en forma correcta la
posición registrada por el usuario (a) puede tener un error de varios cientos de metros.
Uno de los aspectos fundamentales al utilizar un receptor del SPG es la selección del datum; el
cual es un conjunto de parámetros que definen un sistema de coordenadas; así como un conjunto de
puntos de control cuyas relaciones geométricas son conocidas; ya sea mediante mediciones o
cálculos matemáticos (Barrantes, 1966; ESRI, 1994; Inter-American Geodetic Survey, 1950). Dado
que la superficie de la Tierra es un esferoide y no una esfera, la exactitud de un datum local está
restringida a un área particular alrededor de su origen. La cartografía nacional creada por el Instituto
Geográfico Nacional utiliza el datum Ocotepeque; el cual está ubicado en el departamento del
mismo nombre en Honduras, a una altura media sobre el nivel del mar es de 807 metros y fue
establecido por V. J. Hanraham en diciembre de 1934. El elipsoide de referencia es el Clarke de
1866 (Inter-American Geodetic Survy, 1950).
Cuando se utiliza un geoposicionador, es frecuente que las coordenadas se registren utilizando el
datum WGS84; el cual utiliza el centro de masa de la Tierra como referencia (o sea es geocéntrico).
Aun cuando algunos receptores manuales del SPG permiten seleccionar y/o configurarlo para utilizar
un datum local (Ej. Ocotepeque), el usuario(a) debe conocer los parámetros de transformación entre
el WGS84 (World Geodetic Datum de 1984) y el datum local para poder realizar esta tarea.
La transformación más simple es la de Molodensky, la cual utiliza tres parámetros de traslación:
?X, ?Y y ?Z. Otra transformación utilizada con frecuencia es la de Bursa-Wolfe, la cual considera
siete parámetros: esta transformación incluye además de las traslaciones de ejes ya citadas, tres
rotaciones, una en X (Rx), otra en Y (Ry) y otra en Z (Rz); así como un factor de escala (?e); los
primeros tres parámetros de miden en metros, los segundos tres en segundos decimales y el tercero
en partes por millón (ppm) o sea en milímetros por kilómetro (Dörries y Roldán, 1999; ESRI, 1994).
2
Al utilizar la cartografía nacional y los datos obtenidas con un geoposicionador satelital, el
usuario(a) debe ser conciente de que la concordancia espacial entre ambas fuentes de datos
dependerá en gran parte de la selección correcta del datum. El efecto de no utilizar el mismo datum
será un desplazamiento (no coincidencia espacial) entre las datos del mapa y los del
geoposicionador. Por ejemplo, si usted registra la posición de una esquina en latitud y longitud en
Costa Rica utilizando el datum WGS84 y luego lo proyecta a Lambert Norte sin realizar la
transformación de datum, su objetó se desplazado 264.77 m hacia el noroeste. Este desplazamiento
sería aceptable para una cartografía a escala 1:600.000; sin embargo es muy superior a la exactitud
nominal de los geoposicionadores actuales (10-15m).
En el país se han realizado algunos trabajos de investigación sobre el tema (Catastro Nacional,
2001; Dörries y Roldán, 1999; Rodríguez, 2001) sin embargo hasta la fecha no existe un set de
parámetros oficiales de Molodensky ni un programa de computo que le permita a los profesionales
en medio ambiente y recursos naturales realizar la transformación del datum WGS84 al Ocotepeque
ó configurar sus georeceptores. Dado este vacío de información sobre el tema se planteó el presente
trabajo, el cual tiene como objetivo evaluar el error esperado al utilizar diferentes parámetros
disponibles para la transformación del datum WGS84 al Ocotepeque y a la vez implementar los
resultados mediante una extensión para ArcView GIS 3.x.
Metodología
El estudio se realizó utilizando 28 puntos de posición conocida tanto en el datum Ocotepeque
como en el WGS84 (Dörries y Roldán, 1999); su distribución a nivel nacional es adecuada con la
excepción del flanco sureste de la Cordillera de Talamanca y la porción extrema de la Zona Norte
(Figura 1).
Figura 1: Ubicación de los puntos de referencia utilizados en la evaluación del error radial asociado a
la transformación del datum WGS84 a Ocotepeque.
3
Los métodos de transformación y los parámetros de Molodensky (?X, ?Y,?Z) evaluados fueron
los siguientes:
1. Datos sin transformar
En este caso los datos en WGS84 fueron proyectados directamente a Lambert Norte. Esto equivale a
utilizar: ?X: 0 ?Y: 0 y ?Z: 0 y muestra el error esperado si los datos se utilizan sin realizar ninguna
transformación.
2.NAD27 América Central
Este set de parámetros corresponde a los valores promedios sugeridos por la Agencia de Imágenes y
Mapas de los Estados Unidos de Norteamérica (NIMA-USA) para la transformación de datum
WGS84 a NAD27 América Central (Belice a Costa Rica) ?X:0 m, ?Y:125 m y ?Z: 194 m. Dichos
valores son también sugeridos en los manuales proveídos por compañías tales como ESRI (1994) y
la mayoría de las compañías que venden geoposicionadores de bajo costo a nivel mundial. En este
caso dichas publicaciones asumen que el datum utilizado en la cartografía de Costa Rica es el
NAD27 América Central, lo cual es incorrecto.
3. Parámetros sugeridos por el señor G. Stentz del NIMA-USA.
Mediante un correo electrónico el señor Stentz, funcionario del NIMA, sugirió los siguientes
parámetros de transformación: ?X: 205 m, ?Y: 96 m y ?Z: -98 m (email: stentzG@[Link]).
Según el señor Stentz el error esperado de la transformación es de 10m.
4. Traslación de ejes sugerida por el Ing. Rody Rodríguez, ICE
A partir de los levantamientos realizado por el Ing. Rody Rodríguez, funcionario del departamento
de Topografía del Instituto Costarricense de Electricidad; el ha propuesto las siguientes traslaciones
de ejes: ?X: -7.072781" y ?Y: 5.03437". En este caso dichos valores se adicionaron a los
respectivos valores de longitud y latitud en WGS84. A partir de los valores del Ing. Rodríguez, el
autor evaluó otras translaciones con el objetivo de determinar si era posible derivar un nuevo set de
valores con un menor error.
5. Ajuste utilizando función de transformación lineal
En este caso se ajustó una ecuación de regresión lineal simple a los valores de latitud y longitud; los
valores en WGS84 se utilizaron como “X” y los valores en Ocotepeque como “Y”.
5. Parámetros derivados por el autor para OSA y Manzanillo
Mediante un proceso de ajuste interactivo, el autor probó los siguientes parámetros:
JF_osa: ?X: 206.75 ?Y: 156.25 y ?Z: -95.25
JF_manzanillo: ?X: 191.5 ?Y: 137 y ?Z: -93.5
Osa y Manzanillo son un caso especial, ya que el error observado en dichos puntos siempre fue
superior al observado en los otros puntos del país; por esta razón se ajustó un set de parámetros
independientes para dichos puntos.
6. Parámetros sugeridos por Dörries y Roldán (1999):
Dörries y Roldán (1999) recomiendan utilizar los siguientes parámetros de traslación ?X: 213.12
(Sx=0.428 m), ?Y: 9.36 (Sy=0.428 m) y ?Z: -74.95 (Sz=0.428 m). En su informe también brindan
los siguientes valores para las rotaciones alrededor de los ejes X, Y y Z: -2.35”, +0.06”y –6.39”;
respectivamente; así como un factor de escala de 5.2 ppp (milímetros por kilómetro). Estos valores
4
conforman la transformación de 7 parámetros recomendada por dichos autores. Estos parámetros de
transformación permiten pasar del datum WGS84 (ITRF98,35; datum propuesto por los autores).
Una vez elegidos los parámetros de transformación se procedió de la siguiente manera:
1. Se creó una tabla que contenía los puntos en latitud-longitud (grados decimales) tanto en
WGS84 como en Ocotepeque; así como los datos proyectados a Lambert Norte y con el
datum Ocotepeque. Estos valores se tomaron del informe de Dörries y Roldán (1999).
2. Se crearon tres temas de eventos; uno en WGS84, otro en Ocoteqeque y otro en Lambert
Norte_ocotepeque.
3. Se realizó la transformación de datum de WGS84 a Ocotepeque utilizando el programa
ArcView GIS ([Link]) y la extensión CR_proyecciones_datum.avx creada por el
autor como parte del presente trabajo.
4. Se evaluó el error para cada uno de los métodos y set de parámetros de Molodensky
utilizados en el estudio. El error se definió como el desplazamiento radial entre la posición
de cada uno de los puntos en Lambert Norte, Ocotepeque y su posición estimada en Lambert
Norte utilizando los diferentes parámetros de transformación de datum (Acharya y Bell,
1992; Giovachino, 1992). La distancia fue estimada utilizando la función de unión espacial
de ArcView GIS. Por ejemplo, un punto con un error de 10m tanto en X como en Y tendrá
un error radial de 14.1m (Figura 2).
Figura 2 : Estimación del error radial
5. El error medio en posicionamiento para cada una de las transformaciones se calculó con la
siguiente ecuación (raíz del error medio cuadrático):
REMC= [ ? (Exy)2 / n ] 0.5; en donde: Exy representa el error total o radial de posicionamiento y
“n” el número de puntos evaluados.
5
Resultados y discusión
El datum es un conjunto de parámetros que definen un sistema de coordenadas; así como un
conjunto de puntos de control cuyas relaciones geométricas son conocidas; ya sea mediante
mediciones o cálculos matemáticos (Barrantes, 1966; ESRI, 1994; Inter-American Geodetic Survey,
1950). Para cada datum existe un elipsoide de referencia, una posición y orientación del acimut
inicial (lat, long) así como una distancia entre el geoide y el elipsoide en el punto de ubicación
inicial. Dado que la superficie de la Tierra es un esferoide y no una esfera, la exactitud de un datum
local está restringida a un área relativamente pequeña alrededor de su origen (Fig.3).
El datum utilizado en la cartografía oficial de Costa Rica se denomina Ocotepeque y está ubicado
en el departamento del mismo nombre en Honduras, a una altura media sobre el nivel del mar de 807
metros y fue establecido por V. J. Hanraham en diciembre de 1934. El elipsoide de referencia es el
de Clarke de 1866 (Inter-American Geodetic Survey, 1950). El datum WGS84 utiliza el centro de
masa de la Tierra como referencia (es geocéntrico) y por ende no tiene un punto de inicio local como
el Ocotepeque; esto hace que su exactitud puede extenderse a todo el planeta Tierra. Dada esta
discrepancia entre del datum local y el WGS84 es crítico la selección correcta del mismo cuando se
desee integrar datos registrados con un geoposicionador con la cartografía del Instituto Geográfico
Nacional.
Figura 3: Conceptos de geoide, elipsoide y superficie terrestre.
Por ejemplo, si usted registra la posición de cualquier objeto en latitud y longitud en Costa Rica
utilizando el datum WGS84 y luego proyecta los datos a Lambert Norte sin realizar una
transformación de datum, incurrirá en un error radial medio de 264.77 m (Fig. 4). Este
6
desplazamiento sería aceptable para una cartografía a escala 1:600.000 pero está muy por encima de
la exactitud (5-10 m) que brindan los geoposicionadores actuales.
Figura 4: Error medio (m) esperado al utilizar el datum WGS84 con la cartografía de Costa Rica
(Lambert Norte).
Error de los parámetros disponibles para realizar la transformación de datum
El WGS84 es el datum que por omisión utilizan los receptores del Sistemas de Posicionamiento
Global (GPS) para registrar posiciones (coordenadas) sobre la Tierra. Las coordenadas referidas al
sistema WGS 84 pueden transformarse al datum local Ocotepeque utilizando tanto una
transformación de siete parámetros como una transformación de tres parámetros (Molodensky)
(Dörries y Roldan,1999; ESRI, 1994). Para realizar dicha transformación es necesario conocer los
parámetros de los elipsoides en uso. El elipsoide oficial de la cartografía de Costa Rica es el Clarke
de 1866; el cual tiene los siguientes parámetros:
?? Semieje mayor (a): 6378206,4 m
?? Semieje menor(b): 6356583,8 m
?? Factor de achatamiento: Df: (a-b)/a = 0,00339000 y 1/f = 294,978698
?? Factor de excentricidad : e2 = 0,006768657997291
El datum WGS84; basado en la red del Catastro Nacional de 1990; es utilizado en la proyección
Costa Rica Transversal de Mercator (CRTM) en uso por dicha institución y por la cartografía escala
1:25000 del proyecto TERRA. Los parámetros del Sistema Geodésico Mundial de 1984 (WGS84)
son los siguientes:
?? Semieje mayor (a): 6378137.00 m
7
?? Semieje menor (b): 6356752.31 m
?? Factor de Achatamiento: Df: (a-b)/a = 0,0033528 y 1/f= 298,25721
Los valores DX, DY, DZ son las diferencias en los ejes X, Y y Z desde el centro de la tierra en el
modelo WGS84 con respecto al modelo local. La diferencia entre el semieje mayor del eje X del
modelo WGS84 y el modelo local es igual a –69.4 y se denomina Da.
Aun cuando el tema de la transformación del datum WGS84 al Ocotepeque es de interés nacional,
hasta la fecha no existe un set oficial de parámetros de Molodensky para realizar dicha
transformación. Con frecuencia, la literatura que acompaña a los geoposicionadores y a los
programas en Sistemas de Información Geográfica recomiendan utilizar para Costa Rica el datum
NAD27 Central América (ESRI, 1994; NIMA); sin embargo dicho datum no debe utilizarse para
Costa Rica ya que el usuario(a) incurriría en un error radial de 351.53m (Cuadro 1 y figura 5). Al
igual que el caso anterior (datos sin transformar) los desplazamientos serían aceptables para una
cartografía a escala 1:750.000 pero estarían muy por encima de la exactitud (5-10 m) que brindan los
geoposicionadores actuales.
Figura 5: Error medio (m) esperado cuando se utilizan datos registrados con el datum WGS84, luego
se transforman a NAD27-CA y luego se proyectan a Lambert Norte.
Los parámetros suministrados por el señor G. Stentz de la Agencia de Imágenes y Mapas de los
Estados Unidos de Norteamérica (NIMA-USA, email: stentzG@[Link]) mostraron un error radial
medio de 8.32m (S=2.27m) (Cuadro1y Fig. 6); lo cual concuerda con la afirmación del señor Stentz;
en cuanto a que el error esperado de la transformación era de 10m. Si se eliminan los puntos Osa y
Manzanillo el error se reduce a 7.69m (S=1.45m). Este error es apropiado para trabajos cartográficos
a escala 1:17000 o inferiores.
8
Las traslaciones (DX, DY y DZ) derivados por Dörries y Roldan (1999) como parte de su
transformación de siete parámetros mostraron el menor error radial de todas las opciones evaluadas;
la raíz del error medio cuadrático fue de 2.81 m (S=2.24m) y se excluyen los puntos Osa y
Manzanillo dicho valor se reduce a 1.58m (S=1.01m) (Cuadro1y Fig.6). Estos errores son
apropiados para realizar trabajos cartográficos a escalas 1:6000 y 1: 4000 o inferiores;
respectivamente.
Figura 6: Error medio (m) esperado cuando se utilizan datos registrados con el datum WGS84, se
transforman aplicando los parámetros de Stentz (NIMA) y luego se proyectan a Lambert Norte.
9
Figura 7: Raíz del error medio cuadrático (m) esperado cuando se utilizan datos registrados con el
datum WGS84, se transforman aplicando los parámetros DX, DY y DZ derivados por Dörries y
Roldan (1999) y luego se proyectan a Lambert Norte. Los errores de Osa y Manzanillo son muy
superiores a los otros y por esta razón se sugiere eliminarlos del análisis.
Ajustes empíricos
En esta sección del documento se discuten los resultados obtenidos por métodos empíricos, o sea
en este caso no se asume ninguna relación entre elipsoides solo se considera la medición realizada
en los dos datum en puntos homólogos. Este método de transformación podría ser útil para aquellos
usuarios que no disponen de ArcView GIS.
Traslación de ejes de Rody Rodríguez
La traslación de ejes sugerida por el Ing. Rody Rodríguez, funcionario de Área de Topografía del
departamento de Estudios Básicos, UEN PySA del Instituto Costarricense de Electricidad mostró un
error radial promedio de 6.6m (S=2.89m); el cual se redujo a 6.1m (S=2.56m) al eliminar los puntos
Osa y Manzanillo (Cuadro 1 y Fig 7). El error puede considerarse como relativamente pequeño
cuando se utilizan georeceptores con errores en el orden 10 a 15 m; sin embargo sería superior al
error esperado para los nuevos receptores con WAAS (Wide Area Augmentation System). En
general, puede considerarse que el error observado es apropiado para trabajos cartográficos a escala
1:14000 o inferiores.
10
Figura 7: Error medio (m) esperado cuando se utilizan datos registrados con el datum WGS84, se les
aplica la traslación de ejes del Ing. Rody Rodríguez y luego se proyectan a Lambert Norte.
Ajuste ecuación lineal
La ecuación de regresión lineal ajustada a los 28 puntos homólogos mostró un error medio de
4.38m (S=1.48m); el cual se redujo a 4.01m (S=1.12m) al eliminar los puntos Osa y Manzanillo
(Fig. 8 y cuadro 1). El error de Osa (8.5m) y Manzanillo (6.8m) fue inferior al registrado para la
traslación de ejes sugerida por Rody Rodríguez y por tanto puede considerarse que los parámetros
del modelo pueden considerarse como una alternativa simple para usuarios que no disponen de
ArcView.
Además del modelo nacional, también se ajustaron modelos para el Caribe y el Pacífico de Costa
Rica. El error medio para el Caribe se redujo a 2.17m (S=0.68m); sin embargo para el Pacífico
aumentó a 4.13m (S=1.22m). Los errores permiten concluir que el ajuste de una ecuación lineal
podría ser una alternativa válida para aquellas personas que no poseen ArcView GIS y que por tanto
no pueden utilizar la extensión CR_proyección_datum. Según la compañía ESRI (1994) el error
esperado para la conversión de datum utilizando un ajuste lineal es de aproximadamente 5m; lo cual
es coincidente con el error obtenido para el set de datos de Costa Rica. El error sería apropiado para
trabajos cartográficos a escala 1:10.000 o inferior.
11
Figura 8: Error medio (m) esperado al ajustar una ecuación de transformación lineal a 28 puntos
homólogos.
Cuadro 1: Error esperado de varios métodos y parámetros de transformación del datum WGS84 a
Ocotepeque.
Método/ Parámetros Error radial (m)
REMC * Media Mín. Máx. Des. Estándar Rango
Dörries y Roldán( 1999) 2.81 1.76 0.23 10.70 2.89 10.47
Stentz 8.34 8.04 4.66 14.65 2.27 9.99
**
Ecuación lineal 4.01 3.86 0.68 5.81 1.12 5.13
Ecuación lineal (Pac) 4.13 3.96 2.36 5.31 1.22 5.40
Ecuación lineal (Car) 2.17 2.08 0.92 3.10 0.68 2.18
Rody Rodríguez 6.57 5.92 2.24 13.12 7.83 10.87
Sin transformar 264.77 264.75 256.60 270.57 3.33 13.97
NAD27 CA 351.53 351.52 343.28 354.98 2.48 11.70
*
Raíz del error medio cuadrático ** Excluye Osa y Manzanillo Pac: Pacífico Car: Caribe
12
Figura 9: Error radial (m) esperado para la transformación de datum utilizando la traslación de ejes
sugerida por Rody Rodríguez, una ecuación lineal a nivel nacional y una ecuación ajustada al Caribe
de Costa Rica.
Error en Osa y Manzanillo
Los puntos transformados de Osa (Fig.10) y Manzanillo (Fig.11) mostraron un error muy
superior con respecto a los otros puntos para todos los métodos y parámetros de Molodensky
probados en el presente estudio. Este es un aspecto que requiere de mayor investigacion.
Figura 10: Error observado en el punto Osa. Los parámetros de Dörries y Roldan (1999) que
mostraron el mejor ajuste de todos los evaluados, en este caso produjeron un error de 10.7 metros.
13
Figura 11: Error observado en el punto Manzanillo. Los parámetros de Dörries y Roldan (1999) que
mostraron el mejor ajuste de todos los evaluados, en este caso produjeron un error de 6.5 metros.
La extensión CR_proyecciones_datum
Una vez seleccionado el mejor set de parámetros de transformación procedí a automatizar el
proceso de transformación de datum mediante una extensión que operara con el programa ArcView
GIS 3.x (Fallas, 2003). La extensión se compone de tres botones los cuales realizan las siguientes
operaciones:
?? Transformación de datum. Este botón le permite transformar sus datos del datum WGS84 a
Ocotepeque y viceversa.
?? Cambio de proyección. El usuario(a) puede elegir entre las siguientes proyecciones Lambert
Norte, Lambert Sur, CRTM y UTM.
Error esperado al proyectar datos de CRTM a Lambert y viceversa
La extensión CR_proyecciones_datum también permite realizar una transformación proyectiva
entre las proyecciones CRTM, Lambert Norte, Lambert Sur y UTM. Utilizando los 28 puntos cuya
posición se conoce tanto en CRTM98 como en Lambert Norte se obtuvo un error (REMC) de 0.41m
para la transformación de Lambert Norte a CRTM y de 2.89m para la transformación de Lambert
Sur a CRTM; excluyendo los puntos de Osa y Manzanillo, el error se reduciría a 1.21m. Se esperan
errores similares para la transformación de CRTM a Lambert Sur y viceversa.
14
Conclusiones
Los resultados indican que los parámetros seleccionados de Dörries y Roldan (1999) permiten
realizar una conversión de datum de WGS84 a Ocotepeque con un error medio cuadrático de 2.81m
ó de 1.58m se si excluyen los puntos Osa y Manzanillo. El error es homogéneo y consistente a lo
largo y ancho del país. A los usuarios de ArcView GIS se les recomienda utilizar la extensión
CR_proyecciones_datum para realizar la transformación de datum. En caso de no tener acceso a
ArcView GIS se le recomienda utilizar las siguientes ecuaciones para realizar la transformación de
datum:
Cobertura espacial Ecuación REMC (m)**
Nacional* Lat_ocote = 0.000684+1.000074*lat_wgs84 4.01
N= 26 Long_ocote = -0.0013*1.000008*long_wgs84
Pacífico Lat_ocote = 0.000777+1.000064* long_wgs84 4.13
N=17 Long_ocote = -0.001488+1.000005*lat_wgs84
Caribe Lat_ocote = 0.000406+1.000101*lat_wgs84 2.17
N=8 Long_ocote = 0.000436*1.000028*long_wgs84
*
No incluye Osa ni Manzanillo
**
Raíz del error medio cuadrático
Referencias
Acharya, B. y Bell, W. 1992. Universal accuracy standards for a GIS. In. ASPRS/ACSM/RT 92
Technical Papers. Vol 3. GIS and Cartography. pp. 368-379. American Society for Photogrammetry
and Remote Sensing and American Congress on Surveying and Mapping. Bethesda, MD. USA.
Ashtech. 1997. Reliance under-canopy performance evaluation. USDA Forest Service Clackamas
GPS test Network. USA. 7p. [Link]
Barrantes Ferrero, Mario. 1966. El mapa básico de Costa Rica. Instituto Geográfico Nacional,
Ministerio de Transportes. 18p.
Bennett, V. 1990. GPS program status report. In. Protecting Natural Resources with remote sensing.
Proceedings of the Third Forest Service Remote Sensing Applications Conference-April 9-13, 1990.
pp.334-338 American Society for Phtogrammetry and Remote Sensing. Bethesda, Maryland USA.
Bolstad, P. V. 1993. GPS basics: Forestry applicationes. The Compiler, Vol.11(3): 4-8.
Kruczynski, L. y Jasumback, A. 1993. Forestry management applications. Forest Service experience
with GPS. Jour. of Forestry. Vol. 91(8): 20-24.
Deckert, C. y Bolstad, P.V. 1996. Forest canopy, terrain, and distance effects on Global Positioning
System point accuracy. Photogrammetric Engineering and Remote Sensing. Vol. 62(3): 317-321.
Dörries y Roldán, 1999. Estudio comparativo del datum geodésico de Ocotepeque y el datum
satelitario del sistema WGS84. Informe final del proyecto de investigación. Escuela de Topografía,
Catastro y Geodesia. Facultad de Ciencias Exactas y Naturales. Universidad Nacional. Heredia,
Costa Rica. 172p.
15
ESRI, 1994. Map projections. Georeferencing spatial data. Environmental System Research
Institute. United States of America.
Fallas, J. 2003. Extensión CR_proyecciones_datum.avx para ArcView GIS. Laboratorio de
Teledetección y Sistemas de Información Geográfica. Universidad Nacional. Costa Rica.
jfallas@[Link]
Giovachino, D. 1992. Graphic digitizing accuracy. In. ASPRS/ACSM/RT 92 Technical Papers. Vol
3. GIS and Cartography. pp. 380-389. American Society for Photogrammetry and Remote Sensing
and American Congress on Surveying and Mapping. Bethesda, MD. USA.
Hurn, J. 1989. GPS a Guide to the next utility. Trimble [Link]. 76p.
Inter-American Geodetic Survey. 1950. Proyección Lambert Norte para Costa Rica. Army Map
Service, Washington, D.C. USA. 25p.
Oregon State University, Department of Forestry.2000. Extensión Xtools para ArcView Gis.
[Link]/sfgis/
Registro Nacional. Costa Rica. Transformación del antiguo (Datum de Ocotepeque) al nuevo
(Datum WGS84) Tres traslaciones, un factor de escala y tres giros, entre el Elipsoide de Clarke de
1866 y Datum de Ocotepeque y el Elipsoide asociado al Sistema WGS84, calculados con datos de
1990. [Link] registronac@[Link]
Rodríguez, Rody. 2003. Valores para la traslación de ejes para la trasformación de datum WGS84 a
Ocotepeque. Comunicación persona Área Topografía, Estudios Básicos, UEN PySA. Instituto
Costarricense de Electricidad RRodriguezM@[Link] Tel 506 2205123.
Stentz, G. 1997. Parámetros para transformación de datum. Ocotepeque a WGS. NIMA, USA.
stentz@[Link]
Wilkie, D. S. 1989. Performance of a backpack GPS in a tropical rain forest. Photogrammetric
Engineering and Remote Sensing. 55: 1747-1749