logo       

Re: problems loading spproj library for transform(): msg#00057

lang.r.geo

Subject: Re: problems loading spproj library for transform()

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"


<Prev in Thread] Current Thread [Next in Thread>
Google Custom Search

News | FAQ | advertise