Why do I have to use proj4string() instead of projection() in this example?
> class(Montseny20090409sh)
[1] "SpatialPointsDataFrame"
attr(,"package")
[1] "sp"
> projection(Montseny20090409sh)
[1] "+proj=utm +zone=31 +ellps=intl +units=m +no_defs"
> a<- SpatialPoints(coordinates(Montseny20090409sh))
> class(a)
[1] "SpatialPoints"
attr(,"package")
[1] "sp"
> projection(a) <- projection(Montseny20090409sh)
Error in checkSlotAssignment(object, name, value) :
"crs" is not a slot in class "SpatialPoints"
Calls: projection<- -> @<- -> slot<- -> checkSlotAssignment
> proj4string(a) <- projection(Montseny20090409sh)
> geogWGS84 <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
> aGW <- spTransform(a, geogWGS84)
> Montseny20090409sh2 <- Montseny20090409sh
> montseny20090409...@data <- cbind(montseny20090409...@data,coordinates(aGW))
Thanks for any clarification,
Agus
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
