[R] Re-projecting geotiff

[MacQueen, Don]

r-sig-geo would be a better place to ask this question.

-Don

-- 
Don MacQueen

Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062





On 10/18/12 6:52 AM, "Filoche" <pmassicotte at hotmail.com> wrote:

>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 at 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.