Dear R users,
I'm currently trying to re-project a geotiff in another coordinate system.
For instance, I have a tif image in UTM 19 zone which I would like to
reproject into UTM 18. I was wondering if it was possible in R.
Furthermore, I looked into 'rgdal' package, but I can't really find out if
I'm doing the right thing. So far, here is what I'm doing:
library(rgdal)
myTiff = 'L5014028_02819840617_B10.TIF'
a = readGDAL(myTiff)
aa = spTransform(a, CRS("+proj=utm +zone=18 +datum=WGS84"))
However, after few seconds, the memory is full and R crashes. I used other
proprietary tools to reproject and it takes only like 10 seconds to
reproject. Therefore, I suspect I'm not using the 'rgdal' package correctly.
Regards,
Phil
--
View this message in context:
http://r.789695.n4.nabble.com/Re-projecting-geotiff-tp4646616.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
