<div dir="ltr"><div class="gmail_quote"><div class="gmail_attr">(Forwarding to PG User lisdt thread as well)</div><div class="gmail_attr"><br></div><div dir="ltr"><div>Good news - overlaying the full dataset works using JTS OverlayNG in Snap-Rounding mode.  It completes in 195 s using a grid size of 1e-8 (which is quite fine).  </div><div><br></div><div>Interestingly it didn't complete when using OverlayNG with Snapping noding (which is  heuristic-based noding strategy intended to preserve more precision).  I'll have to look into that further.</div><div><br></div><div>Here's some instructions on how to run this in JTS:</div><div><br></div><div>- Clone the branch <a href="https://github.com/dr-jts/jts/tree/overlay-sr" target="_blank">https://github.com/dr-jts/jts/tree/overlay-sr</a></div><div>- Build using maven clean install</div><div>- Run the TestBuilder (<a href="https://github.com/locationtech/jts/blob/master/doc/JTSTestBuilder.md" target="_blank">https://github.com/locationtech/jts/blob/master/doc/JTSTestBuilder.md</a>)</div><div>- Convert the SQL file to a WKB file by removing the SQL header and footer, leaving only the WKB text</div><div>- Drop the WKB file onto the TestBuilder A input box to load it</div><div>- Run the geometry function OverlayNG/union with Grid Scale = 100000000</div><div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jul 7, 2020 at 10:49 PM Jonathan Schultz <<a href="mailto:jonathan@schultz.la" target="_blank">jonathan@schultz.la</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Oh sorry about that. Here is a link to the file on my own VPS: <a href="http://barraqda.org/dbca_dump_116_35_117_34.zip" rel="noreferrer" target="_blank">http://barraqda.org/dbca_dump_116_35_117_34.zip</a><br>
<br>
I've been playing with JTS but it's a bit of a steep learning curve. Getting stuck on Java NoClassDefFoundError messages, I presume related to incompatible versions and such-like. I'd very much appreciate advice on running OverlayNG.<br>
<br>
Cheers,<br>
Jonathan<br>
<br>
On 8/7/20 1:41 pm, Martin Davis wrote:<br>
> That's interesting that ST_SnapToGrid doesn't help for the full dataset - but not totally surprising, since it's a pretty blunt instrument.<br>
> <br>
> I'd like to try out the larger dataset.  That free download website seems to have a bunch of unpleasant strings attached though - is there any other way you can make the data available?<br>
> <br>
> If you are able to run Java, by all means try OverlayNG.  Happy to provide advice on how to run it if needed.<br>
> <br>
> cheers - Martin<br>
> <br>
> On Tue, Jul 7, 2020 at 8:10 PM Jonathan Schultz <<a href="mailto:jonathan@schultz.la" target="_blank">jonathan@schultz.la</a> <mailto:<a href="mailto:jonathan@schultz.la" target="_blank">jonathan@schultz.la</a>>> wrote:<br>
> <br>
>     Thank you for the prompt feedback.<br>
> <br>
>     As you suggested I tried using ST_SnapToGrid and that worked very well even with a very fine grid - st_snaptogrid(geometry,0.000000000001). However when it came to a collection of lines more useful for real work (13540 records compared with 6 in the minimal dataset) I found that snapping even to a grid so coarse as to render the data useless - st_snaptogrid(geometry,0.01) - still produced an incorrect result.<br>
> <br>
>     Here is the dump of my larger dataset: <a href="https://anonfiles.com/rdYcG1Edo2/dbca_dump_116_35_117_34_sql" rel="noreferrer" target="_blank">https://anonfiles.com/rdYcG1Edo2/dbca_dump_116_35_117_34_sql</a><br>
> <br>
>     I'll have a go myself with OverlayNG.<br>
> <br>
>     And I will file a PostGIS issue depending on the conclusions from this discussion.<br>
> <br>
> <br>
> <br>
<br>
</blockquote></div></div>
</div></div>