<!DOCTYPE html>
<html><head>
<meta charset="UTF-8">
</head><body><p>Frank,</p><p>Thank you for this. I had a suspicion feat1 had a purpose at some point. It'd probably be interesting to hear if you ever do remember what that might have been.</p><p>I've made a couple other changes to my version of this since contacting you last week. Precision of the LAT/LONG fields is now 7, as per S-52 Appendix B.2. But then, to add flexibility for the user of the KML file, I ditched my additional second argument to specify metres or feet in favor of creating and filling both a DEPTHM and a DEPTHF field. This gives the user the choice of switching between the two without generating a separate KML file.</p><p>Thanks again,</p><p>Dave</p><blockquote><p>On July 29, 2017 at 10:00 PM Frank Warmerdam <warmerdam@pobox.com> wrote:</p><p>Dave,</p><p>Thanks for the notice. I am not updating the upstream get_soundg.py<br>but I've taken the liberty of cc:ing this reply to the gdal-dev<br>mailing list so more potentially interested folks will be aware of<br>your changes, and your blog post:</p><p><a href="http://micuisine.com/code-projects/?page_id=10">http://micuisine.com/code-projects/?page_id=10</a></p><p>I was particularly amused by change entry "Changed feat2 to feat1 to<br>prevent searching for unused feat1.". Indeed looking back at the code<br>there must have been some history behind my use of feat2 that was lost<br>and now is just confusing.</p><p>Best regards,<br>Frank</p><p>On Sat, Jul 29, 2017 at 7:46 AM, Dave Liske <dave@micuisine.com> wrote:</p><blockquote><p>Mr. Warmerdam,</p><p>In recently looking for a method to split the depth values from NOAA S-57<br>files for a personal art project, I came across your 2003 get_soundg.py file<br>for accomplishing the task via Python. In running your file with<br>us3fl28m.000 I found that it did, in fact, accomplish what I was trying to<br>do. But it also generated a number of errors, the reasons of which are not<br>clear as it likely didn't generate those same errors in 2003 when you<br>originally released the script.</p><p>Sir I know this is out-of-the-blue but I wanted you to know that, with some<br>changes to correct current errors, your 2003 code still does what you<br>intended. I then took it a few more steps to hopefully make it more useful<br>for some as well.</p><p>I've made the following changes, saving the new version as<br>s57_kml_extract_soundg.py:</p><ul><li>Added error handling if SOUNDG layer or source S-57 file is not found.</li><li>Deleted second usage argument: Added code to create target file in same<br>directory as source S-57 file.</li><li>Added new second usage argument: Select soundings as provided as 'metre',<br>or calculate new sounding values as 'feet'</li><li>Target File Type: Changed to KML. Errors are subsequently handled<br>correctly, i.e. the errors are suppressed while the FFPT_RIND and LNAM_REFS<br>fields are also not included in the resulting KML file.</li><li>Added "_metre" or "_feet" to the new filename based on second usage<br>argument</li><li>Added WGS84 CSR to the target KML file.</li><li>Changed created layer name from 'out' to 'KML_SOUNDG'. When imported into<br>QGIS the resulting layer is displayed as <kmlfile>_KML_SOUNDG_Point in the<br>Browser Panel.</li><li>Changed target field 'ELEV' to 'DEPTH', which is accurate for soundings.</li><li>Changed precision for 'DEPTH' field from 4 to 0 to reflect chart usage.</li><li>Added target fields 'LAT' and 'LONG' for xcoord and ycoord.</li><li>Changed feat2 to feat1 to prevent searching for unused feat1.</li><li>Added completion messages.</li></ul><p>I've posted the new code at <a href="http://micuisine.com/code-projects/?page_id=10">http://micuisine.com/code-projects/?page_id=10</a><br>along with further descriptions, as well as a link to a ZIP file containing<br>s57_kml_extract_soundg.py.</p><p>I hope this email finds you well, and that these updates are useful and not<br>offensive in some manner. Please feel free to contact me with any comments<br>or questions about this updated work.</p><p>Sincerely,<br>Dave Liske<br>Luna Pier, Michigan</p></blockquote><p>--<br>---------------------------------------+--------------------------------------<br>I set the clouds in motion - turn up | Frank Warmerdam, warmerdam@pobox.com<br>light and sound - activate the windows |<br>and watch the world go round - Rush | Geospatial Software Developer</p></blockquote></body></html>