[R-br] sobre gdal_translate e gdalwarp

[Antonio Silva]

Olá

Tenho algumas dúvidas sobre o uso gdal no R.

Depois de converter um arquivo HDF, com dados de temperatura de superfície
do mar, para geotiff eu quero reprojetar para WGS84 ou SIRGAS2000 e contar
entre as lats de 22 a 29 S e longs de 40 a 50 W.

o HDF para exemplo está em
https://app.box.com/s/16cf7qv6af6gsz1v66staori2mtneu0r
Estou seguindo basicamente
https://scottishsnow.wordpress.com/2014/08/24/many-rastered-beast/

Bem, converti o HDF com:

gdal_translate("A20080012008031.L3m_MO_SST_4","georef.tif",sd_index=1,a_ullr=c(0,4320,8640,0),
a_srs="+proj=eqc +lat_ts=0 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +a=6371007
+b=6371007 +units=m +no_defs")

Sem a_ullr e a_srs quando utilizava o gdalwarp eu recebia mensage: "ERROR
1: Unable to compute a transformation between pixel/line and georeferenced
coordinates"

Copiei os valores de a_ullr and a_srs values de GDALinfo("georef.tif")

map <- raster("georef.tif")
plot(map)

Meu problema agora é reprojetar e cortar a imagem georef.tif:

gdalwarp("georef.tif", "georef2.tif",
s_srs="+proj=eqc +lat_ts=0 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +a=6371007
+b=6371007 +units=m +no_defs",
t_srs="+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs")

map2 <- raster("georef2.tif")
plot(map2)

Map2 não está em SIRGAS2000 e a opção de corte  te=c(-50,-29,-40,-22) não
funciona.

Onde será que estou errando? Talvez em  a_ullr e a_srs, não sei.

Agradeço muito qualquer ajuda.

Sds

Antonio Olinto
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150921/83e4cd56/attachment.html>