[R-sig-Geo] Re-projection with project raster results in rotated plot

Swagath Navin Manohar swagathnavin82 at gmail.com
Tue Sep 30 08:51:41 CEST 2014


Hello all,

I have a problem concerning re-projecting a raster from EPSG:3035 to 
sr-org: 29, wrf-lambert-conformal-conic which results in a rotated plot.

Here is the code I used:

raster1=rasterFromXYZ(radon[, c("x", "y", "rn")])

####http://spatialreference.org/ref/epsg/3035/proj4/

crs(raster1)='+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 
+ellps=GRS80 +units=m +no_defs'

####to this new projection
########http://spatialreference.org/ref/sr-org/wrf-lambert-conformal-conic/

newproj <- "+proj=lcc +lat_1=33 +lat_2=45 +lat_0=40 +lon_0=-97 +x_0=0 
+y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"

raster2 <- projectRaster(raster2, crs=newproj, res=5000, method='bilinear')

Here is the output: http://s7.postimg.org/g8adezsdn/Page_1_rasters.jpg

I hope someone can give me a hint about what I am doing wrong here.
Thanks a lot for your time.
Navin



More information about the R-sig-Geo mailing list