[GRASS-dev] [GRASS GIS] #4009: 'where' option in v.to.rast can select wrong feature for raster attribute, in areas where features overlap

[GRASS GIS]

#4009: 'where' option in v.to.rast can select wrong feature for raster attribute,
in areas where features overlap
------------------------+---------------------------------
  Reporter:  florisvdh  |      Owner:  grass-dev@…
      Type:  defect     |     Status:  new
  Priority:  normal     |  Milestone:  7.8.3
 Component:  Vector     |    Version:  svn-releasebranch76
Resolution:             |   Keywords:  v.to.rast sql where
       CPU:  x86-64     |   Platform:  Linux
------------------------+---------------------------------

Comment (by florisvdh):

 **@mmetz**: see the attachment.

 A few examples of sets of polygons which have the same polygon but
 different attributes:

 - cat 21936 and 21937
 - cat 21829 and 21830
 - cat 22114 and 22115

 The original data set from which it was derived (involved more
 [https://github.com/inbo/n2khab-mne-
 design/tree/bd4819099f149ab71346deca1727f5942a88e6c3/030_preparations/010_explore_targetpop
 steps] apart from clipping), is at https://doi.org/10.5281/zenodo.3540740
 .

 Some additional output to show region and a few steps:


 {{{
 $ g.region -p
 projection: 99 (Belge 1972 / Belgian Lambert 72)
 zone:       0
 datum:      bel72
 ellipsoid:  international
 north:      244030.1
 south:      153054.1
 west:       22029.6
 east:       258861.6
 nsres:      32
 ewres:      32
 rows:       2843
 cols:       7401
 cells:      21041043

 $ g.region w=196065 e=199009 s=188877 n=190989 save=debug

 $ g.region -p
 projection: 99 (Belge 1972 / Belgian Lambert 72)
 zone:       0
 datum:      bel72
 ellipsoid:  international
 north:      190989
 south:      188877
 west:       196065
 east:       199009
 nsres:      32
 ewres:      32
 rows:       66
 cols:       92
 cells:      6072

 $ v.clip -r input=polygon_layer output=polygon_layer_debug

 $ v.out.ogr input=polygon_layer_debug output=polygon_layer_debug.gpkg

 }}}

-- 
Ticket URL: <https://trac.osgeo.org/grass/ticket/4009#comment:5>
GRASS GIS <https://grass.osgeo.org>