[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