[R-sig-Geo] spTransform error

[Roger Bivand]

On Wed, 22 May 2013, Agustin Lobo wrote:

> Just got the very same problem on osx:
> CRS("+init=epsg:3857")
> Error in CRS("+init=epsg:3857") : no system list, errno: 2

I have given the diagnosis in detail - the OSX CRAN rgdal binary has been 
built without the required configure argument --with-data-copy=yes. I am 
awaiting a response from the OSX maintainer (of all of R, not just rgdal). 
When he gets back to me, we'll know more about when the binary package can 
be rebuilt. While waiting, consider other OS, or consider installing 
William Kyngesburye's framework-based rgdal 0.8-8 for R 3.0 with GDAL 1.10 
(64-bit only).

Roger

>
> I've reinstalled rgdal (using RStudio) but the error persists.
> I also have the autodetected in the path:
>> require(rgdal)
> Loading required package: rgdal
> Loading required package: sp
> rgdal: version: 0.8-9, (SVN revision 470)
> Geospatial Data Abstraction Library extensions to R successfully loaded
> Loaded GDAL runtime: GDAL 1.9.2, released 2012/10/08
> Path to GDAL shared files:
> Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480]
> Path to PROJ.4 shared files: (autodetected)
>
> Is there anything I can do or should I wait for a newer version of
> rgdal to be on cran?
> Not a big problem if have to wait, can use a linux machine.
>
> Thanks
>
> Agus
>
>> sessionInfo()
> R version 3.0.0 (2013-04-03)
> Platform: x86_64-apple-darwin10.8.0 (64-bit)
>
> locale:
> [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
>
> attached base packages:
> [1] stats     graphics  grDevices utils     datasets  methods   base
>
> other attached packages:
> [1] raster_2.1-16 rgdal_0.8-9   sp_1.0-9
>
> loaded via a namespace (and not attached):
> [1] grid_3.0.0      lattice_0.20-15 tools_3.0.0
>
> On Wed, May 15, 2013 at 10:23 PM, Roger Bivand <Roger.Bivand at nhh.no> wrote:
>> On Wed, 15 May 2013, Ross Ahmed wrote:
>>
>>> I¹ve recently reinstalled QGIS - might this have anything to do with it?
>>
>>
>> Probably not.
>>
>>
>>>
>>>
>>> # spTransform() works fine in the other example I posted to the group
>>> yesterday:
>>>
>>> library(sp); library(rgeos); library(rgdal)
>>> polygonCoords <- matrix(c(-1.9450, -1.9075, -1.9075, -1.9450, -1.9450,
>>> 55.72476, 55.72476, 55.70267, 55.70267, 55.72476), ncol=2)
>>> p = Polygon(polygonCoords)
>>> myPolygon = Polygons(list(p),1)
>>> myPolygonSpatial = SpatialPolygons(list(myPolygon))
>>> proj4string(myPolygonSpatial) <- CRS("+proj=longlat +datum=WGS84
>>> +ellps=WGS84 +towgs84=0,0,0")
>>> gArea(spTransform(myPolygonSpatial, CRS("+proj=laea +lon_0=-1.92
>>> +lat_0=55.7 +datum=WGS84")))/1e6
>>> #[1] 5.796854
>>>
>>>
>>>
>>>
>>> # But spTransform() fails in all other instances, for example:
>>>
>>> Library(maptools)
>>> data(wrld_simpl)
>>> spTransform(wrld_simpl[1,], CRS("+init=epsg:4326"))
>>> # Error in spTransform(wrld_simpl[1, ], CRS("+init=epsg:4326")) :
>>> # error in evaluating the argument 'CRSobj' in selecting a method for
>>> function 'spTransform': Error in CRS("+init=epsg:4326") : no system list,
>>> errno: 2
>>>
>>>
>>
>> The failure is caused by no proj (or gdal) data files being available.
>>
>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> # rgdal() start up message
>>>
>>> library(rgdal)
>>> Loading required package: sp
>>> rgdal: version: 0.8-9, (SVN revision 470)
>>> Geospatial Data Abstraction Library extensions to R successfully loaded
>>> Loaded GDAL runtime: GDAL 1.9.2, released 2012/10/08
>>> Path to GDAL shared files:
>>> Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480]
>>> Path to PROJ.4 shared files: (autodetected)
>>>
>>
>> These are the clues:
>>
>>
>> Path to GDAL shared files:
>> Path to PROJ.4 shared files: (autodetected)
>>
>> which should not look like this if you are on:
>>
>> Platform: x86_64-apple-darwin10.8.0 (64-bit)
>>
>> and installed rgdal as an OSX binary from CRAN. The current OSX binary is
>> packaged without the proj/ and gdal/ subdirectories. I'll try to get them
>> reinstated.
>>
>> Roger
>>
>>
>>>
>>>
>>>
>>>
>>>
>>> # sessionInfo() output
>>>
>>> sessionInfo()
>>>
>>> R version 3.0.0 (2013-04-03)
>>> Platform: x86_64-apple-darwin10.8.0 (64-bit)
>>>
>>> locale:
>>> [1] en_GB.UTF-8/en_GB.UTF-8/en_GB.UTF-8/C/en_GB.UTF-8/en_GB.UTF-8
>>>
>>> attached base packages:
>>> [1] tcltk stats4 splines grid stats graphics grDevices utils datasets
>>> methods
>>> [11] base
>>>
>>> other attached packages:
>>> [1] Hmisc_3.10-1.1 adehabitat_1.8.12 shapefiles_0.7 tkrplot_0.0-23
>>> [5] ggmap_2.3 formatR_0.7 lme4_0.999999-2 Matrix_1.0-12
>>> [9] adehabitatHR_0.4.7 adehabitatLT_0.3.11 CircStats_0.2-4 MASS_7.3-26
>>> [13] adehabitatMA_0.3.6 ade4_1.5-2 spatialsegregation_2.38 dismo_0.8-11
>>> [17] raster_2.1-25 geosphere_1.2-28 knitr_1.2 lubridate_1.3.0
>>> [21] stringr_0.6.2 boot_1.3-9 coin_1.0-22 modeltools_0.2-19
>>> [25] mvtnorm_0.9-9994 survival_2.37-4 gmodels_2.15.4 ggplot2_0.9.3.1
>>> [29] spatstat_1.31-2 deldir_0.0-22 mgcv_1.7-22 bibtex_0.3-5
>>> [33] plyr_1.8 reshape2_1.2.2 ProjectTemplate_0.4-2 testthat_0.7.1
>>> [37] maptools_0.8-23 lattice_0.20-15 foreign_0.8-53 rgdal_0.8-9
>>> [41] rgeos_0.2-17 sp_1.0-9
>>>
>>> loaded via a namespace (and not attached):
>>> [1] cluster_1.14.4 colorspace_1.2-2 dichromat_2.0-0 digest_0.6.3
>>> evaluate_0.4.3
>>> [6] gdata_2.12.0.2 gtable_0.1.2 gtools_2.7.1 labeling_0.1 mapproj_1.2-1
>>> [11] maps_2.3-2 munsell_0.4 nlme_3.1-109 png_0.1-4 proto_0.3-10
>>> [16] RColorBrewer_1.0-5 RgoogleMaps_1.2.0.3 rjson_0.2.12 scales_0.2.3
>>> tools_3.0.0
>>>
>>>
>>>
>>>
>>>
>>>
>>> On 14/05/2013 22:33, "Roger Bivand" <Roger.Bivand at nhh.no> wrote:
>>>
>>>> On Tue, 14 May 2013, Ross Ahmed wrote:
>>>>
>>>>> I get the following error when running SpTransform, and error is
>>>>> occurring
>>>>> not matter what object I try to transform:
>>>>>
>>>>> Error in spTransform(columbus, CRS("+init=epsg:27700 +datum=WGS84")) :
>>>>>  error in evaluating the argument 'CRSobj' in selecting a method for
>>>>> function 'spTransform': Error in CRS("+init=epsg:27700 +datum=WGS84") :
>>>>> no
>>>>> system list, errno: 2
>>>>>
>>>>
>>>> You must include all of the important information: all of sessionInfo()
>>>> output, all of the startup messages given by rgdal on loading, and a
>>>> reasonable example (this one isn't, the actual projection is unknown). I
>>>> think you have installed rgdal from source on a system without the
>>>> necessary extra files, but without knowledge of how you installed rgdal,
>>>> it is not possible to tell.
>>>>
>>>> Roger
>>>>
>>>>>
>>>>>
>>>>> I?ve tried to replicate my code using the columbus
>>>>> SpatialPolygonsDataFrame
>>>>> from maptools package:
>>>>>
>>>>> setwd(system.file("shapes", package="maptools"))
>>>>> columbus <- readShapeSpatial(?columbus.shp?)
>>>>> proj4string(columbus) <- CRS("+proj=longlat +datum=WGS84 +ellps=WGS84
>>>>> +towgs84=0,0,0")
>>>>> spTransform(columbus, CRS("+init=epsg:4326?))
>>>>>
>>>>> Error in spTransform(columbus, CRS("+init=epsg:27700 +datum=WGS84")) :
>>>>>  error in evaluating the argument 'CRSobj' in selecting a method for
>>>>> function 'spTransform': Error in CRS("+init=epsg:27700 +datum=WGS84") :
>>>>> no
>>>>> system list, errno: 2
>>>>>
>>>>>
>>>>>
>>>>> I?m using R 3.0.0
>>>>>
>>>>> Ross
>>>>>
>>>>>
>>>>>
>>>>>         [[alternative HTML version deleted]]
>>>>>
>>>>>
>>>>
>>>> --
>>>> Roger Bivand
>>>> Department of Economics, NHH Norwegian School of Economics,
>>>> Helleveien 30, N-5045 Bergen, Norway.
>>>> voice: +47 55 95 93 55; fax +47 55 95 95 43
>>>> e-mail: Roger.Bivand at nhh.no
>>>>
>>>
>>>
>>>
>>
>> --
>> Roger Bivand
>> Department of Economics, NHH Norwegian School of Economics,
>> Helleveien 30, N-5045 Bergen, Norway.
>> voice: +47 55 95 93 55; fax +47 55 95 95 43
>> e-mail: Roger.Bivand at nhh.no
>>
>> _______________________________________________
>> R-sig-Geo mailing list
>> R-sig-Geo at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
>

-- 
Roger Bivand
Department of Economics, NHH Norwegian School of Economics,
Helleveien 30, N-5045 Bergen, Norway.
voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: Roger.Bivand at nhh.no