sharpetronics/QGIS
1
0
Fork 0
mirror of https://github.com/qgis/QGIS.git synced 2025-02-25 00:58:06 -05:00
Code Issues Packages Projects Releases Wiki Activity

Test if polygon rings are still inside the outer boundary after reshape and fix a memory leak. Fixes bug #2220

Browse Source 
git-svn-id: http://svn.osgeo.org/qgis/trunk@12585 c8812cc2-4d05-0410-92ff-de0c093fc19c
This commit is contained in:
mhugent 2009-12-22 21:07:32 +00:00
parent a4689579d0
commit f9c0682aba
1 changed files with 37 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

46
src/core/qgsgeometry.cpp
View File

@ -531,6 +531,7 @@ QgsGeometry & QgsGeometry::operator=( QgsGeometry const & rhs )
mGeometrySize = rhs.mGeometrySize;
// deep-copy the GEOS Geometry if appropriate
GEOSGeom_destroy( mGeos );
mGeos = rhs.mGeos ? GEOSGeom_clone( rhs.mGeos ) : 0;
mDirtyGeos = rhs.mDirtyGeos;
@ -5145,22 +5146,49 @@ GEOSGeometry* QgsGeometry::reshapePolygon( const GEOSGeometry* polygon, const GE
newOuterRing = GEOSGeom_clone( outerRing );
}
GEOSGeometry** newInnerRings = new GEOSGeometry*[nRings];
for ( int i = 0; i < nRings; ++i )
//check if all the rings are still inside the outer boundary
QList<GEOSGeometry*> ringList;
if ( nRings > 0 )
{
if ( lastIntersectingRing == i )
GEOSGeometry* outerRingPoly = GEOSGeom_createPolygon( GEOSGeom_clone( newOuterRing ), 0, 0 );
if ( outerRingPoly )
{
newInnerRings[i] = newRing;
}
else
{
newInnerRings[i] = GEOSGeom_clone( innerRings[i] );
GEOSGeometry* currentRing = 0;
for ( int i = 0; i < nRings; ++i )
{
if ( lastIntersectingRing == i )
{
currentRing = newRing;
}
else
{
currentRing = GEOSGeom_clone( innerRings[i] );
}
//possibly a ring is no longer contained in the result polygon after reshape
if ( GEOSContains( outerRingPoly, currentRing ) == 1 )
{
ringList.push_back( currentRing );
}
else
{
GEOSGeom_destroy( currentRing );
}
}
}
GEOSGeom_destroy( outerRingPoly );
}
GEOSGeometry** newInnerRings = new GEOSGeometry*[ringList.size()];
QList<GEOSGeometry*>::const_iterator it = ringList.constBegin();
for ( int i = 0; i < ringList.size(); ++i )
{
newInnerRings[i] = ringList.at( i );
}
delete [] innerRings;
GEOSGeometry* reshapedPolygon = GEOSGeom_createPolygon( newOuterRing, newInnerRings, nRings );
GEOSGeometry* reshapedPolygon = GEOSGeom_createPolygon( newOuterRing, newInnerRings, ringList.size() );
delete[] newInnerRings;
if ( !reshapedPolygon )
{