|
Re: problems loading spproj library for transform(): msg#00057lang.r.geo
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> |
|---|---|---|
| Previous by Date: | Re: problems loading spproj library for transform(), Roger Bivand |
|---|---|
| Next by Date: | Re: problems loading spproj library for transform(), karl . sommer |
| Previous by Thread: | Re: problems loading spproj library for transform(), Roger Bivand |
| Next by Thread: | Re: problems loading spproj library for transform(), Roger Bivand |
| Indexes: | [Date] [Thread] [Top] [All Lists] |
| News | FAQ | advertise |