[GRASS-dev] [GRASS GIS] #636: r.thin creates triangle at acute intersection of two lines

#636: r.thin creates triangle at acute intersection of two lines
-----------------------+----------------------------------------------------
Reporter: needelsd | Owner: grass-dev@lists.osgeo.org
     Type: defect | Status: new
Priority: normal | Milestone: 6.4.0
Component: Raster | Version: 6.4.0 RCs
Keywords: r.thin | Platform: MacOSX
      Cpu: OSX/Intel |
-----------------------+----------------------------------------------------
I have rasters that were generated using v.rast followed by r.buffer.
Using r.thin on such a raster occasionally creates a triangle at an
intersection instead of the expected two intersecting lines. All of the
examples I have seen involve relatively acute angles (less than 25
degrees??).

This result is easily seen by using r.to.vect to create a vector, which
then needs v.clean to remove dangles. I have attached a screenshot that
shows the original (fat) raster in magenta, the vector with dangles in
white, and the cleaned vector as (a thinner) green. There are numerous
acute angles that generate the vectors I would expect, with one exception
in the center.

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/636&gt;
GRASS GIS <http://grass.osgeo.org>

#636: r.thin creates triangle at acute intersection of two lines
-----------------------+----------------------------------------------------
  Reporter: needelsd | Owner: grass-dev@lists.osgeo.org
      Type: defect | Status: new
  Priority: normal | Milestone: 6.4.0
Component: Raster | Version: 6.4.0 RCs
Resolution: | Keywords: r.thin, r.to.vect
  Platform: MacOSX | Cpu: OSX/Intel
-----------------------+----------------------------------------------------
Changes (by hamish):

  * keywords: r.thin => r.thin, r.to.vect

Comment:

I think this is really a r.to.vect issue and not a r.thin one.

I am not sure if there is a correct solution though. i.e. this could be a
(desired in some cases) feature, not a bug.

to remove the undangles maybe v.to.db upload line length then v.extract
where line length > some threshold?

Hamish

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/636#comment:1&gt;
GRASS GIS <http://grass.osgeo.org>