img2grdgmt - Online na nuvem

PROGRAMA:

NOME

img2grd - Extrair subconjunto do arquivo img no formato Mercator ou Geográfico

SINOPSE

img2grd arquivo img arquivo grd região tipo [] [[minlat / maxlat]] [] [ minutos ] [] [ navegação
] [[escada]] [[nível]] [ maxlon ] [ -n]

Observação: Nenhum espaço é permitido entre o sinalizador de opção e os argumentos associados.

DESCRIÇÃO

img2grd lê um arquivo de formato img, extrai um subconjunto e o grava em um arquivo de grade. o -M
opção determina se a projeção Spherical Mercator do arquivo img é ou não
preservada ou se uma grade geográfica deve ser escrita desfazendo a projeção de Mercator. Se
a grade geográfica é selecionada, você também pode solicitar uma reamostragem na exata -R dado.

É REQUERIDO ARGUMENTOS

arquivo img
Um arquivo de formato Mercator img, como a gravidade marinha ou campos de topografia do fundo do mar
estimado a partir de dados de altímetro de satélite por Sandwell e Smith. Se o usuário tiver definido
uma variável de ambiente $ GMT_DATADIR, Em seguida img2grd tentará encontrar arquivo img in
$ GMT_DATADIR; senão ele tentará abrir arquivo img diretamente.

-Garquivo grd
arquivo grd é o nome do arquivo de grade de saída.

-R [unidade]oeste/leste/sul/norte[/zmin/zmax] [r]
oeste, leste, sul e norte especifique a região de interesse, e você pode especificar
em graus decimais ou no formato [+ -] dd: mm [: ss.xxx] [W | E | S | N]. Acrescentar r se mais baixo
as coordenadas do mapa à esquerda e à direita superior são fornecidas em vez de w / e / s / n. Os dois
taquigrafias -Rg e -Rd significa domínio global (0/360 e -180 / + 180 de longitude
respectivamente, com -90 / + 90 na latitude). Alternativamente para a criação da grade, dê
Rcódigo/lat/nx/ny, Onde código é uma combinação de 2 caracteres de L, C, R (para a esquerda,
centro ou direita) e T, M, B para superior, meio ou inferior. por exemplo, BL para canto esquerdo inferior.
Isso indica qual ponto em uma região retangular o lon/lat coordenar refere-se
para, e as dimensões da grade nx e ny com espaçamentos de grade via -I é usado para criar
a região correspondente. Alternativamente, especifique o nome de um arquivo de grade existente
e os votos de -R as configurações (e espaçamento da grade, se aplicável) são copiadas da grade.
utilização -Runidade espera coordenadas projetadas (cartesianas) compatíveis com as escolhidas -J
e nós projetamos inversamente para determinar a região geográfica retangular real. Para
Visão de perspectiva (-p), opcionalmente anexar /zmin/zmax. Em caso de visão em perspectiva
(-p), um intervalo z (zmin, zmax) pode ser anexado para indicar a terceira dimensão. Esse
precisa ser feito apenas ao usar o -Jz opção, não ao usar apenas a -p opção.
Neste último caso, uma vista em perspectiva do plano é traçada, sem terceiros
dimensão.

OPCIONAL ARGUMENTOS

-C Defina as coordenadas de Mercator xey em relação ao centro de projeção [o padrão é
em relação ao canto esquerdo inferior da grade]. Requer -M.

-D [minlat / maxlat]
Use o intervalo de latitude estendido -80.738 / + 80.738. Alternativamente, anexe
minlat / maxlat como a extensão de latitude do arquivo img de entrada. [O padrão é
-72.006 / 72.006]. Normalmente não é necessário, pois podemos determinar a extensão de
inspeção do tamanho do arquivo.

-E Pode ser usado quando -M não está definido para forçar a grade final a ter exatamente o mesmo
região conforme solicitado com -R. Por padrão, a região final é uma projeção direta de
a região Mercator original e normalmente se estende um pouco além da
intervalo de latitude solicitado e, além disso, o incremento da grade na latitude não
coincidir com o incremento de longitude. No entanto, a reamostragem extra introduz pequenos
erros de interpolação e só devem ser usados ​​se a grade de saída deve corresponder ao
região solicitada e tem x_inc = y_inc. Neste caso, a região definida por -R devemos ser
dado em múltiplos do incremento (.eg, -R0/45/45/72).

-I Indicar minutos como a largura de um pixel img de entrada em minutos de longitude.
[O padrão é 2.0]. Normalmente não é necessário, pois podemos determinar o tamanho do pixel a partir de
inspeção do tamanho.

-M Gerar uma grade esférica de Mercator [O padrão é uma grade geográfica lon / lat]. o
A projeção esférica de Mercator do arquivo img é preservada, de modo que a região -R
definido pelo usuário é ligeiramente modificado; a região modificada corresponde às bordas
de pixels [ou grupos de navegação píxeis]. O cabeçalho do arquivo de grade é definido de forma que x e
os comprimentos do eixo y representam a distância das bordas oeste e sul da imagem,
medido em unidades padrão do usuário, com -Jm1 e o ajustado -R. Definindo o
omissão PROJ_ ELIPSOIDE = Esfera, o usuário pode fazer sobreposições com o -R
para que correspondam. Ver EXEMPLOS abaixo. O ajustado -R também está escrito na grade
observação do cabeçalho, para que possa ser encontrada mais tarde. Ver -C para definir as coordenadas em relação a
centro de projeção.

-Nnavegação Faça a média dos valores nos pixels img de entrada em navegação by navegação quadrados, e criar
um pixel de saída para cada quadrado. Se usado com -T3 vai relatar uma média
restrição entre 0 e 1. Se usado com -T2 a saída será o valor médio dos dados
ou NaN de acordo com se a restrição média é> 0.5. navegação deve dividir uniformemente
nas dimensões do imgfile em pixels. [Predefinição 1 não faz média].

-S [escada]
Multiplique os valores do arquivo img por escada antes de armazenar no arquivo de grade. [O padrão é
1.0]. Para arquivos img recentes: os arquivos img topo são armazenados em medidores (corrigidos) [-S1];
arquivos de gravidade ao ar livre em mGal * 10 [-S0.1 para obter mGal]; arquivos de deflexão vertical em
microrradianos * 10 [-S0.1 para obter micro-radianos], arquivos gradientes de gravidade vertical em
Eotvos * 50 [-S0.02 para obter Eotvos, ou -S0.002 para obter mGal / km]). Se não escada é dada
tentamos determinar a escala examinando o nome do arquivo em busca de pistas.

-Ttipo tipo lida com a codificação de informações de restrição. tipo = 0 indica que não
tais informações são codificadas no arquivo img (usado para versões anteriores a 1995 do
dados de gravidade) e obtém todos os dados. tipo > 0 indica que as informações de restrição são
codificado (1995 e versões posteriores (atuais) dos arquivos img) para que se possa
produza um arquivo de grade da seguinte maneira: -T1 obtém valores de dados em todos os pontos, -T2 obtém dados
valores em pontos restritos e NaN em pontos interpolados; -T3 ganha 1 em
pontos restritos e 0 nos pontos interpolados [o padrão é 1].

-V [nível] (mais ...)
Selecione o nível de verbosidade [c]. Particularmente recomendado aqui, pois é útil ver
como as coordenadas são ajustadas.

-Wmaxlon
Indicar maxlon como a extensão máxima de longitude do arquivo img de entrada. Versões
desde 1995 teve maxlon = 360.0, enquanto alguns arquivos anteriores tinham maxlon = 390.0.
[O padrão é 360.0].

-^ or apenas por -
Imprime uma mensagem curta sobre a sintaxe do comando e sai (NOTA: no Windows
use apenas -).

-+ or apenas por +
Imprima uma mensagem de uso extensivo (ajuda), incluindo a explicação de qualquer
opção específica do módulo (mas não as opções comuns GMT) e, em seguida, sai.

-? or não argumentos
Imprima uma mensagem completa de uso (ajuda), incluindo a explicação das opções e, em seguida,
saídas.

--versão
Imprima a versão GMT e saia.

--show-datadir
Imprima o caminho completo para o diretório de compartilhamento GMT e saia.

GEOGRÁFICO EXEMPLOS

A -M opção deve ser excluída se você precisar que a grade de saída esteja geograficamente
coordenadas. Para extrair dados na região -R-40 / 40 / -70 / -30 de mundo_grav.img.7.2 e
reprojetar para produzir coordenadas geográficas, você pode tentar

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

Como o espaçamento de latitude no arquivo img é equidistante em unidades de Mercator, o
a grade resultante não corresponderá ao especificado -R exatamente, e o espaçamento de latitude não
igual ao espaçamento da longitude. Se você precisa de uma correspondência exata com o seu -R e o mesmo espaçamento
em longitude e latitude, use o -E opção:

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

MERCADOR EXEMPLOS

Uma vez que os arquivos img estão em uma projeção de Mercator, você NÃO deve extrair uma grade geográfica
se o seu plano é fazer um mapa de Mercator. Se você fez isso, você termina de projetar e
reprojetar a grade, perdendo detalhes de comprimento de onda curto. Melhor usar -M e plotar a grade
usando uma projeção linear com a mesma escala da projeção de Mercator desejada (ver GMT
Exemplo 29). Para extrair dados na região -R-40 / 40 / -70 / -30 de mundo_grav.img.7.2, corre

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

Observe que o -V opção nos diz que o intervalo foi ajustado para
-R-40 / 40 / -70.0004681551 / -29.9945810754. Nós também podemos usar Grdinfo para descobrir que o arquivo de grade
cabeçalho mostra que sua região é -R0/80/0 / 67.9666667 Este é o intervalo de x, y que obteremos
de uma projeção esférica de Mercator usando -R-40 / 40 / -70.0004681551 / -29.9945810754 e
-Jm1. Assim, para pegar ship.lonlatgrav e usá-lo para fazer uma amostra do merc_grav.nc, podemos fazer o seguinte:

gmt set PROJ_ELLIPSOID Esfera

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

Recomenda-se usar o método acima para projetar e desprojetar os dados em tais
um aplicativo, porque então há apenas uma etapa de interpolação (em faixa de grd) Se um
primeiro tenta converter o arquivo de grade para lon, lat e depois faz uma amostra, há dois
etapas de interpolação (na conversão e na amostragem).

Para fazer uma grade lon, lat da grade acima, podemos usar

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

Em alguns casos, isso não será fácil, pois o -R nos dois sistemas de coordenadas pode não alinhar
Nós vamos. Quando isso acontece, também podemos usar (na verdade, pode ser sempre melhor usar)

gmt grd2xyz merc_grav.nc | projeto de mapa gmt \
-R-40/40 / -70.0004681551 / -29.994581075 -Jm1i -I | \
superfície gmt -R-40/40 / -70 / 70 -I2m -Ggrav.nc

Para fazer um mapa de Mercator da região acima, suponha que nosso valor gmt.conf para
PROJ_LENGTH_UNIT é polegadas. Então, como o arquivo merc_grav.nc acima é projetado com -Jm1i
tem 80 polegadas de largura. Podemos fazer um mapa de 8 polegadas de largura usando -Jx0.1i em qualquer programa de mapa
aplicado a esta grade (por exemplo, grdcontour, imagem, vista de grd) e, em seguida, para sobreposições que
trabalho em lon, lat (por exemplo, psxy, costa) podemos usar o acima ajustado -R e -Jm0.1 para obter
os dois sistemas combinam.

No entanto, podemos ser mais espertos do que isso. Percebendo que o arquivo img de entrada tinha pixels 2.0
minutos de largura (ou verificando nx e ny com grdinfo merc_grav.nc), percebemos que
merc_grav.nc usou a resolução total do arquivo img e tem 2400 por 2039 pixels, e
com 8 polegadas de largura são 300 pixels por polegada. Nós decidimos que não precisamos de tantos e nós
ficará satisfeito com 100 pixels por polegada, portanto, queremos calcular a média dos dados em 3 por 3
praças. (Se quisermos um gráfico de contorno, provavelmente escolheremos calcular a média dos dados muito mais
(por exemplo, 6 por 6) para obter contornos suaves.) Como 2039 não é divisível por 3, obteremos um
diferente ajustado -R desta vez:

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

Desta vez, descobrimos que a região ajustada é -R-40 / 40 / -70.023256525 / -29.9368261101 e o
a saída é de 800 por 601 pixels, um tamanho melhor para nós. Agora podemos criar um artificial
arquivo de iluminação para isso usando gradiente:

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

e se também tivermos um arquivo CPT chamado "grav.cpt", podemos criar um mapa de relevo colorido sombreado
como isso:

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

Suponha que você deseja obter apenas os valores de dados restritos de um arquivo img, em lat / lon
coordenadas. Então corra img2grd com o -T2 opções, use grd2xyz para despejar os valores, canalize
por meio de grep -v NaN para eliminar NaNs e canalizar projeto de mapa com o inverso
projeção como acima.

Use img2grdgmt online usando serviços onworks.net