Hi !

I am currently working on a wombling (boundary detection) function and I have come across a nag. Hopefully you will be able to give me a hand on this.

Following is a small example to explain what I want to do using the spdep package.

xy<-expand.grid(1:4,1:2)
nbobj<-cell2nb(4,2,"rook")

links<-listw2sn(nb2listw(nbobj))
links.new<-links[-c(4,6,15,17),]

nbobj.new<-sn2listw(links.new)$neighbours

plot(nbobj.new,xy)

"nbobj.new" has two separate groups of points links between each other. Is there a way to automatically counts the number of independent (non-linked) groups of sites ?

Thanks in advance !

Guillaume Blanchet

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

Reply via email to