Manipulando imagens GeoCover com gdal – Parte II

Conversão dos arquivos MrSID para Geotiff

O formato MrSID foi desenvolvido pela LizardTech para codificar imagens de alta-resolução, numa fração do seu tamanho, mantendo-se a qualidade original das imagens, ou seja, o MrSID é um formato para a compressão de imagens com pouca perda de qualidade.

MrSID is a unique raster file format developed by LizardTech to encode large, high-resolution images to a fraction of their original size while maintaining the original image quality.
(link: http://www.mcaggis.com/glossary.htm)

Por ser um formato proprietário, utilizou-se para a conversão de MrSID em Geotiff o pacote MrSID Decode da LizardTech (link: http://www.lizardtech.com). Este arquivo é um executável e seguiu-se a seguinte configuração para a conversão:

./mrsiddecode -i input.sid -o geotiff.tif -of tifg -s 1

O parâmetro -of significa que o formato exportado é um arquivo do tipo tiff geográfico e o parâmetro -s define a resolução de saída do arquivo. Neste caso o -s 1 faz uma reamostragem espacial alterando a resolução do arquivo de 14.25 m para 28.5 m.

Mudança do sistema de referência espacial de UTM para geográfica

Após a compilação do pacote gdal alguns programas são instalados, a descrição destes programas está disponível na página sobre o gdal utilities. Para a mudança do sistema de referência espacial utilizou-se o comando gdalwarp.

A imagem que serviu de input.sid é a S-24-10-2000.sid. Após a execução do comando acima, a imagem de saída (geotiff.tif) é S-24-10-200.tif, com resolução espacial de 28.5 m. O sistema de referência espacial do arquivo original sid/tiff é: 1) projeção: UTM; 2) zona: 24; 3) falso N: 0 m; 4) falso E: 500.000 m; e 5) datum: WGS84.

A adoção do falso norte como 0 m não é comum para o hemisfério sul, então, a seguinte conversão foi executada para que o falso norte da nova imagem fosse de 10.000.000 m.

gdalwarp -s_srs "+proj=utm +zone=24 +datum=WGS84" -t_srs "+proj=utm +zone=24 +datum=WGS84 +south" input.tif output.tif

Os parâmetros -s_srs e -t_srs são da biblioteca PROJ.4, eles indicam o sistema de coordenadas da imagem de entrada e o sistema de coordenadas da transformação ou imagem de saída. O parâmetro +south corrige o problema relacionado ao valor do falso norte para hemisférios norte ou sul. Após se aplicar a função acima, a imagem gerada S-24-10-2000-UTM.tif, possui o seguinte sistema de referência espacial: 1) projeção: UTM; 2) zona: 24; 3) falso N: 10.000.000 m; 4) falso E: 500.000 m; e 5) datum: WGS84.

Como resultado desta parte foram geradas as seguintes imagens: S-24-10-2000-UTM.tif e S-24-15-2000-UTM.tif.