quinta-feira, fevereiro 2

Manipulando imagens GeoCover com gdal – Parte III

Corte das imagens compatível ao padrão Ao Milionésimo

Antes da conversão para coordenadas geográficas, os valores das coordenadas do retângulo envolvente de cada folha 1:1.000.000 foram obtidos. Cada folha é composta por duas imagens Geocover, ou seja, a folha SD24 é composta por duas partes, uma parte da imagem S-24-10-2000-UTM.tif e outra parte da imagem S-24-15-2000-UTM.tif. Utilizou-se novamente o programa gdalwarp para converter os pedaços referentes à composição da folha SD24, do sistema de coordenadas UTM para o sistema geográfico e para cortar o limite da folha. As sintaxes abaixo correspondem a dois passos:

Passo 1 - Conversão para o sistema geográfico.
gdalwarp -t_srs "+proj=latlong +datum=WGS84" input.tif output.tif
Passo 2 - Corte da imagem que está inserida no retângulo envolvente.
gdalwarp -te Xmin Ymin Xmax Ymax input.tif output.tif
O passo 1 permitiu que as imagens S-24-10-2000-UTM.tif e S-24-15-2000-UTM.tif fossem convertidas para as coordenadas geográficas. As novas imagens foram nomeadas de acordo com a seguinte nomenclatura: S-24-10-2000-GEO.tif.

O passo 2 serviu para cortar as partes de cada imagem que compõem a folha. O parâmetro -te corresponde aos limites de cada imagem que estarão dentro do limite da folha SD24 (retângulo envolvente: Xmin = -42º, Ymin = -16º, Xmax = -36º e Ymax = -12º).

As funções acima foram aplicadas e geraram duas novas imagens. Estas devem ser somadas para compor uma única cena que represente a folha. Para que estas duas imagens seja moisaicadas utilizou-se o programa gdal_merge.py para executar esta tarefa.
gdal_merge.py -o output.tif -of GTiff -v input1.tif input2.tif
O parâmetro -o defina a imagem de saída que será gerada, -of define que o formato da imagem de saída será GTiff (geotiff), -v significa verbose e mostra o andamento do processamento de mosaicagem. As imagens que darão origem ao mosaico vêm no final.

No caso da folha SD24 os dois inputs, S-24-10-2000-GEO-S.tif e S-24-10-2000-GEO-I.tif, foram convertidos no arquivo SD24.tif. A imagem S-24-10-2000-GEO-S.tif corresponde ao pedaço superior da folha, enquanto que a imagem S-24-15-2000-GEO-I.tif faz parte do pedaço inferior.

Marcadores: , ,