[R-sig-Geo] problems loading spproj library for transform()

[Michael Sumner]

Thanks Roger, it certainly does work using spTransform.

I had not noticed this 8)  Below is my revised example.

Cheers, Mike.

Note:  "transform" has been renamed to "spTransform" in rgdal 0.4-4

## using data from ?project example

library(rgdal)

## project works OK
 data(state)
 res <- project(cbind(state.center$x, state.center$y), "+proj=lcc 
+lat_1=48 +lat_2=33 +lon_0=-100")
 res1 <- project(res, "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100", 
inv=TRUE)

## and "transform" has been renamed to "spTransform"

res <- SpatialPoints(res, CRS("+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100"))

spTransform(res, CRS("+proj=longlat"))

sessionInfo()
#Version 2.3.0 (2006-04-24)
#i386-pc-mingw32

#attached base packages:
#[1] "methods"   "stats"     "graphics"  "grDevices" "utils"     "datasets"
#[7] "base"   
#other attached packages:
#   rgdal   pixmap    abind       sp
# "0.4-4"  "0.4-4"  "1.1-0" "0.8-14"