[GRASS-user] How to find intersection of two lines
Moritz Lennert
mlennert at club.worldonline.be
Thu Mar 20 08:56:39 EDT 2008
On 20/03/08 11:17, Jarosław Jasiewicz wrote:
> Hi I have lots of raster, thinned lines and second raster with with one
> line. I must find all points where second line intersect with lines on
> first raster
> it is very simple:
>
> r.mapcalc points=if(raster_two==1&&raster_one==1,1,null())
>
> but there is one exception:
>
> AB
> BA
>
> when lines intersect with 45 grades. There is intersection but without
> common point. Naither buffor nor grow is not a solution because I need
> exactly one point lying exactly on line from raster one
>
> Any Ideas
Either ressample your rasters to coarser resolution to make them
coincide, or convert to vector and use a combination of v.patch +
v.clean (e.g. v.clean input=lines1 output=lines2 err=points tool=break)
Moritz
More information about the grass-user
mailing list