[Gvsig_usuarios] Transformação radiométrica de imagens Landsat 8 (16 bits para 8 bits)

Óscar Martínez omartinez en gvsig.com
Sab Nov 5 11:47:26 CET 2016


Por mi parte voy a seguir investigando un poco sobre este tema.. Voy a 
seguir mirando las herramientas que se suelen usar de gdal para 
transformar imagenes, voy a leer un poco sobre las correcciones que se 
suelen hacer a estas imágenes, y poco a poco a ver si vamos sacando algo 
realmente útil para la comunidad.

Solo añadir que con mi script anterior hago un escalado de valores, 
desde el raster de 16 bits cogiendo sus valores de pixel min y max, los 
escalo a la imagen de 8 bits para que entren entre 0 y 255, ocupando 
todo ese rango de valores. Si no me equivoco, estaría haciendo un realce 
de la imagen. No se si esto es correcto o afecta el resultado, a la hora 
de usarlo en el cálculo de indices, pero era como aparecía en el video 
de referencia del inicio de este post usando -scale in gdal_translate:

http://www.gdal.org/gdal_translate.html

Gracias Gilberto por la aportación de los pasos a seguir, lo utilizaré 
de referencia para comprobar valores.



El 03/11/16 a las 21:12, Gilberto Cugler escribió:
> Olá  amigo Eliazer
> Fiz uma prova para ler, calcular NDVI e classificação não 
> supervisionada com a imagem (LC82200772015161LGN00 ) do Landsat 8 do 
> município de Cajati/SP fuso 22-UTM coordenadas E=791.000 e N=7.263.000.
> Usei a versão 2.3.1., e este foi o roteiro que segui:
> 1-baixei as bandas 2,3,4,5,6 pelo endereço https://lv.eosda.com
> 2-adicionei as imagens no projeto com EPSG31982. O gvSIG reprojetou 
> automaticamente da projeção32622(WGS 84/UTM 22N) para 31982.
> 3-Reprojetei a camada para EPSG 31982, usando a ferramenta de 
> reprojetar camada.
> 4-Criei um Shape para recorte de 30kmx30km sobre a sede do 
> município.Recortei com a ferramenta de recorte do sextante.
> 5-Calculei o NDVI, com a ferramenta do sextante.
> 6-Fiz uma composição falsa cor 6-5-4
> 7-Fiz uma classificação não supervisionada com 6 classes
> Tudo sem precisar transformar para 8 Bit.
> Seguindo este roteiro, tudo foi perfeito.
> Vou fazer uma reclassificação e vetorização para comparar melhor os 
> resultados e depois comento o resultado.
>
> A versão 2.3.1 trouxe muitas novidades e esta é uma delas(a biblioteca 
> GDAL)
> Tem um inconveniente, que é preciso testar melhor... me parece que em 
> alguns casos a camada não pode estar reprojetada na vista.
> Você que é um especialista em tutorias, poderia testar melhor este 
> roteiro e ver se não cometi erros e nos blindar com um.
> Também testei o script que o Óscar fiz para transformação de 16 para 8 
> bits, muito bom pois mostra as possibilidades de trabalhar mais 
> facilmente com a nova versão do gvSIG
>
> forte abraço
>
>
> 2016-11-01 21:53 GMT-02:00 Óscar Martínez <omartinez en gvsig.com 
> <mailto:omartinez en gvsig.com>>:
>
>     Se me olvidaba comentar. Haré alguna prueba extra para asegurarme
>     y comprobar que el escalado de valores es acorde con el resultado
>     de otras herramientas.
>
>
>     El 02/11/16 a las 00:50, Óscar Martínez escribió:
>
>         Genial!
>
>         Voy a probar a hacer algo más completo, darle una interfaz
>         visual, y una barra de proceso. Ya estuve haciendo pruebas
>         durante el día de hoy. Volveré a escribir cuando tenga algo
>         con un mejor acabado.
>
>         Saludos
>
>
>         El 01/11/16 a las 19:36, Eliazer Kosciuk escribió:
>
>
>             Testei o script e funcionou perfeitamente. A única coisa
>             que eu acrescentaria seria uma barra de progresso da
>             operação, pois ela pode ser bem demorada e alguém pode
>             achar que não está funcionando. No mais, foi possível
>             gerar corretamente o procedimento de NDVI com as imagens
>             geradas.
>
>
>         _______________________________________________
>         gvSIG_usuarios mailing list
>         gvSIG_usuarios en listserv.gva.es
>         <mailto:gvSIG_usuarios en listserv.gva.es>
>
>         Para ver histórico de mensajes, editar sus preferencias de
>         usuario o darse de baja en esta lista, acuda a la siguiente
>         dirección:
>
>         https://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_usuarios
>         <https://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_usuarios>
>
>
>     _______________________________________________
>     gvSIG_usuarios mailing list
>     gvSIG_usuarios en listserv.gva.es <mailto:gvSIG_usuarios en listserv.gva.es>
>
>     Para ver histórico de mensajes, editar sus preferencias de usuario
>     o darse de baja en esta lista, acuda a la siguiente dirección:
>
>     https://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_usuarios
>     <https://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_usuarios>
>
>
>
>
> -- 
>
>
> Gilberto Cugler .´.
>
>
> _______________________________________________
> gvSIG_usuarios mailing list
> gvSIG_usuarios en listserv.gva.es
>
> Para ver histórico de mensajes, editar sus preferencias de usuario o darse de baja en esta lista, acuda a la siguiente dirección:
>
> https://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_usuarios

------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <http://listserv.gva.es/pipermail/gvsig_usuarios/attachments/20161105/e2e6ddd8/attachment.html>


Más información sobre la lista de distribución gvSIG_usuarios