[R-sig-Geo] rgdal error in pj_transform: failed to load datum shift file

[Olivia Burge]

Dear list,

I am having some trouble transforming data (ie reprojecting) after
installing rgdal, etc on a new mac.

The error message is:
> wgs84 <- CRS("+proj=longlat +datum=WGS84 +no_defs ")
> df <- spTransform(locs_sp, CRS = wgs84)
Error in spTransform(xSP, CRSobj, ...) :
  error in pj_transform: failed to load datum shift file.

locs_sp dput output:

new("SpatialPointsDataFrame"
    , data = structure(list(Plot = c(1L, 2L, 3L, 4L, 5L, 6L, 1L, 1L, 4L,
9L,
19L, 20L), Wetland = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 3L,
4L, 4L, 4L, 4L, 4L), .Label = c("", "Awarua", "Bayswater", "Dunearn"
), class = "factor")), .Names = c("Plot", "Wetland"), row.names = c(NA,
12L), class = "data.frame")
    , coords.nrs = numeric(0)
    , coords = structure(c(2159130, 2159243, 2159242, 2159240, 2159504,
2157624,
2126982, 2138135, 2138087, 2137919, 2137472, 2137450, 5397279,
5397304, 5397304, 5392166, 5397215, 5395541, 5439946, 5456416,
5456657, 5456657, 5456659, 5456662), .Dim = c(12L, 2L), .Dimnames = list(
    c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11",
    "12"), c("NZMGE", "NZMGN")))
    , bbox = structure(c(2126982, 5392166, 2159504, 5456662), .Dim = c(2L,
2L), .Dimnames = list(c("NZMGE", "NZMGN"), c("min", "max")))
    , proj4string = new("CRS"
    , projargs = "+proj=nzmg +lat_0=-41 +lon_0=173 +x_0=2510000
+y_0=6023150 +ellps=intl +units=m +nadgrids=nzgd2kgrid0005.gsb"
)
)

I have a windows computer with rgdal & R on and it works fine on there.

When I load rgdal, I get:
> require(rgdal)
Loading required package: rgdal
Loading required package: sp
rgdal: version: 1.2-5, (SVN revision 648)
 Geospatial Data Abstraction Library extensions to R successfully loaded
 Loaded GDAL runtime: GDAL 2.1.2, released 2016/10/24
 Path to GDAL shared files:
 Loaded PROJ.4 runtime: Rel. 4.9.1, 04 March 2015, [PJ_VERSION: 491]
 Path to PROJ.4 shared files: (autodetected)
 Linking to sp version: 1.2-3

I have seen reference to the output of this:
> system.file("proj", package= "rgdal")
elsewhere. The output is simply ""

SessionInfo:
R version 3.3.2 (2016-10-31)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: macOS Sierra 10.12.1

locale:
[1] en_NZ.UTF-8/en_NZ.UTF-8/en_NZ.UTF-8/C/en_NZ.UTF-8/en_NZ.UTF-8

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

other attached packages:
[1] ggmap_2.7          ggplot2_2.2.1.9000 rgdal_1.2-5        sp_1.2-4


loaded via a namespace (and not attached):
 [1] Rcpp_0.12.9       magrittr_1.5      maps_3.1.1        munsell_0.4.3
  colorspace_1.3-2  geosphere_1.5-5
 [7] lattice_0.20-34   R6_2.2.0          rjson_0.2.15      jpeg_0.1-8
 dplyr_0.5.0       stringr_1.1.0
[13] plyr_1.8.4        tools_3.3.2       grid_3.3.2        gtable_0.2.0
 png_0.1-7         DBI_0.5-1
[19] lazyeval_0.2.0    assertthat_0.1    tibble_1.2        mapproj_1.2-4
  reshape2_1.4.2    bitops_1.0-6
[25] stringi_1.1.2     RgoogleMaps_1.4.1 scales_0.4.1      proto_1.0.0

Any help would be greatly appreciated. I will be out the office for a week
on fieldwork in a couple of hours, after which I won't be able to test
anything until I return.

Kind regards,

Olivia

	[[alternative HTML version deleted]]