<div class="gmail_quote">On Thu, Dec 29, 2011 at 12:42 PM, Adam Stylinski <span dir="ltr"><<a href="mailto:stylinae@mail.uc.edu">stylinae@mail.uc.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello,
<div><br></div><div>I'm having issues with LADAR files which do not use the WGS84 coordinate system. I am translating coordinates into longitude and latitude and testing with google earth. Whenever a LAS file has NAD83 coordinate it seems to have a tendency to be completely wrong (most of which end up being somewhere in the ocean). Is there something I'm doing that doesn't work?</div>
<div><br></div><div>I use the GetProj4() to get the proj4 string. I then do a pj_init_plus(proj4str) to get the represented coordinate system, and then do a pj_latlong_from_proj() on that coordinate system to get the lat/long coordinate system. Then to translate to lat/long I do a pj_transform with the src being the original coordinate system, the destination being the lat/long coordinate system, the next two parameters being 1,1, and the x and y being from a point pulled from GetX(), GetY(). Am I doing this wrong? </div>
</blockquote></div><br><div>I should probably mention that I'm failing on this one (though any NAD83 coordinate seems to fail):</div><div><a href="http://liblas.org/samples/IowaDNR-CloudPeakSoft-1.0-UTM15N.las">http://liblas.org/samples/IowaDNR-CloudPeakSoft-1.0-UTM15N.las</a></div>
<div>The GetProj4() returns this: </div><div><div>GetProj4() = +proj=tmerc +lat_0=0 +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=NAD83 +units=m +geoidgrids=g2003conus.gtx,g2003alaska.gtx,g2003h01.gtx,g2003p01.gtx +no_defs</div></div>
<div>And GetWKT():</div><div><div>GetWKT() = COMPD_CS["unknown",</div><div> PROJCS["NAD83 / UTM zone 15N",</div><div> GEOGCS["NAD83",</div><div> DATUM["North_American_Datum_1983",</div>
<div> SPHEROID["GRS 1980",6378137,298.2572221010002,</div><div> AUTHORITY["EPSG","7019"]],</div><div> AUTHORITY["EPSG","6269"]],</div>
<div> PRIMEM["Greenwich",0],</div><div> UNIT["degree",0.0174532925199433],</div><div> AUTHORITY["EPSG","4269"]],</div><div> PROJECTION["Transverse_Mercator"],</div>
<div> PARAMETER["latitude_of_origin",0],</div><div> PARAMETER["central_meridian",0],</div><div> PARAMETER["scale_factor",1],</div><div> PARAMETER["false_easting",0],</div>
<div> PARAMETER["false_northing",0],</div><div> UNIT["metre",1,</div><div> AUTHORITY["EPSG","9001"]],</div><div> AUTHORITY["EPSG","26915"]],</div>
<div> VERT_CS["NAVD88 height",</div><div> VERT_DATUM["North American Vertical Datum 1988",2005,</div><div> AUTHORITY["EPSG","5103"],</div><div> EXTENSION["PROJ4_GRIDS","g2003conus.gtx,g2003alaska.gtx,g2003h01.gtx,g2003p01.gtx"]],</div>
<div> UNIT["metre",1,</div><div> AUTHORITY["EPSG","9001"]],</div><div> AXIS["Up",UP],</div><div> AUTHORITY["EPSG","5703"]]]</div>
</div>