<div dir="ltr"><div dir="ltr"><div dir="ltr">Hi Aruna,<div><br></div><div>This would be a good question to post on <a href="https://gis.stackexchange.com">https://gis.stackexchange.com</a> , there are good technical GIS experts there from around the world, and there's points to be earned so good incentives for properly resolving. Also check if a similar question has already been posted there.</div><div><br></div><div>Taking a shot:</div><div>I'm not sure if the "right hand rule" validation/rectification for geojson polygons would work here or not, but it did help me in a calculation function using python's geojson_utils module. The area was computing to be negative, and reversing the order of the co-ordinates of the polygon set it right.</div><div><br></div><div>More theory on Right Hand Rule and rectifier tool available here: <a href="https://mapster.me/right-hand-rule-geojson-fixer/">https://mapster.me/right-hand-rule-geojson-fixer/</a></div><div>Validation site that flags if your polygons are adhering to the right hand rule or not: <a href="http://geojsonlint.com">http://geojsonlint.com</a></div><div><br></div><div>My brief understanding: It's a validation for geojson that has been made non-compulsory in recent times, but still causes problems at some places. When traversing through a polygon, the math formulae assume one side to be 'in' and the other side 'out'.</div><div>A quick rectification without having to alter the polygon data : if area computes to be negative (<0), then multiply by -1 to make it positive. But I don't know how this may apply to your specific function calculation.</div><div><br></div><div>I ran your geojson through the tool and it has reversed the order of your polygons' co-ordinates to solve for right hand rule. So try the reversed version in your program and let us know if you were able to resolve your problem.</div><div><br></div><div>If this works, then please post this question on <a href="http://gis.stackexchange.com">gis.stackexchange.com</a> and share the link so I can post the right hand rule answer there! </div><div>If this doesn't work, then definitely post this on <a href="http://gis.stackexchange.com">gis.stackexchange.com</a> all the same :)</div><div><br clear="all"><div><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>--<br>Cheers,<br>Nikhil VJ<br>+91-966-583-1250<br>Pune, India<br><a href="http://nikhilvj.co.in" target="_blank">http://nikhilvj.co.in</a></div></div></div></div></div></div></div></div></div></div><br></div></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Thu, Nov 8, 2018 at 8:21 AM Aruna S <<a href="mailto:safincrazy@gmail.com">safincrazy@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Hi everyone!</div><div><br></div><div>I am working with the mbtiles-extracts library: <a href="https://github.com/mapbox/mbtiles-extracts" target="_blank">https://github.com/mapbox/mbtiles-extracts</a>. Given an mbtiles file and a geojson, the library extracts mbtiles cutouts of the polygons in the geojson.</div><div><br></div><div>Except, it does not work for some geojsons.</div><div><br></div><div>As documented in <a href="https://github.com/mapbox/mbtiles-extracts/issues/10" target="_blank">https://github.com/mapbox/mbtiles-extracts/issues/10</a>, it does not work for any polygon present in the following geojson for the mbtiles listed below:<br></div><div><ul><li>MBTILES: Pretty much every tileset on <a href="https://osmlab.github.io/osm-qa-tiles/historic.html" rel="nofollow" target="_blank">https://osmlab.github.io/osm-qa-tiles/historic.html</a>. Can definitely replicate it using for <a href="https://s3.amazonaws.com/mapbox/osm-qa-tiles/historic/2008-Q1-qa-tiles.mbtiles" rel="nofollow" target="_blank">https://s3.amazonaws.com/mapbox/osm-qa-tiles/historic/2008-Q1-qa-tiles.mbtiles</a></li><li>GEOJSON: <a href="https://gist.github.com/arunasank/c98ef28abd51b54579c7609f8ebc1de1" target="_blank">https://gist.github.com/arunasank/c98ef28abd51b54579c7609f8ebc1de1</a></li></ul><div>When I debug, I see that:</div><div>* This function: <a href="https://github.com/mapbox/mbtiles-extracts/blob/master/index.js#L175" target="_blank">https://github.com/mapbox/mbtiles-extracts/blob/master/index.js#L175</a></div><div>returns a falsy value for all the features in the above geojson, resulting in a falsy value here: <a href="https://github.com/mapbox/mbtiles-extracts/blob/master/index.js#L58-L62" target="_blank">https://github.com/mapbox/mbtiles-extracts/blob/master/index.js#L58-L62</a></div><div><br></div><div>Any help in understanding why these polygons are different from other polygons, or if there's a pattern I am missing, would be really helpful!</div><div><br></div><div>If there are other mailing lists where this could be answered, please point me to them.</div><div><br></div><div>Warmly,</div><div>Aruna<br></div></div></div></div></div></div></div>
_______________________________________________<br>
Talk-in mailing list<br>
<a href="mailto:Talk-in@openstreetmap.org" target="_blank">Talk-in@openstreetmap.org</a><br>
<a href="https://lists.openstreetmap.org/listinfo/talk-in" rel="noreferrer" target="_blank">https://lists.openstreetmap.org/listinfo/talk-in</a><br>
</blockquote></div>