grdgradientgmt - Online na nuvem

PROGRAMA:

NOME

grdgradient - Calcula a derivada direcional ou gradiente de uma grade

SINOPSE

gradiente in_grdfile out_grdfile [ Azim[/azim2]] [[a][c][o][n]] [
[s | p]azim / elev[/ambiente/difundir/especular/brilhar]] [ bandeira ] [[e][t][amp] [/sigma[/compensar]]
] [ região ] [ arquivo de inclinação ] [[nível]] [ -fg ] [ -n]

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

DESCRIÇÃO

gradiente pode ser usado para calcular a derivada direcional em uma determinada direção (-A),
ou a direção (-S) [e a magnitude (-D)] do gradiente vetorial dos dados.

Os valores estimados na primeira / última linha / coluna de saída dependem das condições de limite (ver
-L).

É REQUERIDO ARGUMENTOS

in_grdfile
Arquivo de grade 2-D a partir do qual calcular a derivada direcional. (Veja FORMATOS DE ARQUIVOS DE GRID
abaixo).

-Gout_grdfile
Nome do arquivo de grade de saída para a derivada direcional. (Veja FORMATOS DE ARQUIVOS DE GRID
abaixo).

OPCIONAL ARGUMENTOS

-AAzim[/azim2]
Direção azimutal para uma derivada direcional; Azim é o ângulo em x, y
plano medido em graus positivos no sentido horário do norte (a direção + y) em direção
leste (a direção + x). O negativo da derivada direcional,
- [dz / dx * sin (Azim) + dz / dy * cos (Azim)], seja encontrado; negação produz valores positivos
quando a inclinação de z (x, y) é descendente no Azim direção, o sentido correto para
sombreando a iluminação de uma imagem (ver imagem e vista de grd) por uma fonte de luz
acima do plano x, y brilhando do Azim direção. Opcionalmente, forneça dois
azimutes, -AAzim/azim2, caso em que os gradientes em cada uma dessas direções são
calculado e aquele de maior magnitude é mantido; isso é útil para
iluminar dados com duas direções de estruturas lineares, por exemplo, -A0/270
ilumina do norte (topo) e oeste (esquerda).

-D [a] [c] [o] [n]
Encontre a direção do gradiente positivo (inclinação para cima) dos dados. Para ao invés
encontre o aspecto (a direção do declive), use -Sim. Por padrão, as direções são
medido no sentido horário do norte, como Azim in -A acima de. Acrescentar c usar convencional
Ângulos cartesianos medidos no sentido anti-horário a partir da direção x positiva (leste).
Acrescentar o para relatar orientações (0-180) em vez de direções (0-360). Acrescentar n
para adicionar 90 graus a todos os ângulos (por exemplo, para dar golpes locais na superfície).

-E [s | p]azim / elev[/ambiente/difundir/especular/brilhar]
Calcule a radiância lambertiana apropriada para usar com imagem e vista de grd. O
A reflexão lambertiana assume uma superfície ideal que reflete toda a luz que
atinge-o e a superfície parece igualmente brilhante em todas as direções de visualização. Azim
e criado são o azimute e a elevação do vetor de luz. Opcionalmente, forneça ambiente
difundir especular brilhar que são parâmetros que controlam as propriedades de refletância
da superfície. Os valores padrão são: 0.55/0.6/0.4/10 Para deixar alguns dos valores
intocado, especifique = como o novo valor. Por exemplo -E60/30/=/0.5 define o Azim criado
e difundir para 60, 30 e 0.5 e deixa os outros parâmetros de refletância
intocado. Acrescentar s para usar um algoritmo lambertiano mais simples. Observe que com este formulário
você só precisa fornecer os parâmetros de azimute e elevação. Acrescentar p usar o
Aproximação linear por partes de Peucker (algoritmo mais simples, mas mais rápido; neste caso
que o Azim e criado são conectados a 315 e 45 graus. Isso significa que mesmo se você
forneça outros valores, eles serão ignorados.)

-Lbandeira Condição de contorno bandeira pode ser x or y or xy indicando que os dados são periódicos na faixa de
x ou y ou ambos, ou bandeira pode ser g indicando as condições geográficas (x e y são
lon e lat). [O padrão usa condições "naturais" (segunda derivada parcial normal
até a borda é zero).]

-Internet][amp] [/sigma[/compensar]]
Normalização. [Padrão: sem normalização.] Os gradientes reais g são compensados ​​e
dimensionado para produzir gradientes normalizados gn com uma magnitude de saída máxima de amp.
If amp não é fornecido, padrão amp = 1. Se compensar não é fornecido, é definido para o
média de g. -N rendimentos gn = amp * (g - compensar) / máx (abs (g - compensar)). -Não
normaliza usando uma distribuição de Laplace cumulativa que produz gn = amp * (1.0 -
exp (sqrt(2) * (g - compensar)/ sigma)) Onde sigma é estimado usando a norma L1 de
(g - compensar) se não for fornecido. -Não normaliza usando um Cauchy cumulativo
distribuição rendendo gn = (2* amp / PI) * atan ((g - compensar)/ sigma) Onde sigma
é estimado usando a norma L2 de (g - compensar) se não for fornecido.

-R [unidade]xmin/xmax/ymin/ymax[r] (mais ...)
Especifique a região de interesse. Usando o -R opção irá selecionar uma subseção de
in_grdfile grade. Se esta subseção ultrapassar os limites da grade, apenas o
região comum será extraída.

-Sarquivo de inclinação
Nome do arquivo de grade de saída com magnitudes escalares de vetores de gradiente. Requer -D
mas faz -G opcional.

-V [nível] (mais ...)
Selecione o nível de verbosidade [c].

-fg Grades geográficas (dimensões de longitude, latitude) serão convertidas em metros
por meio de uma aproximação de "Terra plana" usando os parâmetros elipsoides atuais.

-n [b | c | l | n] [+ a] [+ bBC] [+ c] [+ tlimiar] (mais ...)
Selecione o modo de interpolação para grades.

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

GRID DISTÂNCIA UNIDADES

Se a grade não tiver medidor como unidade horizontal, anexe +uunidade para o arquivo de entrada
nome para converter da unidade especificada em medidor. Se sua grade é geográfica, converta
distâncias em metros, fornecendo -fg ao invés.

HINTS

Se você não sabe o que -N opções para usar para fazer um arquivo de intensidade para imagem or
vista de grd, uma boa primeira tentativa é -Não0.6.

Normalmente, 255 tons são mais do que suficientes para fins de visualização. Você pode economizar 75% do disco
espaço anexando = nb / a ao nome do arquivo de saída out_grdfile.

Se você quiser fazer vários mapas iluminados de sub-regiões de um grande conjunto de dados, e você
precisa que os efeitos de iluminação sejam consistentes em todos os mapas, use o -N opção e
fornecer o mesmo valor de sigma e compensar para gradiente para cada mapa. Um bom palpite é
compensar = 0 e sigma encontrado por Grdinfo -L2 or -L1 aplicado a um gradiente grd não normalizado.

Se você simplesmente precisa do x- ou y- derivados da rede, use Grdmath.

GRID ARQUIVO FORMATOS

Por padrão, o GMT grava a grade como flutuadores de precisão simples em um netCDF de reclamação COARDS
formato de arquivo. No entanto, o GMT é capaz de produzir arquivos de grade em muitas outras grades comumente usadas
formatos de arquivo e também facilita o chamado "empacotamento" de grades, escrevendo em ponto flutuante
dados como números inteiros de 1 ou 2 bytes. Para especificar a precisão, escala e deslocamento, o usuário deve
adicione o sufixo =id[/escada/compensar[/nan]], Onde id é um identificador de duas letras da grade
tipo e precisão, e escada e compensar são fator de escala opcional e deslocamento para ser
aplicado a todos os valores da grade, e nan é o valor usado para indicar dados ausentes. Em caso
os dois personagens id não é fornecido, como em =/escada do que um id=nf é assumido. Quando
leitura de grades, o formato é geralmente reconhecido automaticamente. Se não, o mesmo sufixo
pode ser adicionado para inserir nomes de arquivo de grade. Ver grdconverter e Seção de grade-formato de arquivo do
Referência técnica e livro de receitas da GMT para obter mais informações.

Ao ler um arquivo netCDF que contém várias grades, o GMT lerá, por padrão, o
primeira grade bidimensional que pode ser encontrada naquele arquivo. Para persuadir o GMT a ler outro
variável multidimensional no arquivo de grade, anexar ?nome da var ao nome do arquivo, onde
nome da var é o nome da variável. Observe que você pode precisar escapar do significado especial
of ? em seu programa de shell, colocando uma barra invertida na frente dele, ou colocando o
nome do arquivo e sufixo entre aspas ou aspas duplas. o ?nome da var sufixo também pode ser usado
para grades de saída para especificar um nome de variável diferente do padrão: "z". Ver
grdconverter e Seções modificadores-para-CF e formato de arquivo de grade do GMT Técnico
Referência e livro de receitas para obter mais informações, particularmente sobre como ler emendas de 3,
Grades de 4 ou 5 dimensões.

EXEMPLOS

Para fazer um arquivo para iluminar os dados em geoid.nc usando gradientes exp normalizados em
o intervalo [-0.6,0.6] imitando fontes de luz nas direções norte e oeste:

gmt grdgradient geoid.nc -A0 / 270 -Ggradients.nc = nb / a -Ne0.6 -V

Para encontrar as orientações de azimute do tecido do fundo do mar no arquivo topo.nc:

gmt grdgradiente topo.nc -Dno -Gazimuths.nc -V

REFERÊNCIAS

Horn, BKP, Hill-Shading and the Reflectance Map, Proceedings of the IEEE, Vol. 69, No.
1, janeiro de 1981, pp. 14-47. (http://people.csail.mit.edu/bkph/papers/Hill-Shading.pdf)

Use grdgradientgmt online usando serviços onworks.net