<div class="gmail_quote">On Tue, Jan 10, 2012 at 1:25 PM, Howard Butler <span dir="ltr"><<a href="mailto:hobu.inc@gmail.com">hobu.inc@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
and LAStools now does vertical datum transformations too? If so, even more reason to not use libLAS :)<br>
<div><div></div><div class="h5"><br>
On Jan 10, 2012, at 12:10 PM, Martin Isenburg wrote:<br>
<br>
> Hi,<br>
><br>
> when debugging GDAL it may help to look at what LASlib's<br>
> geoprojectionconverter does, which finds that the file<br>
> "IowaDNR-CloudPeakSoft-1.0-UTM15N.las" has correct GeoTIFFs . When I<br>
> create an intensity grid with<br>
><br>
> lasgrid -i IowaDNR-CloudPeakSoft-1.0-UTM15N.las -intensity -opng -step 3<br>
><br>
> the generated intensity raster is automatically geo-referenced and a<br>
> corresponding KML file is generated that places the PNG onto Google<br>
> Earth at the exact right spot as you can see here:<br>
><br>
> <a href="http://lastools.org/download/data/IowaDNR-CloudPeakSoft-1.0-UTM15N.kmz" target="_blank">http://lastools.org/download/data/IowaDNR-CloudPeakSoft-1.0-UTM15N.kmz</a><br>
><br>
> Cheers,<br>
><br>
> Martin @lastools<br>
><br>
> PS: below is the report of lasinfo -i IowaDNR-CloudPeakSoft-1.0-UTM15N.las<br>
><br>
> reporting all LAS header entries:<br>
> file signature: 'LASF'<br>
> file source ID: 0<br>
> global_encoding: 0<br>
> project ID GUID data 1-4: 0 0 0 ''<br>
> version major.minor: 1.0<br>
> system identifier: ''<br>
> generating software: 'CPS/RTW LAS Lib v1.08'<br>
> file creation day/year: 0/0<br>
> header size 227<br>
> offset to point data 639<br>
> number var. length records 3<br>
> point data format 1<br>
> point data record length 28<br>
> number of point records 5847380<br>
> number of points by return 2920986 2926394 0 0 0<br>
> scale factor x y z 0.001 0.001 0.001<br>
> offset x y z 0 4000000 0<br>
> min x y z 614000.000 4791999.999 311.676<br>
> max x y z 615999.990 4793999.989 398.959<br>
> variable length header record 1 of 3:<br>
> reserved 43707<br>
> user ID 'LASF_Projection'<br>
> record ID 34735<br>
> length after header 88<br>
> description 'Projection Info'<br>
> GeoKeyDirectoryTag version 1.1.0 number of keys 10<br>
> key 1024 tiff_tag_location 0 count 1 value_offset 1 -<br>
> GTModelTypeGeoKey: ModelTypeProjected<br>
> key 1026 tiff_tag_location 34737 count 13 value_offset 0 -<br>
> GTCitationGeoKey: UTM zone 15N<br>
> key 2050 tiff_tag_location 0 count 1 value_offset 6269 -<br>
> GeogGeodeticDatumGeoKey: Datum_North_American_Datum_1983<br>
> key 2056 tiff_tag_location 0 count 1 value_offset 7019 -<br>
> GeogEllipsoidGeoKey: Ellipse_GRS_1980<br>
> key 3072 tiff_tag_location 0 count 1 value_offset 26915 -<br>
> ProjectedCSTypeGeoKey: PCS_NAD83_UTM_zone_15N<br>
> key 3075 tiff_tag_location 0 count 1 value_offset 1 -<br>
> ProjCoordTransGeoKey: CT_TransverseMercator<br>
> key 3076 tiff_tag_location 0 count 1 value_offset 9001 -<br>
> ProjLinearUnitsGeoKey: Linear_Meter<br>
> key 4096 tiff_tag_location 0 count 1 value_offset 5103 -<br>
> VerticalCSTypeGeoKey: VertCS_North_American_Vertical_Datum_1988<br>
> key 4097 tiff_tag_location 34737 count 19 value_offset 13 -<br>
> VerticalCitationGeoKey: GEOID 2003 Applied<br>
> key 4099 tiff_tag_location 0 count 1 value_offset 9001 -<br>
> VerticalUnitsGeoKey: Linear_Meter<br>
> variable length header record 2 of 3:<br>
> reserved 43707<br>
> user ID 'LASF_Projection'<br>
> record ID 34737<br>
> length after header 32<br>
> description 'Projection Info'<br>
> GeoAsciiParamsTag (number of characters 32)<br>
> UTM zone 15N GEOID 2003 Applied<br>
> variable length header record 3 of 3:<br>
> reserved 43707<br>
> user ID 'CloudPeakSoft'<br>
> record ID 0<br>
> length after header 128<br>
> description 'CPSInfo'<br>
> the header is followed by 2 user-defined bytes<br>
> reporting minimum and maximum for all LAS point record entries ...<br>
> x 614000000 615999990<br>
> y 791999999 793999989<br>
> z 311676 398959<br>
> intensity 1 3240<br>
> edge_of_flight_line 0 0<br>
> scan_direction_flag 0 0<br>
> number_of_returns_of_given_pulse 2 2<br>
> return_number 1 2<br>
> classification 1 12<br>
> scan_angle_rank 0 0<br>
> user_data 0 0<br>
> point_source_ID 14 18<br>
> gps_time 597987.305116 600731.680834<br>
> overview over number of returns of given pulse: 0 5847380 0 0 0 0 0<br>
> histogram of classification of points:<br>
> 2762810 Unclassified (1)<br>
> 2156929 Ground (2)<br>
> 6318 Low Point (noise) (7)<br>
> 921323 Overlap Points (12)<br>
<br>
_______________________________________________<br>
Liblas-devel mailing list<br>
<a href="mailto:Liblas-devel@lists.osgeo.org">Liblas-devel@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/liblas-devel" target="_blank">http://lists.osgeo.org/mailman/listinfo/liblas-devel</a><br>
<br>
</div></div></blockquote></div><div><br></div><div>I should probably mention that my proj.4 is built with datum grids as well. Liblas is built with GDAL and geotiff as well.</div>