[gdal-dev] How to iterate over each point of an OGRMultiPoint

[sepideh]

I have this code and I don't know how to iterate over each part ( point ) of
the feature of type wkbMultiPoint:  

   //Multipoint Shapefile
   if ( wkbFlatten ( LayerGeometryType ) == wkbMultiPoint )
   {
	   OGRFeature *poFeature;
	   MultipointFeature MultiPoint;
	   for ( int i = 0; i < NumberOfFeatures; i++ )
	   {
		   poFeature = poLayer ->GetNextFeature();
		   OGRGeometry *poGeometry;
		   poGeometry = poFeature ->GetGeometryRef();
		   if ( poGeometry != NULL && wkbFlatten ( poGeometry ->getGeometryType()
) == wkbMultiPoint )
		   {
			   OGRMultiPoint *poMultipoint = ( OGRMultiPoint * )poGeometry;
			   int NumberOfGeometries = poMultipoint ->getNumGeometries();
			   MultiPoint.PointsOfFeature.resize( NumberOfGeometries );
			   for ( int j = 0; j < NumberOfGeometries; j++ )
			   {
				   poMultipoint ->
			   }
		   }
	   }
   }



--
View this message in context: http://osgeo-org.1560.x6.nabble.com/How-to-iterate-over-each-point-of-an-OGRMultiPoint-tp5077787.html
Sent from the GDAL - Dev mailing list archive at Nabble.com.