<div dir="ltr"><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">
Hi List,<br><br>I have a large shapefile (land use) with a lot of small shapes (392k).<br>I also have a smaller shapefile (catchments) with larger shapes (2.5k).<br>The smaller shapefile overlaps the larger shapefile in the center.<br><br>I need to create a new shapefile (expected number of shapes is several 100k) that have the values of the land use and the catchments combined. <br>Currently, I'm doing a GEOSClip() of the land use by the catchment first and then a GEOSUnion of the clipped result and the catchment.<br><br>This is a bit slow. The clipping already takes 40 minutes.<br><br>Am I using the right approach or should I use a different one?<br><br>Extra info. I'm using MapWinGIS (C++) which includes GDAL and GEOS.<br>This is the implementation of GEOSClip:<br><a href="https://github.com/MapWindow/MapWinGIS/blob/develop/src/COM%20classes/Shapefile_Geoprocessing.cpp#L1936">https://github.com/MapWindow/MapWinGIS/blob/develop/src/COM%20classes/Shapefile_Geoprocessing.cpp#L1936</a><br>It was created years ago. Perhaps the implementation can/should be optimized as well.<br><br>Thanks,<br></div><div dir="ltr"><br></div><div>Paul Meems</div><div><br></div></div></div></div></div></div></div>