Centro de Transferencia de Tecnología en Transportación Departamento de Ingeniería Civil y Agrimensura UPR-‐Recinto Universitario de Mayagüez Call Box 9000 * Mayagüez, PR 00681 Tel. 787-‐834-‐6385 * Fax: 787-‐265-‐5695 * www.uprm.edu/prt2 29 Años de Excelencia en el Adiestramiento de Oficiales de Transportación a Nivel Municipal, Estatal, y Federal en Puerto Rico e Islas Vírgenes Instructora Prof. Linda L. Vélez Rodríguez Catedrá/ca Departamento de Ingeniería Civil y Agrimensura UPR – Recinto Universitario de Mayagüez 8 de mayo de 2015 Índice Módulo 1 – Notas de la Instructora Módulo 2 – Definiciones de Datums Módulo 3 – Base de Datos Geográficos Módulo 4 – Geo-‐Based Informa/on Matrix Módulo 5 – Ley Núm. 264 del año 2002: Ley para adoptar el Sistema de Coordenadas Planas Estatales de Puerto Rico Módulo 6 – Sección 5.2.3 -‐ Sistema de Coordenadas Planas Estatales Módulo 7 – Modern Terrestrial Reference Systems – Part 1 Módulo 8 – Modern Terrestrial Reference Systems – Part 2: The Evolu/on of NAD 83 Módulo 9 – Modern Terrestrial Reference Systems – PART 3: WGS 84 and ITRS Módulo 10 – Modern Terrestrial Reference Systems – Part 4: Prac/cal Considera/ons for Accurate Posi/oning Módulo 11 – Notas sobre Agrimensura Módulo 1 Notas de la Instructora CONCEPTOS BÁSICOS DE SISTEMAS DE COORDENADAS EN LA RELOCALIZACIÓN DE UTILIDADES EN PROYECTOS DE CARRETERAS por Prof. Linda L. Vélez-Rodríguez, MS, PE, PLS Catedrática Univ. de PR-Recinto Universitario de Mayaguez e-mail: [email protected] Tel. 787-265-5405 (Ofic) 787-313-4740 (Cel) Temas - AM • Introducción • Requerimientos por Ley y Reglamento • Conceptos Básicos de: – Gedésia y sus superficies – Marcos de Referencia y la Ley 264 del 2002 – Agrimensura y como se obtienen las Coordenadas Planas X,Y • Referencia “NOTAS SOBRE AGRIMENSURA” – Cartografía y sus proyecciones cartográficas • Referencia “SISTEMA LAMBERT” en www.revistatp.com Temas - PM • National Geodetic Survey (NGS): sus Datasheets y sistemas de coordenadas – Referencia: www.ngs.noaa.gov/ • Coordenadas Geo-Referenciadas N,E • Utilidades y sus Atributos • Caso de carretera PR-106 Km 1.6 Mayaguez • Carretera PR-3108 conecta la PR-2 y PR108 Realidad vs. transformación Realidad vs. transformación Estructura VS. Datos Geo-Espaciales Que nos requieren • Ley 264 del 16 de noviembre del 2002 - Para adoptar el Sistema de Coordenadas Planas Estatales utilizando el método de la Proyección Conforme Cónica Lambert y el North American Datum of 1983, o su versión más reciente Que nos requieren • Reglamento Conjunto para Obras de Construcción y Usos de Terrenos vigende desde el 29 de noviembre del 2010 en su Sección 5.2.3 – Sistema de Coordenadas Planas Estatales Geodesia • Se define como la rama de la matemática aplicada que tiene que ver con el estudio de la forma y tamaño de la Tierra al igual que su campo gravitacional. Es una de las Geo-Ciencias o Ciencias Terrestres, junto con: Geografia, Geologia, Geofisica y Geomorfologia Eratostenes: Padre de la Geodesia Tres Superficies Una Realidad • Topografía - superficie de nuestro planeta donde nosotros hacemos las observaciones • Elipsoide - superficie matemática que usamos para describir la tierra • Geoide- superficie física que denota el potencial de gravedad de nuestro planeta Geoide • nosotros habitamos un planeta cuya forma es única y puede ser descrita por el término geoide. Esta palabra procede del griego “geoeides” que quiere decir “perteneciente a la tierra” (GE de tierra y EIDOS de forma). • El geoide es la forma que tiene la superficie de la tierra allí donde el nivel del mar es uniforme y siempre es perpendicular a la dirección de la gravedad. EL GEOIDE ¿Que es ? TEORICAMENTE: G UNA SUPERFICIE A NIVEL EN ESPECIFICO DEL CAMPO GRAVITACIONAL EN LA PRACTICA: • REFERENCIA CERO PARA ELEVACIONES Y PROFUNDIDAD • “MEAN SEA LEVEL” CONTINUO EN LOS CONTINENTES • LA SURPERFICIE PERPENDICULAR A LA DIRECCION DE LA GRAVEDAD 5000 ft 3000 ft G G G G En busca del Geoide… ONDULACION GEOIDE PERPENDICULAR al GEOIDE (linea de la plomada) PERPENDICULAR al ELIPSOIDE DEFLECCION DE LA VERTICAL RELACION GEOIDE-ELIPSOIDE Relación de alturas (height) y las Tres superficies Gravity Recovery And Climate Experiment (GRACE) DATUMS • Un DATUM es un conjunto de parámetros que definen un sistema de coordenadas y un conjunto de puntos de control cuya relación geométrica es conocida ya sea por medidas o cálculos. Dewhurst, 1990 • Todos los DATUMS se fundamentan en un elipsoide, el cual aproxima la forma de la Tierra. EL ELIPSOIDE MODELO MATEMATICO DE APROXIMACION DE LA TIERRA a = Semi Eje Mayor b = Semi Eje Menor f = a-b = Achatamiento a N b a S Datum Horizontales y sus Elipsoides • • • • PR Datum NAD27 NAD83 WGS84 Clarke 1866 Clarke 1866 GRS80 WGS84 NORTH AMERICAN DATUMS NAD 27 & NAD 83 EL GEOIDE Y DOS ELIPSOIDES GEOIDE CLARK 1866 ELLIPSOID CENTRO DE MASA GRS 80 ELLIPSOID Mezcla de Datums Mapa en PR Datum y Datos en WGS84 NSRS Coordinate Systems Latitud & Longitud Coordenadas Planas Estatales Coordenadas UTM “Earth-Centered Earth-Fixed” NAD 83 NAD 27 NGVD 29 NAVD 88 PRVD 02 ITRF00 AGRIMENSURA: localización = coordenadas • Se define como el arte, ciencia y tecnologia de determinar la posicion de un punto en, sobre o bajo la superficie de la tierra, o establecer la posicion de dichos puntos. Envuelve todo lo que tenga que ver con nuestro ambiente pues TODO tiene una definicion en el espacio que habitamos, y de eso es de lo que se trata el NEGOCIO de la AGRIMENSURA, de la localizacion. La superficie se considera plana, asi que se geometria plana. Al agrimensor se le conoce como geometra. Medidas Básicas en Agrimensura • • • • • Distancias Horizontales Distancias Verticales Distancias en Pendiente Angulos Horizontales Angulos Verticales Total Station, Prisma y Baston con Bipode Distancias Ángulos Poligono de Control que se usa en agrimensura; existen poligonos cerrados, o poligonos abiertos, o lineas bases o un punto. Los mas optimos son los poligonos cerrados por la redundancia que ofrecen. Se denominan de control pues es un control o referencia para obtener los datos necesarios para hacer el levantamiento de los datos necesarios para hacer nuestro diseño En los vertices del Poligono de Control debe existir una marca Permanente que se denomina PUNTO DE CONTROL, muchos Agrimensores usan varillas de acero. Se procede a medir los angulos interiores del poligono cerrado y Las distancias entre dichos vertices, y con estos datos se procede a ajustar esos datos (angulos interiores y distancias) Datos sin ajustar nos dan una diferencia en la ubicación de A A’ Las coordenadas se calculan usando los "Balanced Latitude" y "Balance Departure". Estos son los componentes en Y y X ajustados por eso “Balanced” La coordenada X de un punto es igual a la coordenada X del punto anterior más el "Balanced Departure" entre esos puntos. La coordenada Y de un punto es igual a la coordenada Y del punto anterior más el "Balanced Latitude" entre esos puntos. Esto es: Xn = Xn-1 + "Balanced Departure”(n @ n-1) Yn = Yn-1 + "Balanced Latitude”(n @ n-1) Si usamos la analogia del vector y sus componentes en X & Y; el vector es uno del lados de la poligonal y su componente en X es del Departure y el componente en el eje de Y es el Latitude, refierase al poligono ABCD ilustrado a continuacion: El Ajuste de los datos se puede hacer con uno de varios metodos disponibles, que cada uno de ellos pondera una u otras de los observaciones que se hagan. Sistemas de Coordenadas Planas Sistemas de Coordenadas Esfericas COORDENADAS Por definición son cantidades lineales o ángulares, o ambas, que designa la posición de un punto relativo a un marco de referencia. Se dividen en agrimensura generalmente en dos: (1) polares y (2) rectangulares. Se subdividen en tres clases: coordenadas; (1) planas, (2) esfericas y (3) Espaciales. Distancia Lineal DISTANCIA LINEAL • Entre dos puntos con coordenadas conocidas tenemos que la distancia D1-2 es igual a la raiz cuadrada del cuadrado de las diferencias de las Y mas el cuadrado de las diferencias de las X. • Es decir: D1-2 = ((Y2 - Y1)2 + (X2 - X1) 2)1/2 Rosa de los vientos ANGULOS/DIRECCIONES • Para rumbos y acimuts tenemos que la tangente del rumbo o acimut es igual a las diferencias de las X dividido por la diferencia de la Y. • Es decir: tan B1-2 = X2 - X1/ Y2 - Y1 AREA CON COORDENADAS • Dos veces el área de un polígono cerrado usando coordenadas se calcula de la siguiente forma: • 2A = X1(Y2 - Yn) + X2(Y3 - Y1)+. . .+ Xn-1(Yn - Yn-2) + Xn(Y1 - Yn-1) ó • 2A = Y1(X2 - Xn) + Y2(X3 - X1)+. . .+ Yn-1(Xn - Xn-2) + Yn(X1 - Xn-1) Se usa el valor absoluto, pues una sera el negativo de la otra, y no tenemos áreas negativas. Que nos requieren: •Ley 264 del 16 de noviembre del 2002 •Sección 5.2.3 – Sistema de Coordenadas Planas Estatales del Reglamento Conjunto para Obras de Construcción y Usos de Terrenos vigende desde 29 de noviembre del 2010 TRANSFORMACION Se define como una función que relaciona coordenadas en un sistema de coordenadas con otro sistema de coordenadas. Usualmente envuelve cambios en traslación, rotación y/o escala. - “affine transformation” - “conformal transformation” LAS PROYECCIONES CARTOGRAFICAS SON UNA ABSTRACCION MATEMATICA DE LA REALIDAD Matriz del GIS de Puerto Rico Sistema Nacional de Referencia Espacial “National Spatial Reference System” (NSRS) De esto consiste el Sistema Nacional de Coordenadas • Latitud • Longitud • Altura • Escala • Gravedad • Orientacioón Y como estos valores cambian con el tiempo Variedad de fuentes de datos geo-espaciales – OJO con GI-GO Datos de Infraestructura – – – – – – – – – – Localización Geográfica Carreteras Aceras Servicios de Agua Servicios de Alcantarillado Sanitarios Luminarias en Calles Servicios Electricos Tuberías (pipelines) Facilidades de Comunicación Servicios de Control de Trafico – Facilidades de Parques/Recreación – Carriles (lanes) – Escuelas – Inventario de Arboles/Arbustos – Estacionamiento – Puentes – Utilidades extranjeras – Red/Facilidades de Carreteras (Viales) Infraestructura y Atributos – Carreteras • • • • • • • • Tipo Superficie Fecha de Construcción Propietario Fecha del Ultimo Mantenimiento Número de Carriles Ancho Profundidad – Aceras • • • • • • Tipo Superficie Fecha de Construcción Propietario Fecha del Ultimo Mantenimiento Ancho/Profundidad Infraestructura y Atributos – Luminarias en Calles • • • • • Postes Transformadores Alambrado Edad Profundidad – Servicios de Control de Trafico • • • • Tipo Activo/Pasivo Fecha de Instalación Fecha del Ultimo Mantenimiento Infraestructura y Atributos – Puentes • • • • • – Tipo Superficie Fecha de Construcción Número de Carriles Fecha del Ultimo Mantenimiento Utilidades extranjeras • • • • • • • • • Cablevision Propietario Localizacion de Planta Tuberías Propietario Capacidad Localizacion de Planta Caracteristicas Especiales Profundidad – Redes Viales • • • • Tipo Ruta Rasgos Especiales Facilidades – Carriles (lanes) • • • • • Localización Tipo Superficie Propietario Dimensiones Fisicas Infraestructura y Atributos – Servicios de Agua • • • • • • • • • Hidrantes Capacidad Válvulas Líneas Presión/Edad Profundidad Fecha de Instalación Fecha de Ultimo Mantenimiento Estaciones de Bombeo – Servicios de Alcantarillado Sanitarios • • • • • Líneas Válvulas Profundidad Estaciones de Bombeo Tipos Utilidades Soterradas ¿Como las ubicamos? Cartografía • Definición: – Ciencia, Arte y Tecnología de hacer mapas junto con su estudio como documento científico y trabajo de arte. Los mapas son una representación de la tierra o parte de ella a escala. – La escala es la relación entre el tamaño real del área y la porción reducida que aparece en el mapa. Son productos cartográficos, los planos, secciones transversales, modelos en tercera dimensión, entre otros. Escala Escala como fracción representativa: Una razón, tal como 1:1,000,000. Si se asume que el mapa esta en centímetros, significa que un centímetro en el mapa es igual a un millón de centímetros en la Tierra. Esta escala se podría también escribir como fracción proporcional (1/1,000,000). Superficies Desarrollables • son superficies que se desarrollan en un plano, fundamentales para las Proyecciones Cartograficas. Ellas son: – el cono – el cilindro y – el plano Sistema de Coordenadas Planas Estatales • El gobierno federal desarrollo el Sistema de Coordenadas Planas Estatales (SPCS)para cada estado y sus territorios. • En el caso en particular de Puerto Rico e Islas Virgenes se usa la proyección conforme cónica Lambert con dos paralelos standard, teniendo los siguientes parametros para el NAD 83: φN = 18°-26’ N; φS = 18°-02’ N; φ0 = 17°-50’ N; λ0 = 66°- 26’ W; Nb = 200,000.m; E0 = 200,000.m Sistema de Coordenadas Planas Estatales Limites de las zonas son 158 millas Sistema de Coordenadas Planas Estatales Sistema de Coordenadas Planas Estatales = f (, ) = f (, ) DESIGNATION - VÉLEZ PID - DE5545 COUNTRY - PUERTO RICO NAD 83(2002)- 18 26 41.28060(N) SPC PRVI UTM 19 - 067 08 48.93357(W) ADJUSTED North East Units Scale Converg. 267,825.241 124,618.857 MT 1.00000071 -0 13 23.8 2,040,399.173 695,701.034 MT 1.00007349 +0 35 11.3 Universal Transversa Mercator • Para el sistema Universal Transversa Mercator, mejor conocido por sus siglas en ingles UTM, el globo esta dividido en 60 zonas. – Se asume que el esferoide del DATUM es una esfera • Cada zona se extiende seis grados de longitud, con un solape de 30 minutos con las zonas adyacentes. • Cada zona tiene su propio meridiano central del cual se extiende 3 grados al este y 3 grados al oeste Universal Transversa Mercator • Las coordenadas de este sistemas están en metros. • El origen de cada zona es en el ecuador y en su meridiano central. • El valor dado al meridiano central es un "falso este" igual a 500,000 metros. Universal Transversa Mercator • La cuadrícula del sistema UTM se extiende alrededor del globo desde la latitud 80° N a la 80° S. • Cada zona en el Ecuador tiene dos falsos nortes 0 metros para la mitad norte y 10,000,000 metros en la mitad sur. – Lo cual significa que cada Iocalización en el Ecuador tiene dos conjuntos de coordenadas UTM 60 Zonas del Sistema UTM Estas zonas dependen de la Longitud # de Zona = (180° - λW)/6 para Longitudes al Oeste (λW ) del Meridiano de Greenwich # de Zona = (180° + λE)/6 para Longitudes al Este (λE) del Meridiano de Greenwich Si el resultado de la ecuación no es un entero se designa el entero proximo como el número de zona Puerto Rico esta en las zonas 19 y 20 del sistema UTM Red de Triangulación Red del CRIM del 1995 (Centro de Recaudación de Ingresos Municipales) Red de Control establecida en 1995 por el CRIM Las Marias 2 Control Geodésico establecido en 1995 por el CRIM Posiciones • Horizontal • Vertical • Temporal , o , , ,h o ,, t Datum Verticales • Los Datums Verticales son locales – Puerto Rico Vertical Datum 2002 (PRVD02) – Se usa en PR en el epoch del 1982-2001 de datos del mareógrafo de La Puntilla – El cero esta en las Oficinas de la Guardia Costanera en La Puntilla en el Viejo San Juan – Se corrió una nivelación en Mayo del 2002 desde La Puntilla hasta Aguadilla bajo el auspicio del National Geodetic Survey (NGS) – Este proyecto del PRVD02 lo están realizando dos firmas de agrimensores de Puerto Rico, siguiendo los “standards” del NGS pues es una nivelación de primer orden. Rutas propuestas en el proyecto de control vertical de Puerto Rico-PRVD02 para un total de 925 kilómetros incluyendo Vieques y Culebra PRVD02 –Importancia • Como parte de este proyecto se realizaron observaciones gravimétricas en Mayagüez y Aguadilla, y se efectuó un vuelo en enero del 2009 con un gravímetro en el avión como parte del proyecto denominado GRAV-D (Gravity for the Re-definition of the American Vertical Datum). • Se integrara esta tecnología para tener un solo datum uniendo el Horizontal y el Vertical, obteniendo asi alturas ortometricas derivadas. Daniel Winester, Geodesta del NGS monumentando estación en CI-019-RUM MAYAGUEZ A A 2008 MAYAGUEZ A A 2008 Observaciones Gravimetricas por 48 horas PRVD02 –Importancia • El gobierno del Estado Libre Asociado de Puerto Rico aporto $3,000,000 por conducto de tres de sus agencias ellas son la Oficina de Presupuesto y Gerencia; el Departamento de Transportación y Obras Públicas y la Autoridad de Energía Eléctrica. • El Colegio de Ingenieros y Agrimensores de Puerto Rico ha tenido un rol muy importante en este proyecto, pues se esta muy conciente de la importancia del mismo y lo que significa para el desarrollo de la infraestructura del país, junto con los cambios que ocasionara en el desempeño de la profesión de la agrimensura en lo referente a obtener elevaciones. Controles Verticales del USGS Vélez 2002 Vélez 2002 PID - DE5545 • HORZ ORDER - B • NAD 83(2002)– φ=18° 26’ 41.28060” N λ=67° 08’ 48.93357” W – SPC PRVI – N=267,825.241m E=124,618.857m – UTM 19 – N=2,040,399.173m E=695,701.034m • VERT ORDER - FIRST CLASS II • PRVD02 - 134.320meters 440.68feet National Geodetic Survey, Retrieval Date = APRIL 26, 2007 *********************************************************************** DESIGNATION - VELEZ PID - DE5545 STATE/COUNTY- PR/AGUADILLA H=h-N USGS QUAD - DE5545 DE5545 DE5545 DE5545 DE5545 DE5545 134.32 = 90.42 - (- 44.20) DE5545 *CURRENT SURVEY CONTROL 134.32 134.62 DE5545 ___________________________________________________________________ DE5545* NAD 83(2002)- 18 26 41.28060(N) 067 08 48.93357(W) ADJUSTED DE5545* PRVD02 134.320 (meters) 440.68 (feet) ADJUSTED DE5545 ___________________________________________________________________ DE5545 X 2,350,644.082 (meters) COMP DE5545 Y - -5,577,493.184 (meters) COMP DE5545 Z 2,005,175.775 (meters) COMP DE5545 LAPLACE CORR3.52 (seconds) DEFLEC99 DE5545 ELLIP HEIGHT90.42 (meters) (08/30/05) GPS OBS DE5545 GEOID HEIGHT-44.20 (meters) GEOID03 DE5545 DE5545 HORZ ORDER - B DE5545 VERT ORDER - FIRST CLASS II DE5545 ELLP ORDER - THIRD CLASS I NAD 83 (NSRS 2007) = 90.388 DE5545 DE5545.The horizontal coordinates were established by GPS observations DE5545.and adjusted by the National Geodetic Survey in August 2005.. DE5545.This is a SPECIAL STATUS position. See SPECIAL STATUS under the DE5545.DATUM ITEM on the data sheet items page. DE5545 DE5545.The orthometric height was determined by differential leveling DE5545.and adjusted by the National Geodetic Survey in April 2003 Los GPS (Sistemas de Posicionamiento Global) Segmentos que componen el sistema GPS Civil GPS Use Satellite Ops -Ephemeris, Timing Power Grid Interfaces Personal Navigation Surveying & Mapping Trucking & Shipping Aviation Recreation Communications -Network Synchronization and Timing Railroads Fishing & Boating Off shore Drilling Un Satelite Dos Satelites Tres Satelites Cuatro Satelites X23 y23 z23 d23 x14 y14 z14 x19 y19 z19 d19 x21 y21 z21 d14 d21 Measured: x y z Non-Differential GPS (Autonomous or Stand-alone) X23 y23 z23 x19 y19 z19 Measured: x y z Delta: x y z _________ Corrections applied after survey True: x y z Differential GPS x14 y14 z14 x21 y21 z21 True: x y z Measured: x y z ______________ Delta: x y z X23 y23 z23 x19 y19 z19 x14 y14 z14 x21 y21 z21 Measured: x y z Delta: x y z _________ True: x y z Corrections applied during the survey True: x y z Measured: x y z _____________ Delta: x y z Real-Time Differential GPS GPS • • Los datos obtenidos con los receptores de GPS se pueden incorporar a los mapas topográficos en formato digital, si estos estan en el mismo datum y sistemas de coordenadas. Como aplicacion, si obtenemos las coordenadas de los límites municipales de los cuadrángulos estos límites se pueden establecer sobre el terreno usando GPS Transformación de Coordenadas • El procedimiento para transformar unas coordenadas consiste de tres partes: – Escala – Rotación – Traslación • Se tiene un listado de coordenadas planas y se conocen las coordenadas planas del sistema al que se transformará en dos o más puntos. Transformación de Coordenadas • Frecuentemente se trasforman coordenadas de sistemas arbitrarios a sistemas estales o globales, por lo general con la observación de GPS/RTK en dos de los puntos de controles (Linea Base). • Al completar las observaciones de GPS/RTK se obtienen coordenadas en el sistema de coordenadas planas estatales (SPCS) o coordenadas geodésicas. • De obtener coordenadas geodésicas existen herramientas (Geodetic Tools) del National Geodetic Survey (NGS) para la transformación de estas al SPCS. Transformación de Coordenadas • El procedimiento discutido a continuación considera que se conocen las coordenadas en dos de los puntos de control en el sistema al que se transformará. • De conocer las coordenadas en más de dos controles puedes completar el procedimiento con dos de los puntos y utilizar los demás para cotejar y comparar los resultados obtenidos de la transformación. Transformación de Coordenadas • De desear utilizar más de dos puntos para establecer la solución de la transformación, deberá completar un procedimiento de mínimos cuadrados. • Esto pudiera ser recomendable en proyectos de grandes extensiones, donde el cambio en escala en distintas localizaciones del proyecto y los efectos de curvatura deben ser considerados. Transformación de Coordenadas • Procedimiento: Para diferenciar ambos sistemas de coordenadas, se identificaran: X, Y sistema de original E, N sistema al que se transformará. de tener algún subscrito éste identificará el punto. Transformación de Coordenadas • Aplicación de escala E = SX N = SY • Aplicación de traslación E = X + ΔE N = Y + ΔN • Aplicación de escala, traslación y rotación E = (S cos Θ) X – (S sin Θ) Y + ΔE N = (S sin Θ) X + (S cos Θ) Y + ΔN Transformación de Coordenadas • Luego de considerado los efectos de escala, rotación y traslación, obtenemos las siguientes ecuaciones para la obtención de E y N: (S cos Θ) X – (S sin Θ) Y + ΔE = E (S cos Θ) Y + (S sin Θ) X + ΔN = N Transformación de Coordenadas • Con el objetivo de simplificar las pasadas ecuaciones, se reescribirán considerando: a = S cos Θ c = ΔE b = S sin Θ d = ΔN aX–bY+c=E aY+bX+d=N Transformación de Coordenadas • Para continuar se utilizará un listado de coordenadas que contiene los controles 1 y 2, de los cuales se obtuvieron coordenadas planas estales mediante el uso de GPS. • Al tener datos de dos de los puntos, obtendríamos cuatros ecuaciones con cuatros desconocidas (a, b, c, d): a X1 – b Y1 + c = E1 a Y1 + b X1 + d = N1 a X2 – b Y2 + c = E2 a Y2 + b X2 + d = N2 Transformación de Coordenadas • Reescribiremos el sistema de ecuaciones de forma matricial para su solución. a X1 – b Y1 + c = E1 a Y1 + b X1 + d = N1 a X2 – b Y2 + c = E2 a Y2 + b X2 + d = N2 X1 Y 1 X2 Y2 Y1 X1 Y2 X2 1 0 0 1 1 0 0 1 E1 a N b 1 E2 c d N2 Transformación de Coordenadas • Recordando que: X1 a Y b 1 X2 c d Y2 Y1 X1 Y2 X2 1 0 0 1 1 0 0 1 1 E1 N 1 E2 N2 podemos así obtener los valores de las variables a, b, c y d. Transformación de Coordenadas • Considerando los siguientes, obtenemos como resultado de las variables a, b, c y d: • • • • Punto X Y E N 1 1000.000 1000.00 133880.807 241647.576 2 1032.454 1018.109 133910.982 241625.858 a = 0.424274607 b = -0.905934788 c = 132550.5976 d = 242129.2362 Transformación de Coordenadas • Translación ΔE = c = 132550.5976 ΔN = d = 242129.2362 • Rotación b a - 0.905934788 tan 1 1.132805821rad 64.90499255 0.424274607 Escala S a cos 0.424274607 S 1.000363325 cos(1.132805821) tan 1 • Transformación de Coordenadas • Una vez obtenidos todos los valores de a, b, c y d, se puede computar las coordenadas de cualquier punto utilizando las ecuaciones antes mencionadas: aX–bY+c=E E=aX–bY+c aY+bX+d=N N =aY+bX+d • Ejemplo: (Punto 12) E12 = a X12 – b Y12 + c E12 = a 1001.019 – b 1026.355 + c = 133905.115 N12 = a Y12 + b X12 + d N12 = a 1026.355 + b 1001.019 + d = 241657.835 Transformación de Coordenadas • Unidades – Al momento de efectuar la transformación, las coordenadas originales pueden encontrarse en cualquier sistema de unidades. – En caso que el sistema de coordenadas al cual se transformará sea distinto, éste se armonizará incluyendo la conversión en la escala. – Las coordenadas resultantes contendrán las unidades utilizadas por el sistema de coordenadas al cual se transformó. Transformación de Coordenadas • Escala – La escala incluirá el factor para el cambio de unidades entre sistemas en caso de ser unidades distintas. – Incluirá el factor de escala combinado perteneciente a la proyección cartografica utilizada de las coordenadas transformadas. • En caso de que las coordenada originales y las que se transformarán tengan las misma unidades, el factor escala debe ser similar al que obtendría del reporte de las observaciones de GPS en la localización. Transformación de Coordenadas • Escala – Al computar la escala en la transformación se está utilizando la distancia entre los puntos observados con GPS. Por tal razón al utilizar resultados de la observación de GPS erróneos o de pobre calidad, no solo se está localizando (trasladando y rotando) incorrectamente, sino que también se están modificando todas las dimensiones (escala) erróneamente. Transformación de Coordenadas • Escala – Recordemos que al trasformar las coordenadas a un sistema que proviene de una proyección cartografica, estas coordenadas son propias de la proyección y las dimensiones y áreas entre estas son en la proyección. Carretera PR-106 cerrada en el Km 1.6 por Puente Roto desde comienzos de marzo de 2012, luego de estar en ese lugar cerrado un carril, por mas de 12 meses, al comenzar los trabajos colapsa el carril que servia. Desvio Oficial por PR-351 a la PR-108 PR-106 Km 1.6 CERRADA Puente Colapsado en PR-106 Km.1.6 desde Marzo del 2012-Desvio Oficial PR-351 www.ngs.noaa.gov www.revistatp.com –un “Dominio” al servicio de la comunidad geo-espacial www.gitpr.org Módulo 2 Definiciones de Datums GEODETIC DATUMS Classical HORIZONTAL – 2 D (Latitude and Longitude) (e.g. NAD 27, NAD 83 (1986)) VERTICAL – 1 D (Orthometric Height) (e.g. NGVD 29, NAVD 88, IGLD 85) Contemporary PRACTICAL – 3 D (Latitude, Longitude and Ellipsoid Height) Fixed and Stable – Coordinates seldom change (e.g. NAD 83 (1992)) SCIENTIFIC – 4 D (Latitude, Longitude, Ellipsoid Height, Velocity) – Coordinates change with time (e.g. ITRF00, ITRF05) VERTICAL DATUMS A set of fundamental elevations to which other elevations are referred. Datum Types Tidal – Defined by observation of tidal variations over a specified epoch of time Geodetic – Typically based on Mean Sea Level at one or more tidal stations for a specified epoch of time U.S. VERTICAL DATUMS NAVD 88 Conterminous US Tied to primary CHS BM (60-78) at Father’s Point – Quebec, CD NAVD 88 Alaska GPS Ellipsoid Height + Geoid06 Puerto Rico Vertical Datum 2002 (PRVD 02) Defined by MSL (83-01) at 975 5371 A American Samoa Vertical Datum of 2002 Defined by MSL (83-01) at 177 0000 S Guam Vertical Datum of 2003 Defined by MSL (83-01) at 163 0000 TIDAL 4 Marianas Vertical Datum of 2003 Defined by MSL (83-01) at 163 3227 UH-2C Hawaii Vertical Datums ???? Módulo 3 Base de Datos Geográficos BASE DE DATOS GEOGRAFICOS FUENTE: Robert Liley & Associates Management Consultants Exhibit 33 1- Records de Propiedad a. Información Maestra • Localización Geografica • Dirección Civil • Nombre del dueño • Dirección/teléfono del dueño (si es diferente al anterior) • Numero de residentes • Descripción Legal (including stata title units) • Valor de la ultima transacción de la venta • Fecha de la ultima venta • “Covenants” b. Tasación • Clasificación • Exención contributiva/razones • Avalúo (assessment) • Impuestos actuales • Taza de tasación (TaxRate) • Balance Actual de impuestos • Status • Mejoras Locales • Historia c. Propiedades Gubernamentales • Fecha de adquisición/método de adquisición • Código de localización • Costo original • Valor actual del mercado • Seguros 1. Cobertura 2. Carrier Policy Cross Referente • Hipotecas 1. Cantidad original 2. Balance Pendiente 3. Datos de Interés 4. Términos 5. Tenedor de la Hipoteca • Datos de Alquiler (Rental Data) 1. Renter 2. Covenants 3. Terms 4. Fecha de Pago 5. Renovaciones 6. Pago 7. Balance Corriente 8. Insurance Cross Referente (Third Party) • Datos de Renta 1. Renter 2. Covenants 3. Terms 4. Fecha de Pago 5. Renovaciones 6. Pago 7. Deposito de seguridad 8. Balance Corriente 9. Insurance Cross Referente (Third Party) • Disposición 1. Precio de Venta 2. Costos 3. Dólares que no se venden 4. “Grants” d. Características Físicas • Indicadores de agua en lote • Dimensiones del lote • Huella del edificio • Fecha de la Primera construcción • Fecha de la ultima modificación • Referencias de permisos “cross” • Zonificación actual • Uso de Terreno actual • Indicador de no-conformidad legalmente • Conexión de utilidades • Servidumbres • Conexión de Utilidades extranjeras e. Información de Incendios (FIRE) • Características funcionales del Edificio • Peligros (Hazards) • Necesidades Especiales/Limitaciones • Incidentes con Incendios • Inspección 1. Itinerario 2. Hallazgos 3. Asuntos por hacer (Action Ítems) 4. Fecha de la Ultima Inspección • Perfil del Edificio • Tipo de equipo en expediente/localización f. Información Policiaca • Perfil del Edificio (Planificación Táctica) • querellas/Incidentes • Cursos de Acción • Resoluciones g. Asistencia Social • Educación • PAN • Desempleo • Otras Formas de Asistencia Social 2- Datos Ambientales • Tipo de Suelo • Drenaje • Cuenca Hidrográfica • Hidrológia • Vertedero de sustancias Toxicas • Incidentes Ambientales • Topografía • Recursos Naturales 1. Tipo 2. Localización de Referencia • Identificación de Áreas Inundables • Huracanes • Áreas de Deslizamientos • Estudios de Impacto Ambiental 3- Datos Estadísticos a. Crímenes • Fecha del Evento(s) • Tipo • Localización/Dirección • Historia b. Incendios • Fecha del Evento(s) • Tipo • Localización/Dirección • Historia c. Salud/Ambiente • Fecha del Evento(s) • Tipo • Localización/Dirección • Historia d. Datos Censales/Demográficos • Tracto Censal • Nacionalidad • Numero de Personas • Distribución por Edad • Niveles de Ingresos • Idioma • Nivel de Educación • Status Civil • Sexo • Afiliación Política e. Datos Económicos • Importar/Exportar • Turismo • Tipos de Industrias • Producción de Cemento • Gastos/Ingresos • Gastos Capitales • Salario Bruto Pagado • Indicadores del Componentes de Actividad Económica f. Cuentas Nacionales • Balance de Pagos • GNP • Componentes de la Tasación • Datos de Deposito del Banco • Ingresos Gubernamentales • Gastos Gubernamentales • Actividad del Sector Privado • Datos de Importación/Exportación g. Accidentes • Tipos • Fecha del ultimo evento • Localización/Dirección • Historia Previa 4- Datos de Infraestructura a. Localización Geográfica b. Carreteras • Tipo • Superficie • Fecha de Construcción • Propietario • Fecha del Ultimo Mantenimiento • Número de Carriles • Ancho • Profundidad c. Aceras • Tipo • Superficie d. e. f. g. h. • Fecha de Construcción • Propietario • Fecha del Ultimo Mantenimiento • Ancho/Profundidad Servicios de Agua • Hidrantes • Capacidad • Válvulas • Líneas • Presión/Edad • Profundidad • Fecha de Instalación • Fecha de Ultimo Mantenimiento • Estaciones de Bombeo Servicios de Alcantarillado Sanitarios • Líneas • Válvulas • Profundidad • Estaciones de Bombeo • Tipos Luminarias en Calles • Postes • Transformadores • Alambrado • Edad • Profundidad Servicios Eléctricos • Edad • Estaciones de Potencia • Sub-Estaciones • Transformadores • Embobinado (Wiring) • Postes • Cuadricula (Área de Cobertura) • voltaje • “Wattage” • Profundidad Tuberías (pipelines) • Tipo • Construcción • Edad • Profundidad • Tamaño • Capacidad i. Facilidades de Comunicación • Red/Planta de Teléfonos • Postes • Repetidores • Tamaño del Cable/Capacidad • Tipo de Cable • Fecha de Instalación • Fecha de la Ultima Modificación • Propietario • Facilidades de microondas/satélite • Área de Servicio/mantenimiento j. Servicios de Control de Trafico • Tipo • Activo/Pasivo • Fecha de Instalación • Fecha del Ultimo Mantenimiento k. Facilidades de Parques/Recreación • Tipo • Tamaño/Dimensiones • Fecha de Construcción • Rasgos • Código de Mantenimiento • Capacidad • Características Únicas l. Carriles (lanes) • Localización • Tipo • Superficie • Propietario • Dimensiones Físicas m. Escuelas • Tipo • Número de Salones de Clases • Capacidad • Fecha de Edificación • Características Físicas • Rasgos Especiales • Itinerario de usos n. Inventario de Árboles/Arbustos • Tipo • Numero • Fecha de Sembrado • Fecha de Remoción • Fecha de Transferencia o. Estacionamiento • Capacidad • Medido/Sin Medidor • “Signed” p. Puentes • Tipo • Superficie • Fecha de Construcción • Número de Carriles • Fecha del Ultimo Mantenimiento q. Utilidades extranjeras • Cablevisión • Propietario • Localización de Planta • Tuberías • Propietario • Capacidad • Localización de Planta • Características Especiales • Profundidad r. Red/Facilidades de Carreteras (Viales) • Tipo • Ruta • Rasgos Especiales • Facilidades Módulo 4 Geo-‐Based Informa2on Matrix Geo-Based Information Matrix Mapa Base “Assessments” Facilidades de Comunicaciones Estadísticas criminales Demografía Estadísticas Económicas Plantas Eléctricas Datos Ambientales Datos Financieros Edificios Gubernamentales Peligros (Hazards) Vertederos de Sustancias Peligrosas Hidrológica Datos de Inmigración Datos de Importación/Exportación Parcelas de Terrenos Datos de Usos de Terrenos Licencias Recursos Naturales Datos de Tenencia/Ocupación Permisos Utilidades Privadas Facilidades de Recreación Residencial, Comercial, Industrial Carreteras/Aceras Facilidades Agua/Alcantarillados Datos de Servicios Sociales Suelos Luminarias de Calles Impuestos Topografía Facilidades de Transportación o Viales Drenajes Pluviales y Cuencas Zonificación Lista de Agencias Junta de Planificación Autoridad de Edificios Públicos ARPE Defensa Civil Departamento de Agricultura Administración de Desarrollo Económico Departamento de Comercio Federal Administración de Servicios Generales Banco Gubernamental de Fomento Municipios Guardia Nacional Autoridad de Acueductos y Alcantarillados Autoridad de Comunicaciones Departamento de Comercio Departamento de Recursos Naturales y Ambientales Departamento de Transportación y Obras Públicas Autoridad de Energía Eléctrica Compañía de Fomento Industrial (PRIDCO) Departamento de Salud Compañía de Teléfonos (PRTC) Junta de Calidad Ambiental Departamento de Servicios Sociales o Familia Policía de Puerto Rico Compañía de Turismo Departamento de Haciendas Universidad de Puerto Rico Módulo 5 Ley Núm. 264 del año 2002: Ley para adoptar el Sistema de Coordenadas Planas Estatales de Puerto Rico Ley Núm. 264 del año 2002 (P. del S. 990), 2002, ley 264 Ley para adoptar el Sistema de Coordenadas Planas Estatales de Puerto Rico Ley Núm. 264 aprobada en 16 de noviembre de 2002 Para adoptar el Sistema de Coordenadas Planas Estatales utilizando el método de la Proyección Conforme Cónica Lambert y el North American Datum of 1983, o su versión más reciente, como el sistema de mapa oficial de las agencias, instrumentalidades y municipios del Estado Libre Asociado de Puerto Rico, a los fines de uniformar el uso de proyecciones cartográficas por los organismos gubernamentales, disponer lo relativo a la certificación de ubicaciones o posiciones geográficas, y para otros fines. EXPOSICION DE MOTIVOS La Junta de Planificación del Estado Libre Asociado de Puerto Rico es el organismo gubernamental responsable de guiar el desarrollo integral de Puerto Rico de modo coordinado, adecuado y económico, para atender necesidades actuales y futuras. La Junta también sirve como centro de información para el Gobierno y para la ciudadanía en general. Una de las herramientas de mayor importancia, para la evaluación y estudio de casos que se someten ante la Junta de Planificación, es el Archivo Gráfico. En este Archivo se recopila información gráfica en mapas, relacionada con las consultas que se someten ante la Junta, asi como la información relacionada con mapas de zonificación, zonas inundables y planes de usos de terrenos y expansión urbana. Este recurso de la Junta ofrece un marco de referencia abarcador y preciso, el desarrollo físico y estructural de Puerto Rico. Las agencias gubernamentales, los municipios y las empresas privadas utilizan este banco de datos continuamente para planificar sus proyectos Más eficientemente. Por ello, la Junta de Planificación mecanizó sus sistemas de información a fin de poder curnplir con la continua demanda de información del sector público y privado y para atender la cantidad de casos radicados. Para mantener actualizada esta información, se ha requerido, que toda consulta radicada ante la Junta presente la ubicación exacta del centro del proyecto propuesto, utilizando para ello la proyección de coordenadas geográficas conocida como Coordenadas Lambert. El Servicio Geodésico Nacional ha diseñado para los estados y territorios de los Estados Unidos de América, una red de sistemas de coordenadas conocido como Sistemas de Coordenadas Planas Estatales. El propósito de dichos sistemas es lograr la mejor representación cartográfica de las características geográficas de cada uno de los estados o territorios. Las proyecciones que utilizan dichos Sistemas toman en consideración diferentes criterios. Así por ejemplo, se utiliza la extensión del área cuando se trata de un hemisferio, un estado o un continente; la localización del área, esto es, si está en los polos, latitudes medias o en el Ecuador; o la forma de la extensión predominante, si el área discurre de norte a sur o de este a oeste. La Proyección Conforme Cónica Lambert, con dos paralelos estindares, es la base para el Sistema de Coordenadas Planas Estatales diseñado para zonas de dimensiones limitadas de norte a sur y para extensiones de terrenos indefinidas que corren de este a oeste. Varias agencias e instrumentalidades gubernamentales, como el Departamento de Recursos Naturales y Ambientales, la Autoridad de Carreteras, el Departamento de la Vivienda, la Autoridad de Tierras, la Administración de Terrenos y la propia Junta de Planificación, utilizan como referencia los mapas oficiales del Servicio Geológico de los Estados Unidos. Estos mapas de Puerto Rico están representados y utilizan la Proyección Conforme Cónica Lambert. La adopción en la jurisdicción del Estado Libre Asociado de Puerto Rico de una proyección cartográfica uniforme como el Sistema de Coordenadas Lambert y el North American Datum of 1983, o su versión más reciente, contribuirá a que el intercambio de información geográfica se realice de forma rápida, clara y eficiente. Por ello, conviene al interés público que se establezca por Ley que el Sistema de Coordenadas Planas Estatales con el método de Proyección Cónica Lambert serd el sistema de mapa oficial para toda la Isla, de forma tal que se utilice una sola proyeeción cartográfica para determinar la ubicación de un lugar específico en la jurisdicción de Puerto Rico. También se dispone que el gobierno central adoptará un modelo matemático que represente el esferoide terrestre, conocido como datum, de modo que sea estindar o uniforme para todas las agencias e instrumentalidades públicas. DECRETASE POR LA ASAMBLEA LEGISLATIVA DE PUERTO RICO: Artículo 1.- Adopción del Sistema de Coordenadas Planas Estatales.Se adopta el Sistema de Coordenadas Planas Estatales con el método de Proyección Conforme Cónica Lambert como el sisterna de mapa oficial para todas las agencias, instrumentalidades y municipios del Estado Libre Asociado de Puerto Rico, a los fines de uniformar el uso de proyecciones cartográficas por los organismos gubernamentales. Se dispone, además, que el gobierno central adoptará dentro de un (1) año de la aprobación de esta Ley, el North American Datum of 1983, o su versión más reciente, como el modelo matemático para representar el elipsoide terrestre de modo que sea estándar o uniforme para todas las agencias e instrumentalidades públicas. Artículo 2.- Certificación de Ubicaciones o Posiciones Geográficas.Toda ubicación o posición geográfica, que se presente ante las agencias, instrumentalidades y municipios del Estado Libre Asociado de Puerto Rico deberá cumplir con una de las siguientes alternativas: (a) En el caso de la ubicación de estructuras, edificaciones o construcciones deberá estar certificada por un profesional autorizado, de acuerdo con lo establecido en la Ley Núm. 185 de 26 de diciernbre de 1997, según enmendada, conocida como "Ley de la Junta Examinadora de Ingenieros, Arquitectos, Agrimensores y Arquitectos Paisajistas de Puerto Rico." (b) En el caso de estudios, investigaciones, o evaluaciones de naturaleza técnica o científica, deberá estar acompañada de un documento debidamente juramentado por la persona que lo presenta, en el cual se incluya por lo menos la siguiente información: (i) Equipo de mensura utilizado, o sea si es, un equipo del sistema de posicionamiento global o un equipo convencional, asi como la marca, modelo y precisión máxima del equipo. (ii) Modelo matemático o datum utilizado al momento de hacer la medición. (iii) Tiempo requerido en tomar la medición, expresado en horas:minutos:segundos o en la forma establecida por reglamentación. Todas las agencias, instrumentalidades y municipios del Estado Libre Asociado de Puerto Rico ante las cuales haya que presentarse una certificación sobre la ubicación o posición geográfica de una estructura, propiedad, proyecto, desarrollo, estudio científico o técnico, adoptarán un documento o formulario uniforme para tales propósitos. Mediante dicho documento o formulario uniforme para todas las entidades públicas se dispondrá referente a los métodos de certificación de ubicación o posición establecidos en los incisos (a) y (b) de este Artículo. Este documento o formulario uniforme se adoptará por reglamento dentro de los noventa (90) días contados a partir de la aprobación de esta Ley, luego de consultar y contar con las recomendaciones de la Oficina de Gerencia y Presupuesto, del Departamento de Recursos Naturales y Ambientales, de la Junta de Planificación y de la Administración de Reglamentos de Permisos, del Instituto de Cultura Puertorriqueña, así como de la Junta de Calidad Ambiental. Artículo 3.- Esta Ley comenzará regir inmediatamente después de su aprobación. ___________________ Presidente del Senado ____________________ Presidnete de la Cámara Nota Importante: Esta ley es copia de la ley original cuando fue aprobada, no incluye enmiendas posteriores. Módulo 6 Sección 5.2.3 -‐ Sistema de Coordenadas Planas Estatales Sección 5.2.3 Sistema de Coordenadas Planas Estatales a. A tenor con las disposiciones de la Ley 264 de 16 de noviembre de 2002, se adopta un Sistema de Coordenadas Planas Estatales homogéneo que permita el establecimiento de un sistema cartográfico preciso y confiable. b. Toda solicitud ante la OGPe que implique la creación de nuevas parcelas, que contenga medidas agrarias, que requiera la delimitación para su ubicación precisa, que conlleve cambios en la topografía existente, nuevas construcciones o funciones similares, deberá estar geo-referenciado al Sistema de Coordenadas Oficial del Gobierno de Puerto Rico y contar con la Certificación de un Agrimensor Profesional o autorizado a ejercer la profesión de la agrimensura, la cual debe cumplir con lo siguiente: 1. El Agrimensor certificará el Datum Horizontal y su revisión en o posterior al 1997. Esta debe contener un tracto que demuestre que la determinación procede de un punto de control oficial del National Geodetic Survey (NGS) y evidenciar el método de Global Positioning System (GPS) utilizado para la determinación. 2. Evidenciar el equipo y los monumentos de control establecidos en el proyecto. Los monumentos deben estar representados en el plano con las coordenadas tanto planas como geodésicas. 3. En aquellos proyectos que contengan elevaciones, se especificará el Datum vertical y se evidenciará mediante un tracto que su determinación proviene de un Bench Mark (BM) oficial de las entidades relacionadas de gobierno, como pueden ser: NGS y FEMA. 4. Se requerirá que los planos contengan la información de los Datums horizontal y/o vertical, junto a la certificación del Agrimensor responsable de esa determinación. 5. La posición del polígono debe garantizar un error estándar que no exceda el cero punto cero cinco (±0.05) metros equivalentes a dos (2) pulgadas dentro del Datum utilizado. 6. El Archivo deberá estar en formato DXF y georeferenciado en coordenadas State Plane Nad 83 que incluya parcela, huella de los edificios y/o estructura, aceras y calles. Cada uno de estos deben estar hechos en capas (layers) individuales y en polígonos cerrados. 7. Para aquellos proyectos que no impacten toda la parcela deben entregar Archivo en formato DXF y georeferenciado en coordenadas State Plane NAD83 que incluya parcela, límite del proyecto, huella de los edificios y/o estructura, aceras y calles. Cada uno de estos deben estar hechos en capas (layers) individuales y en polígonos cerrados. Módulo 7 Modern Terrestrial Reference Systems – Part 1 REFERENCE SYSTEMS Modern Terrestrial Reference Systems (Part 1) Dr. Richard A. Snay and Dr. Tomás Soler S urveyors, GIS/LIS professionals, engineers, cartographers, and others who work in North America face the challenge of dealing with at least three different 3D terrestrial reference systems. For many legal activities, these people express positional coordinates in the reference system known as the North American Datum of 1983 (NAD 83). Alternatively, they often favor using the World Geodetic System of 1984 (WGS 84) for various practical positioning activities involving the Global Positioning System (GPS), or they find the International Terrestrial Reference System (ITRS) more suitable for achieving superior positional accuracy. While these three reference systems differ from one another only slightly in concept, they differ significantly in how they have been realized, where the realization of a particular reference system is called a “reference frame”. A particular reference frame is usually established by designating positions and velocities for several identifiable points. To date there have been several realizations of each of these three reference systems, as institutions have systematically revised positions and velocities from time to time to keep pace with how evolving technology has improved positioning accuracy. Here, we review the evolution of these reference frames, and we discuss transforming positions between different reference frames. Finally, we address some practical considerations for accurate positioning and discuss plans for a new NAD 83 realization. Defining a Reference System The modern approach to defining a 3D terrestrial reference system may be divided into four steps. The first step links the axes of a 3D cartesian coordinate system to a configuration of physically measurable locations on or within the earth. As a result, the location and orientation of the three coordinate axes are defined. The second step relates the concept of distance to physically measurable quantities whereby a unit of length is introduced. The third step introduces an auxiliary geometric surface that approximates the size and shape of the earth. Finally, Figure 1 the fourth step addresses the question of how Earth’s gravity field contributes to the notion of position, and especially that of height. We shall be concerned here with only the first three steps, thus focusing on the geometric aspects involved in defining a reference system. For the first step, most scientists involved in defining modern reference systems agree that the origin of the 3D cartesian system should be located at Earth’s center of mass (geocenter); also that the cartesian system’s z-axis should pass through the conventional definition of the North Pole, or more precisely, the International Reference Pole (IRP) as defined by the International Earth Rotation Service (IERS), an international organization established in 1988 and headquartered in Paris, France. The x-axis should go through the point of zero longitude located on the plane of the conventional equator, which is also defined by the IERS. The meridian going through this point is located very close to the meridian of Greenwich although the two are not coincident. The y-axis forms a righthanded coordinate frame with the x- and z-axes. Indeed, each of the three refer- ence systems—NAD 83, WGS 84, and ITRS—has been accordingly defined in concept. They differ, however, as we shall soon discuss, in their realizations; that is, in how the location and orientation of their respective cartesian axes have been physically materialized as well as their respective concepts of distance. Unfortunately, what initially appears to be a simple geometric procedure is complicated by Earth’s dynamic behavior. For example, Earth’s center of mass is moving relative to Earth’s surface. Also, there are variations of Earth’s rotation rate as well as motions of Earth’s rotation axis both with respect to space (precession and nutation) and to Earth’s surface (polar motion). Moreover, points on the earth’s crust are moving relative to one another as a result of plate tectonics, earthquakes, volcanic/magmatic activity, postglacial rebound, people’s extraction of underground fluids, solid Earth tides, ocean loading, and several other geophysical phenomena. Modern terrestrial reference systems, hence, need to account for these motions. One option is to relate the cartesian axes to the locations of selected points measured at a particu- COPYRIGHT PROFESSIONAL SURVEYOR • December 1999 • All Rights Reserved 1 REFERENCE SYSTEMS lar instant of time (epoch). This alternative is generally used when dealing with the motion of the earth’s rotational axis and with the motions associated with plate tectonics. Other types of motion (for example, subsidence) are accounted for by fixing the cartesian axes to some temporal average of the locations for selected points. As we shall see, a fundamental difference among the various reference frames involves how they address the motion associated with plate tectonics. For the second step, scientists concerned with defining up-to-date terrestrial reference systems agree that the unit of length, the meter or “metre”, corresponds to the length of path traveled by light in a vacuum during a time interval of exactly 1/299,792,458 seconds. This solves the problem of relating the concept of distance to a physically measurable quantity in theory, but not in realization. Each of the various reference frames associated with NAD 83, WGS 84, and ITRS relies on a distinct set of measurements that were performed using one or more of several widely different types of instruments and techniques, among the most representative: GPS, electrooptical distance measuring instrumentation, Doppler satellite positioning, very long baseline interferometry (VLBI), and satellite laser ranging (SLR). While each measurement type had been calibrated to fit the definition of a meter as best as possible, the observations, nevertheless, contain uncertainties. Consequently, the “scale” of any particular reference frame is somewhat less than perfect. In particular, when old classical terrestrial frames are compared with modern “space-age” frames, scale errors at the part-per-million (ppm) level may often be detected. Because of recent technological advances in the measurement of time and, consequently, distance, scale differences between modern frames are now approaching the part-per-billion (ppb) level. For the third step, the earth’s surface is approximated in size and shape by the geometric surface that is formed by rotating an ellipse about its smaller axis (Figure 1). The generated surface is termed an “ellipsoid of revolution” or simply ellipsoid. The ellipsoid’s geometric center should be located at the origin of the 3D cartesian system, and its axis of radial symmetry (semi-minor axis) should coincide with the cartesian z-axis of the selected terrestrial reference frame. The size and shape of the rotated ellipse may be completely specified using two parameters: the length of its semi-major axis, usually denoted a, which approximates the distance from the geocenter to a point on the equator (approximately 6,378 km); and the length of the semi-minor axis, de2 noted b, which approximates the distance from the geocenter to the North Pole. The value of b is about 0.3% shorter than a. The fact that a is longer than b is a consequence of the force imparted by Earth’s rotation causing our home planet to bulge outward around its equator. Often, instead of b, the ellipsoid’s flattening, f = (a - b)/a, is used. D i ff e rent ellipsoids have been adopted for the different reference systems. NAD 83 uses the same ellipsoid as the Geodetic Reference System of 1980 which was adopted by the International Association of Geodesy. WGS 84 uses an ellipsoid adopted by the National Imagery and Mapping Agency (NIMA, formerly the Defense Mapping Agency), and ITRS uses an ellip- soid adopted by the IERS. Corresponding values of a and f are presented in the table below. Given values for a and b (or f), people can convert 3D cartesian coordinates—x, y, z— into the curvilinear coordinates—latitude, longitude, and ellipsoidal height—and vice versa. These curvilinear coordinates embody a certain intuitive appeal in specifying locations on and near the earth’s surface, as they relate to our innate sense of the horizontal and vertical dimensions. DR. RICHARD A. SNAY is Manager of the National Continuously Operating Reference Station (CORS) program and a geodesist with the National Geodetic Survey. DR. TOMÁS SOLER, is Chief, Global Positioning System Branch, Spatial Reference Systems Division, National Geodetic Survey. Reference System Semi-major axis, m Flattening, unitless NAD 83 6,378,137.0 1/298.257222101 WGS 84 6,378,137.0 1/298.257223563 ITRS 6,378,136.49 1/298.25645 COPYRIGHT PROFESSIONAL SURVEYOR • December 1999 • All Rights Reserved Módulo 8 Modern Terrestrial Reference Systems – Part 2: The Evolu2on of NAD 83 Modern Terrestrial Reference Systems PART 2: The Evolution of NAD 83 Dr. Richard A. Snay and Dr. Tomás Soler he first realization of NAD 83 was introduced in 1986 by a group of institutions representing the various North American countries to upgrade the previous reference system; that is, the North American Datum of 1927 or NAD 27. In particular, the National Geodetic Survey (NGS) represented the United States, and this federal institution officially refers to the first NAD 83 realization as NAD 83 (1986). For this realization, the group of institutions relied heavily on Doppler satellite observations collected at a few hundred sites to estimate the lo cation of the Earth’s center of mass and the orientation of the 3D cartesian axes. They also relied on these same Doppler observations to provide scale for NAD 83 (1986). More precisely, the group of institutions relied on 3D Doppler-derived positions that had been transformed by: • a translation of 4.5 m along the z-axis • a clockwise rotation of 0.814 arc sec onds about the z-axis • a scale change of -0.6 ppm The Doppler-derived positions were so transformed to make them more consistent with the very long baseline interferometry (VLBI), satellite laser ranging (SLR), and terrestrial azimuth measurements that were available in the early 1980s. While NAD 83 (1986) is 3D in scope, NGS adopted only horizontal coordinates (latitude and longitude) for over 99% of the approximately 250,000 U.S. control points that were involved in defining this reference frame. Unfortunately, this first realization of NAD 83 occurred a few years before GPS technology made the vertical dimension economically accessible. spirit leveling, and not geometric heights relative to an abstract mathematical surface (the ellipsoid), as obtained with GPS. The attitude towards using GPS to measure heights gradually evolved, however, as NGS and other institutions developed improved geoidal models for determining the spatial separation between mean sea level and the ellipsoid. These improvements enabled people to convert ellipsoidal heights into orthometric heights with greater and greater accuracy. Moreover, practitioners can measure heights much more economically with GPS than with spirit leveling. As GPS matured, so did other spaceage geodetic technologies; in particular, SLR and VLBI. Within a few years after 1986, both GPS and SLR measurements had allowed geodesists to locate Earth’s center of mass with a precision of a few centimeters. In doing so, these technologies revealed that the center of mass that was adopted for NAD 83 (1986) is displaced by about 2 m from the true geocenter. Similarly, GPS, SLR, and VLBI revealed that the orientation of the NAD 83 (1986) cartesian axes is misaligned by over 0.03 arc seconds relative to their true orientation, and that the NAD 83 (1986) scale differs by about 0.0871 ppm from the true definition of a meter. These discrepancies caused significant concern as the use of highly accurate GPS measurements proliferated. In particular, starting with Tennessee in 1989, each state–in collaboration with NGS and various other institutions–used GPS technology to establish regional reference frames that were to be consistent with NAD 83. The corresponding networks of GPS control points were originally called High Precision Geodetic Networks (HPGN). Currently, they are referred to as High Accuracy Reference Networks (HARN). This latter name reflects the fact that relative accuracies among HARN control points are better than 1 ppm, whereas relative accuracies among pre-existing control points were nominally only 10 ppm. For defining these regional reference frames, NGS retained the location of the geocenter and the orientation of the 3D cartesian axes which had been derived in 1986 from the transformed Doppler observations. This agency, however, opted to introduce a new scale that would be consistent with the scale of the then cur- GPS Changed Everything Around the same time that NGS adopted NAD 83 (1986), the agency had begun using GPS technology, instead of triangulation and/or trilateration, for horizontal positioning. The fact that GPS technology also provided accurate ellipsoidal heights was somewhat overlooked in the 1980s because surveyors, hydrologists, and other users of vertical positions required orthometric heights relative to mean sea level, as obtained with tide gauges and COPYRIGHT PROFESSIONAL SURVEYOR • February 2000 • All Rights Reserved 1 A P P L I C AT I O N rent ITRS realization which is known as the International Terrestrial Reference Frame of 1989 (ITRF89). As we shall discuss later, the ITRF89 scale was based on a combination of GPS, SLR, VLBI, and lunar-laser-ranging (LLR) measurements. The resulting scale change, equal to -0.0871 ppm, altered existing NAD 83 latitudes and longitudes insignificantly, but it systematically decreased all ellipsoidal heights by about 0.6 m ( = 0.0871•10-6 • R where R is the radius of the Earth). Nevertheless, this change to a more accurate scale facilitated the migration toward using GPS technology for deriving accurate heights. Let us call this second realization NAD 83 (HPGN) or NAD 83 (HARN); but keep in mind that this second realization is actually a collection of regional realizations that were formulated over a period of several years (1989-1997) with each new regional realization being “adjusted” to fit with those that preceded it. Third Realization of NAD 83 Incorporated CORS In late 1994, NGS introduced a third realization of NAD 83 when the agency organized a network of continuously operating reference stations (CORS). Each CORS includes a GPS receiver whose data NGS collects, processes, and disseminates for public use. Surveyors and other professionals can apply CORS data to position points at which other GPS data have been collected with accuracies that approach a few centimeters, both horizontally and vertically. The CORS system started with about a dozen sites in December 1994, and it has grown at a rate of about three sites per month. Figure 2 depicts the recent status of the National CORS network. Positional coordinates of the early CORS sites were first computed in the ITRS realization known as ITRF93. Equivalent NAD 83 coordinates were then computed by applying a Helmert transformation; that is, a transformation of the form xNAD83 = Tx + (1+s) • xITRF + Rz • yITRF - Ry • zITRF (1a) yNAD83 = Ty - Rz • xITRF + (1+s) • yITRF + Rx • zITRF (1b) zNAD83 = Tz + Ry • xITRF - Rx • yITRF + (1+s) • zITRF (1c) where Tx, Ty, and Tz represent three translation along the xaxis, y-axis, and z-axis, respectively, which will bring the origin of the two frames into coincidence. Rx, Ry, and Rz represent three rotations about the x-axis, y-axis, and z-axis, respectively, which, in combination, will bring the three axes of one frame into parallel alignment with their corresponding axes in the other frame. Finally, s represents the difference in scale between the two frames. The values of Tx, Ty, Tz, Rx, Ry, and Rz had been estimated so that the ITRF93 positional coordinates of nine VLBI sites in the United States would transform as best as possible (in a least squares sense) to their adopted NAD 83 (HARN) positional coordinates. The scale difference, s, was set to equal zero. These VLBI sites were used because they had highly accurate positions (cm-level) in both ITRF93 and NAD 83 (HARN). Let us use the label NAD 83 (CORS93) to identify the reference frame obtained by applying this transformation to convert CORS positions from ITFR93 to NAD 83. In the spring of 1996, NGS computed positional coordinates for all the then existing CORS in yet another ITRS realization, known as ITRF94. Similarly, the agency developed a Helmert 2 transformation from ITRF94 to NAD 83 using eight of the same VLBI sites (a VLBI site in California was not used here because of crustal motion concerns). Again, the scale difference was set to equal zero. NGS applied this new transformation to convert ITRF94 coordinates for the CORS sites to a fourth NAD 83 realization which we call NAD 83 (CORS94). Most recently in the fall of 1998, NGS computed positional coordinates for all the then existing CORS in the ITRS realization known as ITRF96. This time, however, NGS collaborated with representatives from Canada’s Geodetic Survey Division to derive a Helmert transformation based on eight VLBI sites in the United States and four VLBI sites in Canada. Again, a scale difference equal to zero was enforced. Other adopted parameters for this transformation are presented in the sidebar box. NGS applied the resulting transformation to convert ITRF96 positional coordinates for the CORS sites to a fifth NAD 83 realization which we call NAD 83 (CORS96). NGS, however, continues to use NAD 83 (CORS94) positions for all CORS sites, except those that have come online since the fall of 1998 and those whose NAD 83 (CORS96) position differs from their corresponding NAD 83 (CORS94) position by more than 2 cm horizontally or 4 cm vertically. In summary, surveyors and others have witnessed five realizations of NAD 83 in the United States. A similar evolution occurred in Canada, but the two countries at least agree on their first and last realizations. The five U.S. realizations are consistent in their choice of origin and orientation; they differ, however, in their choice of scale. While the scale difference between NAD 83 (1986) and NAD 83 (HARN) equals -0.0871 ppm, the scale difference between NAD 83 (HARN) and any NAD 83 (CORSxx) is smaller than 0.005 ppm in magnitude. It should be noted that the NAD 83 (HARN) latitude and/or longitude of a given control point may differ by up to a meter from its corresponding NAD 83 (1986) coordinate. Fortunately, the horizontal discrepancy between the NAD 83 (CORS93) and NAD 83 (HARN) positions for a control point is almost always less than 10 cm, and the horizontal discrepancy between any two NAD 83 (CORSxx) positions for a control point is almost always less than 2 cm. In addition, as NAD 83 has evolved from mostly a horizontal reference system to a full 3D reference system, the number of control points with measured ellipsoidal heights has grown dramatically. Finally, for the present discussion, we intentionally omitted the role of crustal motion on the evolution of NAD 83. We will address this topic in some detail, however, after discussing the evolution of ITRS, as this international reference system provides the global perspective that will help us understand certain concepts associated with crustal motion. DR. RICHARD A. SNAY is Manger of the National Continuously Operating Reference Station (CORS) program and a geodesist with the National Geodetic Survey. DR. TOMÁS SOLER is Chief, Global Positioning System Branch, Spatial Reference Systems Division, National Geodetic Survey. COPYRIGHT PROFESSIONAL SURVEYOR • February 2000 • All Rights Reserved Módulo 9 Modern Terrestrial Reference Systems – Part 3: WGS 84 and ITRS Modern Terrestrial Reference Systems PART 3: WGS 84 and ITRS Dr. Richard A. Snay and Dr. Tomás Soler he Department of Defense (DoD) realization of the International Terrestrial (starting at 0h UTC, 2 January 1994) when developed the WGS 84 reference Reference System (ITRS). NIMA started expressing their derived GPS system to support global activities involvThe original WGS 84 realization essenorbits in this frame. The latest WGS 84 reing mapping, charting, positioning, and tially agrees with NAD 83 (1986). Subsealization, called WGS 84 (G873), is also navigation. More specifically, DoD introquent WGS 84 realizations, however, apbased completely on GPS observations. duced WGS 84 to express satellite orbits; proximate certain ITRS realizations. BeAgain, the letter G reflects this fact, and that is, satellite positions as a function of cause GPS satellites broadcast the predict“873” refers to the GPS week number starttime. Accordingly, WGS 84 is widely used ed WGS 84 orbits, people who use this ing at 0h UTC, 29 September 1996. Alfor “absolute” positioning activities wherebroadcast information for positioning though NIMA started computing GPS orby people assume that satellite orbits are points automatically obtain coordinates bits in this frame on this date, the GPS Opsufficiently accurate to serve as the sole that are consistent with WGS 84. Hence, erational Control Segment did not adopt source of control for positioning points of the popularity of using GPS for real-time WGS 84 (G873) until 29 January 1997. interest. In particular, absolute positioning positioning has promoted greater use of The origin, orientation, and scale of does not rely on using positional coordiWGS 84. Despite its popularity, people WGS 84 (G873) are determined relative to nates for pre-existing terrestrial points for generally do not use WGS 84 for high-preadopted positional coordinates for 15 GPS control, except indirectly in that orbits are cision positioning activities, because such tracking stations: five of them are mainderived from adopted positions for a small activities require the use of highly accurate tained by the Air Force and ten by NIMA set of tracking stations (Fig. 3). The generpositions on pre-existing terrestrial points (see Fig. 3). NIMA chose their sites to comal user, however, never needs to know the for control. For example, various differenplement the somewhat equatorial distribupositions of these tracking stations. tial GPS techniques use known positions tion of the Air Force sites and to optimize DoD provides both “predicted” and for one or more pre-existing terrestrial multiple station visibility from each GPS “postfit” orbits in the WGS 84 reference points to remove certain systematic errors satellite. People may anticipate further imsystem. As implied by the name, predicted in computing highly precise positions for provements of WGS 84 in the future, as orbits are calculated ahead of time by apnew points. Consequently, before WGS 84 new GPS tracking sites may be added or plying physical principles to extrapolate can support high-precision positioning acexisting antennas may be relocated or recurrently observed satellite positions. On tivities, a rather extensive network of acplaced. NIMA is dedicated to take approthe other hand, postfit orbits are calculated curately positioned WGS 84 terrestrial conpriate measures to guarantee the highest from previously observed satellite positrol points would have to be established. possible degree of quality and to perpetutions. Postfit orbits are more precise than DoD established the original WGS 84 ate the accuracy of WGS 84. As mentioned predicted orbits both because they do not reference frame in 1987 using Doppler obearlier, however, most regions lack a netinvolve predicting the future and because servations from the Navy Navigation Satelwork of accessible reference points that they are usually derived using a larger lite System (NNSS) or TRANSIT. The WGS might serve as control points from which number of tracking stations. GPS predicted 84 frames have evolved significantly since highly accurate WGS 84 coordinates may orbits and satellite clock parameters are the mid-1980s. In 1994, DoD introduced a be propagated using an appropriate static generated by the Air Force at the GPS Oprealization of WGS 84 that is based comdifferential GPS technique involving carrierational Control Segment, located at pletely on GPS observations, instead of er phase observables. Another minor Schriever AFB, Colorado. The Air Force Doppler observations. This new realization drawback affecting accurate GPS work is then uploads these predicted quantities to is officially known as WGS 84 (G730) the unavailability to the general GPS user the GPS satellites so that this information where the letter G stands for “GPS” and of the crustal velocities at the WGS 84 may be included in the radio signal trans“730” denotes the GPS week number tracking stations. More information about mitted by these satellites. These predicted WGS 84 may be obtained via the Internet orbits support all real-time positionby accessing: http:// ing and navigation activities in164.214.2.59/GandG/tr8350_2.html volving GPS. Postfit GPS orbits and satellite clock parameters The Evolution of ITRS are generated by the National In the late 1980s, the InternaImagery and Mapping tional Earth Rotation Service Agency (NIMA), who current(IERS) introduced ITRS to suply makes this information port those scientific activities that available on its Geodesy and require highly accurate positionGeophysics World Wide Web al coordinates; for example, pages. A number of other ormonitoring crustal motion and ganizations also generate the motion of Earth’s rotational postfit GPS orbits which they axis. The initial ITRS realization Figure 3. Combined DoD tracking network that defines WGS 84. usually express in a particular was called the International Ter- T COPYRIGHT PROFESSIONAL SURVEYOR • March 2000 • All Rights Reserved 1 A P P L I C AT I O N restrial Reference Frame of 1988 (ITRF88). Accordingly, IERS published positions and velocities for a worldwide network of several hundred stations. The IERS, with the help of several cooperating institutions, derived these positions and velocities using various highly precise geodetic techniques including GPS, VLBI, SLR, LLR, and DORIS (Doppler orbitography and radiopositioning integrated by satellite). Every year or so since introducing ITRF88, the IERS has developed a new ITRS realization –ITRF89, ITRF90, ..., ITRF97– whereby they have published revised positions and velocities for previously existing sites, as well as new positions and velocities for those sites that had been established since after earlier realizations had been developed. Each new realization not only incorporated at least an additional year of data, but also the most current understanding of Earth’s dynamic behavior. The ITRF96 frame is defined by the positions and velocities of 508 stations dispersed among 290 globally distributed sites (Fig. 4). Recall that a particular site may involve one or more co-located instruments employing various space-related techniques (e.g., GPS, VLBI, SLR, LLR, and DORIS). The accuracy and rigor of ITRS has proven contagious, and its popularity is steadily growing among those who engage in positioning activities. Furthermore, ITRS is the first major international reference system to directly address plate tectonics and other forms of crustal motion by publishing velocities as well as positions for its control points. To appreciate the need for velocities, consider the theory of plate tectonics. According to this theory, Earth’s outer shell consists of about 20 plates that are essentially rigid, and these plates move mostly laterally relative to one another like several large sheets of ice on a body of water. The relative motion between points on different plates are, in some cases, as large as 150 mm/yr, which is easily detectable using GPS and other modern day positioning techniques. Given the fact that each tectonic plate is moving relative to the others, one may ask how crustal velocities may be expressed in “absolute” terms. The people responsible for ITRS currently address this 2 Figure 4. Sites defining ITRF96. dilemma by assuming that the Earth’s surface, as a whole, does not move “on average” relative to Earth’s interior. Said differently, the ITRS developers assume that the total angular momentum of Earth’s outer shell is zero. Hence, the angular momentum associated with the motion of any one plate is compensated by the combined angular momentum associated with the motions of the remaining plates. Consequently, points on the North American plate generally move horizontally at measurable rates according to the ITRS definition of absolute motion. In particular, horizontal ITRF96 velocities have magnitudes between 10 and 20 mm/yr in the coterminous 48 states. Moreover, horizontal ITRF96 velocities have even greater magnitudes in Alaska and Hawaii. In contrast, the NAD 83 reference system addresses plate motion under the assumption that the North American plate, as a whole, does not move “on average” relative to Earth’s interior. Hence, points on the North American plate generally have no horizontal velocity relative to NAD 83 unless they are located near the plate’s margin (California, Oregon, Washington, and Alaska) and/or they are affected by some other deformational process (volcanic/magmatic activity, postglacial rebound, etc.). The NAD 83 reference system, however, does make special accommodations for certain U.S. regions that are located completely on another plate. In Hawaii, for example, NAD 83 positional coordinates are defined as if the Pacific plate is not moving. This approach is convenient for people who are involved with positioning activities solely in Hawaii. This approach, however, introduces a layer of complexity for people who are involved in positioning points in Hawaii relative to points in North America. In the realm of crustal motion, it is inappropriate to specify positional coordinates without specifying the “epoch date” for these coordinates; that is, the date to which these coordinates correspond. Accordingly, ITRF96 positions are usually specified for the epoch date of 1 January 1997 (often denoted in units of years as 1997.0). To obtain positions for another time, t, people need to apply the formula x(t) = x(1997.0) + vx • ( t - 1997.0 ) (2) and similar formulas for y(t) and z(t). Here, x(t) denotes the point’s x-coordinate at time t, x(1997.0) denotes the point’s xcoordinate on 1 January 1997, and vx denotes the x-component of the point’s velocity. NGS furnishes easy access to the ITRF96 reference frame through a set of over 170 stations belonging to the National CORS network (recall Fig. 2). Positions, velocities, and other pertinent information for these stations are available via the Internet by accessing: ftp://www.ngs.noaa. gov/cors/coord/coord_96. Transforming Between Reference Frames In 1998, U.S. and Canadian officials jointly adopted a Helmert transformation to convert positional coordinates between ITRF96 and NAD 83 (CORS96). The IERS has also adopted appropriate Helmert transformations for converting between ITRF96 and other ITRS realizations. NGS has encoded all these transformations into a software package, called HTDP (Horizontal Time-Dependent Positioning), which is freely available via the Internet: http://www.ngs.noaa.gov/TOOLS/ Htdp/Htdp.html This software enables people to transform individual positions entered interactively or a collection of positions entered as a formatted file. Also, if people expect to transform only a few positions, then they may run HTDP interactively from this web page. While Helmert transformations, as encoded into HTDP, are appropriate for transforming positions between any two ITRS realizations or between any ITRS realization and NAD 83 (CORS96), more complicated transformations are required COPYRIGHT PROFESSIONAL SURVEYOR • March 2000 • All Rights Reserved A P P L I C AT I O N for conversions that involve NAD 27, NAD 83 (1986), or NAD 83 (HARN). These complications arise because these frames contain large local and regional distortions that can not be quantified by a simple Helmert transformation. For instance, NAD 27 contains distortions at the 10 m level. That is, if one applied the best possible Helmert transformation from NAD 27 to NAD 83 (CORS96), then the converted NAD 27 positions may still be in error by as much as 10 m. In a similar manner, NAD 83 (1986) contains distortions at the 1 m level, and NAD 83 (HARN) contains distortions at the 0.1 m level. NGS has developed a software package, called NADCON (), that embodies rather intricate transformations to convert positional coordinates between any pair of the following reference frames: NAD 27, NAD 83 (1986), and NAD 83 (HARN). Referring to a pair of 2D grids that span the United States, NADCON contains appropriate values for each grid node to transform its positional coordinates from one reference frame to another. Furthermore, NADCON interpolates these gridded values to transform points located within the grid’s span. It should be noted that NADCON may be used only to transform horizontal coordinates (latitude and longitude), because ellipsoidal heights— relative to NAD 27 or NAD 83 (1986)— have never been adopted for most control points. While HTDP may be used with pairs of certain reference frames (NAD 83 (CORS96), ITRF88, ITRF89,. . . , and ITRF97) and NADCON with pairs of other reference frames (NAD 27, NAD 83 (1986), and NAD 83 (HARN)), no NGSsanctioned software exists for transforming coordinates from any member of one set to any member of the other. Also, no NGS-sanctioned software exists for transforming NAD 83 (CORS93) and/or NAD 83 (CORS94) positions to other reference frames. Regarding the WGS 84 reference system, it is generally assumed that WGS 84 (original) is identical to NAD 83 (1986), that WGS 84 (G730) is identical to ITRF92, and that WGS 84 (G873) is identical to ITRF96. Other transformations between pairs of the WGS 84 realizations, however, have also appeared in the literature. DR. RICHARD A. SNAY is Manager of the National Continuously Operating Reference Station (CORS) program and a geodesist with the National Geodetic Survey. DR. TOMÁS SOLER is Chief, Global Positioning System Branch, Spatial Reference Systems Division, National Geodetic Survey. COPYRIGHT PROFESSIONAL SURVEYOR • March 2000 • All Rights Reserved 3 Módulo 10 Modern Terrestrial Reference Systems – Part 4: Prac2cal Considera2ons for Accurate Posi2oning Modern Terrestrial Reference Systems—Part 4: Practical Considerations for Accurate Positioning Dr. Richard A. Snay and Dr. Tomás Soler P rior to the space-age, people primarily used terrestrial techniques like triangulation and/or trilateration to position points. In those days, they had essentially only one choice for a nationwide reference frame; namely, NAD 27. To position new points, a person would measure distances and angles that related the unknown positional coordinates of the new points to the known coordinates of some pre-existing terrestrial control points. More often than not, these known coordinates were referred to NAD 27. Hence, calculated positions for the new points were automatically referred to NAD 27. In the 1970s, people started using Earth-orbiting satellites for positioning points. In particular, before GPS became operational, a number of people used the Doppler satellite system, known as TRANSIT, to position points with an accuracy of several meters relative to the geocenter. Now people use GPS to posi- tion points even more accurately. Hence, satellites have become flying control points whose time-dependent positions are well known thanks to diligent work by several institutions. Currently, GPS satellite orbits are available in appropriate realizations of both WGS 84 and ITRS, but not in any NAD 83 realization. Confusion Can Be Introduced Techniques that position new points by using only satellites as control automatically yield positional coordinates that are referred to the same reference frame as that used for the orbits, either some realization of WGS 84 or some realization of ITRS. Satellite-only positioning techniques, however, do not currently enable the centimeter-level positioning accuracy that is possible with techniques, like the static differential carrier phase GPS technique, that involve using both satellite and pre-existing terrestrial points as control. The use of an integrat- Figure 5: GPS stations defining the National Spatial Reference System: NAD 83 (NSRS) ed approach, however, introduces confusion when the satellite orbits are referred to one reference frame and the terrestrial control points are referred to another, as the newly computed positions will be referred to some ambiguous hybrid reference frame. Although this hybrid reference frame will approximate the reference frame of the terrestrial control points, the error of this approximation at a given point depends on many factors including the distances and directions from this point to the various terrestrial control points being used as well as the magnitude and orientation of the vector connecting the origins of the two reference frames. For the case of a single terrestrial control point referred to NAD 83 with the GPS orbits referred to ITRF96, the error of this approximation at a point will grow less than 1 cm for every 100 km of distance to the terrestrial control point. Hence, for accurate positioning, people should use the same reference frame for both the orbits and the terrestrial control points, especially if the distance to the nearest terrestrial control point is large. If for some reason the available orbits are expressed in one reference frame while the terrestrial control is expressed in another, then the positions of the terrestrial control points should ideally be transformed into the given reference frame of the orbits, as may be accomplished, for example, by using the HTDP software. Alternatively, one could transform the orbits to the reference frame of the terrestrial control points, but this option is likely to involve more computations. Moreover, the software being used to process GPS data may require the orbits to be expressed in a realization of a specific reference system. Thus, it is usually better to transform the terrestrial control points rather than to transform the orbits. Now, once the orbits and the terrestrial control are expressed in a common reference frame, the computed coordi- COPYRIGHT PROFESSIONAL SURVEYOR • April 2000 • All Rights Reserved 1 REFERENCE SYSTEMS nates of the newly positioned points will also be referred to this reference frame. Hence, if people need these positional coordinates to be expressed relative to some other reference frame, then they will need to transform these coordinates. Again, the HTDP software will serve this purpose. For rigorous computations, the epoch date of the positional coordinates of the terrestrial control should agree with the date that the GPS observations were performed. This agreement is especially important when computations are performed in an ITRS realization because terrestrial points are moving relative to this reference system. As previously mentioned, ITRF96 velocities have magnitudes ranging between 10 and 20 mm/yr in the coterminous United States, and these magnitudes are even higher in Alaska and Hawaii. If velocities at the control points are known, then equation (2) may be applied to convert the positional coordinates of these points from their given epoch date to the date that the GPS observations were performed. If velocities of these control points are unknown, then these velocities may be predicted using the HTDP software. sites have been positioned with great accuracy, as they form the foundation of the National Spatial Reference System. Each site contains a geodetic-quality GPS receiver whose carrier phase and code range measurements are transmitted to NGS headquarters in Maryland. NGS then makes these data freely available via the Internet. In particular, people may download CORS data either via Use of Precise Orbits Recommended As previously implied, we are recommending the use of highly precise postfit GPS orbits expressed in ITRS coordinates for accurate GPS positioning. The International GPS Service (IGS) freely distributes such orbits through their website http://igscb.jpl.nasa.gov/. The IGS orbits are the result of computational efforts involving seven different institutions located in four different countries: Canada, Germany, Switzerland, and the United States. Three of these seven institutions are U.S. based: the Jet Propulsion Laboratory, Scripps Institute of Oceanography, and NGS. Each of the seven institutions computes GPS orbits essentially independently, and then IGS rigorously combines the seven solutions to produce the orbits that they distribute. Also, the NGS-derived postfit orbits are available at www.ngs.noaa.gov/GPS/GPS.html Finally for accurate positioning, sites from the National CORS system provide good terrestrial control. These CORS 2 COPYRIGHT PROFESSIONAL SURVEYOR • April 2000 • All Rights Reserved REFERENCE SYSTEMS anonymous FTP (file transfer protocol) using the address ftp://www.ngs. noaa.gov/cors/ or through the World Wide Web using the address http://www.ngs.noaa.gov/CORS/. Furthermore, people may pose CORS-related questions to NGS via the email address [email protected]. Whether using CORS or other sites as terrestrial control, it is usually best to include several points for the sake of redundancy and to suppress those systematic errors that grow as a function of distance. The Species of "Origins" In 1859, Charles Darwin published his revolutionary book, The Origin of Species, in which he documented evidence for the theory of how various forms of life evolve to better adapt to their environment. Inasmuch as this theory of evolution applies to humans, it may also apply to the products that humans develop. In particular, the theory of evolution may apply in large measure to terrestrial reference systems, as these human artifacts evolve to adapt to our knowledge of our environment and to become commensurate with our ability to position points with ever increasing accuracy. In this series of articles, we have discussed how three particular species of reference systems—NAD 83, WGS 84, and ITRS—have evolved over the past couple of decades. Surely, they will continue to evolve into the future. Indeed, U.S. and Canadian representatives are already planning to introduce a new realization of NAD 83. The U.S. geospatial community will come to know this new realization as NAD 83 (NSRS) where NSRS is an acronym for National Spatial Reference System. NAD 83 (NSRS) will serve to supercede the conglomeration of regional reference frames that comprise NAD 83 (HARN). In particular, NAD 83 (NSRS) will incorporate newly computed positional coordinates for all HARN control points, as well as other control points, so that these new positions will be consistent with adopted CORS positions (Figure 5). Currently, horizontal discrepancies as large as 7 cm exist between NAD 83 (HARN) positions of control points and their idealized NAD 83 (CORS96) positions. Significantly greater discrepancies exist in the ver- tical dimension, because the accuracy of ellipsoidal heights measured during the earlier HARN surveys compares poorly relative to today’s height-measuring capability. Consequently to obtain more accurate heights, NGS is cooperating with several organizations to resurvey the HARN in each State. These HARN resurveys reflect NGS’s Height Modernization Initiative to support the public’s growing use of GPS to measure accurate heights. NGS expects to complete the HARN resurveys around 2002. Moreover, NGS will wait until these resurveys have been completed before releasing NAD 83 (NSRS), as this agency plans to combine these new observations rigorously with existing observations to compute accurate ellipsoidal heights that are con- sistent across the country. NGS is also developing more and more accurate geoid models to convert such ellipsoidal heights to appropriate orthometric heights. Additional details about this forthcoming NAD 83 realization may be found at the Internet address: http://www.ngs.noaa.gov/initiatives/new _reference.shtml. DR. RICHARD A. SNAY is Manager of the National Continuously Operating Reference Station (CORS) program and a geodesist with the National Geodetic Survey. DR. TOMÁS SOLER is Chief, Global Positioning System Branch, Spatial Reference Systems Division, National Geodetic Survey. Transformation Parameters from ITRF96 to NAD83 (CORS96) In the summer of 1998, NGS and Canada’s Geodetic Survey Division adopted the following parameters for transforming positions in North America from ITRF96 to NAD 83 (CORS96) via equation (1) : For the three shifts: Tx = 0.9910 m Ty = -1.9072 m Tz = -0.5129 m. For the scale factor: s = 0.0 (unitless). For the three rotations (counterclockwise sense positive): Rx = [25.79 + 0.0532 • (t - 1997.0)]• mr radians Ry = [ 9.65 - 0.7423 • (t -1997.0)]• mr radians Rz = [11.66 - 0.0316 • (t - 1997.0)]• mr radians. In the above equations, mr = 4848.13681• (10**-12) is the conversion factor from milli-arcseconds to radians; t is the “epoch date” in years (e.g., 1999.3096 0h UTC, 23 April 1999). Recall that the epoch date is the date to which the given positions correspond. Note that the rotations, Rx, Ry, and Rz, are time dependent because the North American tectonic plate moves relative to ITRF96 while this plate is essentially stable relative to NAD 83 (CORS96). Because ITRF94 and ITRF97 have the same origin, orientation, and scale as ITRF96, these same parameters are applicable for transforming these other two ITRS realizations to NAD83 (CORS96). On the other hand, these parameters are not applicable for transforming ITRF93 and earlier ITRS realizations to NAD 83 (CORS96). COPYRIGHT PROFESSIONAL SURVEYOR • April 2000 • All Rights Reserved 3 Módulo 11 Notas sobre Agrimensura
© Copyright 2024