OGRBuildPolygonFromEdges() segfaults on MULTIPOLYGON inputs

[dron]

OGRBuildPolygonFromEdges() can't properly handle a case when the MULTIPOLYGON geometry supplied as an input. The test case follows:

#!/usr/bin/env python

from osgeo import ogr

g=ogr.Geometry(wkt="MULTIPOLYGON (((20 15,15 15,15 20,20 20,20 15)),((17 10,17 12,20 12, 20 10,17 10)))")
ogr.BuildPolygonFromEdges(g, bAutoClose=1)

This script results in segmentation fault. I think that OGRBuildPolygonFromEdges() should either refuse to work with such a geometry or return a multipolygon as a result.

Best regards, Andrey

[warmerdam]

Agreed.

Chaitanya, could you please modify the function (in gdal/ogr/ograssemblepolygon.cpp to error out in an informative way for non-linestring (or linearring) geometries in the passed in collection? Also, please extend the test suite with this change. It would be sufficient to do this in trunk.

[Even Rouault]

Fixed in trunk (r17556). Tests added in r17557