Eu tenho dois rasters: um modelo de elevação digital (DEM) e uma banda Landsat. A resolução da grade de ambos é de 30 m, mas as grades não são coincidentes (veja a imagem: a varredura escura no topo é o DEM; a varredura clara abaixo é a banda Landsat).
Eu quero realizar alguma transformação da banda Landsat para que sua grade se alinhe perfeitamente com a grade DEM. Tentei redimensionar a banda do Landsat para 30 m no ArcGIS usando o DEM como Snap Raster nas configurações do ambiente, mas isso não tem nenhum efeito. Também tentei recortar a banda do Landsat na extensão de um polígono de vetor (no qual o DEM também é recortado), enquanto o definia como Snap Raster, mas sem sucesso.
Eu suspeito que isso esteja acontecendo porque os dois rasters já estão na mesma resolução (apenas em grades diferentes), mas eu realmente preciso que suas grades estejam alinhadas para poder fazer qualquer análise adicional.
Posso fazer isso no ArcGIS Desktop?
fonte
Respostas:
Após mais pesquisas, descobri que a configuração incorreta do Snap Raster foi realmente causada por um bug no ArcGIS 10. A solução parece ser desativar o Geoprocessamento em Segundo Plano ( Opções de Geoprocessamento \ Geoprocessamento ... e desmarcar a opção Ativar na seção Processamento em Segundo Plano ), e execute a operação de recorte, com o Snap Raster especificado, no ArcCatalog.
Há uma discussão nos fóruns do ArcGIS aqui: http://forums.arcgis.com/threads/50808-Another-Snap-to-Raster-problem
Aparentemente, o bug foi corrigido no ArcGIS SP 4. Instalei o service pack, mas ainda não testei.
Curiosamente (ou melhor, irritante), se eu executar o clipe no meu polígono da área de interesse com o Snap Raster especificado (Snap Raster é DEM) conforme a solução na postagem do fórum acima, o Snap Raster funciona, mas o Landsat A banda possui mais uma coluna e mais uma linha que a varredura DEM, que também foi cortada no ArcGIS para o polígono da área de interesse. Isso é irritante porque, para uma análise mais aprofundada, preciso que os dois rasters tenham exatamente o mesmo número de linhas e colunas. Uma solução simples é executar um cálculo com a
Con()
função que cria uma nova varredura com o mesmo número de linhas e colunas que o DEM, mas que possui os valores da banda Landsat.fonte
Na verdade, lutei com esse problema por um longo tempo e descobri o que considero uma solução muito eficaz e bastante fácil. Escolha um dos seus rasters, aquele que você acredita estar no lugar "correto", aumente o zoom o máximo que o ArcGIS permitir e crie um novo shapefile de ponto diretamente no topo da interseção de quatro células raster. O único truque aqui é garantir que você aumente o zoom o máximo possível e crie o novo ponto diretamente em cima da 'cruz' formada por quatro células raster. Em seguida, registre as coordenadas do ponto que você acabou de criar. Finalmente, reprojete sua segunda varredura (mas mantenha inalterado o tamanho da célula de projeção e saída) e use as coordenadas do seu ponto recém-criado como o "Ponto de Referência". Isso alinhará perfeitamente seus dois rasters. Sem erros,
fonte
Você pode tentar usar a ferramenta de mascaramento para sobrepor o DEM com o Landsat.
Pode ser necessário criar uma caixa de área de estudo para capturar pixels que caem completamente um sobre o outro.
fonte
Meça o espaço entre uma varredura e a outra e, em seguida, mude uma varredura de acordo.
Você pode usar a ferramenta shift na caixa de ferramentas de gerenciamento de dados.
fonte
Eu tive o mesmo problema. Tentei a solução Louis Longchamps , que funcionou maravilhosamente em uma pequena amostra de teste, mas meu conjunto de dados era grande demais para usar praticamente esse método. Após um longo período de tempo e várias maneiras de tentar fazer isso, finalmente, por um capricho, usei o "Extrair por retângulo" após uma nova amostra, que funcionava exatamente como a função snap raster na ferramenta de clipe deveria ter funcionado. Eu configurei e tentei usar um construtor de modelos, mas isso me causou um erro; portanto, acabei tendo que fazer isso uma etapa de cada vez, mas, eventualmente, obtive todas as correspondências da minha varredura.
fonte
Eu não estava confortável com a ideia de mudar a varredura, porque, de acordo com mim, você altera as informações. Para evitar isso, usei uma sequência de conversões.
Aqui estão os passos:
Properties/Source
e anoteCellsize
em um pedaço de papel com todas as casas decimais disponíveis.Conversion Tools/Raster to Point
.Spatial Analyst Tools/Extraction/Extract Multi Values to Points
, especificando o recurso Ponto que você acabou de criar (na etapa 2) como o recurso de ponto de entrada e todas as rasters que deseja provar nas rasters de Entrada.Conversion Tools/To Raster/Point to Raster
, especificando qual camada você está criando no campo Valor e usando um nome de arquivo correspondente E oCellsize
valor que você anotou no pedaço de papel.Eles devem se alinhar perfeitamente e você pode transferi-los para R, por exemplo, e ter matrizes e vetores das mesmas dimensões exatas.
fonte
Estávamos tão preocupados com a nossa imagem DEM que ela também não estava alinhada com a imagem rasterizada classificada supervisionada. Também usamos o método de "Matthew Long" .
adicionamos os instantâneos de imagens não alinhadas e alinhadas.
Após o método de Matthew Long, nossa imagem DEM ficou totalmente alinhada com a nossa imagem rasterizada classificada supervisionada.
Esta é a imagem da ferramenta que usamos para alinhar nossa imagem DEM, no ARCGIS 10.2.1 nas ferramentas de gerenciamento de dados:
Esta é a imagem antes do alinhamento, a imagem preta é DEM e a imagem colorida é a imagem raster classificada supervisionada.
Esta é a imagem após o alinhamento, a imagem preta é DEM e a imagem colorida é supervisionada pela imagem raster classificada.
fonte
Eu resolvo meu problema com o método "Matthew Long". No entanto, para tornar ainda mais fácil e 100% de precisão. Sem qualquer zoom, basta usar a coordenada de extensão superior esquerda da varredura alinhada.
Para fazer isso: Clique com o botão direito do mouse, propriedade, painel Origem na propriedade Extensão
Além disso, usando a coordenada de extensão superior esquerda, é possível criar uma solução python sem intervenções manuais.
fonte
Eu tive um problema semelhante ao ter dois arquivos de fontes diferentes. Primeiro, um arquivo rasterizado para o nível de referência do desmatamento no nível do país (Peru). Segundo, um arquivo consolidado na plataforma GEE foi baixado e reprojetado no sistema do primeiro elemento. Em conclusão, a sobreposição de ambos me mostrou que havia um deslocamento entre as células (30 * 30 células) não superior a 5 metros. Para alinhar e ajustar o segundo arquivo ao primeiro, ocorreu-me estabelecer uma * base estendida do primeiro arquivo que possa cobrir ambos os arquivos * e manter o alinhamento de pixels do primeiro arquivo com seu sistema de projeção e outros recursos . Assim, no momento de reprojetar meu segundo arquivo ou qualquer subproduto com um deslocamento semelhante,
fonte