Dear List,
For finding the centroid of a catchment, I tried with the function:
get.Pcent(gallego_catchment)
but I got the following error:
Error en get.Pcent(gallego_catchment) : not a Map object
However, I checked the class:
> class(gallego_catchment)
[1] "SpatialPolygonsDataFrame"
attr(,"package")
[1] "sp"
and I read the catchment boundary with:
p4s <- CRS("+proj=utm +zone=30 +ellps=intl +units=m +no_defs")
# reading the boundary of the gallego catchment
library(maptools)
gallego_catchment <- readShapePoly("only_gallego_catchment.shp",
proj4string=p4s)
Could somebody tell me if I'm doing something wrong?
Thanks in advance
Mauricio
--
hzambran
Linux user #454569 -- Ubuntu user #17469
_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
