sharpetronics/QGIS
1
0
Fork 0
mirror of https://github.com/qgis/QGIS.git synced 2025-04-24 00:04:42 -04:00
Code Issues Packages Projects Releases Wiki Activity

* fix QgsGeometry::moveVertex for multipolygons()

Browse Source 
* add test for QgsGeometry::addPart()
* remove expectedFailure decoration from PyQgsGeometry.testSimplifyIssue4189
* modify move and translate tests by 1,1 to 1,2
This commit is contained in:
Juergen E. Fischer 2014-01-29 21:28:24 +01:00
parent d63104ae50
commit bc757ffe50
3 changed files with 81 additions and 34 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
python/core/qgsgeometry.sip
View File

@ -200,7 +200,7 @@ class QgsGeometry
* Searches for the closest vertex in this geometry to the given point. * Searches for the closest vertex in this geometry to the given point.
* @param point Specifiest the point for search * @param point Specifiest the point for search
* @param atVertex Receives index of the closest vertex * @param atVertex Receives index of the closest vertex
* @return The squared cartesian distance, negative number on error * @return The squared cartesian distance is also returned in sqrDist, negative number on error
*/ */
double closestVertexWithContext( const QgsPoint& point, int& atVertex /Out/ ); double closestVertexWithContext( const QgsPoint& point, int& atVertex /Out/ );
@ -415,6 +415,7 @@ class QgsGeometry
class Error class Error
{ {
public: public:
Error();
Error( QString m ); Error( QString m );
Error( QString m, QgsPoint p ); Error( QString m, QgsPoint p );

15
src/core/qgsgeometry.cpp
View File

@ -1198,13 +1198,16 @@ bool QgsGeometry::moveVertex( QgsWkbPtr &wkbPtr, const double &x, const double &
const int ps = ( hasZValue ? 3 : 2 ) * sizeof( double ); const int ps = ( hasZValue ? 3 : 2 ) * sizeof( double );
// Not this linestring/ring? // Not this linestring/ring?
if ( atVertex > pointIndex + nPoints ) if ( atVertex >= pointIndex + nPoints )
{ {
wkbPtr += ps * nPoints; wkbPtr += ps * nPoints;
pointIndex += nPoints; pointIndex += nPoints;
return false; return false;
} }
if ( isRing && atVertex == pointIndex + nPoints - 1 )
atVertex = pointIndex;
// Goto point in this linestring/ring // Goto point in this linestring/ring
wkbPtr += ps * ( atVertex - pointIndex ); wkbPtr += ps * ( atVertex - pointIndex );
wkbPtr << x << y; wkbPtr << x << y;
@ -2449,7 +2452,7 @@ int QgsGeometry::addPart( const QList<QgsPoint> &points, QGis::GeometryType geom
break; break;
case QGis::Line: case QGis::Line:
// Line needs to have at least two points and must be closed // line needs to have at least two points
if ( points.size() < 2 ) if ( points.size() < 2 )
{ {
QgsDebugMsg( "line must at least have two points: " + QString::number( points.size() ) ); QgsDebugMsg( "line must at least have two points: " + QString::number( points.size() ) );
@ -2458,14 +2461,14 @@ int QgsGeometry::addPart( const QList<QgsPoint> &points, QGis::GeometryType geom
break; break;
case QGis::Polygon: case QGis::Polygon:
// Polygon needs to have at least three points and must be closed // polygon needs to have at least three distinct points and must be closed
if ( points.size() < 3 ) if ( points.size() < 4 )
{ {
QgsDebugMsg( "polygon must at least have three points: " + QString::number( points.size() ) ); QgsDebugMsg( "polygon must at least have three distinct points and must be closed: " + QString::number( points.size() ) );
return 2; return 2;
} }
// polygon must be closed // Polygon must be closed
if ( points.first() != points.last() ) if ( points.first() != points.last() )
{ {
QgsDebugMsg( "polygon not closed" ); QgsDebugMsg( "polygon not closed" );

77
tests/src/python/test_qgsgeometry.py
View File

@ -216,7 +216,6 @@ class TestQgsGeometry(TestCase):
("True", crossesGeom)) ("True", crossesGeom))
assert crossesGeom == True, myMessage assert crossesGeom == True, myMessage
@expectedFailure
def testSimplifyIssue4189(self): def testSimplifyIssue4189(self):
"""Test we can simplify a complex geometry. """Test we can simplify a complex geometry.
@ -802,8 +801,8 @@ class TestQgsGeometry(TestCase):
def testMoveVertex(self): def testMoveVertex(self):
multipoint = QgsGeometry.fromWkt( "MULTIPOINT(5 0,0 0,0 4,5 4,5 1,1 1,1 3,4 3,4 2,2 2)" ) multipoint = QgsGeometry.fromWkt( "MULTIPOINT(5 0,0 0,0 4,5 4,5 1,1 1,1 3,4 3,4 2,2 2)" )
for i in range(0,10): for i in range(0,10):
assert multipoint.moveVertex( i+1, i+1, i ), "move vertex %d failed" % i assert multipoint.moveVertex( i+1, -1-i, i ), "move vertex %d failed" % i
expwkt = "MULTIPOINT(1 1, 2 2, 3 3, 4 4, 5 5, 6 6, 7 7, 8 8, 9 9, 10 10)" expwkt = "MULTIPOINT(1 -1, 2 -2, 3 -3, 4 -4, 5 -5, 6 -6, 7 -7, 8 -8, 9 -9, 10 -10)"
wkt = multipoint.exportToWkt() wkt = multipoint.exportToWkt()
assert compareWkt( expwkt, wkt ), "Expected:\n%s\nGot:\n%s\n" % (expwkt, wkt ) assert compareWkt( expwkt, wkt ), "Expected:\n%s\nGot:\n%s\n" % (expwkt, wkt )
@ -833,12 +832,15 @@ class TestQgsGeometry(TestCase):
wkt = polygon.exportToWkt() wkt = polygon.exportToWkt()
assert compareWkt( expwkt, wkt ), "Expected:\n%s\nGot:\n%s\n" % (expwkt, wkt ) assert compareWkt( expwkt, wkt ), "Expected:\n%s\nGot:\n%s\n" % (expwkt, wkt )
assert polygon.moveVertex( 1, 1, 0 ), "move vertex failed" assert polygon.moveVertex( 1, 2, 0 ), "move vertex failed"
expwkt = "MULTIPOLYGON(((1 1,1 0,1 1,2 1,2 2,0 2,1 1)),((4 0,5 0,6 2,3 2,3 1,4 1,4 0)))" expwkt = "MULTIPOLYGON(((1 2,1 0,1 1,2 1,2 2,0 2,1 2)),((4 0,5 0,6 2,3 2,3 1,4 1,4 0)))"
wkt = polygon.exportToWkt()
assert polygon.moveVertex( 1, 1, 8 ), "move vertex failed" assert compareWkt( expwkt, wkt ), "Expected:\n%s\nGot:\n%s\n" % (expwkt, wkt )
expwkt = "MULTIPOLYGON(((1 1,1 0,1 1,2 1,2 2,0 2,1 1)),((1 1,5 0,6 2,3 2,3 1,4 1,1 1)))"
assert polygon.moveVertex( 2, 1, 7 ), "move vertex failed"
expwkt = "MULTIPOLYGON(((1 2,1 0,1 1,2 1,2 2,0 2,1 2)),((2 1,5 0,6 2,3 2,3 1,4 1,2 1)))"
wkt = polygon.exportToWkt()
assert compareWkt( expwkt, wkt ), "Expected:\n%s\nGot:\n%s\n" % (expwkt, wkt )
def testDeleteVertex(self): def testDeleteVertex(self):
# 2-+-+-+-+-3 # 2-+-+-+-+-3
@ -957,27 +959,26 @@ class TestQgsGeometry(TestCase):
def testTranslate(self): def testTranslate(self):
point = QgsGeometry.fromWkt( "POINT(1 1)" ) point = QgsGeometry.fromWkt( "POINT(1 1)" )
assert point.translate( 1, 1 )==0, "Translate failed" assert point.translate( 1, 2 )==0, "Translate failed"
expwkt = "POINT(2 2)" expwkt = "POINT(2 3)"
wkt = point.exportToWkt() wkt = point.exportToWkt()
assert compareWkt( expwkt, wkt ), "Expected:\n%s\nGot:\n%s\n" % (expwkt, wkt ) assert compareWkt( expwkt, wkt ), "Expected:\n%s\nGot:\n%s\n" % (expwkt, wkt )
point = QgsGeometry.fromWkt( "MULTIPOINT(1 1,2 2,3 3)" ) point = QgsGeometry.fromWkt( "MULTIPOINT(1 1,2 2,3 3)" )
assert point.translate( 1, 1 )==0, "Translate failed" assert point.translate( 1, 2 )==0, "Translate failed"
expwkt = "MULTIPOINT(2 2, 3 3, 4 4)" expwkt = "MULTIPOINT(2 3, 3 4, 4 5)"
wkt = point.exportToWkt() wkt = point.exportToWkt()
assert compareWkt( expwkt, wkt ), "Expected:\n%s\nGot:\n%s\n" % (expwkt, wkt ) assert compareWkt( expwkt, wkt ), "Expected:\n%s\nGot:\n%s\n" % (expwkt, wkt )
linestring = QgsGeometry.fromWkt( "LINESTRING(1 0,2 0)" ) linestring = QgsGeometry.fromWkt( "LINESTRING(1 0,2 0)" )
assert linestring.translate( 1, 2 )==0, "Translate failed"
assert linestring.translate( 1, 1 )==0, "Translate failed" expwkt = "LINESTRING(2 2, 3 2)"
expwkt = "LINESTRING(2 1, 3 1)"
wkt = linestring.exportToWkt() wkt = linestring.exportToWkt()
assert compareWkt( expwkt, wkt ), "Expected:\n%s\nGot:\n%s\n" % (expwkt, wkt ) assert compareWkt( expwkt, wkt ), "Expected:\n%s\nGot:\n%s\n" % (expwkt, wkt )
polygon = QgsGeometry.fromWkt( "MULTIPOLYGON(((0 0,1 0,1 1,2 1,2 2,0 2,0 0)),((4 0,5 0,5 2,3 2,3 1,4 1,4 0)))" ) polygon = QgsGeometry.fromWkt( "MULTIPOLYGON(((0 0,1 0,1 1,2 1,2 2,0 2,0 0)),((4 0,5 0,5 2,3 2,3 1,4 1,4 0)))" )
assert polygon.translate( 1, 1 )==0, "Translate failed" assert polygon.translate( 1, 2 )==0, "Translate failed"
expwkt = "MULTIPOLYGON(((1 1,2 1,2 2,3 2,3 3,1 3,1 1)),((5 1,6 1,6 1,4 3,4 2,5 2,5 1)))" expwkt = "MULTIPOLYGON(((1 2,2 2,2 3,3 3,3 4,1 4,1 2)),((5 2,6 2,6 2,4 4,4 3,5 3,5 2)))"
wkt = polygon.exportToWkt() wkt = polygon.exportToWkt()
ct = QgsCoordinateTransform() ct = QgsCoordinateTransform()
@ -1077,5 +1078,47 @@ class TestQgsGeometry(TestCase):
bb = polygon.boundingBox() bb = polygon.boundingBox()
assert expbb == bb, "Expected:\n%s\nGot:\n%s\n" % (expbb.toString(), bb.toString()) assert expbb == bb, "Expected:\n%s\nGot:\n%s\n" % (expbb.toString(), bb.toString())
def testAddPart(self):
# 2-3 6-+-7
# | | | |
# 0-1 4 5 8-9
points = [
[ QgsPoint(0,0), QgsPoint(1,0), QgsPoint(1,1), QgsPoint(2,1), QgsPoint(2,0), ],
[ QgsPoint(3,0), QgsPoint(3,1), QgsPoint(5,1), QgsPoint(5,0), QgsPoint(6,0), ]
]
polyline = QgsGeometry.fromPolyline( points[0] )
assert polyline.addPart( points[1][0:1] ) == 2, "addPart with one point line unexpectedly succeeded."
assert polyline.addPart( points[1][0:2] ) == 0, "addPart with two point line failed."
expwkt = "MULTILINESTRING((0 0, 1 0, 1 1, 2 1, 2 0), (3 0, 3 1))"
wkt = polyline.exportToWkt()
assert compareWkt( expwkt, wkt ), "Expected:\n%s\nGot:\n%s\n" % (expwkt, wkt )
polyline = QgsGeometry.fromPolyline( points[0] )
assert polyline.addPart( points[1] ) == 0, "addPart with %d point line failed." % len(points[1])
expwkt = "MULTILINESTRING((0 0, 1 0, 1 1, 2 1, 2 0), (3 0, 3 1, 5 1, 5 0, 6 0))"
# 5-+-4 0-+-9
# | | | |
# | 2-3 1-2 |
# | | | |
# 0-1 7-8
points = [
[ [ QgsPoint(0,0), QgsPoint(1,0), QgsPoint(1,1), QgsPoint(2,1), QgsPoint(2,2), QgsPoint(0,2), QgsPoint(0,0), ] ],
[ [ QgsPoint(4,0), QgsPoint(5,0), QgsPoint(5,2), QgsPoint(3,2), QgsPoint(3,1), QgsPoint(4,1), QgsPoint(4,0), ] ]
]
polygon = QgsGeometry.fromPolygon( points[0] )
assert polygon.addPart( points[1][0][0:1] ) == 2, "addPart with one point ring unexpectedly succeeded."
assert polygon.addPart( points[1][0][0:2] ) == 2, "addPart with two point ring unexpectedly succeeded."
assert polygon.addPart( points[1][0][0:3] ) == 2, "addPart with unclosed three point ring unexpectedly succeeded."
assert polygon.addPart( [ QgsPoint(4,0), QgsPoint(5,0), QgsPoint(4,0) ] ) == 2, "addPart with 'closed' three point ring unexpectedly succeeded."
assert polygon.addPart( points[1][0] ) == 0, "addPart failed"
expwkt = "MULTIPOLYGON(((0 0,1 0,1 1,2 1,2 2,0 2,0 0)),((4 0,5 0,5 2,3 2,3 1,4 1,4 0)))"
wkt = polygon.exportToWkt()
assert compareWkt( expwkt, wkt ), "Expected:\n%s\nGot:\n%s\n" % (expwkt, wkt )
if __name__ == '__main__': if __name__ == '__main__':
unittest.main() unittest.main()