Neojp, Try AddGeometry with a capital A.
Best regards, Frank On Aug 2, 2012 2:33 PM, "neojp" <[email protected]> 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 > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev >
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
