Data from intersecting Lines of shp files
On Thu, 11 Jul 2013, Barry Rowlingson wrote:
On Thu, Jul 11, 2013 at 1:27 PM, Tscheikner-Gratl, Franz <Franz.Tscheikner-Gratl at uibk.ac.at> wrote:
Thanks for your reply. The problem is that I have data of a pipe network with missing attributes (construction year). So I would like to see what are the constructions years of the intersecting (neighbouring) pipes. Then I would fill the missing attributes with the oldest construction age. I could make it by hand with ArcGis but it would be lots of work - so I thought I could make it with R.
That's a bit clearer - if you could make a simple example shapefile that shows the problem (you don't have to share all your data) I'm pretty sure its just a case of selecting the pipe features and looping over one of the intersection functions, which should return the The problem is if you have more than one labelled pipe crossing an unlabelled pipe, or if no pipes cross - or in that case do you want the nearest pipe? Maybe nearest is a better solution in any case (since crossing pipes have distance zero)...
Maybe gTouches() rather than gIntersects() as maybe the contiguous pipe
segments to segments with NA attribute values can by listed, while I'm not
sure whether they intersect:
library(rgeos)
l1 = readWKT("LINESTRING(0 0,2 2)")
l2 = readWKT("LINESTRING(1 1,2 0)")
gTouches(l1,l2)
gIntersects(l1,l2)
so they do touch and intersect. I would suggest identifying the NA
attribute Lines objects first, then loop over them to find their
contiguous neighbours. If there are two neighbbours, and the value of
interest is identical, it may be copied across, if different, or for three
or more with different values, manual intervention may be needed, unless
there is some other non-NA variable indicating which bit of the jigsaw
puzzle belongs where.
Roger
Worst case scenario is a double loop over features with missing data over features with data... Barry
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Roger Bivand Department of Economics, NHH Norwegian School of Economics, Helleveien 30, N-5045 Bergen, Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43 e-mail: Roger.Bivand at nhh.no