GDAL 1.11 Conversão Radiométrica

Jorge Santos2014

GDAL 1.10: Conversão Radiométrica

Conteúdo

SumárioCapítulo 1.................................................................................................................................. 3

Sobre a Conversão Radiométrica.......................................................................................31.1 Transformação Radiométrica..................................................................................31.2 Sobre o GDAL............................................................................................................. 31.3 Procedimento para Execução do GDAL.................................................................31.4 Primeiros Passos no GDAL.......................................................................................3

1.4.1 Comando para Mudança de Unidade de Disco...........................................41.4.2 Comando cd..................................................................................................... 41.4.3 Comando cls..................................................................................................... 51.4.4 Modo de Edição Rápida..................................................................................5

1.5 O programa gdal_translate (16 Bit para 8 Bit).....................................................61.5.1 Conversão Radiométrica para um único arquivo.........................................61.5.2 Conversão Radiométrica para Múltiplos Arquivos Raster..........................7

1.6 [VIDEO] GDAL 1.10: Conversão Radiométrica para Múltiplos Arquivos.............8Capítulo 2.................................................................................................................................. 9

Contatos................................................................................................................................ 9Autoria deste artigo......................................................................................................... 9Processamento Digital – Geotecnologias e Software Livre.......................................9Competências................................................................................................................... 9

2

Capítulo 1

Sobre a Conversão Radiométrica

1.1 Transformação Radiométrica

É o processo utilizado para modificar o tipo de dado nos arquivos raster. É frequentemente utilizada para conversão de formatos do tipo Inteiro não Sinalizado de 16 Bit para 8 Bit.

A Transformação Radiométrica é o processo que deve ser executado antes da Equalização de Imagens ou Correspondência de Histograma, que são processos para modificação do histograma do raster através do balanço de cores.

A Transformação Radiométrica nunca deve ser executada antes da Classificação de Imagens.

1.2 Sobre o GDAL

GDAL <http://www.gdal.org/> é uma biblioteca geoespacial para tradução de formatos de dados raster. Apesar de não possuir interface gráfica de usuário, o programa possui muitas rotinas interessantes. É preciso possuir pequenas noções de MS-DOS para utilizar o programa em ambiente Windows. GDAL é um software livre.

1.3 Procedimento para Execução do GDAL

Para os usuários que possuem o aplicativo SIG QGIS <http://www.qgis.org>, não é necessário instalar o aplicativo GDAL. Após o download e instalação do QGIS, encontre o GDAL no Menu Iniciar do Windows através do aplicativo OSGeo4W:

Figura 01. Localização do GDAL no grupo de programas do QGIS 2.4 Chugiak.

1.4 Primeiros Passos no GDAL

Os comandos GDAL devem ser digitados diretamente no Prompt de Comando, porém, temos algumas dicas que podem auxiliar iniciantes no programa.

3

GDAL 1.10: Conversão Radiométrica

1.4.1 Comando para Mudança de Unidade de Disco

Em ambiente Windows, a principal unidade de disco é representada pela letra C. Outras letras podem representar outros tipos de mídia como pendrives, drives de CD/DVD-ROM, discos externos, etc. Assumindo que a Unidade F representa outra unidade de disco do computador que contém dados, digite o comando abaixo na janela do GDAL:

F:

Pressione ENTER. A nova unidade de disco será montada e estará pronta para receber novos comandos.

Figura 02. Comando para mudança de unidade de disco no MS-DOS.

Para retornar ao drive padrão, basta digitar C: no Prompt da janela.

1.4.2 Comando cd

Utilizado para acessar um diretório ou subdiretórios de arquivos. Por exemplo, temos uma imagem armazenada na pasta IMAGEM. Esta pasta está na unidade C: e queremos ter acesso a ela para executar comandos.

Digite o comando abaixo no Prompt de Comando:

cd C:\imagem

Figura 03. Comando cd. utilizado para acessar um subdiretório.

O diretório será modificado. Digite cd\ para regressar até o diretório raiz.

4

GDAL 1.10: Conversão Radiométrica

1.4.3 Comando cls

Comando utilizado para limpar a tela. Digite o comando cls para limpar a tela:

Figura 04. Comando cls. Utilizado para limpar a tela.

1.4.4 Modo de Edição Rápida

Permite acessar rapidamente qualquer comando copiado para a janela do Prompt.

Na janela principal, clique no botão de controle e acesse as Propriedades:

Figura 05. Propriedades do Prompt de Comando.

Marque a opção Modo de Edição Rápida e clique no botão OK:

Figura 06. Opção para acessar o Modo de Edição Rápida.

5

GDAL 1.10: Conversão Radiométrica

1.5 O programa gdal_translate (16 Bit para 8 Bit)

Após tomar nota dos principais comandos do MS-DOS, chegou o momento de executar o programa gdal_translate para modificar o tipo de dado dos arquivos raster.

1.5.1 Conversão Radiométrica para um único arquivo

gdal_translate -of GTiff -ot Byte -scale -co TFW=YES input.tif output.tif

gdal_translate: Programa para tradução entre diferentes formatos de arquivo raster.-of GTiff: Parâmetro que modifica o formato do arquivo raster de saída para GeoTIFF.-ot Byte: Parâmetro que define o Tipo de Dado para Inteiro Não Sinalizado de 8 Bit.-scale: Parâmetro que transforma a escala do dado, modificando valores.-co TFW=YES: Opções de Criação. Gera um arquivo externo TFW para o raster de saída.-input.tif/output.tif : Nome para as imagens de entrada e saída.

Acesse a pasta onde está a imagem Landsat-8. Primeiro, realizamos a Composição Colorida RGB 6/5/4 e agora vamos transformar a imagem para 8 Bit.

Figura 07. Local do computador onde encontra-se a imagem Landsat-8.

Selecione e copie o diretório das imagens:

Figura 08. Cópia do diretório onde encontra-se a imagem.

Na janela OSGeo4W, digite cd + [ESPAÇO] e clique com o botão direito sobre o Prompt:

Figura 09. Cópia do diretório onde encontra-se a imagem.

6

GDAL 1.10: Conversão Radiométrica

Pressione [ENTER]. Com a edição rápida, você será capaz de colar rapidamente longos comandos e diretórios do Windows diretamente para a janela do Prompt de Comando. Vamos utilizar esses procedimentos para executar todos os comandos do GDAL.

Copie e cole o comando original no Bloco de Notas. troque os parâmetros input.tif e output.tif pelo nome da sua imagem (na saída, além do nome original da imagem, eu adicionei o sufixo _8bits). No final, para nossa imagem de teste, o comando ficou assim:

gdal_translate -of GTiff -ot Byte -scale -co TFW=YES LC82220742013121LGN01_654.tif LC82220742013121LGN01_654_8bits.tif

Cole este comando na janela do Prompt e pressione [ENTER]. O processo será executado:

Figura 10. Processo para criação de um raster em 8 bits no GDAL.

O novo arquivo raster será adicionado na pasta. Imagem foi convertida para 8 bits.

Figura 11. Processo para criação de um raster em 8 bits no GDAL.

1.5.2 Conversão Radiométrica para Múltiplos Arquivos Raster

Para converter inúmeras imagens na mesma pasta, utilize o comando abaixo:

for %i in (*.tif ) do gdal_translate -of GTiff -ot Byte -scale -co TFW=YES %i %~ni_.tif

for %i in (*.tif ) do: Repetição ou Looping (executa o processo em todos os arquivos da pasta).

gdal_translate: Programa para tradução entre diferentes formatos de arquivo raster.

-of GTiff : Parâmetro que modifica o formato do arquivo raster de saída para GeoTIFF.

-ot Byte: Parâmetro que define o Tipo de Dado para Inteiro Não Sinalizado de 8 Bit.

-scale: Parâmetro que transforma a escala do dado, modificando valores.

-co TFW=YES: Opções de Criação. Gera um arquivo externo TFW para o raster de saída.

%i: Variável que representa todos os arquivos de entrada.

%~ni_.tif: Variável que representa todos os arquivos de saída.

7

GDAL 1.10: Conversão Radiométrica

1.6 [VIDEO] GDAL 1.10: Conversão Radiométrica para Múltiplos Arquivos

Acompanhe todo o procedimento através deste vídeo de cerca de dois minutos.

• GDAL 1.10: Conversão Radiométrica

8

Capítulo 2

Contatos

Autoria deste artigo

Jorge Santos

Lattes: http://lattes.cnpq.br/1910845468254276

E-mail : jorgepsantos@outlook.com

Skype: jorgepsantos2002

Processamento Digital – Geotecnologias e Software Livre

Endereço: http://www.processamentodigital.com.br/

Twitter: http://twitter.com/jpsantos2002

Facebook: http://www.facebook.com/ProcessamentoDigital

LinkedIn: http://br.linkedin.com/pub/jorge-santos/10/38b/8a4/

Competências

• Consultoria em SIG/PDI;

• Treinamento em Geotecnologias / Processamento Digital de Imagens;

• Gerência de Projetos.

9