[R-sig-Geo] spTransform problem

[Mathieu Basille]

Hi Stephen,

Did you *actually* try to load rgdal in that R session? It doesn't seem
to appear in the list of attached packages...

HTH,
Mathieu.



Le 2010-07-22 14:55, stephen sefick a écrit :
> I am having problems with spTransform in packages rgdal.  If this is
> not enought information please let me know and I will send whatever
> you need.  It doesn't seem like it is avaliable.
> 
> library(rgdal)
> w <- (new("SpatialPoints"
>     , coords = structure(c(33.31665, 33.31834, 33.3239, 33.37435, 33.36018,
> 33.23301, 33.23423, 33.17796, 33.18121, 33.1949, 33.31868, 33.30408,
> 33.29747, 33.17787, 33.29818, 33.29995, 33.22684, 35.1132, 35.17292,
> 35.17181, 35.06247, 35.11168, 35.13462, 33.29927, 33.29815, 33.29678,
> 33.22417, 35.18823, 35.18498, 35.11643, 35.01896, 32.50986, 32.49736,
> 33.31731, 33.3194, 33.3244, 33.37333, 33.36077, 33.2325, 33.23349,
> 33.1781, 33.18051, 33.19357, 33.31922, 33.30454, 33.29836, 33.17613,
> 33.2981, 33.36369, 33.30074, 33.22576, 35.1132, 35.17424, 35.17208,
> 35.06282, 35.11102, 35.13354, 33.30036, 33.29896, 33.29789, 33.22353,
> 35.18748, 35.1838, 35.11509, 35.01986, 32.51054, 32.49748, -81.59802,
> -81.59618, -81.60159, -81.55044, -81.55657, -81.6235, -81.63911,
> -81.56433, -81.56165, -81.57855, -81.57896, -81.62487, -81.63039,
> -81.58159, -81.62553, -81.58714, -81.63564, -79.33718, -79.1805,
> -79.22391, -79.30669, -79.35294, -79.34676, -81.66629, -81.67698,
> -81.68295, -81.5105, -79.0738, -79.09837, -79.32681, -79.51664,
> -84.64703, -84.70618, -81.59759, -81.59619, -81.60101, -81.54981,
> -81.55763, -81.62424, -81.63809, -81.5657, -81.56335, -81.5788,
> -81.58011, -81.6263, -81.62959, -81.58174, -81.62629, -81.55806,
> -81.5868, -81.6357, -79.33718, -79.18026, -79.23305, -79.30483,
> -79.35168, -79.34534, -81.66618, -81.6779, -81.68343, -81.50894,
> -79.07469, -79.09892, -79.32548, -79.51553, -84.64841, -84.70769
> ), .Dim = c(67L, 2L), .Dimnames = list(NULL, c("lat", "long")))
>     , bbox = structure(c(32.49736, -84.70769, 35.18823, -79.0738), .Dim = c(2L,
> 2L), .Dimnames = list(c("lat", "long"), c("min", "max")))
>     , proj4string = new("CRS"
>     , projargs = "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"
> )
> ))
> 
> 
> str(w)
> Formal class 'SpatialPoints' [package "sp"] with 3 slots
>   ..@ coords     : num [1:67, 1:2] 33.3 33.3 33.3 33.4 33.4 ...
>   .. ..- attr(*, "dimnames")=List of 2
>   .. .. ..$ : NULL
>   .. .. ..$ : chr [1:2] "lat" "long"
>   ..@ bbox       : num [1:2, 1:2] 32.5 -84.7 35.2 -79.1
>   .. ..- attr(*, "dimnames")=List of 2
>   .. .. ..$ : chr [1:2] "lat" "long"
>   .. .. ..$ : chr [1:2] "min" "max"
>   ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots
>   .. .. ..@ projargs: chr "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"
> 
> 
> w.utm <- spTransform(y, CRS=CRS("+proj=utm +zone=17 +ellps=GRS80
> +datum=NAD83 +units=m +no_defs"))
> Error: could not find function "spTransform"
> 
> sessionInfo()
> R version 2.11.1 (2010-05-31)
> x86_64-unknown-linux-gnu
> 
> locale:
>  [1] LC_CTYPE=en_US.utf8       LC_NUMERIC=C
>  [3] LC_TIME=en_US.utf8        LC_COLLATE=en_US.utf8
>  [5] LC_MONETARY=C             LC_MESSAGES=en_US.utf8
>  [7] LC_PAPER=en_US.utf8       LC_NAME=C
>  [9] LC_ADDRESS=C              LC_TELEPHONE=C
> [11] LC_MEASUREMENT=en_US.utf8 LC_IDENTIFICATION=C
> 
> attached base packages:
> [1] graphics  grDevices utils     datasets  stats     grid      methods
> [8] base
> 
> other attached packages:
>  [1] gpclib_1.5-1            RPostgreSQL_0.1-6       DBI_0.2-5
>  [4] StreamMetabolism_0.03-3 maptools_0.7-34         lattice_0.18-8
>  [7] sp_0.9-65               foreign_0.8-40          chron_2.3-35
> [10] zoo_1.6-4               vegan_1.17-3            ggplot2_0.8.8
> [13] proto_0.3-8             reshape_0.8.3           plyr_1.0.3
> 
> loaded via a namespace (and not attached):
> [1] tools_2.11.1
> 
> 

-- 

~$ whoami
Mathieu Basille, Post-Doc

~$ locate
Laboratoire d'Écologie Comportementale et de Conservation de la Faune
+ Centre d'Étude de la Forêt
Département de Biologie
Université Laval, Québec

~$ info
http://ase-research.org/basille

~$ fortune
``If you can't win by reason, go for volume.''
Calvin, by Bill Watterson.