img2grdgmt - Online en la nube

PROGRAMA:

NOMBRE

img2grd - Extrae un subconjunto del archivo img en formato Mercator o geográfico

SINOPSIS

img2grd archivo img archivo grd región tipo [] [[minlat / maxlat]] [] [ minutos ] [] [ navegación
] [[escala]] [[nivel]] [ Maxlon ] [ -n]

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

DESCRIPCIÓN

img2grd lee un archivo de formato img, extrae un subconjunto y lo escribe en un archivo de cuadrícula. los -M
La opción dicta si la proyección esférica de Mercator del archivo img es
conserva o si se debe escribir una cuadrícula geográfica deshaciendo la proyección de Mercator. Si
se selecciona la cuadrícula geográfica, también puede solicitar un remuestreo en el exacto -R dado.

REQUERIDO ARGUMENTOS

archivo img
Un archivo de formato img de Mercator, como los campos de gravedad marina o topografía del fondo marino
estimado a partir de datos de altímetro satelital de Sandwell y Smith. Si el usuario ha configurado
una variable de entorno $ GMT_DATADIR, entonces img2grd trataré de encontrar archivo img in
$ GMT_DATADIR; de lo contrario, intentará abrir archivo img .

-Garchivo grd
archivo grd es el nombre del archivo de cuadrícula de salida.

-R [unidad]oeste/este/sur/norte[/zmín/zmáx] [r]
oeste, este, sury norte especificar la región de interés, y puede especificar
en grados decimales o en formato [+ -] dd: mm [: ss.xxx] [W | E | S | N]. Adjuntar r si es menor
Las coordenadas del mapa izquierda y superior derecha se dan en lugar de w / e / s / n. Los dos
taquigrafía -Rg y -Rd significa dominio global (0/360 y -180 / + 180 de longitud
respectivamente, con -90 / + 90 en latitud). Alternativamente, para la creación de cuadrículas, dé
Rcodelón/lat/nx/ny, Donde código es una combinación de 2 caracteres de L, C, R (para la izquierda,
centro o derecha) y T, M, B para arriba, medio o abajo. por ejemplo, BL para la parte inferior izquierda.
Esto indica en qué punto de una región rectangular lon/lat coordinar se refiere
a, y las dimensiones de la cuadrícula nx y ny con espaciamiento de cuadrícula vía -I se usa para crear
la región correspondiente. Alternativamente, especifique el nombre de un archivo de cuadrícula existente
y del -R la configuración (y el espaciado de la cuadrícula, si corresponde) se copian de la cuadrícula.
Usar -Runidad espera coordenadas proyectadas (cartesianas) compatibles con las elegidas -J
y proyectamos inversamente para determinar la región geográfica rectangular real. Para
vista de perspectiva (-p), opcionalmente anexar /zmín/zmáx. En caso de vista en perspectiva
(-p), un rango z (zmín, zmáx) se puede agregar para indicar la tercera dimensión. Esta
debe hacerse solo cuando se utiliza el -Jz opción, no cuando se usa solo la -p .
En el último caso, se traza una vista en perspectiva del plano, sin una tercera
dimensión.

OPCIONAL ARGUMENTOS

-C Establezca las coordenadas de Mercator xey relativas al centro de proyección [el valor predeterminado es
relativo a la esquina inferior izquierda de la cuadrícula]. Requiere -M.

-D[minlat / maxlat]
Utilice el rango de latitud extendido -80.738 / + 80.738. Alternativamente, agregue
minlat / maxlat como la extensión de latitud del archivo img de entrada. [El valor predeterminado es
-72.006 / 72.006]. Por lo general, no es necesario, ya que podemos determinar el alcance de
inspección del tamaño del archivo.

-E Puede usarse cuando -M no está configurado para forzar que la cuadrícula final tenga exactamente el mismo
región según lo solicitado con -R. Por defecto, la región final es una proyección directa de
la región original de Mercator y típicamente se extenderá un poco más allá de la
rango de latitud solicitado, y además el incremento de la cuadrícula en latitud no
coincidir con el incremento de longitud. Sin embargo, el remuestreo adicional introduce pequeños
errores de interpolación y solo debe usarse si la cuadrícula de salida debe coincidir
región solicitada y tienen x_inc = y_inc. En este caso, la región establecida por -R debe ser
dado en múltiplos del incremento (.eg, -R0/45/45/72).

-I Indicar minutos como el ancho de un píxel de imagen de entrada en minutos de longitud.
[El valor predeterminado es 2.0]. No suele ser necesario, ya que podemos determinar el tamaño de píxel a partir de
inspección del tamaño.

-M Genere una cuadrícula esférica de Mercator [El valor predeterminado es una cuadrícula geográfica lon / lat]. los
Se conserva la proyección esférica de Mercator del archivo img, de modo que la región -R
establecido por el usuario se modifica ligeramente; la región modificada corresponde a los bordes
de píxeles [o grupos de navegación píxeles]. El encabezado del archivo de cuadrícula se establece de modo que la xy
Las longitudes del eje y representan la distancia desde los bordes oeste y sur de la imagen,
medido en unidades predeterminadas por el usuario, con -Jm1 y el ajustado -R. Al establecer el
tu préstamo estudiantil PROY_ ELIPSOIDE = Esfera, el usuario puede hacer superposiciones con el ajuste -R
para que coincidan. Ver EJEMPLOS debajo. El ajustado -R también está escrito en la cuadrícula
comentario del encabezado, por lo que se puede encontrar más tarde. Ver -C para establecer coordenadas relativas a
centro de proyección.

-Nnavegación Promedio de los valores en los píxeles img de entrada en navegación by navegación cuadrados y crear
un píxel de salida para cada uno de esos cuadrados. Si se usa con -T3 informará un promedio
restricción entre 0 y 1. Si se usa con -T2 la salida será el valor de datos promedio
o NaN según si la restricción media es> 0.5. navegación debe dividir uniformemente
en las dimensiones del archivo img en píxeles. [Defecto 1 no promedia].

-S[escala]
Multiplica los valores del archivo img por escala antes de almacenarlo en un archivo de cuadrícula. [El valor predeterminado es
1.0]. Para archivos img recientes: los archivos img topo se almacenan en metros (corregidos) [-S1];
archivos de gravedad al aire libre en mGal * 10 [-S0.1 para obtener mGal]; archivos de deflexión vertical en
micro-radianes * 10 [-S0.1 para obtener micro-radianes], archivos de gradiente de gravedad vertical en
Eotvos * 50 [-S0.02 para obtener Eotvos, o -S0.002 para obtener mGal / km]). Si no escala es dado
tratamos de determinar la escala examinando el nombre del archivo en busca de pistas.

-Ttipo tipo maneja la codificación de la información de restricción. tipo = 0 indica que no
dicha información está codificada en el archivo img (utilizado para las versiones anteriores a 1995 del
datos de gravedad) y obtiene todos los datos. tipo > 0 indica que la información de restricción es
codificado (1995 y versiones posteriores (actuales) de los archivos img) para que uno pueda
producir un archivo de cuadrícula de la siguiente manera: -T1 obtiene valores de datos en todos los puntos, -T2 obtiene datos
valores en puntos restringidos y NaN en puntos interpolados; -T3 obtiene 1 en
puntos restringidos y 0 en puntos interpolados [El valor predeterminado es 1].

-V [nivel] (más ...)
Seleccione el nivel de verbosidad [c]. Particularmente recomendado aquí, ya que es útil ver
cómo se ajustan las coordenadas.

-WMaxlon
Indicar Maxlon como la extensión de longitud máxima del archivo img de entrada. Versiones
desde 1995 han tenido Maxlon = 360.0, mientras que algunos archivos anteriores tenían Maxlon = 390.0.
[El valor predeterminado es 360.0].

-^ 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.

GEOGRÁFICO EJEMPLOS

El -M La opción debe excluirse si necesita que la cuadrícula de salida esté en
coordenadas. Para extraer datos en la región -R-40 / 40 / -70 / -30 desde mundo_grav.img.7.2 y
reproyectar para obtener coordenadas geográficas, puede intentar

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40 / -70 / -30 -V

Debido a que el espaciado de latitud en el archivo img es equidistante en unidades de Mercator, el
la cuadrícula resultante no coincidirá con la especificada -R exactamente, y el espaciado de latitud no
igual al espaciado de longitud. Si necesita una coincidencia exacta con su -R y el mismo espaciado
en longitud y latitud, utilice el -E opción:

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40 / -70 / -30 -E -V

MERCATOR EJEMPLOS

Dado que los archivos img están en una proyección de Mercator, NO debe extraer una cuadrícula geográfica
si su plan es hacer un mapa de Mercator. Si hiciste eso terminaste de proyectar y
reproyecta la rejilla, perdiendo detalles de longitud de onda corta. Mejor usar -M y trazar la cuadrícula
utilizando una proyección lineal con la misma escala que la proyección de Mercator deseada (ver GMT
Ejemplo 29). Para extraer datos en la región -R-40 / 40 / -70 / -30 desde mundo_grav.img.7.2, correr

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav.nc -R-40/40 / -70 / -30 -V

Tenga en cuenta que -V La opción nos dice que el rango se ajustó a
-R-40 / 40 / -70.0004681551 / -29.9945810754. También podemos usar información grd para encontrar que el archivo de cuadrícula
El encabezado muestra que su región es -R0/80/0 / 67.9666667 Este es el rango de x, y que obtendremos
de una proyección esférica de Mercator usando -R-40 / 40 / -70.0004681551 / -29.9945810754 y
-Jm1. Por lo tanto, para tomar ship.lonlatgrav y usarlo para probar el merc_grav.nc, podemos hacer esto:

gmt set PROJ_ELLIPSOID Esfera

gmt mapproject -R-40/40 / -70.0004681551 / -29.9945810754 -Jm1i ship.lonlatgrav | \
gmt grdtrack -Gmerc_grav.nc | proyecto de mapa gmt \
-R-40/40 / -70.0004681551 / -29.9945810754 -Jm1i -I> ship.lonlatgravsat

Se recomienda utilizar el método anterior para proyectar y desproyectar los datos en tales
una aplicación, porque entonces solo hay un paso de interpolación (en pista de tierra). Si uno
primero intenta convertir el archivo de cuadrícula a lon, lat y luego probarlo, hay dos
Pasos de interpolación (en conversión y en muestreo).

Para hacer una cuadrícula lon, lat de la cuadrícula anterior, podemos usar

gmt grdproject merc_grav.nc -R-40/40 / -70.0004681551 / -29.9945810754 -Jm1i -I -D2m -Ggrav.nc

En algunos casos, esto no será fácil ya que -R en los dos sistemas de coordenadas puede no alinearse
bien. Cuando esto sucede, también podemos usar (de hecho, siempre puede ser mejor usar)

gmt grd2xyz merc_grav.nc | proyecto de mapa gmt \
-R-40/40 / -70.0004681551 / -29.994581075 -Jm1i -I | \
superficie gmt -R-40/40 / -70 / 70 -I2m -Ggrav.nc

Para hacer un mapa de Mercator de la región anterior, suponga nuestro valor gmt.conf para
PROJ_LENGTH_UNIT son pulgadas. Entonces, dado que el archivo merc_grav.nc anterior se proyecta con -Jm1i
tiene 80 pulgadas de ancho. Podemos hacer un mapa de 8 pulgadas de ancho usando -Jx0.1i en cualquier programa de mapas
aplicado a esta cuadrícula (p. ej., grdcontorno, grdimagen, grdview), y luego para superposiciones que
trabajar en lon, lat (p. ej., psxy, costa rica) podemos usar el ajuste anterior -R y -Jm0.1 para obtener
los dos sistemas para que coincidan.

Sin embargo, podemos ser más inteligentes que esto. Al darse cuenta de que el archivo img de entrada tenía píxeles 2.0
minutos de ancho (o comprobando nx y ny con grdinfo merc_grav.nc) nos damos cuenta de que
merc_grav.nc usó la resolución completa del archivo img y tiene 2400 por 2039 píxeles, y
a 8 pulgadas de ancho, esto es 300 píxeles por pulgada. Decidimos que no necesitamos tantos y
estará satisfecho con 100 píxeles por pulgada, por lo que queremos promediar los datos en 3 por 3
cuadrícula. (Si queremos un gráfico de contorno, probablemente elegiremos promediar los datos mucho más
(por ejemplo, 6 por 6) para obtener contornos suaves.) Dado que 2039 no es divisible por 3, obtendremos un
diferente ajustado -R esta vez:

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav_2.nc -R-40/40 / -70 / -30 -N3 -V

Esta vez encontramos que la región ajustada es -R-40 / 40 / -70.023256525 / -29.9368261101 y el
La salida es de 800 por 601 píxeles, un tamaño mejor para nosotros. Ahora podemos crear un artificial
archivo de iluminación para este uso degradado:

gmt grdgradient merc_grav_2.nc -Gillum.nc -A0 / 270 -Ne0.6

y si también tenemos un archivo CPT llamado "grav.cpt" podemos crear un mapa de relieve sombreado por colores
Me gusta esto:

gmt grdimage merc_grav_2.nc -Iillum.nc -Cgrav.cpt -Jx0.1i -K> map.ps
gmt psbasemap -R-40/40 / -70.023256525 / -29.9368261101 -Jm0.1i -Ba10 -O >> map.ps

Suponga que desea obtener solo los valores de datos restringidos de un archivo img, en lat / lon
coordenadas. Entonces corre img2grd con el -T2 opción, uso grd2xyz para volcar los valores, pipe
a través de grep -v NaN para eliminar los NaN y canalizar maproject con la inversa
proyección como arriba.

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