[GRASS-SVN] r66460 - in grass/trunk/vector/v.distance: . testsuite
svn_grass at osgeo.org
svn_grass at osgeo.org
Fri Oct 9 17:02:55 PDT 2015
Author: wenzeslaus
Date: 2015-10-09 17:02:55 -0700 (Fri, 09 Oct 2015)
New Revision: 66460
Added:
grass/trunk/vector/v.distance/testsuite/
grass/trunk/vector/v.distance/testsuite/test_areas_points.py
Modified:
grass/trunk/vector/v.distance/distance.c
grass/trunk/vector/v.distance/main.c
Log:
v.distance: use bbox in 2D and 3D accordingly (see #2734)
Check if we are considering z and use appropriate point in box function. This should
fix #2734 where point inside an area is classified as outsie of that area when point
3D but the area is 2D. The expected behavior is reduction to 2D.
The test works with this commit but its 3D part fails with code prior to this commit.
This commit depens on r66459 which adds Vect_point_in_box_2d() function.
Modified: grass/trunk/vector/v.distance/distance.c
===================================================================
--- grass/trunk/vector/v.distance/distance.c 2015-10-09 23:51:11 UTC (rev 66459)
+++ grass/trunk/vector/v.distance/distance.c 2015-10-10 00:02:55 UTC (rev 66460)
@@ -289,9 +289,14 @@
all_inside_outer = all_outside_outer = 1;
all_outside_inner = 1;
+ int in_box;
for (i = 0; i < Points->n_points; i++) {
- if (Vect_point_in_box(Points->x[i], Points->y[i],
- Points->z[i], abox)) {
+ if (with_z)
+ in_box = Vect_point_in_box(Points->x[i], Points->y[i],
+ Points->z[i], abox);
+ else
+ in_box = Vect_point_in_box_2d(Points->x[i], Points->y[i], abox);
+ if (in_box) {
int poly;
@@ -323,8 +328,13 @@
int inside_isle = 0;
for (j = 0; j < nisles; j++) {
- if (Vect_point_in_box(Points->x[i], Points->y[i], Points->z[i],
- &ibox[j])) {
+ if (with_z)
+ in_box = Vect_point_in_box(Points->x[i], Points->y[i],
+ Points->z[i], &ibox[j]);
+ else
+ in_box = Vect_point_in_box_2d(Points->x[i],
+ Points->y[i], &ibox[j]);
+ if (in_box) {
poly = Vect_point_in_poly(Points->x[i], Points->y[i], iPoints[j]);
Modified: grass/trunk/vector/v.distance/main.c
===================================================================
--- grass/trunk/vector/v.distance/main.c 2015-10-09 23:51:11 UTC (rev 66459)
+++ grass/trunk/vector/v.distance/main.c 2015-10-10 00:02:55 UTC (rev 66460)
@@ -13,7 +13,7 @@
*
* PURPOSE: Calculates distance from a point to nearest feature in vector layer.
*
- * COPYRIGHT: (C) 2002-2012 by the GRASS Development Team
+ * COPYRIGHT: (C) 2002-2015 by the GRASS Development Team
*
* This program is free software under the
* GNU General Public License (>=v2).
Added: grass/trunk/vector/v.distance/testsuite/test_areas_points.py
===================================================================
--- grass/trunk/vector/v.distance/testsuite/test_areas_points.py (rev 0)
+++ grass/trunk/vector/v.distance/testsuite/test_areas_points.py 2015-10-10 00:02:55 UTC (rev 66460)
@@ -0,0 +1,150 @@
+"""
+TEST: v.distance
+
+AUTHOR(S): Vaclav Petras
+
+PURPOSE: Test v.distance 2D and 3D points with areas
+
+COPYRIGHT: (C) 2015 Vaclav Petras, and by the GRASS Development Team
+
+ This program is free software under the GNU General Public
+ License (>=v2). Read the file COPYING that comes with GRASS
+ for details.
+"""
+
+from grass.gunittest.case import TestCase
+from grass.gunittest.main import test
+from grass.gunittest.gmodules import call_module
+
+areas = """\
+VERTI:
+B 6
+ 11.6 29.55
+ 12.01 27.01
+ 15.25 26.68
+ 16.32 29.42
+ 13.49 31.06
+ 11.6 29.55
+C 1 1
+ 13.85 28.69
+ 1 1
+B 10
+ 18.49 34.1
+ 21.73 28.85
+ 28.24 33.28
+ 26.6 37.17
+ 23.9 34.92
+ 22.75 36.92
+ 18.86 37.17
+ 14.1 36.19
+ 14.19 32.54
+ 18.49 34.1
+C 1 1
+ 21.11 34.03
+ 1 2
+"""
+
+points = """\
+13.45|29.96
+14.64|27.5
+18.04|35.2
+21.56|34.22
+23.2|32.13
+26.48|35.98
+14.76|33.44
+21.85|29.38
+"""
+
+points_3d = """\
+13.45|29.96|200
+14.64|27.5|250
+18.04|35.2|893
+21.56|34.22|350
+23.2|32.13|296
+26.48|35.98|250
+14.76|33.44|258
+21.85|29.38|257
+"""
+
+table_ref = """\
+1|South-west area
+2|South-west area
+3|North-east area
+4|North-east area
+5|North-east area
+6|North-east area
+7|North-east area
+8|North-east area
+"""
+
+
+class TestPointsAndAreas(TestCase):
+ """Test how points get attributes from areas
+
+ Created for #2734 (3D point inside area is classified as outside)
+ https://trac.osgeo.org/grass/ticket/2734
+ """
+ # TODO: replace by unified handing of maps
+ to_remove = []
+ areas = "test_vdistance_areas"
+ points = "test_vdistance_points"
+ points_3d = "test_vdistance_points_3d"
+
+ @classmethod
+ def setUpClass(cls):
+ """Create vector maps (areas with attributes)"""
+ cls.runModule('v.in.ascii', input='-', output=cls.areas,
+ format='standard', stdin=areas)
+ cls.to_remove.append(cls.areas)
+ cls.runModule('v.db.addtable', map=cls.areas,
+ columns="number INTEGER, label VARCHAR(250)")
+ # TODO: this should be done in more effective way
+ cls.runModule('v.db.update', map=cls.areas, column='label',
+ value="South-west area", where="cat=1")
+ cls.runModule('v.db.update', map=cls.areas, column='label',
+ value="North-east area", where="cat=2")
+
+ cls.runModule('v.in.ascii', input='-', output=cls.points,
+ format='point', separator='pipe', flags='t',
+ stdin=points)
+ cls.to_remove.append(cls.points)
+ cls.runModule('v.db.addtable', map=cls.points,
+ columns="area_label VARCHAR(250)")
+
+ cls.runModule('v.in.ascii', input='-', output=cls.points_3d,
+ format='point', separator='pipe', flags='zt', z=3,
+ stdin=points_3d)
+ cls.to_remove.append(cls.points_3d)
+ cls.runModule('v.db.addtable', map=cls.points_3d,
+ columns="area_label VARCHAR(250)")
+
+ @classmethod
+ def tearDownClass(cls):
+ """Remove vector maps"""
+ if cls.to_remove:
+ cls.runModule('g.remove', flags='f', type='vector',
+ name=','.join(cls.to_remove), verbose=True)
+
+ def test_area_attrs_to_2d_points(self):
+ """Check that values are uploaded to 2D points in areas (dmax=0)"""
+ # using call_module because PyGRASS doen't accept form_
+ call_module('v.distance', from_=self.points, to=self.areas,
+ dmax=0, upload='to_attr',
+ column='area_label', to_column='label')
+ # using call_module because it is easier
+ table = call_module('v.db.select', map=self.points,
+ separator='pipe', flags='c')
+ self.assertMultiLineEqual(table, table_ref)
+
+ def test_area_attrs_to_3d_points(self):
+ """Check that values are uploaded to 3D points in areas (dmax=0)"""
+ call_module('v.distance', from_=self.points_3d, to=self.areas,
+ dmax=0, upload='to_attr',
+ column='area_label', to_column='label')
+ table = call_module('v.db.select', map=self.points_3d,
+ separator='pipe', flags='c')
+ self.assertMultiLineEqual(table, table_ref)
+
+
+if __name__ == '__main__':
+ test()
Property changes on: grass/trunk/vector/v.distance/testsuite/test_areas_points.py
___________________________________________________________________
Added: svn:mime-type
+ text/x-python
Added: svn:eol-style
+ native
More information about the grass-commit
mailing list