Message-ID: <alpine.LRH.2.03.1406031119330.2723@reclus.nhh.no>
Date: 2014-06-03T09:22:45Z
From: Roger Bivand
Subject: Cut SpatialLines into two Lines when intersecting with SpatialPolyon
In-Reply-To: <E9C18A23CC90B64F91C6A03B094B1B7E19820D@MAIL01.arthur.nscr.nl>
On Tue, 3 Jun 2014, Wouter Steenbeek wrote:
> [My apologies if this is a double post. I wrote the first post in HTML,
> and I think it may not have gotten through.]
Thanks for respecting our guidelines - it is much harder to infect plain
text, which is why HTML is not welcome.
>
> I have a SpatialLines(DataFrame) (e.g. streets) and
> SpatialPolygons(DataFrame) (e.g. neighborhoods). When a Line intersects
> with a polygon boundary, I want to split this Line into two new Lines.
> So: if a street begins in one neighborhood and ends in another
> neighborhood, I want to remove the original SpatialLine and replace that
> by two SpatialLines (streets): one street from the starting point until
> the intersection with the polygon boundary, and another street from the
> intersection with the polygon boundary to the original end point.
>
> Can someone give me advice on how to do this?
As in most cases, you will increase your chances of a reply if you provide
a toy example, for instance using readWKT() in rgeos to generate the
objects. Probably all you need is in rgeos, but without a toy example, it
is hard to give advice that isn't more like hand-waving.
Roger
>
> Cheers,
> Wouter
>
> _______________________________________________
> 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, Norwegian School of Economics,
Helleveien 30, N-5045 Bergen, Norway.
voice: +47 55 95 93 55; fax +47 55 95 91 00
e-mail: Roger.Bivand at nhh.no