Hi Tom, your example works for me; maybe update rgdal? If this is the
problem, than raster DESCRIPTION needs to update its rgdal dependency.
R version 3.1.3 (2015-03-09)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 14.04.2 LTS
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
[5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
[7] LC_PAPER=en_US.UTF-8 LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] raster_2.3-33 sp_1.1-0
loaded via a namespace (and not attached):
[1] grid_3.1.3 lattice_0.20-30 rgdal_0.9-2
On 03/22/2015 09:46 PM, Tomislav Hengl wrote:
Hi Robert,
I think something happened in the raster package v2.3-33 and your
function projectRaster reports now an error (which probably affects
several other packages).
# create a new (not projected) RasterLayer with cellnumbers as values
r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40)
r <- setValues(r, 1:ncell(r))
projection(r)
[1] "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"
# proj.4 projection description
newproj <- "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100 +ellps=WGS84"
pr1 <- projectRaster(r, crs=newproj)
Error: 'rawTransform' is not an exported object from 'namespace:rgdal'
Thank you!
T. Hengl