[gdal-dev] generating multipolygon from polygons

Frank Warmerdam warmerdam at pobox.com
Thu Aug 2 21:14:06 PDT 2012


Neojp,

Try AddGeometry with a capital A.

Best regards,
Frank
On Aug 2, 2012 2:33 PM, "neojp" <bluejp at gmail.com> wrote:

> Dear List,
>
> I need to generate a multipolygon by adding all the polygons in a
> shapefile.
> Based on a previous post I plan to use the addgeometry() to pass each
> feature to the multipolygon layer. I have implemented this in python but I
> cannot figure out what is going wrong.  Python says "'Geometry' object has
> no attribute 'addGeometry'" that i can understand but not really know how
> to
> fix. Below is the code, I would appreciate any help.
>
>
> import ogr, os, sys
> driver=ogr.GetDriverByName('ESRI Shapefile')
>
> #open shape
> shapefileName= 'c:\\test_multi\\polygonize_binary.shp'
> dataset=driver.Open(shapefileName,0)
> fn = 'c:\\test_multi\\multi_out.shp'
> if os.path.exists(fn):
>   driver.DeleteDataSource(fn)
> outDS = driver.CreateDataSource(fn)
> if outDS is None:
>   print 'Could not create file'
>   sys.exit(1)
>
> #Create a layer to receive the multiple geometry
> #create a layer
> outLayer = outDS.CreateLayer('multi_out', geom_type=ogr.wkbMultiPolygon)
> # create new id and fields in the output shapefile
> #add a new feature
> # create an empty multipolygon geometry
> m_poly = ogr.Geometry(ogr.wkbMultiPolygon)
>
> layer=dataset.GetLayer()
> for index in xrange(layer.GetFeatureCount()):
>     feature = layer.GetFeature(index)
>     geometry = feature.GetGeometryRef()
>     #print geometry
>     m_poly.addGeometry(geometry)
>
>
>
> --
> View this message in context:
> http://osgeo-org.1560.n6.nabble.com/generating-multipolygon-from-polygons-tp4992794.html
> Sent from the GDAL - Dev mailing list archive at Nabble.com.
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20120802/bccf2dc6/attachment.html>


More information about the gdal-dev mailing list