Mais

Transformando dados usando 2D Costum Transformation no Arcmap 10.2 for Desktop

Transformando dados usando 2D Costum Transformation no Arcmap 10.2 for Desktop


Eu tenho camadas de ArcMap de estradas, edifícios etc. em um sistema de referência (Israel_TM). Eu preciso transformar todos esses dados (enquanto mantenho a topologia) para o "sistema de referência" do espaço de pixels. Já calculei os parâmetros de transformação (usando Matlab) de uma Transformação Conformada (2D) entre os dois sistemas de referência: Israel_TM para Pixels. É uma transformação de 4 parâmetros que tenho: 1xRotation, 1xScale, 2xTranslations. Como posso criar uma transformação 2D definida pelo usuário e usá-la para transformar todos os meus dados em espaço de pixels? Tudo que consegui encontrar foi Helmert 7-Par. transformação ou outros métodos de transformação 3D em Propriedades da camada -> Sistema de Coordenadas -> Transformação. Tentei usar a caixa de ferramentas para isso, mas não consegui encontrar nada relevante.


Não é possível transformar o sistema de coordenadas do British National Grid para WGS_OSGB_1984 Mercat

Vou começar dizendo que sou muito novo nisso, e embora eu tenha feito o curso introdutório ao ArcGIS Desktop 10.0, ele não está realmente ajudando com o meu problema.

Eu tenho reunido informações de várias partes interessadas e marcadas em um mapa (gráfico raster). Para analisar as informações, importei os mesmos mapas de fundo e os dados coletados para um novo documento de mapa, no entanto, eles não estão alinhados. Parece que a razão para isso é que no software que estamos usando (um questionário especialmente desenvolvido, o CS foi definido como British National Grid, mas o raster cartografado é WGS_OSGB_1984 Mercator - quando importei pela primeira vez os mapas para o pacote do questionário me informou que havia necessidade de transformar as coordenadas que eu fiz.

Pesquisei isso até onde pude e, aparentemente, há uma transformação de BNG para WGS84 disponível para 9.1, mas não consigo encontrar nada para 10.0. Todas as transformações fornecidas no Arc 10 são inúteis, e passei dias tentando calcular uma transformação mylef, mas ela está indo por todo o lugar sem nenhum padrão lógico real.

Eu acrescentaria que trabalho para uma organização sem fins lucrativos e, portanto, não temos orçamento para pagar um desenvolvedor para resolver isso, mas se você puder me ajudar, será um herói absoluto!

por MelitaKennedy

A varredura está em WGS_1984_Web_Mercator ou WGS_1984_Web_Mercator_Auxiliary_Sphere? Isso fará a diferença. O primeiro usa um sistema de coordenadas geográficas personalizado (WGS 1984 modificado, "GCS_WGS_1984_Major_Auxiliary_Sphere) e não há transformação direta entre ele e o OSGB 1936. Você tem que usar duas transformações, o que é. Difícil. No ArcMap, mas pode ser feito convertendo o dados no ArcToolbox utilizando as ferramentas Project ou Project Raster O 2º utiliza o padrão WGS 1984 e deverão estar disponíveis várias transformações.

Obrigdo por sua contribuição. O sistema de coordenadas é descrito como WGS_1984_Mercator, não há menção web_mercator ou web_mercator_auxillary_sphere.

Tentarei usar as transformações que você mencionou, mas qualquer informação adicional que você possa fornecer será bem-vinda.

por MelitaKennedy

Sim, essa é uma definição genérica de Mercator (centrada em 0,0) que adicionamos como um caso de teste. Se eu adicionar dados WGS84 ao ArcMap e definir o sistema de coordenadas do quadro de dados para British National Grid, recebo esta lista de transformações com as áreas de uso e precisões adicionadas:

OSGB_1936_To_WGS_1984_1 "Reino Unido - Grã-Bretanha Ilha de Man" 21.0
OSGB_1936_To_WGS_1984_2 "RU - Inglaterra" 10.0
OSGB_1936_To_WGS_1984_3 "Reino Unido - Inglaterra, País de Gales, Ilha de Man" 21.0
OSGB_1936_To_WGS_1984_4 "RU - Escócia" 18.0
OSGB_1936_To_WGS_1984_5 "Reino Unido - País de Gales" 35.0
OSGB_1936_To_WGS_1984_NGA_7PAR "Reino Unido - Grã-Bretanha Ilha de Man" 21.0
OSGB_1936_To_WGS_1984_Petroleum "UK - Grã-Bretanha Ilha de Man" 2.0

Ao verificar possíveis transformações, apenas os sistemas de coordenadas geográficas importam.

Eu tentei as transformações listadas para transformar meus dados de partes interessadas (reunidos no GCS 1936 / British National Grid) em 1984 (usados ​​para o mapa de fundo) e transformar os dados do mapa em 1936 / BNG GCS. Coloquei um marcador em um ponto distinto para tentar alinhá-los sem muito sucesso, abaixo estão as coordenadas de cada conjunto de transformações para mostrar a escala do problema. Não sei se isso vai permitir que você me ajude mais, mas é tudo em que posso pensar no momento (a lógica não está me levando muito longe).

Transformando dados de 1936 em 1984 (todos os valores em metros)

Ponto de referência -475663.757, 6844471.977
Transformação 1 -283557.103, 4080120.052
Transformação 2 -283551.527, 4080094.476
Transformação 3 -283602.731, 4080103.425
Transformação 4 -283547.172, 4080068.931
Transformação 5 -283547.172, 4080107.315
Transformation_NGA_7PAR -283568.237, 4080103.425
Transformation_Petroleum -283568.237, 4080068.931

Transformando os dados de 1984 em 1936 (todos os valores em metros)

Ponto de referência 245573.383, 283857.525
Transformação 1 215822.76, 2396583.867
Transformação 2 215825.84, 2396594.223
Transformação 3 215825.84, 2396594.223
Transformação 4 215825.84, 2396610.354
Transformação 5 215809.709, 2396578.091
Transformation_NGA_7PAR 215825.84, 2396642.617
Transformation_Petroleum 215825.24, 2396642.617

Eu já havia tentado trabalhar minha própria transformação usando esse ponto de referência e consegui chegar bem perto (em cerca de um km), mas acho que havia um problema de escala ou rotação (estava usando uma transformação xy simples). Quando tento desenvolver uma transformação mais complexa, os dados parecem pular para todos os lados, sem nenhum padrão lógico.


Criação de um arquivo de transformação geográfica personalizado

Um arquivo de transformação geográfica customizado é criado usando a ferramenta Criar Transformação Geográfica Customizada. Este arquivo .gtf é criado em uma das duas pastas, dependendo das permissões de gravação da sua conta de usuário.

  • No Windows, o arquivo normalmente é salvo em C: Users & ltUSERNAME & gt AppData Roaming ESRI Desktop10.X ArcToolbox CustomTransformations. Este caminho pode ser acessado diretamente digitando% appdata% ESRI Desktop10.X ArcToolbox CustomTransformations no Windows Explorer. (Certifique-se de atualizar o X no caminho para corresponder à sua versão do ArcGIS for Desktop.)
  • Se você não tiver permissão de gravação para o caminho acima (ele foi definido como somente leitura), o arquivo será criado na pasta temporária da sua conta de usuário. Isso pode ser acessado digitando% temp% ArcToolbox CustomTransformations no Windows Explorer.
  • Em um sistema Unix, o .gtf será salvo no diretório inicial do usuário em ArcToolbox CustomTransformations, na pasta Application Data ou na pasta tmp, novamente, dependendo de suas permissões de gravação.

Consulte as notas de uso da ferramenta Criar transformação geográfica personalizada para obter informações mais detalhadas sobre onde este arquivo foi salvo.


Nordpil

O melhor conjunto de dados disponível gratuitamente sobre placas tectônicas, orógenos e limites de placas foi publicado em 2003 por Peter Bird (antiga Universidade da Califórnia).

Os dados foram publicados originalmente no jornal Um modelo digital atualizado de limites de placa (Geochemistry Geophysics Geosystems, 4 (3), 1027, doi: 10.1029 / 2001GC000252, 2003). O conjunto de dados está disponível para download gratuito e irrestrito, atualmente no site privado do pesquisador. No entanto, ele está disponível apenas em um formato de texto complicado, o que significa que não é facilmente acessível para uso em aplicativos GIS.

Para trazer este conjunto de dados para a era moderna, os dados originais foram analisados, limpos e verificados usando ArcGIS 10.2 e convertidos em arquivos de forma.

Todos os dados originais foram recuperados de http://peterbird.name/oldFTP/PB2002/ em junho de 2014 e analisados ​​usando Python e Global Mapper 11. A manipulação de dados adicional foi realizada no ArcGIS 10.2. As principais edições consideraram segmentos que abrangem o limite -180/180, que tiveram que ser divididos e movidos manualmente.

Esta coleção está disponível sob o Licença de atribuição Open Data Commons, o que impõe poucas restrições ao uso dos dados.

O conjunto de dados apresenta placas tectônicas e seus limites, bem como orógenos e informações sobre os limites. Os dados são úteis para aplicações geológicas, análise e educação, e devem ser fáceis de usar em qualquer aplicativo de software GIS moderno. Para informações sobre os campos e valores, consulte a documentação original e o artigo científico.


1.3 Tutorial interativo

O resto do workshop foi dedicado a um exercício prático em que vivemos codificados juntos. Os dados que usamos foram dados de fronteira do distrito aldermanic de Chicago de 1986, que podem ser encontrados na página da Coleção de Mapas de UChicago aqui.

O roteiro deste workshop pode ser encontrado aqui.

1.3.1 Funções aprendidas

1.3.2 Importar dados espaciais

A primeira coisa que faremos é importar os dados espaciais. Nossos dados espaciais são armazenados como um arquivo de forma (.shp, mas na verdade 4 arquivos). Eles estão ficando menos na moda (consulte GeoJSON, etc.), mas muitos dados espaciais ainda são armazenados dessa forma.

Primeiro carregue a biblioteca sf. Se você não tiver, instale-o em seu console ou no painel Packages RStudio.

Usaremos a função st_read. Também lê bancos de dados GeoJSON, PostGIS e muito mais.

Você pode usar o atalho de teclado Ctrl-Enter para executar uma linha de código em R.

Pressione tab após digitar st_read (& quot & quot) com o cursor entre as aspas e você obterá um bom recurso de preenchimento automático.

Verifique em qual projeção os dados estão:

Sike, isso não é projetado! Você pode dizer porque proj4string começa com + proj = longlat. Você ainda pode plotar isso, mas as coisas começarão a ficar duvidosas se você tentar fazer cálculos de distância ou área, ou plotar esses limites de ala com outras camadas.

Ainda posso traçar o mapa, mas como um bom analista de dados geográficos, vou precisar projetá-lo.

1.3.3 Projete seus dados

Precisamos projetar os dados: mas como escolhemos qual projeção precisamos? Esta é uma questão filosófica profunda, mas felizmente a biblioteca UChicago já nos disse as melhores projeções para usar para nossos dados: Illinois East State Plane ou Zona 16 UTM.

Eu geralmente faço uma pesquisa online para procurar os códigos EPSG para a projeção que desejo usar ou uso spatialreference.org, que tem um banco de dados para todos os códigos EPSG para projeções.

Depois de fazer algumas investigações, descobri que o código EPSG que desejo para UTM Zona 16 é 32616. Em seguida, uso st_transform () para projetar os dados e salvá-los como um novo dataframe sf.

Eu verifico o CRS - sim, parece certo!

Pergunta: em que dado está isso? Quais são as unidades para esta projeção?

  • Projete os dados da ala de 1986 na projeção do plano do estado de Illinois. Siga as etapas acima para fazer isso. Dica: procure o código EPSG online!

Observe que + unidades = us-ft parte da proj4string significa que quaisquer cálculos de distância que você fizer nos dados ward86_stateplane serão feitos em pés.

Lembre-se de que existem projeções “ruins” para seus dados. Por exemplo, se eu acidentalmente escolher a projeção Alaska Albers para meus dados ...

Obviamente, essa projeção não é boa para Chicago! Mas algo para ter em mente, se, digamos, você estiver fazendo mapas de lugares mais próximos dos pólos. `` `


Editando os Parâmetros do Transformer

Usando um conjunto de opções de menu, os parâmetros do transformador podem ser atribuídos referenciando outros elementos na área de trabalho. Funções mais avançadas, como um editor avançado e um editor aritmético, também estão disponíveis em alguns transformadores. Para acessar um menu dessas opções, clique ao lado do parâmetro aplicável. Para obter mais informações, consulte Opções do menu de parâmetros do transformador.

Definindo Valores

Existem várias maneiras de definir um valor para uso em um Transformer. O mais simples é simplesmente digitar um valor ou string, que pode incluir funções de vários tipos, como referências de atributos, funções matemáticas e de string e parâmetros do espaço de trabalho. Existem várias ferramentas e atalhos que podem ajudar na construção de valores, geralmente disponíveis no menu de contexto suspenso adjacente ao campo de valor.

Usando o Editor de Texto

O Editor de Texto fornece uma maneira conveniente de construir strings de texto (incluindo expressões regulares) a partir de várias fontes de dados, como atributos, parâmetros e constantes, onde o resultado é usado diretamente dentro de um parâmetro.

Usando o Editor Aritmético

O Editor Aritmético fornece uma maneira conveniente de construir expressões matemáticas de várias fontes de dados, como atributos, parâmetros e funções de recurso, onde o resultado é usado diretamente dentro de um parâmetro.

Valores Condicionais

Defina os valores de acordo com uma ou mais condições de teste aprovadas ou reprovadas.

Contente

Expressões e strings podem incluir várias funções, caracteres, parâmetros e muito mais.

Ao definir valores - sejam inseridos diretamente em um parâmetro ou construídos usando um dos editores - strings e expressões contendo String, Math, Date / Time ou FME Feature Functions terão essas funções avaliadas. Portanto, os nomes dessas funções (na forma @ & ltfunction_name& gt) não devem ser usados ​​como valores de string literais.


Métodos de autenticação multifator

Os usuários são obrigados a se autenticar com todos ou alguns destes fatores:

Autenticação de fator único: A autenticação baseada em senha é o tipo mais comumente usado de autenticação de fator único, considerada insegura, pois a maioria das pessoas usa senhas que são facilmente comprometidas ou são fracas.

Autenticação de dois fatores (2FA): O usuário deve fornecer um código por meio de um aplicativo de autenticação em seu smartphone, um token de software (& ldquosoft token & rdquo), um token de hardware (& ldquohard token & rdquo) ou uma mensagem SMS. Este código pode ser uma senha descartável que expira após 30 segundos, ou uma senha HMAC descartável, ou HOTP, que não expira até que seja usada.

Autenticação de três fatores (3FA): Fatores biométricos são levados em consideração com a autenticação de três fatores & ndash, algo que é exclusivo para cada usuário. Isso inclui reconhecimento facial, geometria de mão, geometria de dedo, varredura de impressão palmar e varredura de impressão digital.

Autenticação de quatro fatores (4FA): Na autenticação de quatro fatores, uma localização geográfica do usuário e rsquos, junto com a quantidade de tempo que leva para o usuário chegar ao seu destino, é levada em consideração. Esses tipos de medidas de autenticação normalmente não exigem que os usuários validem essas informações. Em vez disso, ele é executado em segundo plano enquanto o risco de autenticação está sendo determinado. A autenticação de quatro fatores confirma a geolocalização e geo-velocidade do usuário, bem como sua localização no mundo. Por exemplo, haveria um problema de segurança se um usuário se autenticasse em um provedor de autenticação multifator em Chicago e tentasse autenticar novamente em Moscou dez minutos depois.

Autenticação de cinco fatores (5FA): Na autenticação de cinco fatores, um usuário pode gerar padrões de toque ou gestos específicos. Por exemplo, um usuário pode usar o sistema operacional Windows mais recente (8 ou posterior) para criar uma senha de imagem usando uma tela de toque de computador, onde pode tocar em uma imagem ou desenhar linhas retas ou círculos para gerar uma senha de gesto exclusiva.

BAutenticação iométrica: A biometria de autenticação está se tornando mais comum, à medida que os recursos de digitalização de impressão digital e reconhecimento facial nos produtos Android e Apple cresceram em popularidade entre os consumidores. Existem outros métodos biométricos que ainda não foram autorizados para uso comercializado de autenticação multifator, como reconhecimento de digitação, análise de assinatura e caligrafia, padrões de veias, padrões de marcha, identidade de odor, varreduras de íris, varreduras de retina, impressões de voz e formato de orelha reconhecimento.


Essa animação é feita sob medida em Flash, que é uma ferramenta perfeita e completa para fazer animações 2D (e até jogos) como essa, mas tem um custo bastante alto, mesmo para a licença de estudante ($ $ 19,99 $ por mês ou $ $ 199 $ por ano).

GeoGebra é uma ferramenta gratuita para desenhar geometria, com recursos para fazer interações, animações e uma versão beta para gráficos 3D. É uma das melhores ferramentas para fazer geometria, tem uma versão online bastante limitada e uma versão completa instalada (ambas gratuitas).

Cut-the-knot Geoboard é um miniaplicativo fantástico (interações, animações), e tem um custo muito confortável ($ $ 5 $) para uma licença.

Para fazer gráficos 3D arquitetônicos de forma rápida e simples, recomendo a versão gratuita do Sketchup, que também é uma ferramenta fantástica para fazer desgins e também é possível fazer animações em 3D.

Como uma ferramenta de design 3D totalmente desenvolvida, eu realmente recomendo o Blender, um programa gratuito para fazer 3D no qual você tem total capacidade artística, mas por isso não possui ferramentas específicas de geometria.

Para geometria 2D:

Existe http://www.robocompass.com/app - gratuito e fantástico e faz construções geométricas animadas.

Para geometria 3D:

Desenvolvi um programa de desenho 3D gratuito que roda em navegadores modernos (html5). A versão em inglês é chamada Geoservant 3D, a versão original em alemão é chamada Geoknecht 3D.

Lá você pode especificar formas 3D (geometria sólida) e animá-las alterando seus valores. Veja a descrição na página como alterar os valores: mantenha pressionada a tecla SHIFT CTRL ALT e use as teclas do cursor. Também dê uma olhada na galeria na página para ver o que você pode conseguir com este programa usando uma sintaxe simples. Para cada objeto, o programa calcula valores básicos (área, volume, vetores).

Para cores diferentes, basta adicionar o código de cores HEX no final, por exemplo, esfera (3 | 3 | 1 2) . Se você adicionar um valor entre 0 e 1 entre colchetes, você pode definir uma transparência para cada objeto, por exemplo, esfera (3 | 3 | 1 2)

Além disso, as variáveis ​​podem ser definidas. Exemplo:

Um hemisfério pode ser criado como essa esfera (0 | 0 | 4 4 0 6,28 1,55 3,14), veja a galeria no site!

Idiomas: Geoservant está agora disponível em alemão, inglês, espanhol, francês, russo, sueco e chinês - e aceita os objetos especificados nesses idiomas, por exemplo: 长方体 (4 | 1 | 1 2 | 4 | 1)

Alternativo

Mais complexo e poderoso: acho a ferramenta de geometria 3D tinkercard.com bastante incrível.

Outro programa é o Geodrafter 2D (inglês) ou Geozeichner 2D (alemão) para formas bidimensionais.


Sintaxe de script

NetCDFConvert2DVariableToArcGISRaster_GeoEco (inputFile, outputRaster, variableName, xLowerLeftCorner, yLowerLeftCorner, cellSize, nodataValue, transpor, mirror, flip, swapHemispheres, coordinateSystem, projectedCoordinateSystem, geográfico de transferência, recortePontoDeRegistro, registro de dados geográficos, gráficoCoordenado, compilação e transferência de dados geográficos).

No momento em que esta ferramenta foi desenvolvida, duas versões principais do netCDF estavam em desenvolvimento: NetCDF-3, que estava na versão 3.6.2, e NetCDF-4, que estava na versão 4.0 beta 1. A implementação atual desta ferramenta pode processar Arquivos NetCDF-3, mas não arquivos NetCDF-4. O suporte a NetCDF-4 será adicionado em um futuro próximo.

Se você fornecer um arquivo compactado em um formato de compactação compatível, ele será descompactado automaticamente. Se for um arquivo (por exemplo, .zip ou .tar), ele deve conter exatamente um arquivo, que não deve estar em um subdiretório.

Nome de uma variável no arquivo netCDF. A variável deve ter duas dimensões e um tipo de dados inteiro ou de ponto flutuante.

Se você não souber o nome da variável, dê seu melhor palpite. Se não existir uma variável com esse nome, um ValueError será gerado. A mensagem de erro listará os nomes de todas as variáveis ​​2D presentes no arquivo de entrada. Selecione um desses nomes e tente novamente.

Você também pode usar a ferramenta Extrair Cabeçalho NetCDF para despejar o cabeçalho netCDF em um arquivo de texto. Você pode encontrar os nomes das variáveis ​​procurando no arquivo a seção que começa com "variáveis:". Cada variável é listada como:

Onde tipo de dados é o tipo de dados da variável (por exemplo, "float"), nome é o nome da variável, e dimensões é a lista de dimensões da variável. Por exemplo, o trecho abaixo de um netCDF de correntes geostróficas do Aviso DT-MADT lista duas variáveis ​​que têm duas dimensões, Grid_0001 e Grid_0002:

Esta ferramenta extrai a variável netCDF para um ArcInfo ASCII Grid e a converte em um raster usando a ferramenta de geoprocessamento ArcGIS ASCII to Raster. Essa ferramenta e o formato raster do ArcGIS têm várias limitações que restringem os tipos de variáveis ​​netCDF que podem ser convertidas com sucesso para o formato raster do ArcGIS.

As variáveis ​​NetCDF que usam o tipo de dados float ou double não devem conter valores "infinito" (INF) ou "não é um número" (NAN). Um ValueError será gerado se esses valores forem descobertos.

O formato raster do ArcGIS suporta o tipo de dados float de 32 bits, mas não o tipo de dados duplos de 64 bits. Se você fornecer uma variável netCDF com o tipo de dados double, ela será convertida para um raster flutuante de 32 bits usando a ferramenta de geoprocessamento ArcGIS ASCII to Raster. O comportamento da ferramenta nesta situação não está documentado. No ArcGIS 9.1 parece ser:

Valores onde o expoente varia de -38 a +38 são representados corretamente no raster flutuante de 32 bits resultante, embora alguma precisão seja perdida devido à mantissa menor do tipo de dados flutuante de 32 bits.

Valores em que o expoente é menor que -38 (por exemplo, -39, -40 e assim por diante) são convertidos para 0.

Valores onde o expoente é maior que +38 são convertidos em -INF ou + INF, dependendo do sinal do valor (por exemplo, -5.3083635279597874e-212 aparece como -1. # INF na GUI ArcCatalog, enquanto 2.5502286890301497e + 084 aparece como 1. # INF).

A ferramenta ArcGIS 9.1 ASCII para Raster também apresenta algumas peculiaridades ao converter rasters inteiros:

Para dados inteiros de 8 bits, a ferramenta criará raster de 16 bits se o valor -128 aparecer, a menos que -128 seja designado o valor NODATA. Especificar um valor NODATA diferente, como 0, ainda produz uma varredura de 16 bits se -128 aparecer.

Da mesma forma, para dados inteiros de 16 bits, a ferramenta criará raster de 32 bits se o valor -32768 aparecer no arquivo ASCII, a menos que seja designado o valor NODATA.

Pior, para dados inteiros de 32 bits, a ferramenta relatará um erro se o valor -2147483648 aparecer, a menos que seja designado o valor NODATA. Ainda mais estranho, o valor -2147483647 é sempre traduzido para NODATA, não importa o quê.

Para todos os tipos de dados inteiros, a ferramenta produz um comportamento estranho quando você especifica um valor NODATA que não é o menor valor possível para o tipo de dados. Por exemplo, se o intervalo de dados de 0 a 255 e 0 for designado como o valor NODATA, a ferramenta produzirá um raster de saída de 8 bits sem sinal. Mas se 1 for designado o valor NODATA, ele produzirá um raster de saída de 16 bits com sinal e ArcCatalog mostra em Propriedades do conjunto de dados raster que o valor NoData é -32768, embora a ferramenta Identify mostre células que tinham o valor 1 são, na verdade, NODATA. Resultados estranhos semelhantes podem ser obtidos para rasters inteiros de outros tipos de dados, quando você designa um NODATA que não é o menor valor possível.

Coordenada X do canto esquerdo inferior do raster.

A coordenada é para o canto da célula inferior esquerda, não o centro dessa célula. Por exemplo, se o raster é uma projeção geográfica de toda a Terra, a coordenada do canto esquerdo inferior seria -180,0, correspondendo a uma longitude de 180 graus Oeste.

Coordenada Y do canto inferior esquerdo do raster.

A coordenada é para o canto da célula inferior esquerda, não o centro dessa célula. Por exemplo, se o raster for uma projeção geográfica de toda a Terra, a coordenada do canto esquerdo inferior seria -90,0, correspondendo a uma latitude de 90 graus ao sul.

Por exemplo, se o raster é uma projeção geográfica de toda a Terra, com 720 colunas e 360 ​​linhas, ele teria um tamanho de célula de 0,5, correspondendo a 1/2 de um grau geográfico.

O formato de dados subjacente requer que as células sejam quadradas. Não é possível especificar um tamanho de célula para cada dimensão.

Valor que indica que uma célula não possui dados.

Se for True, a imagem será transposta (invertida em torno do eixo diagonal) antes da conversão. Use esta opção para corrigir uma imagem que tem o eixo leste / oeste indo para cima e para baixo em vez de esquerda e direita.

Se for True, a imagem será invertida no eixo vertical antes da conversão. Use esta opção para corrigir uma imagem que é a "imagem espelhada" do que deveria ser.

Se for True, a imagem será invertida no eixo horizontal antes da conversão. Use esta opção para corrigir uma imagem que está de cabeça para baixo.

Se for True, os hemisférios leste e oeste da imagem serão trocados. Use esta opção para alterar a orientação de uma imagem global de uma orientação de 0 a 360 centrada no oceano Pacífico para uma orientação de -180 a +180 centrada no oceano Atlântico, ou vice-versa.

Sistema de coordenadas a definir para o raster de saída. Se um valor não for fornecido, o sistema de coordenadas do raster de saída permanecerá indefinido.

Novo sistema de coordenadas para projetar o raster de saída.

O raster só pode ser projetado para um novo sistema de coordenadas se a projeção original for definida. Um erro será gerado se você especificar um novo sistema de coordenadas sem definir o sistema de coordenadas original.

A ferramenta ArcGIS Project Raster é usada para realizar a projeção. A documentação dessa ferramenta recomenda que você também especifique um tamanho de célula para o novo sistema de coordenadas.

Percebi que, para certos sistemas de coordenadas, a ferramenta ArcGIS 9.2 Project Raster parece cortar o raster projetado em uma extensão arbitrária que é muito pequena. Por exemplo, ao projetar uma imagem global de clorofila MODIS Aqua de 4 km em coordenadas geográficas para Lambert_Azimuthal_Equal_Area com meridiano central de -60 e latitude de origem de -63, a imagem resultante é cortada para mostrar apenas um quarto do planeta. Este problema não ocorre quando o Project Raster é chamado interativamente a partir da interface do usuário do ArcGIS, ele ocorre apenas quando a ferramenta é chamada programaticamente (o método ProjectRaster_management do geoprocessador). Portanto, você pode não ver quando usa o Project Raster, mas pode acontecer quando você usa ferramentas MGET que invocam o Project Raster como parte de suas operações de geoprocessamento.

Se você encontrar esse problema, pode contorná-lo desta forma:

Primeiro, execute esta ferramenta sem especificar um novo sistema de coordenadas, para obter o raster de saída no sistema de coordenadas original.

No ArcCatalog, use a ferramenta Project Raster para projetar o raster para o novo sistema de coordenadas. Verifique se todo o raster está presente, se não foi cortado em uma extensão muito pequena.

No ArcCatalog, procure a extensão do raster projetado clicando com o botão direito do mouse na árvore do catálogo, selecionando Propriedades e rolando para baixo até Extent.

Agora, antes de executar a ferramenta MGET que projeta o raster, defina a configuração do ambiente Extent para os valores que você pesquisou. Se você estiver chamando a ferramenta MGET interativamente de ArcCatalog ou ArcMap, clique no botão Ambientes na caixa de diálogo da ferramenta, abra Configurações Gerais, altere a lista suspensa Extensão para "Conforme Especificado Abaixo" e digite os valores que você pesquisou. Se você o está chamando de um modelo de geoprocessamento, clique com o botão direito do mouse na ferramenta no modelo, selecione Criar variável, Do ambiente, Configurações gerais, Extensão. Isso colocará Extent como uma variável em seu modelo, anexada à ferramenta MGET. Abra a variável Extent, altere-a para "As Specified Below" e digite os valores que você pesquisou. Se você estiver chamando a ferramenta MGET programaticamente, deve definir a propriedade Extent do geoprocessador para os valores pesquisados. Consulte a documentação do ArcGIS para obter mais informações sobre isso e as configurações de ambiente em geral.

Execute a ferramenta MGET. A extensão da varredura de saída agora deve ter o tamanho adequado.

Um método de transformação usado para converter entre o sistema de coordenadas original e o novo sistema de coordenadas.

Este parâmetro é uma nova opção introduzida pelo ArcGIS 9.2. Você deve ter ArcGIS 9.2 para usar este parâmetro.

Este parâmetro é necessário apenas quando você especifica que o raster deve ser projetado para um novo sistema de coordenadas e que o novo sistema usa um datum diferente do sistema de coordenadas original, ou há alguma outra diferença entre os dois sistemas de coordenadas que requer uma transformação. Para determinar se uma transformação é necessária, recomendo o seguinte procedimento:

Primeiro, execute esta ferramenta sem especificar um novo sistema de coordenadas, para obter o raster de saída no sistema de coordenadas original.

Em seguida, use a ferramenta ArcGIS 9.2 Project Raster no raster de saída para projetá-lo no sistema de coordenadas desejado. Se uma transformação geográfica for necessária, essa ferramenta solicitará uma. Anote o nome exato da transformação que você usou.

Finalmente, se uma transformação for necessária, digite o nome exato nesta ferramenta, execute-a novamente e verifique se o raster de saída foi projetado conforme desejado.

O algoritmo de reamostragem a ser usado para projetar o raster original em um novo sistema de coordenadas. A ferramenta ArcGIS Project Raster é usada para realizar a projeção e aceita os seguintes valores:

NEAREST - interpolação do vizinho mais próximo

BILINEAR - interpolação bilinear

Você deve especificar um desses algoritmos para projetar para um novo sistema de coordenadas. Um erro será gerado se você especificar um novo sistema de coordenadas sem selecionar um algoritmo.

O tamanho da célula do sistema de coordenadas projetadas. Embora este parâmetro seja opcional, para receber os melhores resultados, a documentação do ArcGIS recomenda que você sempre o especifique ao projetar para um novo sistema de coordenadas.

As coordenadas xey (no espaço de saída) usadas para alinhamento de pixel.

Este parâmetro é uma nova opção introduzida pelo ArcGIS 9.2. Você deve ter ArcGIS 9.2 para usar este parâmetro. Ele será ignorado se você não especificar que a imagem deve ser projetada em um novo sistema de coordenadas.

Retângulo ao qual o raster deve ser recortado.

Se um sistema de coordenadas projetadas foi especificado, o recorte é executado após a projeção e as coordenadas do retângulo devem ser especificadas no novo sistema de coordenadas. Se nenhum sistema de coordenadas projetadas foi especificado, as coordenadas devem ser especificadas no sistema de coordenadas original.

A ferramenta ArcGIS Clip é usada para executar o clipe. O retângulo de recorte deve ser passado para esta ferramenta como uma string de quatro números separados por espaços. A interface de usuário do ArcGIS formata automaticamente a string de maneira adequada ao chamar esta ferramenta da UI do ArcGIS, você não precisa se preocupar com o formato. Mas, ao invocá-lo programaticamente, tome cuidado para fornecer uma string formatada corretamente. Os números são ordenados ESQUERDO, INFERIOR, DIREITO, TOPO. Por exemplo, se o raster estiver em um sistema de coordenadas geográficas, ele pode ser cortado em 10 W, 15 S, 20 E e 25 N com a string:

Podem ser fornecidos números inteiros ou decimais.

Expressão de álgebra do mapa para executar no raster de saída.

AVISO: O ArcGIS Geoprocessing Model Builder pode excluir aleatoriamente e silenciosamente o valor deste parâmetro. Este é um bug no ArcGIS. Antes de executar um modelo que você salvou, abra esta ferramenta e valide se o valor do parâmetro ainda existe.

A expressão é executada depois que o raster convertido é projetado e cortado (se essas opções forem especificadas). Use a string inputRaster que faz distinção entre maiúsculas e minúsculas para representar o raster sobre o qual você agora deseja executar a álgebra de mapa. Por exemplo, para converter a varredura em uma varredura inteira e adicionar 1 a todas as células, use esta expressão:

A string inputRaster diferencia maiúsculas de minúsculas. Prior to executing the map algebra expression, the string is replaced with the path to a temporary raster that represents the output raster being generated. The final expression must be less than 4000 characters long or ArcGIS will report an error.

The ArcGIS Single Output Map Algebra tool is used to execute the map algebra expression. You must have a license for the ArcGIS Spatial Analyst extension in order to perform map algebra.

Map algebra syntax can be very picky. Here are some tips that will help you succeed with this tool:

Before using this tool, construct and test out your map algebra expression using the ArcGIS Single Output Map Algebra tool. Then paste the expression into this tool and edit it to use the inputRaster variable rather than the test value you used with Single Output Map Algebra.

If you do develop your expression directly in this tool, start with a very simple expression. Verify that it works properly, add a little to it, and verify again. Repeat this process until you have built up the complete expression.

Always separate mathematical operators from raster paths using spaces. In the example above, the / operator contains a space on either side. Follow this pattern. In some circumstances, ArcGIS will fail to process raster algebra expressions that do not separate raster paths from operators using spaces. The reported error message usually does not indicate that this is the problem, and tracking it down can be very frustrating.

If True, pyramids will be built for the output raster, which will improve its display speed in the ArcGIS user interface. This is the last step performed in post-conversion processing.


Can we use tf.spectral fourier functions in keras?

Let us start with an input that is a simple time series and try to build an autoencoder that simply fourier transforms then untransforms our data in keras.

Then the third line throws an error when entered:

You see, the output of tf.spectral.irfft is float32 but it looks like Lambda thinks it is complex64?? (Complex64 is the input x from the previous step)

We can fix that error at model entry time with:

This is accepted at input time but then when we try to build the model:

Which I guess is reasonable and was the reason I didn't want to cast it in the first place.

Main question: how do I successfully wrap the tf.spectral.irfft function which outputs float32 ?

More general question for learning: Let's assume I actually want to do something between the rfft and the irfft, how can I cast those imaginary numbers into absolute values without breaking keras so I can apply various convolutions and the like?


Assista o vídeo: A Complete Beginners Guide to ArcGIS Desktop Part 1