<div dir="ltr">Keith,<div><br></div><div>Hmm. It seems two dashe are required in front of options like --a_srs and --t_srs.</div><div><br></div><div>I get:</div><div><br></div><div><div>las2las srs.las --a_srs "+proj=geocent +datum=WGS84" --t_srs EPSG:2151 out.las </div>
<div>error: could not import coordinate system into OSRSpatialReference SetFromUserInput</div></div><div><br></div><div><br></div><div class="gmail_extra"><div class="gmail_extra">$ las2las | head -10</div><div class="gmail_extra">
Input LAS file not specified!</div><div class="gmail_extra">--------------------------------------------------------------------</div><div class="gmail_extra"> las2las (libLAS 1.7.0 with GeoTIFF 1.3.0 GDAL 1.8.1 LASzip 2.1.0)</div>
<div><br></div><div>It appears the GDAL library (1.8.1) linked with las2las lacks support for the +proj=geocent. I'm not sure how to get around that with this build of las2las. If Howard or someone were to rebuild with the stock GDAL 1.9 on OSGeo4W I think it would be ok since the following works:</div>
<div><div>$ gdalsrsinfo "+proj=geocent +datum=WGS84" </div><div><br></div><div>PROJ.4 : '+proj=geocent +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs '</div><div><br></div><div>OGC WKT :</div><div>
GEOCCS["Geocentric",</div><div> DATUM["WGS_1984",</div><div> SPHEROID["WGS 84",6378137,298.257223563,</div><div> AUTHORITY["EPSG","7030"]],</div><div>
TOWGS84[0,0,0,0,0,0,0],</div><div> AUTHORITY["EPSG","6326"]],</div><div> PRIMEM["Greenwich",0,</div><div> AUTHORITY["EPSG","8901"]]]</div></div><div>
<br></div><div>Best regards,</div><div>Frank</div><br><div class="gmail_quote">On Mon, Sep 9, 2013 at 11:26 AM, Keith Williams <span dir="ltr"><<a href="mailto:keithw@unavco.org" target="_blank">keithw@unavco.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div lang="EN-US" link="blue" vlink="purple"><div><p class="">
<span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Hi Frank,<u></u><u></u></span></p><p class=""><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class=""><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Attempted the script that you provided and ended up with the error: too many positional options.<u></u><u></u></span></p><p class="">
<span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class=""><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">C:\Users\keithw\Desktop\TEST>las2las test.las -a_srs "+proj=geocent +datum=WGS84 +units=m +no_defs" -t_srs EPSG:2151<u></u><u></u></span></p>
<p class=""><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">error: too many positional options<u></u><u></u></span></p><p class=""><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class=""><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">C:\Users\keithw\Desktop\TEST>las2las test.las --a_srs "+proj=geocent +datum=WGS84 +units=m +no_defs" --t_srs EPSG:2151<u></u><u></u></span></p>
<p class=""><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">error: unknown option –û<u></u><u></u></span></p><p class=""><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class=""><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Thanks again,<u></u><u></u></span></p><p class=""><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Keith<u></u><u></u></span></p>
<p class=""><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class=""><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class=""><b><span style="font-size:10pt;font-family:Tahoma,sans-serif">From:</span></b><span style="font-size:10pt;font-family:Tahoma,sans-serif"> <a href="mailto:fwarmerdam@gmail.com" target="_blank">fwarmerdam@gmail.com</a> [mailto:<a href="mailto:fwarmerdam@gmail.com" target="_blank">fwarmerdam@gmail.com</a>] <b>On Behalf Of </b>Frank Warmerdam<br>
<b>Sent:</b> Monday, September 09, 2013 12:01 PM<br><b>To:</b> Keith Williams<br><b>Cc:</b> <a href="mailto:liblas-devel@lists.osgeo.org" target="_blank">liblas-devel@lists.osgeo.org</a><br><b>Subject:</b> Re: [Liblas-devel] ECEF?<u></u><u></u></span></p>
<div><div class="h5"><p class=""><u></u> <u></u></p><div><p class="">Keith,<u></u><u></u></p><div><p class=""><u></u> <u></u></p></div><div><p class="">I suspect the problem is that the EPSG catalogs are out of date in OSGeo4W. Could you try the following and see if it works better?<u></u><u></u></p>
</div><div><p class=""><u></u> <u></u></p></div><div><p class=""> las2las MyFile.las –a_srs "+proj=geocent +datum=WGS84 +units=m +no_defs" –t_srs EPSG:2151<u></u><u></u></p></div><div><p class=""><u></u> <u></u></p>
</div><div><p class="">Best regards,<u></u><u></u></p></div><div><p class="">Frank<u></u><u></u></p></div><div><p class=""><u></u> <u></u></p></div></div><div><p class="" style="margin-bottom:12pt"><u></u> <u></u></p><div>
<p class="">On Mon, Sep 9, 2013 at 10:29 AM, Keith Williams <<a href="mailto:keithw@unavco.org" target="_blank">keithw@unavco.org</a>> wrote:<u></u><u></u></p><div><div><p class="">Hello Frank,<u></u><u></u></p><p class="">
<u></u><u></u></p><p class="">I was hoping to request some assistance on converting a LiDAR point cloud from ECEF geocentric coordinates to UTM coordinates.<u></u><u></u></p><p class=""> <u></u><u></u></p><p class="">I noticed that you have previously discussed this at <a href="http://lists.osgeo.org/pipermail/liblas-devel/2011-March/001283.html" target="_blank">http://lists.osgeo.org/pipermail/liblas-devel/2011-March/001283.html</a> however, I am quite new to Liblas, and command line interfaces in general so am not having much luck!<u></u><u></u></p>
<p class=""> <u></u><u></u></p><p class="">I am using: las2las MyFile.las –a_srs EPSG:4328 –t_srs EPSG:2151 through the OSGeo4W shell and I believe that I have the GDAL and liblas trunk files loaded; when installing the OSGeo4W shell I selected every GDAL and Liblas option available. I get the error that “ERROR 6: EPSG PCS/GCS code 4328 not found in EPSG support files. Is this a valid EPSG coordinate system? error: could not import coordinate system into OSRSpatialReference SetFromUserInput”.<u></u><u></u></p>
<p class=""> <u></u><u></u></p><p class="">Thank you very much for any direction you may be willing to supply!<u></u><u></u></p><p class=""> <u></u><u></u></p><p class="">Keith<u></u><u></u></p><p class=""> <u></u><u></u></p>
<p class=""><i><span style="font-size:10pt">----------------------------<br><span style="color:rgb(109,103,67)">Keith Williams, EIT, LSIT<br>UNAVCO, Inc.<br>6350 Nautilus Dr. <br>Boulder, CO 80301</span></span></i><u></u><u></u></p>
<p class=""><i><span style="font-size:10pt;color:rgb(109,103,67)">Cell: <a href="tel:%28541%29%20224-2924" target="_blank">(541) 224-2924</a></span></i><u></u><u></u></p><p class=""><i><span style="font-size:10pt;color:rgb(109,103,67)"><a href="mailto:keithw@unavco.org" target="_blank">keithw@unavco.org</a></span></i><u></u><u></u></p>
<p class=""> <u></u><u></u></p></div></div></div><p class=""><br><br clear="all"><u></u><u></u></p><div><p class=""><u></u> <u></u></p></div><p class="">-- <br>---------------------------------------+--------------------------------------<br>
I set the clouds in motion - turn up | Frank Warmerdam, <a href="mailto:warmerdam@pobox.com" target="_blank">warmerdam@pobox.com</a><br>light and sound - activate the windows | <a href="http://pobox.com/~warmerdam" target="_blank">http://pobox.com/~warmerdam</a><br>
and watch the world go round - Rush | Geospatial Software Developer<u></u><u></u></p></div></div></div></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>---------------------------------------+--------------------------------------<br>
I set the clouds in motion - turn up | Frank Warmerdam, <a href="mailto:warmerdam@pobox.com" target="_blank">warmerdam@pobox.com</a><br>light and sound - activate the windows | <a href="http://pobox.com/~warmerdam" target="_blank">http://pobox.com/~warmerdam</a><br>
and watch the world go round - Rush | Geospatial Software Developer<br>
</div></div>