mapprojectgmt - Online na nuvem

PROGRAMA:

NOME

mapproject - Faça transformações de mapa para frente e inversa, conversões de datum e geodésia

SINOPSE

projeto de mapa [ tabelas ] parâmetros região [ b|B|f|F|o|O[lon0/lat0]] [[dx/dy]] [ c|i|p ] [
[dado]] [[unidade]] [[x0/y0/] [[+|-]unidade][+|-]] [] [ linha.xy[/ [+|-]unidade][+]] [
[a|c|g|m]] [[d|e ] [] [[h]da[/para]] [[nível]] [[w|h]] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o]]

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

DESCRIÇÃO

projeto de mapa lê (longitude, latitude) posições de tabelas [ou entrada padrão] e
calcula as coordenadas (x, y) usando a projeção do mapa e as escalas especificadas. Opcionalmente,
pode ler (x, y) posições e calcular (longitude, latitude) valores fazendo o inverso
transformação. Isso pode ser usado para transformar pontos lineares (x, y) obtidos por digitalização
um mapa de projeção conhecida para coordenadas geográficas. Também pode calcular distâncias ao longo
trilha, para um ponto fixo, ou abordagem mais próxima de uma linha. Finalmente, pode ser usado para realizar
várias conversões de datum. Campos de dados adicionais são permitidos após as 2 primeiras colunas
que deve ter (longitude, latitude) ou (x, y). Veja a opção -: sobre como ler
(latitude, longitude) arquivos.

É REQUERIDO ARGUMENTOS

-Jparâmetros (mais ...)
Selecione a projeção do mapa.

-R [unidade]xmin/xmax/ymin/ymax[r] (mais ...)
Especifique a região de interesse. Caso especial para a projeção UTM: Se -C é usado
e -R não é fornecido, então a região é definida para coincidir com a zona UTM fornecida, então
quanto a preservar a solução elipsoidal completa (Veja RESTRIÇÕES para mais
em formação).

OPCIONAL ARGUMENTOS

mesa Um ou mais ASCII (ou binário, veja -bi[notas][tipo]) arquivo (s) de tabela de dados contendo um
número de colunas de dados. Se nenhuma tabela for fornecida, então lemos a entrada padrão.

-Ab | B | f | F | o | O [lon0/lat0]
-Af calcula o azimute (para frente) do ponto fixo lon / lat para cada ponto de dados.
Use -Ausente para obter o azimute de volta dos pontos de dados para o ponto fixo. Usar -Ao para obter
orientações (-90/90) em vez de azimutes (0/360). Maiúsculas F, B or O precisarão
converter de latitudes geodésicas para geocêntricas e estimar o azimute de geodésicas
(assumindo que o elipsóide atual não é uma esfera). Se nenhum ponto fixo for dado, então
calculamos o azimute (ou azimute posterior) do ponto anterior.

-C [dx/dy]
Defina o centro das coordenadas projetadas para estar no centro de projeção do mapa [O padrão é
canto inferior esquerdo]. Opcionalmente, adicione deslocamentos nas unidades projetadas a serem adicionadas (ou
subtraído quando -I é definido) para (a partir) das coordenadas projetadas, como falso
leste e norte para zonas de projeção específicas [0/0]. A unidade usada para o
offsets é a unidade de distância do gráfico em vigor (consulte PROJ_LENGTH_UNIT), a menos que -F is
usado, caso em que os deslocamentos são em metros.

-Dc | i | p
Substitua temporariamente PROJ_LENGTH_UNIT e use c (cm), i (polegada), ou p (pontos)
em vez de. Não pode ser usado com -F.

-E [dado]
Converter de geodésico (lon, lat, altura) em Terra centrada na Terra (ECEF)
(x, y, z) coordenadas (adicionar -I para a conversão inversa). Anexar ID de referência (ver -Qd)
ou dar elipsóide:dx,dy,dz onde elipsóide pode ser um ID elipsóide (ver -Qe) ou
dado como a[, * inv_f *], onde a é o semi-eixo maior e inv_f é o inverso
achatamento (0 se omitido). Se dado é - ou não dado assumimos WGS-84.

-F [unidade]
Forçar escala 1: 1, ou seja, saída (ou entrada, consulte -I) os dados estão em projeções reais
metros. Para especificar outras unidades, anexe a unidade desejada (consulte UNIDADES). Sem -F,
a saída (ou entrada, veja -I) estão nas unidades especificadas por PROJ_LENGTH_UNIT (mas
Vejo -D).

-G [x0/y0/] [[+ | -]unidade] [+ | -]
Calcular distâncias ao longo da trilha or para o ponto opcional definido com -Gx0 / y0. Acrescentar
a unidade de distância (ver UNIDADES), incluindo c (Distância cartesiana usando entrada
coordenadas) ou C (Distância cartesiana usando coordenadas projetadas). o C unidade
exige -R e -J a ser definido. Sem nenhum ponto fixo é fornecido, calculamos acumular
distâncias ao longo da trilha. Acrescentar - para obter distância incremental entre sucessivas
pontos. Acrescentar + para especificar o segundo ponto por meio de duas colunas extras no arquivo de entrada.

-I Faça a transformação inversa, ou seja, obtenha (longitude, latitude) dos dados (x, y).

-Llinha.xy[/ [+ | -]unidade] [+]
Determine a distância mais curta dos pontos de dados de entrada à (s) linha (s) dada (s) em
o arquivo ASCII multisegment linha.xy. A distância e as coordenadas do
o ponto mais próximo será anexado à saída como três novas colunas. Anexe o
unidade de distância (ver UNIDADES), incluindo c (Distância cartesiana usando coordenadas de entrada)
or C (Distância cartesiana usando coordenadas projetadas). o C unidade requer -R e
-J a ser definido. Finalmente, acrescente + para relatar a id do segmento de linha e o fracionário
número do ponto em vez de lon / lat do ponto mais próximo.

-N [a | c | g | m]
Converter de latitudes geodésicas (usando o elipsóide atual; consulte PROJ_ELLIPSOID)
para uma das quatro latitudes auxiliares diferentes (as longitudes não são afetadas). Escolher
da autálico, cinformal, geocêntrico, e mlatitudes eridionais [geocêntricas]. Usar -I
para converter de latitudes auxiliares em latitudes geodésicas.

-Q [d | e Liste todos os parâmetros de projeção. Para listar apenas datums, use -Qd. Para apenas listar
elipsóides, use -Qe.

-S Suprima os pontos que estão fora da região.

-º]da[/para]
Coordenar conversões entre datums da e para usando o Molodensky padrão
transformação. Usar -º se a 3ª coluna de entrada tiver altura acima do elipsóide [Padrão
assume altura = 0, ou seja, no elipsóide]. Especifique os datums usando o ID do datum (ver
-Qd) ou dar elipsóide:dx,dy,dz onde elipsóide pode ser um ID elipsóide (ver -Qe) ou
dado como a[, * inv_f *], onde a é o semi-eixo maior e inv_f é o inverso
achatamento (0 se omitido). Se dado é - ou não dado assumimos WGS-84. -T pode ser
usado em conjunto com -R -J para alterar o datum antes da projeção coordenada
(Adicionar, acrescentar -I para aplicar a conversão do datum após a projeção inversa). Certifique-se de que
a configuração PROJ_ELLIPSOID está correta para o seu caso.

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

-W [w | h]
Imprime a largura e a altura do mapa na saída padrão. Nenhum arquivo de entrada é lido. Para apenas
produza a largura ou altura, acrescente w or h, respectivamente. As unidades do
as dimensões podem ser alteradas via -D.

-bi[notas] [t] (mais ...)
Selecione a entrada binária nativa. [O padrão é 2 colunas de entrada].

-bo [notas][tipo] (mais ...)
Selecione a saída binária nativa. [O padrão é o mesmo que a entrada].

-d [i | o]sem dados (mais ...)
Substitua as colunas de entrada iguais sem dados com NaN e faça o inverso na saída.

-f [i | o]colinfo (mais ...)
Especifique os tipos de dados de colunas de entrada e / ou saída.

-g [a] x | y | d | X | Y | D | [col] z [+ | -]lacuna[A] (mais ...)
Determine lacunas de dados e quebras de linha.

-h [i | o] [n] [+ c] [+ d] [+ rcomentário] [+ rtítulo] (mais ...)
Pule ou produza registro (s) de cabeçalho.

-icols[l] [sescada] [ocompensar] [,...] (mais ...)
Selecione as colunas de entrada (0 é a primeira coluna).

-ocols[, ...] (mais ...)
Selecione as colunas de saída (0 é a primeira coluna).

-p [x | y | z]Azim/criado[/nível z] [+ wlon0/lat0[/z0]] [+ vx0/y0] (mais ...)
Selecione a vista em perspectiva.

-s [cols] [a | r] (mais ...)
Definir o manuseio de registros NaN.

-: [i | o] (mais ...)
Troque a 1ª e 2ª coluna na entrada e / ou saída.

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

UNIDADES

Para unidade de distância do mapa, anexe unidade d para grau de arco, m para minuto de arco, e s para arco
segundo, ou e para medidor [padrão], f para o pé, k por km, M para milhas estatutárias, n para náutico
milha, e u para o pé de pesquisa dos EUA. Por padrão, calculamos tais distâncias usando um esférico
aproximação com grandes círculos. Prepend - a uma distância (ou a unidade está sem distância é
dado) para realizar cálculos de "Terra plana" (mais rápido, mas menos preciso) ou prefixar + para
realizar cálculos geodésicos exatos (mais lentos, mas mais precisos).

ASCII FORMATO PRECISÃO

Os formatos de saída ASCII de dados numéricos são controlados por parâmetros em seu gmt.conf
Arquivo. Longitude e latitude são formatadas de acordo com FORMAT_GEO_OUT, enquanto outras
os valores são formatados de acordo com FORMAT_FLOAT_OUT. Esteja ciente de que o formato em vigor pode
levar à perda de precisão na saída, o que pode levar a vários problemas a jusante. Se
você descobrir que a saída não foi escrita com precisão suficiente, considere mudar para binário
saída (-bo se disponível) ou especifique mais decimais usando a configuração FORMAT_FLOAT_OUT.

EXEMPLOS

Para converter coordenadas UTM em metros para localizações geográficas, dado um arquivo utm.txt e
conhecendo a zona UTM (e zona ou hemisfério), tente

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

Para transformar um arquivo com (longitude, latitude) em (x, y) posições em cm em um Mercator
grade para uma determinada escala de 0.5 cm por grau, execute

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

Para transformar vários arquivos binários de dupla precisão de 2 colunas com (latitude, longitude)
em posições (x, y) em polegadas em uma grade transversal de Mercator (longitude central 75W) para
escala = 1: 500000 e suprimir aqueles pontos que cairiam fora da área do mapa, execute

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

Para converter as coordenadas geodésicas (lon, lat, altura) no arquivo old.dat do NAD27
Datum CONUS (Datum ID 131 que usa o elipsóide Clarke-1866) para WGS 84, execute

gmt mapproject old.dat -Th131> new.dat

Para calcular a distância mais próxima (em km) entre cada ponto no arquivo de entrada quakes.dat
e os segmentos de linha fornecidos no arquivo ASCII multisegment lifetime.xy, execute

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

RESTRIÇÕES

A região de entrada retangular definida com -R será, em geral, mapeado em um formato não retangular
grade. A não ser que -C é definido, o ponto mais à esquerda nesta grade tem xvalue = 0.0, e o
o ponto mais inferior terá yvalue = 0.0. Assim, antes de digitalizar um mapa, execute o extremo
coordenadas do mapa através de projeto de mapa usando a escala apropriada e ver quais valores (x, y)
eles são mapeados. Use esses valores ao configurar para digitalização, a fim de ter o
a transformação inversa funciona corretamente ou, alternativamente, use awk para dimensionar e mudar o
(x, y) valores antes da transformação.

Para alguma projeção, uma solução esférica pode ser usada apesar de o usuário ter selecionado um
elipsóide. Isso ocorre quando os usuários -R configuração implica em uma região que excede o domínio
em que as expansões em série elipsoidal são válidas. Estas são as condições: (1)
Cônica conformada de Lambert (-JL) e Albers Equal-Area (-JB) usará a solução esférica
quando a escala do mapa excede 1.0E7. (2) Mercator transversal (-JT) e UTM (-JU) vai vai
usar a solução esférica quando o limite oeste ou leste dado em -R é mais que
10 graus do meridiano central, e (3) o mesmo para a Cassini (-JC), mas com um limite de
apenas 4 graus.

ELIPSOIDES E ESFERÓIDES

O GMT usará fórmulas elipsoidais se forem implementadas e o usuário tiver selecionado um
elipsóide como forma de referência (ver PROJ_ELLIPSOID). O usuário precisa estar ciente de alguns
armadilhas potenciais: (1) Para algumas projeções, como Transverse Mercator, Albers e
Na cônica conformada de Lambert, usamos as expressões elipsoidais quando as áreas mapeadas são
pequeno, e mudar para as expressões esféricas (e substituindo o auxiliar apropriado
latitudes) para mapas maiores. As fórmulas elipsoidais são usadas da seguinte forma: (a) Transversal
Mercator: quando todos os pontos estão dentro de 10 graus do meridiano central, (b) projeções cônicas
quando o alcance longitudinal é inferior a 90 graus, (c) projeção Cassini quando todos os pontos
estão dentro de 4 graus do meridiano central. (2) Quando você está tentando combinar algum histórico
dados (por exemplo, coordenadas obtidas com uma certa projeção e uma certa referência
elipsóide), você pode descobrir que o GMT fornece resultados ligeiramente diferentes. Um provável
A fonte dessa incompatibilidade é que os cálculos mais antigos costumavam usar dígitos menos significativos. Para
exemplo, os exemplos de Snyder costumam usar o elipsóide Clarke 1866 (definido por ele como tendo
a achatamento f = 1 / 294.98). De f obtemos a excentricidade ao quadrado para 0.00676862818
(isso é o que GMT usa), enquanto Snyder arredonda e usa 0.00676866. Essa diferença pode
dar discrepâncias de várias dezenas de cm. Se você precisar reproduzir as coordenadas projetadas
com esta excentricidade ligeiramente diferente, você deve especificar seu próprio elipsóide com o
mesmos parâmetros de Clarke 1866, mas com f = 1 / 294.97861076. Além disso, esteja ciente de que
os dados podem ser referenciados a diferentes datums, e a menos que você saiba qual datum foi usado e
converter todos os dados em um datum comum, você pode experimentar incompatibilidades de dezenas a centenas de
metros. (3) Finalmente, esteja ciente de que PROJ_SCALE_FACTOR tem certos valores padrão para alguns
projeções, então você pode ter que substituir a configuração para corresponder aos resultados produzidos
com outras configurações.

Use mapprojectgmt online usando serviços onworks.net