mapprojectgmt - Online en la nube

PROGRAMA:

NOMBRE

mapproject: realiza transformaciones de mapas hacia adelante e inversas, conversiones de datum y geodesia

SINOPSIS

maproject [ tablas ] parámetros región [ b|B|f|F|o|O[can0/lat0]] [[dx/dy]] [ c|i|p ] [
[fecha]] [[unidad]] [[x0/y0/] [[+|-]unidad][+|-]] [] [ línea.xy[/ [+|-]unidad][+]] [
[a|c|g|m]] [[d|e ] [] [[h]en[/a]] [[nivel]] [[w|h]] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o]]

Nota: No se permite ningún espacio entre el indicador de opción y los argumentos asociados.

DESCRIPCIÓN

maproject lee (longitud, latitud) posiciones de tablas [o entrada estándar] y
calcula las coordenadas (x, y) utilizando la proyección del mapa y las escalas especificadas. Opcionalmente,
puede leer posiciones (x, y) y calcular valores (longitud, latitud) haciendo lo inverso
transformación. Esto se puede utilizar para transformar puntos lineales (x, y) obtenidos al digitalizar
un mapa de proyección conocida a coordenadas geográficas. También puede calcular distancias a lo largo
pista, a un punto fijo, o la aproximación más cercana a una línea. Finalmente, se puede utilizar para realizar
varias conversiones de datos. Se permiten campos de datos adicionales después de las 2 primeras columnas
que debe tener (longitud, latitud) o (x, y). Ver opción -: sobre como leer
(latitud, longitud) archivos.

REQUERIDO ARGUMENTOS

-Jparámetros (más ...)
Seleccione la proyección del mapa.

-R [unidad]xmín/xmax/ymin/ymax[r] (más ...)
Especifique la región de interés. Caso especial para la proyección UTM: Si -C se utiliza
y -R no se da, entonces la región se establece para que coincida con la zona UTM dada, por lo que
como para preservar la solución elipsoidal completa (Ver RESTRICCIONES para más
información).

OPCIONAL ARGUMENTOS

mesa Uno o más ASCII (o binarios, consulte -bi[ncoles][tipo]) archivo (s) de tabla de datos con un
número de columnas de datos. Si no se dan tablas, leemos de la entrada estándar.

-Ab | B | f | F | o | O [can0/lat0]
-Af calcula el acimut (hacia adelante) desde un punto fijo lon / lat a cada punto de datos.
Uso -Lejos para obtener el azimut de retorno de los puntos de datos al punto fijo. Usar -Ao para obtener
orientaciones (-90/90) en lugar de acimutes (0/360). Mayúsculas F, B or O seguirá
convertir de latitudes geodésicas a geocéntricas y estimar el azimut de las geodésicas
(asumiendo que el elipsoide actual no es una esfera). Si no se da un punto fijo, entonces
calculamos el azimut (o back-azimut) desde el punto anterior.

-C[dx/dy]
Establecer el centro de las coordenadas proyectadas para que esté en el centro de proyección del mapa [El valor predeterminado es
esquina inferior izquierda]. Opcionalmente, agregue compensaciones en las unidades proyectadas que se agregarán (o
restado cuando -I está establecido) a (desde) las coordenadas proyectadas, como falso
Este y norte para zonas de proyección particulares [0/0]. La unidad utilizada para el
compensaciones es la unidad de distancia de la parcela en efecto (ver PROJ_LENGTH_UNIT) a menos que -F is
utilizado, en cuyo caso las compensaciones están en metros.

-Dc | i | p
Anula temporalmente PROJ_LENGTH_UNIT y usa c (cm), i (pulgada), o p (puntos)
en lugar de. No se puede utilizar con -F.

-MI[fecha]
Convertir de geodésico (largo, latitud, altura) a Tierra Centrada en la Tierra Fija (ECEF)
(x, y, z) coordenadas (agregar -I para la conversión inversa). Anexar ID de referencia (ver -Qd)
o dar elipsoide:dx,dy,dz donde elipsoide puede ser un ID de elipsoide (ver -Qe) o
dado como a[, * inv_f *], donde a es el semi-eje mayor y inv_f es la inversa
aplanamiento (0 si se omite). Si fecha es - o no dado, asumimos WGS-84.

-F[unidad]
Forzar escala 1: 1, es decir, salida (o entrada, ver -I) los datos están en la proyección real
metros. Para especificar otras unidades, agregue la unidad deseada (consulte UNIDADES). Sin -F,
la salida (o entrada, ver -I) están en las unidades especificadas por PROJ_LENGTH_UNIT (pero
ver -D).

-GRAMO[x0/y0/] [[+ | -]unidad] [+ | -]
Calcular distancias a lo largo de la pista or al punto opcional establecido con -Gx0 / y0. Adjuntar
la unidad de distancia (ver UNIDADES), incluyendo c (Distancia cartesiana usando entrada
coordenadas) o C (Distancia cartesiana utilizando coordenadas proyectadas). los C unidad
requiere -R y -J para ajustar. Sin un punto fijo, calculamos acumular
distancias a lo largo de la vía. Adjuntar - para obtener una distancia incremental entre sucesivos
puntos. Adjuntar + para especificar el segundo punto a través de dos columnas adicionales en el archivo de entrada.

-I Realice la transformación inversa, es decir, obtenga (longitud, latitud) de (x, y) datos.

-Llínea.xy[/ [+ | -]unidad] [+]
Determine la distancia más corta desde los puntos de datos de entrada a las líneas dadas en
el archivo ASCII multisegmento línea.xy. La distancia y las coordenadas del
el punto más cercano se agregará a la salida como tres columnas nuevas. Anexar el
unidad de distancia (ver UNIDADES), incluyendo c (Distancia cartesiana usando coordenadas de entrada)
or C (Distancia cartesiana utilizando coordenadas proyectadas). los C la unidad requiere -R y
-J para ajustar. Finalmente, agregue + para informar la identificación del segmento de línea y la fracción
número de punto en lugar de lon / lat del punto más cercano.

-N [a | c | g | m]
Convertir de latitudes geodésicas (usando el elipsoide actual; ver PROJ_ELLIPSOID)
a una de las cuatro latitudes auxiliares diferentes (las longitudes no se ven afectadas). Escoger
en autálico conformal geocéntrico, y mlatitudes eridionales [geocéntricas]. Usar -I
convertir de latitudes auxiliares a latitudes geodésicas.

-Q [d | e Enumere todos los parámetros de proyección. Para enumerar solo datums, utilice -Qd. Solo enumerar
elipsoides, use -Qe.

-S Suprime los puntos que quedan fuera de la región.

-T [h]en[/a]
Coordinar conversiones entre datums en y a utilizando el estándar Molodensky
transformación. Usar -Désimo si la tercera columna de entrada tiene una altura por encima del elipsoide [predeterminado
asume altura = 0, es decir, en el elipsoide]. Especifique los datums utilizando el ID de datum (consulte
-Qd) o dar elipsoide:dx,dy,dz donde elipsoide puede ser un ID de elipsoide (ver -Qe) o
dado como a[, * inv_f *], donde a es el semi-eje mayor y inv_f es la inversa
aplanamiento (0 si se omite). Si fecha es - o no dado, asumimos WGS-84. -T puede ser
utilizado junto con -R -J para cambiar el datum antes de la proyección de coordenadas
(añadir -I para aplicar la conversión de datum después de la proyección inversa). Asegúrate de eso
la configuración de PROJ_ELLIPSOID es correcta para su caso.

-V [nivel] (más ...)
Seleccione el nivel de verbosidad [c].

-W [w | h]
Imprime el ancho y la altura del mapa en la salida estándar. No se leen archivos de entrada. Solo para
generar el ancho o el alto, agregar w or h, respectivamente. Las unidades del
las dimensiones se pueden cambiar a través de -D.

-bi[ncoles] [t] (más ...)
Seleccione la entrada binaria nativa. [El valor predeterminado es 2 columnas de entrada].

-bo [ncoles][tipo] (más ...)
Seleccione la salida binaria nativa. [El valor predeterminado es el mismo que el de la entrada].

-d [i | o]sin datos (más ...)
Reemplazar columnas de entrada que sean iguales sin datos con NaN y haga lo contrario en la salida.

-f [i | o]colinfo (más ...)
Especifique los tipos de datos de las columnas de entrada y / o salida.

-g [a] x | y | d | X | Y | D | [columna] z [+ | -]brecha[u] (más ...)
Determine las lagunas de datos y los saltos de línea.

-h [i | o] [n] [+ c] [+ d] [+ rcomentario] [+ rtítulo] (más ...)
Omitir o producir registros de encabezado.

-icols[l] [sescala] [ocompensar] [,...] (más ...)
Seleccione las columnas de entrada (0 es la primera columna).

-ocols[, ...] (más ...)
Seleccione las columnas de salida (0 es la primera columna).

-p [x | y | z]azim/estudiante[/nivel z] [+ wcan0/lat0[/z0]] [+ vx0/y0] (más ...)
Seleccione la vista en perspectiva.

-s[cols] [a | r] (más ...)
Establecer el manejo de registros NaN.

-: [i | o] (más ...)
Intercambie la 1ª y la 2ª columna en la entrada y / o salida.

-^ or solo -
Imprima un mensaje corto sobre la sintaxis del comando, luego sale (NOTA: en Windows
usar solo -).

-+ or solo +
Imprima un mensaje de uso extenso (ayuda), incluida la explicación de cualquier
opción específica del módulo (pero no las opciones comunes de GMT), luego sale.

-? or no argumentos
Imprima un mensaje de uso completo (ayuda), incluida la explicación de las opciones, luego
salidas

--versión
Imprime la versión GMT y sal.

--show-datadir
Imprima la ruta completa al directorio compartido GMT y salga.

BODEGAS

Para la unidad de distancia del mapa, agregue unidad d para el grado de arco, m por minuto de arco, y s para arco
segundo, o e para el medidor [predeterminado], f para pie, k por km, M por la milla del estatuto, n para náutica
milla, y u para pie topográfico de EE. UU. Por defecto, calculamos tales distancias utilizando un esférico.
aproximación con grandes círculos. Anteponer - a una distancia (o la unidad no es ninguna distancia es
dado) para realizar cálculos de "Tierra plana" (más rápido pero menos preciso) o anteponer + a
realizar cálculos geodésicos exactos (más lento pero más preciso).

ASCII FORMATO PRECISIÓN

Los formatos de salida ASCII de datos numéricos se controlan mediante parámetros en su gmt.conf
expediente. La longitud y la latitud se formatean de acuerdo con FORMAT_GEO_OUT, mientras que otras
los valores se formatean de acuerdo con FORMAT_FLOAT_OUT. Tenga en cuenta que el formato en vigor puede
conducir a la pérdida de precisión en la salida, lo que puede dar lugar a varios problemas en sentido descendente. Si
encuentra que la salida no está escrita con suficiente precisión, considere cambiar a binario
producción (-bo si está disponible) o especifique más decimales usando la configuración FORMAT_FLOAT_OUT.

EJEMPLOS

Para convertir coordenadas UTM en metros a ubicaciones geográficas, dado un archivo utm.txt y
conociendo la zona UTM (y la zona o hemisferio), intente

gmt mapproject utm.txt -Ju + 11/1: 1 -C -I -F

Para transformar un archivo con (longitud, latitud) en posiciones (x, y) en cm en un Mercator
cuadrícula para una escala dada de 0.5 cm por grado, ejecute

gmt mapproject lonlatfile -R20 / 50/12/25 -Jm0.5c> xyfile

Para transformar varios archivos binarios de doble precisión de 2 columnas con (latitud, longitud)
en posiciones (x, y) en pulgadas en una cuadrícula de Mercator transversal (longitud central 75W) para
scale = 1: 500000 y suprime los puntos que quedarían fuera del área del mapa, ejecuta

gmt mapproject tracks. * -R-80 / -70 / 20/40 -Jt-75/1: 500000 -: -S -Di -bo -bi2> tmfile.b

Para convertir las coordenadas geodésicas (lon, lat, height) en el archivo old.dat del NAD27
Datum CONUS (Datum ID 131 que usa el elipsoide Clarke-1866) a WGS 84, ejecute

gmt mapproject old.dat -Th131> new.dat

Para calcular la distancia más cercana (en km) entre cada punto en el archivo de entrada quakes.dat
y los segmentos de línea dados en el archivo ASCII multisegmento coastline.xy, ejecute

gmt mapproject quakes.dat -Lcoastline.xy / k> quake_dist.dat

RESTRICCIONES

La región de entrada rectangular configurada con -R en general, se mapeará en un formato no rectangular
red. A no ser que -C está establecido, el punto más a la izquierda en esta cuadrícula tiene xvalue = 0.0, y el
el punto más bajo tendrá yvalue = 0.0. Por lo tanto, antes de digitalizar un mapa, ejecute el extremo
coordenadas del mapa a través de maproject usando la escala apropiada y vea qué valores (x, y)
están mapeados. Utilice estos valores al configurar la digitalización para tener la
la transformación inversa funciona correctamente, o alternativamente, usa awk para escalar y cambiar el
(x, y) valores antes de transformar.

Para algunas proyecciones, se puede utilizar una solución esférica a pesar de que el usuario haya seleccionado una
elipsoide. Esto ocurre cuando los usuarios -R la configuración implica una región que excede el dominio
en el que son válidas las expansiones de series elipsoidales. Estas son las condiciones: (1)
Cónica conformada de Lambert (-JL) y Albers Equal-Area (-JB) usará la solución esférica
cuando la escala del mapa excede 1.0E7. (2) Mercator transversal (-JT) y UTM (-JU) lo hará
utilizar la solución esférica cuando el límite oeste o este dado en -R Es mas que
10 grados del meridiano central, y (3) lo mismo para Cassini (-JC) pero con un límite de
solo 4 grados.

ELIPSOIDES Y ESFEROIDES

GMT usará fórmulas elipsoidales si están implementadas y el usuario ha seleccionado una
elipsoide como forma de referencia (ver PROJ_ELLIPSOID). El usuario debe conocer algunos
Escollos potenciales: (1) Para algunas proyecciones, como Transverse Mercator, Albers y
Cónica conforme de Lambert usamos las expresiones elipsoidales cuando las áreas mapeadas son
pequeño, y cambie a las expresiones esféricas (y sustituya el auxiliar apropiado
latitudes) para mapas más grandes. Las fórmulas elipsoidales se utilizan de la siguiente manera: (a) Transversal
Mercator: cuando todos los puntos están dentro de los 10 grados del meridiano central, (b) Proyecciones cónicas
cuando el rango longitudinal es inferior a 90 grados, (c) Proyección de Cassini cuando todos los puntos
están dentro de los 4 grados del meridiano central. (2) Cuando intenta hacer coincidir algunos datos históricos
datos (p. ej., coordenadas obtenidas con una determinada proyección y una determinada referencia
elipsoide), es posible que GMT dé resultados ligeramente diferentes. Uno probable
La fuente de esta discrepancia es que los cálculos más antiguos a menudo usaban dígitos menos significativos. Para
Por ejemplo, los ejemplos de Snyder a menudo usan el elipsoide de Clarke 1866 (definido por él como
un aplanamiento f = 1 / 294.98). De f obtenemos que la excentricidad al cuadrado es 0.00676862818
(esto es lo que usa GMT), mientras que Snyder completa y usa 0.00676866. Esta diferencia puede
dar discrepancias de varias decenas de cm. Si necesita reproducir las coordenadas proyectadas
con esta excentricidad ligeramente diferente, debe especificar su propio elipsoide con el
mismos parámetros que Clarke 1866, pero con f = 1 / 294.97861076. Además, tenga en cuenta que los mayores
Los datos pueden estar referenciados a diferentes datos y, a menos que sepa qué dato se utilizó y
convertir todos los datos a un dato común, puede experimentar desajustes de decenas a cientos de
metros. (3) Finalmente, tenga en cuenta que PROJ_SCALE_FACTOR tiene ciertos valores predeterminados para algunos
proyecciones, por lo que es posible que tenga que anular la configuración para hacer coincidir los resultados producidos
con otras configuraciones.

Utilice mapprojectgmt en línea utilizando los servicios de onworks.net