[GRASS-dev] [GRASS GIS] #3567: v.select: doubtful result with the default overlap operator selecting areas with lines
GRASS GIS
trac at osgeo.org
Tue May 29 06:27:52 PDT 2018
#3567: v.select: doubtful result with the default overlap operator selecting areas
with lines
--------------------------+------------------------------
Reporter: mlennert | Owner: grass-dev@…
Type: defect | Status: closed
Priority: normal | Milestone: 7.4.2
Component: Vector | Version: svn-trunk
Resolution: fixed | Keywords: v.select overlap
CPU: Unspecified | Platform: Unspecified
--------------------------+------------------------------
Comment (by mmetz):
Replying to [comment:4 mlennert]:
> Replying to [comment:3 mmetz]:
> > Replying to [ticket:3567 mlennert]:
> > > Trying to select areas with lines, I get a result using the overlap
operator that I don't really understand, and which I suspect being a bug.
> > >
> > > Using the following command in the attached location, test_vselect
mapset:
> > >
> > > {{{
> > > v.select ain=test_polys bin=test_line atype=area btype=line
out=result_overlap op=overlap
> > > }}}
> > >
> > > I get a result of only three areas selected (see attached map)
where, IMHO, all 41 should be.
> > >
> > The problem was not the overlap operator of v.select, but instead that
the vector test_line is 3D, and one of the underlying library functions
did not ignore z as it should have done. This bug has been fixed in trunk
and relbr74 with r72746,7.
>
> Thanks a lot ! Sorry, should have noticed the 3D nature of the vector
myself. I guess there must have been a message at import which I didn't
see.
>
> This obviously raises the question of whether v.select should always
ignore 3D or whether a flag (-3) should allow enabling it so that if one
wants to check spatial relation in 3D this is possible.
This is not possible because `Vect_segment_intersection()` does not
support 3D coordinates.
Also note that 3D topology does not exist yet in GRASS, e.g. a purely
vertical line can cause problems, because it is marked as degenerate at
the corresponding nodes.
--
Ticket URL: <https://trac.osgeo.org/grass/ticket/3567#comment:5>
GRASS GIS <https://grass.osgeo.org>
More information about the grass-dev
mailing list