[R] rgdal functions: spTransform and CRS not working
Ossenbruggen, Paul
Paul.Ossenbruggen at unh.edu
Fri Jan 13 22:21:11 CET 2017
The following are examples from the rgdal library help file.
Neither functions operate properly.
Help is appreciated.
Paul
Example 1
> EPSG <- make_EPSG()
Error in make_EPSG() : Error opening epsg file
> EPSG[grep("Oslo", EPSG$note), 1:2]
Error: object 'EPSG' not found
> EPSG[1925:1927, 3]
Error: object 'EPSG' not found
Example 2
> data(meuse)
> coordinates(meuse) <- c("x", "y")
> proj4string(meuse) <- CRS(paste("+init=epsg:28992",
+ "+towgs84=565.237,50.0087,465.658,-0.406857,0.350733,-1.87035,4.0812"))
NOTE: rgdal::checkCRSArgs: no proj_defs.dat in PROJ.4 shared files
Error in CRS(paste("+init=epsg:28992", "+towgs84=565.237,50.0087,465.658,-0.406857,0.350733,-1.87035,4.0812")) :
no system list, errno: 2
> # see http://trac.osgeo.org/gdal/ticket/1987
> summary(meuse)
Object of class SpatialPointsDataFrame
Coordinates:
min max
x 178605 181390
y 329714 333611
Is projected: NA
proj4string : [NA]
Number of points: 155
Data attributes:
cadmium copper lead zinc
Min. : 0.200 Min. : 14.00 Min. : 37.0 Min. : 113.0
1st Qu.: 0.800 1st Qu.: 23.00 1st Qu.: 72.5 1st Qu.: 198.0
Median : 2.100 Median : 31.00 Median :123.0 Median : 326.0
Mean : 3.246 Mean : 40.32 Mean :153.4 Mean : 469.7
3rd Qu.: 3.850 3rd Qu.: 49.50 3rd Qu.:207.0 3rd Qu.: 674.5
Max. :18.100 Max. :128.00 Max. :654.0 Max. :1839.0
elev dist om ffreq soil lime
Min. : 5.180 Min. :0.00000 Min. : 1.000 1:84 1:97 0:111
1st Qu.: 7.546 1st Qu.:0.07569 1st Qu.: 5.300 2:48 2:46 1: 44
Median : 8.180 Median :0.21184 Median : 6.900 3:23 3:12
Mean : 8.165 Mean :0.24002 Mean : 7.478
3rd Qu.: 8.955 3rd Qu.:0.36407 3rd Qu.: 9.000
Max. :10.520 Max. :0.88039 Max. :17.000
NA's :2
landuse dist.m
W :50 Min. : 10.0
Ah :39 1st Qu.: 80.0
Am :22 Median : 270.0
Fw :10 Mean : 290.3
Ab : 8 3rd Qu.: 450.0
(Other):25 Max. :1000.0
NA's : 1
> meuse.utm <- spTransform(meuse, CRS("+proj=utm +zone=32 +datum=WGS84"))
NOTE: rgdal::checkCRSArgs: no proj_defs.dat in PROJ.4 shared files
Error in spTransform(xSP, CRSobj, ...) :
No transformation possible from NA reference system
> summary(meuse.utm)
Error in summary(meuse.utm) : object 'meuse.utm' not found
>
More information about the R-help
mailing list