Divisão de Observação da Terra e Geoinformática

Essa página está em modo somente de leitura. Você pode visualizar a fonte, mas não alterá-la. Informe-se com o administrador do Wiki, caso você ache que isso está incorreto.
====== Fusão de Imagens do Sensor WPM CBERS-4A no TerraView 5 ======
=== Autoria ===
Laércio M. Namikawa - laercio.namikawa@inpe.br


----


==== Objetivo ====
O objetivo é descrever o procedimento para colorir a imagem pancromática do sensor WPM com as bandas multiespectrais (Bandas IR, R, G e B).
As imagens do CBERS-4A podem ser obtidas gratuitamente do [[http://www.dgi.inpe.br/catalogo/|catálogo de imagens do INPE]].\\
O procedimento descreve a técnica utilizando a transformação por Componentes Principais. Nas transformações, as cores são obtidas das bandas multiespectrais de 8m, e na execução das transformações inversas, a Primeira Componente é substituída pela imagem pancromática (Banda 0) com 2 m de resolução espacial.

=== Obtendo o TerraView 5 ===
O TerraView Versão 5, nas versões Linux, MacOS e Windows, pode ser [[http://www.dpi.inpe.br/terralib5/wiki/doku.php?id=wiki:downloads|obtido livremente]] ("Free and Open Software").

----


==== Importando Imagens WPM ====
Utilize o menu Projeto -> Adiciona Camada -> Arquivo Raster e selecione os arquivos.

Ou arraste e solte sobre a janela do TerraView as bandas 0 a 4 do sensor WPM.

----

==== Recortando Imagens WMP 8m e e WPM 2m ====
Selecione a área de interesse e recorte as imagens para esta área menor. Eventualmente a área de interesse pode estar em duas cenas. neste caso, efetue o Mosaico das cenas.

  * Utilize o menu Projeto -> Processamento Matricial -> Recorte.
  * Selecione a camada.
  * Defina a área de recorte por Região de Interesse (área definida com cursor), Extensão do Layer (Retângulo Envolvente de uma camada - X e Y podem ser alterados), Camada (recorte por polígono definido na camada), ou Dimensão (Retângulo definido por linha, coluna e tamanho da imagem). Utilize o botão de Aplicar quando a opção de recorte selecionada apresentar este botão).
  * Definir diretório e nome do arquivo de Saída (Arquivo de saída será GeoTIFF).
  * Executar.


=== Opcional - Reamostrando A WPM de 8m para 4m com Interpolador Bicúbico ===
O registro pode ser realizado utilizando as bandas Multiespectrais reamostradas para 4m. O método de interpolação preferido é o bicúbico. Um efeito que pode ser gerado pelo bicúbico é a oscilação quando a transição é muito rápida. Para a interpolação, o TerraView tem a ferramenta apropriada - selecione o interpolador Bicúbico e o pixel de saída igual a 4m. Se a ferramenta não estiver disponível, a ferramenta de registro poder ser utilizada, definido o pixel de saída com 4m e o interpolador bicúbico.

----

==== Registrando a imagem WPM 2m com as imagens das bandas da WPM 8m ====
A imagem WPM 2m (B0) pode não ter uma coincidência ótima com a WPM 8m(B4 a 1). Mesmo com nível de correção L4, o registro entre as bandas e com a pancromática podem ser melhorados para área de interesse do recorte).

No TerraView, selecione uma banda WPM 8m (ou 4m caso tenha reamostrado) (se a WPM estiver em composição de bandas, utilize a função de decompor bandas para separar) como sendo a imagem a ser ajustada e a banda WPM 2m como sendo a imagem de referência. Tente utilizar a ferramenta de busca automática. Caso os erros dos pontos de controle sejam muito altos ou o número de pontos encontrados seja muito pequeno em determinada área, amplie para área e tente a ferramenta novamente. Uma vez que o mapeamento definido tem erro médio baixo (use de preferência a transformação Afim), o registro pode ser realizado, gerando uma nova imagem.
  * Selecione uma banda WPM 8m como sendo a imagem a ser ajustada.
  * Selecione a WPM 2m como sendo a imagem de referência.
  * Utiliza a ferramenta de busca automática.
  * Verifique o erro médio e os erros dos pontos de controle.
  * Caso sejam muito altos ou o número de pontos encontrados seja muito pequeno em determinada área:
  - Amplie para área.
  - Utiliza a ferramenta de busca automática.
  - Repita até que o mapeamento definido tenha erro médio baixo (use de preferência a transformação Afim)
  * Defina o nome da imagem WPM 8m registrada.
  * Execute o registro. (Repita para as outras bandas se não utilizou a composição).

----



==== Executando a Fusão ====
O procedimento de fusão no TerraView exige que as bandas da imagem colorida (WPM 8m neste caso) estejam compostas em uma imagem com várias bandas. A ferramenta Compõe/Decompõe executa esta tarefa. Esta imagem composta de bandas WPM 8m será fusionada com a imagem WPM 2m.

=== Executar a Fusão ===
  * Selecionar a camada com a imagem de Resolução Inferior (WPM 8m).
  * Selecionar a imagem pancromática como sendo a imagem de Resolução Superior (WPM 2m).
  * Definir o tipo de fusão.
  - A opção PCA fornece resultados melhores do que a IHS e é mais flexível que a WiSpeR (requer que a banda pancromática intersecte espectralmente com todas as bandas coloridas).
  * Selecionar a banda 0 da imagem de Resolução Maior (WPM 2m).
  * Selecionar as bandas da imagem de Resolução Menor (WPM 8m).
  * Definir Interpolador Bilinear ou Bicúbica (mais demorada).
  * Definir o diretório e nome do arquivo fusionado.
  * Executar a Fusão.

----