[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