../../_images/logo_proj.png ../../_images/OSGeo_project.png

Inicio Rápido de PROJ

PROJ es una función de filtro estándar Unix que convierte coordenadas geográficas latitud y longitud en coordenadas cartesianas (y viceversa), y es una C API para los desarrolladores de software incluyen la transformación de coordenadas en su propio software.

Esta guía rápida describe cómo utilizar estas herramientas:

PROJ

Nota

¿Qué voy a aprender? en esta sección, usted aprenderá cómo utilizar la herramienta de PROJ.

¿Qué es PROJ?

PROJ y invproj realizan respectiva transformación hacia adelante e inversa de datos cartográficos o de datos cartesianos con una amplia gama de funciones de proyección seleccionable.

Sinopsis de proj

proj [ -bcCeEfiIlmorsStTvVwW [ args ] ] [ +args ] file[s]
invproj [ -bcCeEfiIlmorsStTwW [ args ] ] [ +args ] file[s]

Estos comandos son descritos completamente en la documentación de PROJ para la última versión.

Ejemplo de PROJ

El siguiente script

proj +proj=utm  +lon_0=112w  +ellps=clrk66
-r <<EOF
45d15'33.1"   111.5W
45d15.551666667N   -111d30
+45.25919444444    111d30'000w
EOF

realizará una proyección UTM hacia adelante con un meridiano central UTM estándar más cercano a la longitud 112W. Los valores geográficos de este ejemplo son equivalentes y significan ejemplos de varias formas de entrada de DMS. Los datos de salida x-y aparecerán como tres líneas de:

460769.27     5011648.45

cs2cs

Nota

¿Qué voy a aprender? en esta sección, usted aprenderá cómo utilizar la herramienta de cs2cs.

¿Qué es cs2cs?

cs2cs realiza la transformación entre el sistema de coordenadas cartográfico fuente y destino en un conjunto de puntos de entrada. La transformación del sistema de coordenadas puede incluir traducción entre coordenadas geográficas y proyectadas, así como la aplicación de cambios de referencia.

Sinopsis de cs2cs

cs2cs [ -eEfIlrstvwW [ args ] ] [ +opts[=arg] ] [ +to [+opts[=arg]] ] file[s]

Estos comandos son descritos completamente en la documentación de cs2cs para la última versión.

Ejemplo de cs2cs

El siguiente script

cs2cs +proj=latlong +datum=NAD83
      +to +proj=utm +zone=10  +datum=NAD27
-r <<EOF
45d15'33.1"   111.5W
45d15.551666667N   -111d30
+45.25919444444    111d30'000w
EOF

va a transformar las coordenadas geográficas de NAD83 entradas en NAD27 coordenadas en la proyección UTM zona 10 seleccionados. Los valores geográficos de este ejemplo son equivalentes y significados como ejemplos de diversas formas de entrada DMS. Los datos de salida x y aparecerán como tres líneas de:

1402285.99      5076292.42 0.000

geod

Nota

¿Qué voy a aprender? en esta sección, usted aprenderá cómo utilizar la herramienta de geod.

¿Qué es geod?

geod (directa) y invgeod (inverso) realizan cálculos geodésicos («Gran Círculo») para determinar la latitud, longitud y azimut inverso de un punto de término dado una inicial punto de latitud, longitud, Azimut y distancia (directas) o el delantero y la espalda Acimutes y distancia entre inicial y terminal de punto de latitudes y longitudes (inverso). Los resultados son precisos para culminar para |f| < 1/50, donde f es aplanar.

invgeod puede no estar disponible en todas las plataformas; en este caso convoque geod con la opción -I.

Sinopsis de geod

geod +ellps=<ellipse> [ -afFIlptwW [ args ] ] [ +args ] file[s]

invgeod +ellps=<ellipse> [ -afFIlptwW [ args  ]  ]  [  +args  ] file[s]

Estos comandos son descritos completamente en la documentación de PROJ para la última versión.

Ejemplo de Geod

La siguiente secuencia de comandos determina los azimuts geodésicos y la distancia en millas estatutarias de EE. UU. Desde Boston, MA, a Portland, OR:

geod +ellps=clrk66 <<EOF -I +units=us-mi
42d15'N 71d07'W 45d31'N 123d41'W
EOF

que da los resultados:

-66d31'50.141" 75d39'13.083" 2587.504

donde los dos primeros valores son el azimut de Boston a Portland, el azimut inverso de Portland a Boston seguida de la distancia.

Un ejemplo del uso geodésica hacia adelante es usar la loca-ción de Boston y determinar la ubicación de Portland por azimut y distancia:

geod +ellps=clrk66 <<EOF +units=us-mi
42d15'N 71d07'W -66d31'50.141" 2587.504
EOF

que da:

45d31'0.003"N 123d40'59.985"W 75d39'13.094"

Nota

Falta de precisión en el valor de distancia compromete la precisión de la localización de Portland.

Lectura adicional

  1. GeographicLib

  2. C. F. F. Karney, Algorithms for Geodesics, J. Geodesy 87, 43-55 (2013). Apéndice

  3. La bibliografía en línea de geodésica