<div dir="ltr">Thanks Mike for your reply.<div><br><div>After clipping I'm doing a dissolve. The clipping takes 5 min, which is more or less acceptable for 2.8 million shapes.</div><div>But the dissolve is taking over 30 min. for 680k shapes.</div><div><br></div><div>I've added a spatial index using</div><div><font face="monospace, monospace">dsShapefile.ExecuteSQL("CREATE SPATIAL INDEX ON " + layerName, null, "sqlite");</font><br></div><div>but this didn't help with the clipping.</div><div><br></div><div>Before dissolving I add an attribute index:</div><div><font face="monospace, monospace">dsShapefile.ExecuteSQL($"CREATE INDEX ON {layerName} USING Avg1", null, null);</font><br></div><div>and to dissolve I use </div><div><font face="monospace, monospace">.OpenEx(input, 4, null, null, null);<br></font></div><div><font face="monospace, monospace">.wrapper_GDALVectorTranslateDestName();</font><br></div><div><br></div><div>with these options:</div><div><div><font face="monospace, monospace">var options = new[]</font></div><div><font face="monospace, monospace">{</font></div><div><font face="monospace, monospace">    "-f", "ESRI Shapefile",</font></div><div><font face="monospace, monospace">    "-overwrite",</font></div><div><font face="monospace, monospace">    "-explodecollections",</font></div><div><font face="monospace, monospace">    "-dialect", "sqlite",</font></div><div><font face="monospace, monospace">    "-sql", "SELECT ST_Union(geometry), Avg1 FROM {layerName} GROUP BY Avg1"</font></div><div><font face="monospace, monospace">};</font></div></div><div><br></div><div>When I import the shapefile in PostGIS and do the same query it just takes 6 min.</div><div>So I keep thinking I'm not doing it right with GDAL/OGR.</div><div><br></div><div>As Jan Heckman suggested I tried using the Clipper library. Luckily it is implemented in MapWinGIS as well.</div><div>The dissolve using Clipper is just taking 4 min.</div><div>I prefer not to use yet another library.</div><div><br></div><div>So I'm still open for more suggestions.</div><div><br></div></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><br>
<br>
Paul<br>
<br>
<div style="font-size:13.3px;font-family:Verdana,Arial,Helvetica,sans-serif"><div style="border:1px solid #e5e5e5;background-color:#f8f8f8;padding:0.8em 0.8em 0.3em 0.8em"><div style="margin:0 0 8px 0"><p style="margin:0"><span><b><img style="float:left;margin-left:5px;margin-right:5px" src="http://www.bontepaarden.nl/bontepaarden/images/newButton.png" alt="" height="120" width="160">Paul Meems </b><br><span>Release manager, configuration manager</span><br><span>and forum moderator of MapWindow GIS.</span><br><a href="http://www.mapwindow.org/" target="_blank"><span>www.mapwindow.org</span></a><br><br><span>Owner of MapWindow.nl - Support for</span><br><span>Dutch speaking users.</span><br><a href="http://www.mapwindow.nl/" target="_blank"><span>www.mapwindow.nl</span></a><br></span></p>
<p style="margin:0px"><b><br></b></p><p style="margin:0px"><b>The MapWindow GIS project has moved to <a href="https://github.com/MapWindow" target="_blank">GitHub</a>!</b></p><p style="font-size:13.3px;margin:0px"><br><a href="https://github.com/MapWindow/MapWinGIS/releases" target="_blank">Download the latest MapWinGIS mapping engine.</a><br></p><p style="font-size:13.3px;margin:0px"><span style="font-size:13.3px"><a href="https://github.com/MapWindow/MapWindow5/releases" target="_blank">Download the latest MapWindow 5 open source desktop application.</a></span></p></div></div><img src="http://p1.wisestamp.com/pixel.png?p=mozilla&v=3.9.12&t=1322134056890&u=482737fd0c4ca423" height="1" width="1"></div>
</div></div></div></div></div></div></div></div></div>
<br><div class="gmail_quote">2017-09-11 23:34 GMT+02:00 Mike Toews <span dir="ltr"><<a href="mailto:mwtoews@gmail.com" target="_blank">mwtoews@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On 11 September 2017 at 21:09, Paul Meems <<a href="mailto:bontepaarden@gmail.com">bontepaarden@gmail.com</a>> wrote:<br>
> I have a large shapefile with over 2.8 million shapes (fishnet) and I have a<br>
> border file with only 1 multipolygon.<br>
><br>
> I'm trying to clip the fishnet with the border.<br>
> Using code is takes about 5 min. using command line it takes even longer.<br>
><br>
> This is my command:<br>
> ogr2ogr fishnetClipped.shp fishnet.shp -clipsrc border.shp<br>
</span>> ...<br>
<span class="">> Can I somehow improve the speed?<br>
<br>
</span>If the fishnet shapefile does not have a spatial index, try creating<br>
one. There are a few ways to do this:<br>
<a href="http://www.gdal.org/drv_shapefile.html" rel="noreferrer" target="_blank">http://www.gdal.org/drv_<wbr>shapefile.html</a><br>
<br>
I'm not sure if -clipsrc datasource will use the spatial index, but<br>
it's worth trying.<br>
</blockquote></div><br></div>