Como fusionar uma RGB + PANCROMATICA

[Felipe Dantas]

Olá, amigos.

Utilizo as imagens landsat 8 para um teste e gostaria de saber como faço para fusionar a pancromatica (15m) em uma composição colorida falsa cor (30m), tendo como resultado final um resolução espacial de 15m ?

[Rodrigo Moreira Maia]

nobre.. primeiro vc deve " esticar a de 30m " georreferenciando a de 30 com basenos dados de15, depois vc verá que a imagem " esticou"
depois vc domina faz o rgb e ja era....

Motorola XOOM™ Android 4.0.1

Em 16/06/2013 18:11, "Felipe Dantas" <felipera...@gmail.com> escreveu:

Olá, amigos.

Utilizo as imagens landsat 8 para um teste e gostaria de saber como faço para fusionar a pancromatica (15m) em uma composição colorida falsa cor (30m), tendo como resultado final um resolução espacial de 15m ?

--
--
---
Você recebeu esta mensagem do Grupo "QGISBrasil"
Para postar neste grupo, envie um e-mail para qgisb...@googlegroups.com
Para cancelar a sua inscrição neste grupo, envie um e-mail paragisbrasi...@googlegroups.com
Para ver mais opções, visite este grupo em
http://groups.google.com/group/qgisbrasil?hl=pt?hl=pt-BR
 
---
Você está recebendo esta mensagem porque se inscreveu no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar a inscrição neste grupo e parar de receber seus e-mails, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para obter mais opções, acesse https://groups.google.com/groups/opt_out.
 
 

[Rodrigo Moreira Maia]

Felipe cara me ajuda, agora é minha vez de te perguntar...
to usando o qgis 1.9, e o RGB não funfa..
ta usando o 1.9 também ????
se sim, como vc faz RGB nele ?????

deu certo esse esquema de esticar a imagem de 30metros para a de 15 metros por pixel ??? ?????

Em 16-06-2013 18:11, Felipe Dantas escreveu:

Olá, amigos.

Utilizo as imagens landsat 8 para um teste e gostaria de saber como faço para fusionar a pancromatica (15m) em uma composição colorida falsa cor (30m), tendo como resultado final um resolução espacial de 15m ?

[Felipe Dantas]

Rodrigo ou outros, me corrijam se eu estiver errado, mas esse método que tu falou não é muito certo, isso estar me parecendo mais uma gabiarra (sem despreso). Acredito que no Sextante, deva existir um plugin que automatize isso que estou querendo, mas pela busca que fiz não obtive resultado. Então, se alguém souber, por favor, compartilhem seus conhecimentos.

Cara, uso  qgis 1.8, mas fiz um teste aqui o qgis 1.9 pelo gdal (miscelania>mosaico) e rodou legal, sem problema algum !

--

Felipe Ramos Dantas

Graduando em Tecnologia  em Geoprocessamento - IFPI

Contatos: tim (86) 9921-5742 /  claro (86) 9411-4372

Respeite o ciclista.

Uma bike a mais, um carro a menos.
..... __@
...._ \ >_
...(_)/ (_)___________________________

_________
░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░\__

[Rodrigo Moreira Maia]

    Nem precisa corrigir, vc ta certo, é gambiarra,
por uma razão simples, existe menos informações na imagem de 30m com relação a de 15m,
então a única forma de preservar o "tamanho", não é correto essa forma para essa operação.. é assim mesmo, esticando a imagem...

pois, voce vê mais detalhes na imagem de 15 metros pp,  a de 30 tem menos, 4X menos pra ser mais exato..
porém existem pontos que são comuns e identificáveis nas imagens....

sacou o que acontece com esse tipo de operação ?????
sabe responder minha pergunta nobre amigo..
RGB no qgis 1.9 ja que o plugin zicou ????

[Felipe Dantas]

Hum, entendi, vou tentar aqui e ver os resultados e te falo no que deu !

Rodrigo, teste aqui pelos pluguins do sextante, realmente nao funfou... teste outras ferramentas aqui, também nao funfou, sendo que está configurado corretamente,...

Cria um tópico aí pra ver se aparece alguem que já passou por isso e o resolveu !

[Marcelo Bastos]

Oi Felipe

Fusionar imagens com diferentes resoluções deve ser feito a partir de transformações nos espaços RGB->IHS.

Abraço.

2013/6/16 Felipe Dantas <felipera...@gmail.com>

Olá, amigos.

Utilizo as imagens landsat 8 para um teste e gostaria de saber como faço para fusionar a pancromatica (15m) em uma composição colorida falsa cor (30m), tendo como resultado final um resolução espacial de 15m ?

--

[Rodrigo Moreira Maia]

NUsss, valeu mano, vou tentar...

obrigado..

[Marcelo Bastos]

Oi Felipe

Não uso o QGIS para processamento de imagens (só para SIG) mas esta rotina funciona na maioria dos softwares que trabalham com PDI, incluindo Corel, Photoshop...

Primeiro faz uma composção colorida RGB (321, cores verdadeiras, ou qualquer que de sua escolha: 543, 432, dependendo Landsat/IRS-P6 ou CBERS que não tem a banda 5 mas tem ótimos pares CCD/HRC com pixel de 2,5m com passagem na mesma data durante os anos de 2008 e 2009, além de ser grátis)

Passa a composição RGB escolhida por ti para o espaço IHS (RGB->IHS).

Vais gerar três novas bandas IHS (intensidade-saturação e não me lembro do nome da outra componente...)

Volta para o espaço RGB (IHS->RGB), utilizando as bandas S e H da primeira transformação e coloca na banda I a imagem de melhor definição/resolução espacial.

Pede para a resolução da imagem final (RGB) ficar com o mesmo pixel da banda pancromática ou de melhor resolução que tu vais colocar no lugar a banda I resultante da primeira transformação e que será descartada.

Não sei se fui muito didático. Acho que tu encontras um roteiro melhor no WIKI-INPE.

Abraço.

2013/6/17 Felipe Dantas <felipera...@gmail.com>

Olha, sou graduando em Geoprocessamento, já to quase pra terminar o curso, e sinceramente nunca ouvir falar disso. 

Detalhe( minha professora não era muuuuuuito didática e não frequentava muito a instituição kkkk, ela gostava mesmo era de viajar pros "isteites" kk)

Marcelo, se tu puder, me explica aí como funciona essa transformação de RGB-->IHS ?

como faço pelo qgis ?

[Hebert Guilherme de Azevedo]

Essa função, se não me engano, está disponível no GRASS.

Att

Em 17/06/13, Marcelo Bastos<marcelo.almei...@gmail.com> escreveu:

>>>> *Felipe Ramos Dantas**
>>>> *
>>>> *Graduando em Tecnologia em Geoprocessamento - IFPI**
>>>> *
>>>> *Contatos: tim (86) 9921-5742* */ claro (86) 9411-4372*
>>>> *
>>>> Respeite o ciclista.
>>>> *
>>>> *Uma bike a mais, um carro a menos.
>>>> *..... __@

>>>> ...._ \ >_
>>>> ...(_)/ (_)___________________________
>>>> _________
>>>> ░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░\__

>>> *Felipe Ramos Dantas**
>>> *
>>> *Graduando em Tecnologia em Geoprocessamento - IFPI**
>>> *
>>> *Contatos: tim (86) 9921-5742* */ claro (86) 9411-4372*
>>> *
>>> Respeite o ciclista.
>>> *
>>> *Uma bike a mais, um carro a menos.
>>> *..... __@

>>> ...._ \ >_
>>> ...(_)/ (_)___________________________
>>> _________
>>> ░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░\__

>> *Felipe Ramos Dantas**
>> *
>> *Graduando em Tecnologia em Geoprocessamento - IFPI**
>> *
>> *Contatos: tim (86) 9921-5742* */ claro (86) 9411-4372*
>> *
>> Respeite o ciclista.
>> *
>> *Uma bike a mais, um carro a menos.
>> *..... __@

>> ...._ \ >_
>> ...(_)/ (_)___________________________
>> _________
>> ░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░\__
>>

Hebert Guilherme de Azevedo
(21) 7573-5705 - Claro
(31) 7531-2575 - Tim

Skype: hebert.guilherme.de.azevedo

lattes: http://lattes.cnpq.br/6582724988395217

[Marcelo Bastos]

Em anexo, paper do XIIISBSR com referência para a rotina descrita.

2013/6/17 Hebert Guilherme de Azevedo <azevedo...@gmail.com>

[Felipe Dantas]

Obrigado, Marcelo. 

No artigo que tu mencionou, o autor usou o SPRING, porém quero bater cabeça nos QGIS e seus pluguin para o processamento de imagens..

Como o herbert falou, há sim um pluguin no grass que transforma o rgb > ihs e vice-versa. Na primeira operação a ferramenta me pede  as bandas (r-g-b) separadamente (eu pensava que primeiro teria que gerar o RGB e na ferramente ele iria "funfar") e como saida as imagens, separadamente, I - H - S.

Como tu mencionou, para fazer a voltar e fusionar com a pancromatica, eu teria que informar a camada H e S e no lugar da intensidade eu usaria a PAN. Porem, a ferramente me pede com saída as camdas RGB separadamente, e depois de executado as imagem geradas aparecem toda branca, sem nenhuma informação de pixel.

Concerteza tem algo errado, mas nao consegui descobrir. 

Vlw

Felipe Ramos Dantas

Graduando em Tecnologia  em Geoprocessamento - IFPI

Contatos: tim (86) 9921-5742 /  claro (86) 9411-4372

Respeite o ciclista.

Uma bike a mais, um carro a menos.

..... __@
...._ \ >_
...(_)/ (_)___________________________

_________
░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░\__

[Sidney - Geosaber]

Felipe,

Se você tiver tempo e disposição pode tentar fazê-lo pelo Orfeo Toolbox, no QGIS via SEXTANTE.

http://www.orfeo-toolbox.org/SoftwareGuide/SoftwareGuidech13.html

Aqui um exemplo de alguém que fez o Pansharpening de imagem do Landsat 8 pelo OTB, pena que ele não publicou o procedimento:

http://agmapsonline.com/?p=48

Boa sorte,

Sidney Goveia

http://www.geosaber.com.br

[Sidney - Geosaber]

Aqui tem o procedimento para se fazer pelo Orfeo Toolbox, não é muito detalhado mas acho que ajudará em seus estudos:

http://www.mapbox.com/blog/pansharpening-satellite-imagery-openstreetmap/

[Felipe Dantas]

Sidney,

Dei uma olhada nos arquivos que você postou, e tentei realizar o procedimento.

Porém, quando mandei executar a ferramenta, ele me pediu um arquivo e não sei que arquivo é esse (veja o anexo).

Não sei se preciso, antes, realizar algum pré-processamento, mas só executando a ferramenta (Bundle To Perfect Sensor), ele já vai me gerar a imagem fusionada ?

--

--
---
Você recebeu esta mensagem do Grupo "QGISBrasil"
Para postar neste grupo, envie um e-mail para qgisb...@googlegroups.com
Para cancelar a sua inscrição neste grupo, envie um e-mail paragisbrasi...@googlegroups.com
Para ver mais opções, visite este grupo em
http://groups.google.com/group/qgisbrasil?hl=pt?hl=pt-BR
 
---
Você está recebendo esta mensagem porque se inscreveu no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar a inscrição neste grupo e parar de receber seus e-mails, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para obter mais opções, acesse https://groups.google.com/groups/opt_out.
 
 

[Sidney - Geosaber]

Para maiores detalhes você terá que consultar o CookBook (página 154): http://www.orfeo-toolbox.org/packages/nightly/latest/CookBook-3.17.0.pdf

Geoid File: Use a geoid grid to get the height above the ellipsoid in case there is no DEM available, no coverage for some points or pixels with no data in the DEM tiles.

A version of the geoid can be found on the OTB website (http://hg.orfeo-toolbox.org/OTB-Data/raw-file/404aa6e4b3e0/Input/DEM/egm96.grd).

Você tem bastante trabalho pela frente, terá que ler atentamente o manual e testar para conseguir alcançar seus objetivos.

Infelizmente eu estou com muito pouco tempo para poder testar e ver quais os problemas dessa operação.

Boa sorte com seus estudos.

Sidney Goveia

http://www.geosaber.com.br

[Felipe Dantas]

Valeu, Sidney.

Agora é tentar entender esse negócio.

Quando consegui, posto aqui os resultado. Porém, se alguém conseguir de outra maneira dentro do qgis, posta aqui também.

Vlw, pessoal. 

Sidney Goveia

http://www.geosaber.com.br

--

--
---
Você recebeu esta mensagem do Grupo "QGISBrasil"
Para postar neste grupo, envie um e-mail para qgisb...@googlegroups.com
Para cancelar a sua inscrição neste grupo, envie um e-mail paragisbrasi...@googlegroups.com
Para ver mais opções, visite este grupo em
http://groups.google.com/group/qgisbrasil?hl=pt?hl=pt-BR
 
---
Você está recebendo esta mensagem porque se inscreveu no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar a inscrição neste grupo e parar de receber seus e-mails, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para obter mais opções, acesse https://groups.google.com/groups/opt_out.
 
 

[Sidney - Geosaber]

Fez-se a luz!

Abri a imagem do Landsat 8 no Monteverdi (RGB432 - visível e 8 Pan) banda por banda;

• File > Open Dataset.

Concatenei as bandas multiespectrais (empilhamento de bandas - igual ao mosaico > camada acumulada do QGIS);

• File > Concatenate Images.

Depois usei a opção de superimposição das cenas (multi e pan);

• Geometry > Superimpose two images (aqui, pelo que entendi do manual é para fazer o corregistro entre a multi e a pan, como opcional pode-se usar um DEM para ortorretificação das cenas).

E por fim, fiz a fusão das imagens (Pansharpening).

• Filtering > Pansharpening > Simple RCS pan-sharpening

E consegui acender uma vela de sabedoria para iluminar a escuridão de ignorância*.

O passo 2 será fazer isso funcionar dentro do QGIS, mas aí serão outras noite sem dormir!

Boa sorte,

Sidney Goveia

http://www.geosaber.com.br

*Carl Sagan, 1997. O Mundo Assombrado Pelos Demônios: A Ciência Vista Como Uma Vela No Escuro.

[Felipe Dantas]

Não sei se é meu notebook que não presta, ou se fiz algo de errado, porque toda vez que faço isso no monteverdi (concatenação, orto-registo, e fusão) ele trava na hora que coloco pra visualizar a fusão....  oh my god !

Mas valeu, sidney, continuar tentando aqui, uma hora eu consigo rs !

abraço.

--

--
---
Você recebeu esta mensagem do Grupo "QGISBrasil"
Para postar neste grupo, envie um e-mail para qgisb...@googlegroups.com
Para cancelar a sua inscrição neste grupo, envie um e-mail paragisbrasi...@googlegroups.com
Para ver mais opções, visite este grupo em
http://groups.google.com/group/qgisbrasil?hl=pt?hl=pt-BR
 
---
Você está recebendo esta mensagem porque se inscreveu no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar a inscrição neste grupo e parar de receber seus e-mails, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para obter mais opções, acesse https://groups.google.com/groups/opt_out.
 
 

[Felipe Dantas]

Não sei se é meu notebook que não presta, ou se fiz algo de errado, porque toda vez que faço isso no monteverdi (concatenação, orto-registo, e fusão) ele trava na hora que coloco pra visualizar a fusão....  oh my god !

A imagem que você usou, tu deixaste no sistema de coordenadas original ou transformou pra zona utm sul (caso seja no hemisferio sul) ?

Mas valeu, sidney, continuar tentando aqui, uma hora eu consigo rs !

abraço.

--

--
---
Você recebeu esta mensagem do Grupo "QGISBrasil"
Para postar neste grupo, envie um e-mail para qgisb...@googlegroups.com
Para cancelar a sua inscrição neste grupo, envie um e-mail paragisbrasi...@googlegroups.com
Para ver mais opções, visite este grupo em
http://groups.google.com/group/qgisbrasil?hl=pt?hl=pt-BR
 
---
Você está recebendo esta mensagem porque se inscreveu no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar a inscrição neste grupo e parar de receber seus e-mails, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para obter mais opções, acesse https://groups.google.com/groups/opt_out.
 
 

[Sidney - Geosaber]

A imagem que eu usei é uma das primeiras a ser disponibilizada para download, é de uma área na Califórnia, uma cena completa, sem recortes.

O problema dos crashs é questão do cache, o Monteverdi faz cache (na memória) de todas as imagens e também durante o processamento, eu fiz todo o meu processo sem carregar nada no Viewer, na minha máquina (i5, 4GB, win7pro64) o processamento mais pesado chegava a quase 2 minutos.

Agora quero fazer funcionar no QGIs.

Até,

Sidney Goveia

http://www.geosaber.com.br

[Felipe Dantas]

Reinstalei o monteverdi, e refiz os procedimentos, desta vez deu certo, não travou o programa, porém a imagem fusionda que eu exportei, ficou com o tamanho de 1,25 GB. A sua ficou com este tamanho também, sidney, ou não chegaste a exportar ?

--

--
---
Você recebeu esta mensagem do Grupo "QGISBrasil"
Para postar neste grupo, envie um e-mail para qgisb...@googlegroups.com
Para cancelar a sua inscrição neste grupo, envie um e-mail paragisbrasi...@googlegroups.com
Para ver mais opções, visite este grupo em
http://groups.google.com/group/qgisbrasil?hl=pt?hl=pt-BR
 
---
Você está recebendo esta mensagem porque se inscreveu no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar a inscrição neste grupo e parar de receber seus e-mails, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para obter mais opções, acesse https://groups.google.com/groups/opt_out.
 
 

[Felipe Dantas]

Segue as imagens, já no qgis, porém todo procedimentos foi feito no monteverdi.

[Sidney - Geosaber]

Felipe,

Permaneço trabalhando pelo Monteverdi porque pelo Orfeo Toolbox (OTB/SEXTANTE/QGIS) está difícil de configurar, e tem ocorrido muitos "paus" tanto com o OTB quanto o GDAL na versão Master (1.9).

As cenas do Landsat 8 são grandes (maior número de bandas que o L7 e ainda por cima tem resolução radiométrica de 16bits, que era de 8bits no L7 e L5).

E as imagens resultantes dos processamentos ficam bem grandes mesmo, eu crio tudo em cache (raster virtual) e depois eu salvo (export dataset) com resolução radiométrica de 8bits, assim facilita na hora de comparar com imagens do L7 e L5 no QGIS (ficam com o mesmo nível de contraste).

Está para ser lançada uma versão 2 do Monteverdi, por enquanto tem um beta sem muitos recursos, mas ele vai ficar mais agradável de trabalhar.

Como o desenvolvimento do QGIS, do SEXTANTE e do OTB são desconexos (são grupos separados) é difícil conseguir buscar soluções para os problemas que eles apresentam.

[Tadeu Pontes]

Bom dia a todos,

Estou tentando fazer a fusão RGB+pancromática de uma imagem landsat-8, uso o QGIS 2.6 e o windows 8.1, mas estou tendo algumas dúvidas.

1 - baixei o spectral transformer for landsat-8, mas quando vou abrir a GUI aparece um erro de configuração lado a lado,alguém teve esse problema? Alguma sugestão de como resolver? Achei várias sugestões na internet de atualizar o .NET framework, mas a versão que tenho é a mais atual.

2 - Como faço para adicionar a Orfeo toolbox e outras opções na caixa de ferramentas? 

Att,

Tadeu

[Tadeu Pontes]

Compartilhando a solução pra primeira pergunta, baixei a versão do C++ que consta no link: http://www.microsoft.com/en-au/download/details.aspx?id=26347 e deu tudo certo.

[Felipe Bernardi]

Um ótimo tutorial:

https://www.youtube.com/watch?v=0vg66if5KPE