<div dir="ltr"><div><div class="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">
<div>Hi list,</div><div><br></div><div>I have a very large shapefile, with almost 700k shapes which I want to clip with a border shapefile (in red):</div><div><a href="https://ibb.co/gKth1G">https://ibb.co/gKth1G</a></div><div><br></div><div>I'm using this command:</div><div><font face="monospace, monospace">ogr2ogr -clipsrc "border.shp" -overwrite -explodecollections -f "ESRI Shapefile" ogr_clipped.shp Fishnet.shp</font></div><div>with GDAL 2.1.3, released 2017/20/01</div><div>Eventually, I'll transform the command to a call to ogr.VectorTranslate().</div><div><br></div><div>The result is visually correct:</div><div><a href="https://ibb.co/hOvJab">https://ibb.co/hOvJab</a></div><div><br></div><div>But at the red asterix, I have an invalid shape. The reason is according to GEOS:</div><div><font face="monospace, monospace">Polygon must be clockwise</font></div><div>I have more of these errors.</div><div>They all seem to occur when the shape due to the clipping is triangular shaped:</div><div><a href="https://ibb.co/cKuY8w">https://ibb.co/cKuY8w</a></div><div><a href="https://ibb.co/eXcWvb">https://ibb.co/eXcWvb</a></div><div><br></div><div>Is this an error/warning which can be expected or do I need to add an additional parameter?</div><div>If needed I can try to make a test set with as little data as possible.</div><div><br></div><div>Thanks,</div><div><br></div><div>Paul</div></div></div></div></div></div></div>
</div>