get X and Y form point shapefile

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20100210/86cef62d/attachment.pl>
xy <- coordinates(sample.shp)

will give you a 2-column matrix "xy" of the X and Y coordinates.

For more complete reporting for this list, please provide the packages
you are using and the version. The easiest way is to use sessionInfo()

See here for more information, and the posting guide:
http://www.r-project.org/mail.html

Cheers, Mike.

On Thu, Feb 11, 2010 at 8:54 AM, gianni lavaredo
dear Rearchers

I know this is easy but after 5 hours on R, I am honest and ?I don't
rememeber to get X, Y coordinate from point shapefile

sample.shp <- readShapeSpatial("sample.shp")
proj4string(sample.shp) <- CRS("+init=epsg:26911")
str(sample.shp)
?str(sample.shp)

Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots
?..@ data ? ? ? :'data.frame': 50372 obs. of ?1 variable:
?.. ..$ Z: num [1:50372] 1498 1498 1499 1499 1503 ...
?.. ..- attr(*, "data_types")= chr "N"
?..@ coords.nrs : num(0)
?..@ coords ? ? : num [1:50372, 1:2] 268753 268754 268763 268761 268774 ...
?.. ..- attr(*, "dimnames")=List of 2
?.. .. ..$ : chr [1:50372] "0" "1" "2" "3" ...
?.. .. ..$ : chr [1:2] "coords.x1" "coords.x2"
?..@ bbox ? ? ? : num [1:2, 1:2] 268501 4147540 268881 4147786
?.. ..- attr(*, "dimnames")=List of 2
?.. .. ..$ : chr [1:2] "coords.x1" "coords.x2"
?.. .. ..$ : chr [1:2] "min" "max"
?..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots
?.. .. ..@ projargs: chr "+init=epsg:26911"

thanks
gianni

? ? ? ?[[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo