Coverter imagens 0º 360º a -180 180º com GDAL

Um dia de estes estava a procurar dados para um trabalho que temos entre-mãos e deparei-me com umas imagens, as que necessitava, com os dados de 0º a 360º. Assim à primeira parece que está tudo bem e com lógica, o problema está em que por aqui gostamos de pôr a Europa no centro do mapa com os dados de -180º a 180º. Então seria necessário cortar a imagem pelo meio e a segunda metade mover-la para o lado esquerdo da primeira metade e voltar a juntar as duas partes. Parecia complicado, mas depois de googlear um pouco dei com a solução e deixo-vos aqui por se alguma vez vos passa o mesmo.

Screenshot1b

Supondo que GDAL está instalado vamos usar gdalinfo para conhecer o numero de colunas que a imagem tem.

Agora já sabemos que a imagem tem 144 colunas por 73 filas. Com gdal_translate cortamos e georeferênciamos a parte direita (este) e repetimos o processo para a parte esquerda (oeste) da imagem. Com o parâmetro -srcwin indicamos por onde cortamos imagem (coluna 72=144/2) e com o parâmetro -a_ullr a posição que queremos. No final juntamos as duas imagens com gdal_merge.py.

Este é o resultado.

Screenshot2b