spTransform error

require(rgdal)

sessionInfo()

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