<div dir="ltr"><div><div><div>Olá<br><br></div><div>Tenho algumas dúvidas sobre o uso gdal no R.<br><br></div><div>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 <br></div></div>entre as lats de 22 a 29 S e longs de 40 a 50 W.<br></div><div><br></div>o HDF para exemplo está em  <a href="https://app.box.com/s/16cf7qv6af6gsz1v66staori2mtneu0r" target="_blank">https://app.box.com/s/16cf7qv6af6gsz1v66staori2mtneu0r</a><br><div>Estou seguindo basicamente <a href="https://scottishsnow.wordpress.com/2014/08/24/many-rastered-beast/" target="_blank">https://scottishsnow.wordpress.com/2014/08/24/many-rastered-beast/</a></div><div><br></div><div>Bem, converti o HDF com:<br></div><div><br>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")<br><br></div><div>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"<br></div><div><br>Copiei os valores de a_ullr and a_srs values de GDALinfo("georef.tif")<br><br>map <- raster("georef.tif")<br>plot(map)<br><br></div><div>Meu problema agora é reprojetar e cortar a imagem georef.tif: <br></div><div><br>gdalwarp("georef.tif", "georef2.tif",<br>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",<br>t_srs="+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs")<br><br>map2 <- raster("georef2.tif")<br>plot(map2)<br><br></div><div>Map2 não está em SIRGAS2000 e a opção de corte  te=c(-50,-29,-40,-22) não funciona.<br><br><span lang="en"><span>Onde será que estou errando? </span></span><span lang="en"><span>Talvez em  </span></span><span lang="en"><span>a_ullr e a_srs, não sei.<br><br></span></span></div><div><span lang="en"><span>Agradeço muito qualquer ajuda.<br><br></span></span></div><div><span lang="en"><span>Sds<br><br></span></span></div><span lang="en"><span>Antonio Olinto</span></span></div>