Differences
This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
fusaopan10m5mtv:exemplo [2019/07/29 16:08] laercio [Resultado] |
fusaopan10m5mtv:exemplo [2019/07/29 17:18] tkorting [Fusão por Componentes Principais] |
||
---|---|---|---|
Line 1: | Line 1: | ||
- | ====== Fusão de Imagens do Satélite CBERS-2B no SPRING ====== | + | ====== Fusão de Imagens do Satélite CBERS-4 no TerraView 5 ====== |
=== Autoria === | === Autoria === | ||
Laércio M. Namikawa - laercio.namikawa@inpe.br | Laércio M. Namikawa - laercio.namikawa@inpe.br | ||
Line 19: | Line 19: | ||
==== Importando Imagens PAN10M e PAN5M ==== | ==== Importando Imagens PAN10M e PAN5M ==== | ||
- | Utilize o menu Projeto-> Adiciona Camada - > Arquivo Raster e selecione os arquivos. | + | Utilize o menu Projeto -> Adiciona Camada -> Arquivo Raster e selecione os arquivos. |
Ou arraste e solte sobre a janela do TerraView. | Ou arraste e solte sobre a janela do TerraView. | ||
Line 28: | Line 28: | ||
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. | 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. | + | * Utilize o menu Projeto -> Processamento Matricial -> Recorte. |
* Selecione a camada. | * 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). | * 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). | ||
Line 57: | Line 57: | ||
==== Executando a Fusão ==== | ==== Executando a Fusão ==== | ||
O procedimento de fusão no TerraView exige que as bandas da imagem colorida (PAN10M 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 PAN10M será fusionada com a imagem PAN5M. | O procedimento de fusão no TerraView exige que as bandas da imagem colorida (PAN10M 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 PAN10M será fusionada com a imagem PAN5M. | ||
- | |||
- | === Compor as Bandas === | ||
- | A imagem HRC cobre uma área 25 vezes menor que a da imagem CCD ([[http://www.cbers.inpe.br/pt/programas/cbers1-2_cameras.htm|veja aqui]]). Durante a reamostragem , além de reduzir a resolução do pixel das bandas CCD restauradas de 10 para 2.5 metros, as bandas devem ser recortadas para o mesmo retângulo envolvente da imagem HRC. O interpolador bilinear é utilizado para a reamostragem para garantir uma transição suave de cores entre pixels. | ||
=== Executar a Fusão === | === Executar a Fusão === | ||
- | * Selecionar a banda 0 da imagem de Resolução Maior (PAN5M). | + | * Selecionar a camada com a imagem de Resolução Inferior (PAN10M). |
- | * Selecionar a imagem pancromática (PAN5M). | + | * Selecionar a imagem pancromática como sendo a imagem de Resolução Superior (PAN5M). |
* Definir o tipo de fusão. | * 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). | - 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). | ||
Line 112: | Line 109: | ||
|Processamento->Processamento Matricial->Compõe/Decompõe\\ Selecione as camadas com as bandas PAN10M\\ Clique no botão **Avançar**|{{:fusaopan10m5mtv:compoe1.png?400|}}| | |Processamento->Processamento Matricial->Compõe/Decompõe\\ Selecione as camadas com as bandas PAN10M\\ Clique no botão **Avançar**|{{:fusaopan10m5mtv:compoe1.png?400|}}| | ||
- | |Selecione a aba Compõe Bandas\\ Selecione a imagem da componente 0 (Camada com a PAN10M Banda 4)\\ Selecione a imagem da componente 1 (Camada com a PAN10M Banda 3)\\ Selecione a imagem da componente 2 (Camada com a PAN10M Banda 2)\\ Selecione a banda 0 de cada uma das imagens (as imagens de entrada têm somente uma banda)\\ Clique no botão **Avançar**|{{:fusaopan10m5mtv:compoe2.png?400|}}| | + | |Selecione a aba Compõe Bandas\\ Selecione a imagem da componente 0\\ (Camada com a PAN10M Banda 4)\\ Selecione a imagem da componente 1\\ (Camada com a PAN10M Banda 3)\\ Selecione a imagem da componente 2\\ (Camada com a PAN10M Banda 2)\\ Selecione a banda 0 de cada uma das imagens\\ (as imagens de entrada têm somente uma banda)\\ Clique no botão **Avançar**|{{:fusaopan10m5mtv:compoe2.png?400|}}| |
|Clique no botão para definir nome da saída\\ Definir diretório e nome do arquivo recortado\\ Clique no botão **Terminar**|{{:fusaopan10m5mtv:compoe3.png?400|}}| | |Clique no botão para definir nome da saída\\ Definir diretório e nome do arquivo recortado\\ Clique no botão **Terminar**|{{:fusaopan10m5mtv:compoe3.png?400|}}| | ||
|Visualizando a imagem composta em RGB default|{{:fusaopan10m5mtv:compoe4.png?400|}}| | |Visualizando a imagem composta em RGB default|{{:fusaopan10m5mtv:compoe4.png?400|}}| | ||
Line 120: | Line 117: | ||
==== Fusão por Componentes Principais ==== | ==== Fusão por Componentes Principais ==== | ||
- | O método calcula as bandas Componentes Principais e troca a Primeira Componente pela imagem PAN5M, inserindo a resolução espacial de 5 metros nas bandas espectrais da PAN10M. | + | O método calcula as bandas Componentes Principais e troca a \\ Primeira Componente pela imagem PAN5M, inserindo a resolução espacial de 5 metros \\ nas bandas espectrais da PAN10M. |
- | |Processamento->Processamento Matricial->Fusão\\ Selecione a camada com imagem composta PAN10M como sendo a camada com resolução inferior\\ Clique no botão **Avançar**\\ Selecione a camada com imagem PAN5M como sendo a camada com resolução superior\\ Clique no botão **Avançar**|{{:fusaopan10m5mtv:fusao1.png?400|}}| | + | |Processamento->Processamento Matricial->Fusão\\ Selecione a camada com imagem composta PAN10M \\ como sendo a camada com resolução inferior\\ Clique no botão **Avançar**\\ Selecione a camada com imagem PAN5M\\ como sendo a camada com resolução superior\\ Clique no botão **Avançar**|{{:fusaopan10m5mtv:fusao1.png?400|}}| |
- | |Selecione a aba Resolução Menor\\ Observe que a banda 0 é a disponível e está selecionada|{{:fusaopan10m5mtv:fusao2.png?400|}}| | + | |Selecione a aba Resolução Maior\\ Observe que a banda 0 é a disponível e está selecionada|{{:fusaopan10m5mtv:fusao2.png?400|}}| |
- | |Selecione a aba Resolução Maior\\ Selecione as 3 bandas PAN10M a serem fusionadas |{{:fusaopan10m5mtv:fusao3.png?400|}}| | + | |Selecione a aba Resolução Menor\\ Selecione as 3 bandas PAN10M a serem fusionadas.\\ Defina o interpolador Bilinear |{{:fusaopan10m5mtv:fusao3.png?400|}}| |
|Clique no botão para definir nome da saída\\ Definir diretório e nome do arquivo recortado\\ Clique no botão **Terminar**|{{:fusaopan10m5mtv:fusao4.png?400|}}| | |Clique no botão para definir nome da saída\\ Definir diretório e nome do arquivo recortado\\ Clique no botão **Terminar**|{{:fusaopan10m5mtv:fusao4.png?400|}}| | ||
Line 132: | Line 129: | ||
+ | ==== Resultado ==== | ||
+ | Imagem PAN5M pancromática.\\ Resolução espacial de 5 metros.\\ {{:fusaopan10m5mtv:exemplofusao4.png?800|}} \\ \\ | ||
- | ==== Transformando RGB para IHS ==== | ||
- | As cores RGB são obtidas de uma composição colorida de bandas CCD reamostradas e transformadas nas componentes de **I**ntensidade, **S**aturação e **M**atiz.\\ | ||
- | Para transformar, execute a seqüência: | ||
- | |Visualize uma\\ composição colorida\\ de bandas CCD\\ reamostradas e inicie\\ a ferramenta de\\ transformação através do menu\\ **Imagem** - **Transformação\\ IHS <-> RGB**, |{{fusaohrcccdcbers2b:ihsrgb:tp200955145518.png?400|}}| | ||
- | |Selecione transformação\\ **RGB->IHS**, clique\\ no botão **R** para\\ definir o PI com\\ a imagem em vermelho,\\ selecione a imagem\\ com a Banda 3\\ reamostrada e clique\\ no botão **Executar**|{{fusaohrcccdcbers2b:ihsrgb:tp200956145609.png?200|}}{{fusaohrcccdcbers2b:ihsrgb:tp200956145614.png?200|}}| | ||
- | |Repita para o verde com o botão **R**,\\ selecionando a imagem\\ com a Banda 4\\ reamostrada e cliquando\\ no botão **Executar**|{{fusaohrcccdcbers2b:ihsrgb:tp200956145644.png?400|}}| | ||
- | |Repita para o azul com o botão **B**,\\ selecionando a imagem\\ com a Banda 2\\ reamostrada e cliquando\\ no botão **Executar**|{{fusaohrcccdcbers2b:ihsrgb:tp200957145728.png?400|}}| | ||
- | |Selecione o tamanho\\ do pixel da imagem\\ de saída como sendo\\ igual ao da imagem\\ em vermelho (R) e\\ clique no botão **Executar**|{{fusaohrcccdcbers2b:ihsrgb:tp200958145814.png?200|}}| | ||
+ | Imagem PAN10M em composição de bandas: R4G3B2.\\ Resolução espacial de 10 metros.\\ | ||
- | ---- | ||
+ | {{:fusaopan10m5mtv:exemplofusao1.png?800|}} \\ \\ | ||
+ | Imagem PAN10M fusionada com a PAN5M.\\ Resolução espacial de 5 metros.\\ | ||
- | ==== Resultado ==== | + | {{:fusaopan10m5mtv:exemplofusao2.png?800|}} |
- | |Imagem PAM10M em composição de bandas: R4G3B2.// ResoluçÃo espacial de 10 metros|{{exemplofusao1.png?400|}}| | + | === Resultado === |
- | |Selecione transformação\\ **IHS->RGB**, clique\\ no botão **I** para\\ definir o PI com\\ a imagem intensidade,\\ selecione o PI com a imagem\\ HRC recortada e clique\\ no botão **Executar**|{{fusaohrcccdcbers2b:ihsrgb:tp200901150130.png?200|}}{{fusaohrcccdcbers2b:ihsrgb:tp200901150135.png?200|}}| | + | |
- | |Repita para o matiz\\ com o botão **H**,\\ selecionando a imagem\\ com a componente H\\ e cliquando\\ no botão **Executar**|{{fusaohrcccdcbers2b:ihsrgb:tp200901150154.png?400|}}| | + | |
- | |Repita para a saturação\\ com o botão **B**,\\ selecionando a imagem\\ com a componente S\\ e cliquando\\ no botão **Executar**|{{fusaohrcccdcbers2b:ihsrgb:tp200902150207.png?400|}}| | + | |
- | |Selecione o tamanho\\ do pixel da imagem\\ de saída como sendo\\ igual ao da imagem\\ intensidade (I) e\\ clique no botão **Executar**|{{fusaohrcccdcbers2b:ihsrgb:tp200902150239.png?200|}}| | + | |
- | |Imagem resultante|{{fusaohrcccdcbers2b:ihsrgb:tp200904150419.png?400|}}| | + | |
- | ---- | + | Para um resultado melhor, deve-se restaurar as bandas da PAN10M utilizando o SPRING. Veja sobre a Restauração [[http://wiki.dpi.inpe.br/doku.php?id=spring:restauracbers|aqui]].\\ |
- | ==== Métodos Alternativos ==== | + | Imagem PAN10M restaurada para 5 metros e fusionada com a PAN5M.\\ Resolução espacial de 5 metros.\\ |
- | Outros métodos podem ser usados como, por exemplo, o de [[spring:fusaopc|transformação por componentes principais]]. | + | |
- | ---- | + | {{:fusaopan10m5mtv:exemplofusao5.png?800|}} |
- | + | ||
- | ==== Exemplo da Fusão ==== | + | |
- | === Imagens Originais === | + | |
- | {{fusaohrcccdcbers2b:rgb_ori.jpg|Imagem RGB das bandas CCD 3, 4 e 2}} | + | |
- | {{fusaohrcccdcbers2b:mono.jpg|Imagem Monocromática da HRC}} | + | |
- | === Imagem RGB Resultante da Fusão === | ||
- | {{fusaohrcccdcbers2b:rgb.jpg|Imagem RGB resultante da fusão}} | ||
- | Para obter a imagem resultante inteira, acesse o {{fusaohrcccdcbers2b:20080710_154c_1271_hrc342.zip|link}}. |