Dúvida - bandas termais Landsat 8
434 views
Skip to first unread message
Sabrina Mieko
Aug 13, 2014, 8:25:43 AM8/13/14
to qgisb...@googlegroups.com
Bom dia pessoal,
Gostaria de saber se alguém aqui do grupo já trabalhou com as imagens
das bandas termais do Landsat 8 no QGIS.
Estou tentando, via calculadora raster, converter os DN em radiância,
e a radiância em graus Kelvin, e graus Kelvin em Celsius de uma cena
(conforme as equações da USGS -
http://landsat.usgs.gov/Landsat8_Using_Product.php), mas não está
dando certo, os resultados tem saído bem estranhos!
Abraços,
Sabrina Mieko
Gostaria de saber se alguém aqui do grupo já trabalhou com as imagens
das bandas termais do Landsat 8 no QGIS.
Estou tentando, via calculadora raster, converter os DN em radiância,
e a radiância em graus Kelvin, e graus Kelvin em Celsius de uma cena
(conforme as equações da USGS -
http://landsat.usgs.gov/Landsat8_Using_Product.php), mas não está
dando certo, os resultados tem saído bem estranhos!
Abraços,
Sabrina Mieko
Jorge Santos
Aug 13, 2014, 8:35:54 AM8/13/14
to qgisb...@googlegroups.com
Oi Sabrina,
Eu tenho um artigo de referência que utilizou a calculadora raster do ArcGIS para realizar conversões de DN (Digital Number) para Radiância, Radiância para Reflectância, Cálculo de NDVI e SAVI |
Landsat-5 NDVI, SAVI, Radiância, Reflectância, Distância Terra-Sol e outros assuntos | http://goo.gl/pbPYQL | #PDI #geoespacial
Verifique se as fórmulas são parecidas e compare os resultados.
Jorge Santos
Técnico em Geoprocessamento
http://www.processamentodigital.com.br/
http://www.processamentodigital.com.br/
Skype: jorgepsantos2002
Twitter: http://twitter.com/jpsantos2002
Facebook: http://www.facebook.com/ProcessamentoDigital
Twitter: http://twitter.com/jpsantos2002
Facebook: http://www.facebook.com/ProcessamentoDigital
--
--
---
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 inscrição nesse grupo e parar de receber e-mails dele, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para obter mais opções, acesse https://groups.google.com/d/optout.
Franklin Januário
Aug 13, 2014, 9:34:42 AM8/13/14
to qgisb...@googlegroups.com
Olá Sabrina,
Veja se os artigos do link abaixo lhe são úteis na tarefa, caso obtenha resultados positivos, poste aqui para a comunidade.
http://www.gisagmaps.com/landsat-8-data-tutorial/
Abs.
Veja se os artigos do link abaixo lhe são úteis na tarefa, caso obtenha resultados positivos, poste aqui para a comunidade.
http://www.gisagmaps.com/landsat-8-data-tutorial/
Abs.
Luana de Oliveira
Nov 23, 2014, 8:57:30 PM11/23/14
to qgisb...@googlegroups.com
Oi Sabrina, você conseguiu?
estou tendo o mesmo problema e não encontro solução.
E só tenho a opção de fazer com o Qgis...
Se puder me ajudar, agradeço.
Abraços,
Luana
Sabrina Mieko
Nov 24, 2014, 4:43:30 AM11/24/14
to qgisb...@googlegroups.com
Oi Luana,
Acabou que desde aquele e-mail, não consegui mais mexer nos arquivos. Tive de parar tudo pra atender outras demandas do pós doutorado... Mas o assunto está pendente, porque também ainda não achei solução pra ele, infelizmente. Na próxima semana vou ter uma reunião com uma professora de sig da UFSCar e vou ver se ela me dá uma luz!--
Comunidade QGIS Brasil: http://qgisbrasil.wordpress.com
---
Você recebeu essa mensagem porque está inscrito no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar inscrição nesse grupo e parar de receber e-mails dele, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para postar nesse grupo, envie um e-mail para qgisb...@googlegroups.com.
Acesse esse grupo em http://groups.google.com/group/qgisbrasil.
Herminio Sabino
Aug 23, 2016, 3:54:27 PM8/23/16
to qgisbrasil
boa tarde,
vi que faz um tempo a primeira postagem do tópico,
mas esses dias estava fazendo algumas pesquisas e teste no qgis para fazer o calculo de temperatura de superfície
e acabei encontrando um plugin muito útil.
se chama Geosud.
fiz o comparativo dos resultados da calculadora raster e plugin.
e gostei dos resultados.
Sabrina Mieko
Sep 26, 2016, 11:01:50 AM9/26/16
to qgisbrasil
Bom dia galera!
Venho eu desenterrar este tópico e dar um retorno pra algumas pessoas que me perguntaram sobre esse tema... Entre idas e vindas acabei deixando o trabalho com as bandas termais do Landsat 8 de lado por um tempo pra me dedicar a outras tarefas (nesse meio tempo tive uma filha, que ainda engatinha pela casa e corri com outras tarefas do pós-doc). Vou relatar aqui pra vocês minha experiência e perguntar se mais alguém está passando pelo mesmo perrengue.
Minhas primeiras tentativas na extração das temperaturas de superfície foram com uso da calculadora raster do QGIS. Recortei e converti as imagens que estavam em WGS84 23N, pra WGS84 23S e quando seguia os cálculos recomendados na página da USGS com os parâmetros que estão nos metafiles do pacote das imagens do Landsat8, começaram a sair resultados bizarros, com temperaturas de -120°C (algo improvável pro interior de São Paulo). Tentei juntando as fórmulas, separando as fórmulas, com a imagem recortada, imagem inteira, reprojetada, e até voltando pra original sem nenhuma alteração e nada. Isso tudo no QGis... Ainda não tive a oportunidade (e nem vou ter tempo agora) de testar em outro software.
A fórmula original está aqui: http://landsat.usgs.gov/documents/Landsat8DataUsersHandbook.pdf
A fórmula que transcrevi na calculadora raster foi essa (juntei as fórmulas pra transformar DN pra radiância e radiância pra °K, depois tentei separado tb): (1321.0789/Ln(774.8853/(0.0003342* “raster_B10” + 0.1)+1))
Até que eu resolvi usar o Semi Automatic Classification Plugin (que já tinha instalado pra reclassificar os valores de emissividade dos meus rasters), e tem uma função: Brightness temperature in Celsius e... Funcionou! (Fogos de artifício) O porém é que ele só obtém as temperaturas corretas se for feito com a imagem original. Tentei fazer com a imagem recortada e reprojetada e não deu certo.
https://media.readthedocs.org/pdf/semiautomaticclassificationmanual-v4/latest/semiautomaticclassificationmanual-v4.pdf
Gente, será que fiz algo errado na hora de calcular o raster ou escrever a fórmula? (considerando que não sou uma usuária experiente com manuseio de rasters no SIG e ainda estou aprendendo)
Abraços,
Sabrina Mieko
Venho eu desenterrar este tópico e dar um retorno pra algumas pessoas que me perguntaram sobre esse tema... Entre idas e vindas acabei deixando o trabalho com as bandas termais do Landsat 8 de lado por um tempo pra me dedicar a outras tarefas (nesse meio tempo tive uma filha, que ainda engatinha pela casa e corri com outras tarefas do pós-doc). Vou relatar aqui pra vocês minha experiência e perguntar se mais alguém está passando pelo mesmo perrengue.
Minhas primeiras tentativas na extração das temperaturas de superfície foram com uso da calculadora raster do QGIS. Recortei e converti as imagens que estavam em WGS84 23N, pra WGS84 23S e quando seguia os cálculos recomendados na página da USGS com os parâmetros que estão nos metafiles do pacote das imagens do Landsat8, começaram a sair resultados bizarros, com temperaturas de -120°C (algo improvável pro interior de São Paulo). Tentei juntando as fórmulas, separando as fórmulas, com a imagem recortada, imagem inteira, reprojetada, e até voltando pra original sem nenhuma alteração e nada. Isso tudo no QGis... Ainda não tive a oportunidade (e nem vou ter tempo agora) de testar em outro software.
A fórmula original está aqui: http://landsat.usgs.gov/documents/Landsat8DataUsersHandbook.pdf
A fórmula que transcrevi na calculadora raster foi essa (juntei as fórmulas pra transformar DN pra radiância e radiância pra °K, depois tentei separado tb): (1321.0789/Ln(774.8853/(0.0003342* “raster_B10” + 0.1)+1))
Até que eu resolvi usar o Semi Automatic Classification Plugin (que já tinha instalado pra reclassificar os valores de emissividade dos meus rasters), e tem uma função: Brightness temperature in Celsius e... Funcionou! (Fogos de artifício) O porém é que ele só obtém as temperaturas corretas se for feito com a imagem original. Tentei fazer com a imagem recortada e reprojetada e não deu certo.
https://media.readthedocs.org/pdf/semiautomaticclassificationmanual-v4/latest/semiautomaticclassificationmanual-v4.pdf
Gente, será que fiz algo errado na hora de calcular o raster ou escrever a fórmula? (considerando que não sou uma usuária experiente com manuseio de rasters no SIG e ainda estou aprendendo)
Abraços,
Sabrina Mieko
Em quarta-feira, 13 de agosto de 2014 09:25:43 UTC-3, Sabrina Mieko escreveu:
Marcelo Bastos
Sep 27, 2016, 5:53:36 PM9/27/16
to qgisbrasil
Oi Sabrina, tem que cuidar a resolução radiométrica (número de bits) da imagem de saída (reprojetada).
O cálculo (constantes) só funciona para a resolução da imagem original.
Caso diminua a resolução (como é comum ocorrer e a gente não perceber - raster/propriedades) de 16 para 8 bits, não dá para usar a primeira equação (DN para radiãncia) pois os parâmetros de calibração do sensor são para as imagens brutas, apenas corrigidas geometricamente pelo método do vizinho mais próximo (sem alteração do nível de cinza original) que é o caso das imagens do L8.
Abraço.
Sabrina Mieko
Sep 30, 2016, 10:28:36 AM9/30/16
to qgisb...@googlegroups.com
Oi Marcelo!
Eu pensei nisso quando estava fazendo os testes, por isso que acabei voltando as imagens brutas pra refazer os cálculos e mesmo assim as temperaturas finais acabavam em números estranhos (-120°C no interior de SP em dezembro... Estranho!). Mesmo agora que acertei as conversões, pra não dar erro, acabei deixando pra fazer os recortes no final com as imagens já com as infos de temperatura.Sabrina Mieko Viana
PhD in Forest Resources - ESALQ-USP
Brazil
PhD in Forest Resources - ESALQ-USP
Brazil
--
Comunidade QGIS Brasil: http://qgisbrasil.org/
---
Você recebeu essa mensagem porque está inscrito em um tópico no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar inscrição nesse tópico, acesse https://groups.google.com/d/topic/qgisbrasil/z0vK0jZ1KRU/unsubscribe.
Para cancelar inscrição nesse grupo e todos os seus tópicos, envie um e-mail para qgisbrasil+unsubscribe@googlegroups.com.
Acesse esse grupo em https://groups.google.com/group/qgisbrasil.
Jéssica P dos Santos
Nov 17, 2017, 10:53:27 AM11/17/17
to qgisbrasil
Oi Sabrina,
Você fez o download das bandas termais direto pelo site do INPE ou pelo site do USGS?
Você fez o download das bandas termais direto pelo site do INPE ou pelo site do USGS?
Abraços
"Por um mundo onde sejamos socialmente iguais, humanamente diferentes e totalmente livres"
Rosa Luxemburgo
Jéssica
P. dos Santos
Geógrafa/ Prof. de Geografia
Especialista em Meio Ambiente e Geoprocessamento
Em 30 de setembro de 2016 10:28, Sabrina Mieko <sabrin...@gmail.com> escreveu:
Sabrina MiekoAbraços,Oi Marcelo!Eu pensei nisso quando estava fazendo os testes, por isso que acabei voltando as imagens brutas pra refazer os cálculos e mesmo assim as temperaturas finais acabavam em números estranhos (-120°C no interior de SP em dezembro... Estranho!). Mesmo agora que acertei as conversões, pra não dar erro, acabei deixando pra fazer os recortes no final com as imagens já com as infos de temperatura.
Sabrina Mieko Viana
PhD in Forest Resources - ESALQ-USP
Brazil
+55 16 99750-8849 (Vivo) - Telegram
+55 16 99120-5720 (Claro) - Whatsapp
skype: sabrina.mieko
Você recebeu essa mensagem porque está inscrito no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar inscrição nesse grupo e parar de receber e-mails dele, envie um e-mail para qgisbrasil+unsubscribe@googlegroups.com.
Sabrina Mieko
Nov 20, 2017, 4:22:57 PM11/20/17
to qgisb...@googlegroups.com
Boa noite Jéssica!
Baixei direto pela USGS mesmo, um pacote enorme com todas as bandas.
Abçs,
Sabrina Mieko - via smartphone
Reply all
Reply to author
Forward
0 new messages