<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 14 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:Tahoma;
        panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
p
        {mso-style-priority:99;
        mso-margin-top-alt:auto;
        margin-right:0in;
        mso-margin-bottom-alt:auto;
        margin-left:0in;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";}
p.MsoAcetate, li.MsoAcetate, div.MsoAcetate
        {mso-style-priority:99;
        mso-style-link:"Balloon Text Char";
        margin:0in;
        margin-bottom:.0001pt;
        font-size:8.0pt;
        font-family:"Tahoma","sans-serif";}
span.BalloonTextChar
        {mso-style-name:"Balloon Text Char";
        mso-style-priority:99;
        mso-style-link:"Balloon Text";
        font-family:"Tahoma","sans-serif";}
span.EmailStyle20
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-family:"Calibri","sans-serif";}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Nice &#8211; good work, Martin!<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p>&nbsp;</o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>-mpg<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p>&nbsp;</o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p>&nbsp;</o:p></span></p><div style='border:none;border-left:solid blue 1.5pt;padding:0in 0in 0in 4.0pt'><div><div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in'><p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'>From:</span></b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'> Martin Isenburg [mailto:martin.isenburg@gmail.com] <br><b>Sent:</b> Wednesday, April 27, 2011 8:55 AM<br><b>To:</b> mpg@flaxen.com<br><b>Cc:</b> liblas-devel@lists.osgeo.org<br><b>Subject:</b> Re: [Liblas-devel] LASzip versus LT compressor<o:p></o:p></span></p></div></div><p class=MsoNormal><o:p>&nbsp;</o:p></p><p class=MsoNormal style='margin-bottom:12.0pt'>hi,<o:p></o:p></p><div><p class=MsoNormal>On Wed, Apr 27, 2011 at 8:28 AM, Michael P. Gerlek &lt;<a href="mailto:mpg@flaxen.com">mpg@flaxen.com</a>&gt; wrote:<o:p></o:p></p><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;color:#1F497D'>Martin, can you give the data below in terms of compression ratio against the original file sizes?</span><o:p></o:p></p></div></div><div><p class=MsoNormal style='margin-bottom:12.0pt'><br>sure.<br><br>SID&nbsp; LAZ&nbsp;&nbsp; file_name&nbsp; original_file_size (MB)<br><br>3.1&nbsp;&nbsp;&nbsp; 5.9&nbsp;&nbsp;&nbsp; &nbsp;autzen-colorized-1.2-3.las&nbsp;&nbsp;&nbsp; 345<br>2.8&nbsp;&nbsp;&nbsp; 7.0&nbsp;&nbsp;&nbsp; &nbsp;Dallas.las&nbsp;&nbsp;&nbsp; 100<br>6.9&nbsp;&nbsp;&nbsp; 7.2&nbsp;&nbsp;&nbsp; &nbsp;GrassLakeSmall.las&nbsp;&nbsp;&nbsp; 118<br>3.3&nbsp;&nbsp;&nbsp; 8.3&nbsp;&nbsp;&nbsp; &nbsp;IowaDNR-CloudPeakSoft-1.0-UTM15N.las&nbsp;&nbsp;&nbsp; 156<br>6.5&nbsp;&nbsp;&nbsp; 7.4&nbsp;&nbsp;&nbsp; &nbsp;LAS12_Sample_withRGB_QT_Modeler.las&nbsp;&nbsp;&nbsp; 95<br>4.3&nbsp;&nbsp;&nbsp; 4.6&nbsp;&nbsp;&nbsp; &nbsp;LASFile_1.las&nbsp;&nbsp;&nbsp; 46<br>4.5&nbsp;&nbsp;&nbsp; 4.8&nbsp;&nbsp;&nbsp; &nbsp;LASFile_2.las&nbsp;&nbsp;&nbsp; 42<br>4.2&nbsp;&nbsp;&nbsp; 4.6&nbsp;&nbsp;&nbsp; &nbsp;LASFile_3.las&nbsp;&nbsp;&nbsp; 16<br>4.5&nbsp;&nbsp;&nbsp; 4.9&nbsp;&nbsp;&nbsp; &nbsp;LASFile_4.las&nbsp;&nbsp;&nbsp; 46<br>4.4&nbsp;&nbsp;&nbsp; 4.7&nbsp;&nbsp;&nbsp; &nbsp;LDR030828_212242_0.las&nbsp;&nbsp;&nbsp; 57<br>4.5&nbsp;&nbsp;&nbsp; 4.8&nbsp;&nbsp;&nbsp; &nbsp;LDR030828_213023_0.las&nbsp;&nbsp;&nbsp; 56<br>4.3&nbsp;&nbsp;&nbsp; 4.6&nbsp;&nbsp;&nbsp; &nbsp;LDR030828_213450_0.las&nbsp;&nbsp;&nbsp; 51<br>2.8&nbsp;&nbsp;&nbsp; 5.2&nbsp;&nbsp;&nbsp; &nbsp;LDR091111_181233_1.las&nbsp;&nbsp;&nbsp; 52<br>2.8&nbsp;&nbsp;&nbsp; 5.3&nbsp;&nbsp;&nbsp; &nbsp;LDR091111_182803_1.las&nbsp;&nbsp;&nbsp; 52<br>2.8&nbsp;&nbsp;&nbsp; 5.1&nbsp;&nbsp;&nbsp; &nbsp;Ldr100402_220229_1.las&nbsp;&nbsp;&nbsp; 1781<br>6.1&nbsp;&nbsp;&nbsp; 6.5&nbsp;&nbsp;&nbsp; &nbsp;Lincoln.las&nbsp;&nbsp;&nbsp; 177<br>4.0&nbsp;&nbsp;&nbsp; 3.8&nbsp;&nbsp;&nbsp; &nbsp;line_27007_dd.las&nbsp;&nbsp;&nbsp; 103<br>6.3&nbsp;&nbsp;&nbsp; 8.3&nbsp;&nbsp;&nbsp; &nbsp;MARS_Sample_Filtered_LiDAR.las&nbsp;&nbsp;&nbsp; 156<br>2.8&nbsp;&nbsp;&nbsp; 5.2&nbsp;&nbsp;&nbsp; &nbsp;merrick_vertical_1.2.las&nbsp;&nbsp;&nbsp; 52<br>12.9&nbsp;&nbsp;&nbsp; 12.2&nbsp;&nbsp;&nbsp; &nbsp;MountStHelensNov202004.las&nbsp;&nbsp;&nbsp; 110<br>6.4&nbsp;&nbsp;&nbsp; 6.6&nbsp;&nbsp;&nbsp; &nbsp;MountStHelensOct42004.las&nbsp;&nbsp;&nbsp; 129<br>3.1&nbsp;&nbsp;&nbsp; 3.3&nbsp;&nbsp;&nbsp; &nbsp;ncwc000008.las&nbsp;&nbsp;&nbsp; 60<br>6.5&nbsp;&nbsp;&nbsp; 6.8&nbsp;&nbsp;&nbsp; &nbsp;PalmBeachPreHurricane.las&nbsp;&nbsp;&nbsp; 49<br>8.0&nbsp;&nbsp;&nbsp; 8.6&nbsp;&nbsp;&nbsp; &nbsp;radiohead_data1.las&nbsp;&nbsp;&nbsp; 397<br>8.0&nbsp;&nbsp;&nbsp; 8.7&nbsp;&nbsp;&nbsp; &nbsp;radiohead_data2.las&nbsp;&nbsp;&nbsp; 433<br>3.4&nbsp;&nbsp;&nbsp; 7.9&nbsp;&nbsp;&nbsp; &nbsp;S1C1_strip021.las&nbsp;&nbsp;&nbsp; 75<br>3.3&nbsp;&nbsp;&nbsp; 9.1&nbsp;&nbsp;&nbsp; &nbsp;SerpentMoundModelLASData.las&nbsp;&nbsp;&nbsp; 87<br>2.8&nbsp;&nbsp;&nbsp; 5.8&nbsp;&nbsp;&nbsp; &nbsp;Tetons.las&nbsp;&nbsp;&nbsp; 100<br>2.9&nbsp;&nbsp;&nbsp; 5.3&nbsp;&nbsp;&nbsp; &nbsp;USACE_Merrick_lots_of_VLRs.las&nbsp;&nbsp;&nbsp; 96<br>9.7&nbsp;&nbsp;&nbsp; 10.5&nbsp;&nbsp;&nbsp; &nbsp;xyzrgb_manuscript.las&nbsp;&nbsp;&nbsp; 53<o:p></o:p></p></div><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in'><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;color:#1F497D'>&nbsp;Also, are you willing to report the timing data?</span><o:p></o:p></p></div></div></blockquote><div><p class=MsoNormal><br>for the 1.7 GB file it took LASzip 1:33 min to encode and 1:35 min to decode. it took the LT compressor 18:26 min to encode and 4:46 min to decode. your own mpg will vary depending on disk and compressor speeds ...<br><br>cheers,<br><br>martin @lastools<br><span style='font-size:11.0pt;color:#1F497D'>&nbsp;</span><o:p></o:p></p></div><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in'><div><div><div style='border:none;border-left:solid windowtext 1.5pt;padding:0in 0in 0in 4.0pt;border-color:-moz-use-text-color -moz-use-text-color -moz-use-text-color blue'><div><div style='border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0in 0in 0in;border-color:-moz-use-text-color -moz-use-text-color'><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><b><span style='font-size:10.0pt'>From:</span></b><span style='font-size:10.0pt'> <a href="mailto:liblas-devel-bounces@lists.osgeo.org" target="_blank">liblas-devel-bounces@lists.osgeo.org</a> [mailto:<a href="mailto:liblas-devel-bounces@lists.osgeo.org" target="_blank">liblas-devel-bounces@lists.osgeo.org</a>] <b>On Behalf Of </b>Martin Isenburg<br><b>Sent:</b> Wednesday, April 27, 2011 8:12 AM<br><b>To:</b> <a href="mailto:liblas-devel@lists.osgeo.org" target="_blank">liblas-devel@lists.osgeo.org</a><br><b>Subject:</b> [Liblas-devel] LASzip versus LT compressor</span><o:p></o:p></p></div></div><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>&nbsp;<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;margin-bottom:12.0pt'>Hello,<br><br>People sometimes ask me how LASzip compares to the LIDAR compressor from Lizard Tech and I usually refer them to Michael's email (see below). Because his benchmarking was wrong on one model, namely MG4 does not outperform LASzip on 2398_400.las, i did my own experiments that suggest that LASzip compresses about 35% better and is faster.<br><br>A set of 27 LAS files (see below for a listing) compresses to 403 MB with LASzip and to 610 MB with the LIDAR compressor from Lizard Tech.<br><br>A 1.7GB LAS file of a flight swath (see below for the details) compresses to 352 MB with LASzip and to 648 MB with the LIDAR compressor from Lizard Tech. LASzip encoding is about 10 times faster. LASzip decoding is about 3 times faster. timings measurements included all disk I/O from compressed file to uncompressed file (and vice-versa) using two separate drives. disclaimer: the LASzip compressor is tuned for LAS files that contain LIDAR in acquisition order.<br><br>Cheers,<br><br>martin @lastools<br><br>the list of 27 LAS files. the first number is the compressed file size in bytes for the LIDAR Compressor of Lizardtech. the second number is the compressed file size in bytes for LASzip. all files can be found here <a href="http://liblas.org/samples" target="_blank">http://liblas.org/samples</a> except &quot;Dallas.las&quot; and &quot;Tetons.las&quot; which are here: <a href="http://bin.us.lizardtech.com/lidar/LT_LiDAR_Sample_Data.zip" target="_blank">http://bin.us.lizardtech.com/lidar/LT_LiDAR_Sample_Data.zip</a><br><br>115,857,121&nbsp;&nbsp;&nbsp; 61,809,700&nbsp;&nbsp;&nbsp; &nbsp;autzen-colorized-1.2-3.las<br>37,462,947&nbsp;&nbsp;&nbsp; 14,881,473&nbsp;&nbsp;&nbsp; &nbsp;Dallas.las<br>18,035,893&nbsp;&nbsp;&nbsp; 17,128,065&nbsp;&nbsp;&nbsp; &nbsp;GrassLakeSmall.las<br>48,947,417&nbsp;&nbsp;&nbsp; 19,621,507&nbsp;&nbsp;&nbsp; &nbsp;IowaDNR-CloudPeakSoft-1.0-UTM15N.las<br>15,248,628&nbsp;&nbsp;&nbsp; 13,382,538&nbsp;&nbsp;&nbsp; &nbsp;LAS12_Sample_withRGB_QT_Modeler.las<br>11,222,840&nbsp;&nbsp;&nbsp; 10,444,300&nbsp;&nbsp;&nbsp; &nbsp;LASFile_1.las<br>9,805,767&nbsp;&nbsp;&nbsp; 9,154,780&nbsp;&nbsp;&nbsp; &nbsp;LASFile_2.las<br>3,966,869&nbsp;&nbsp;&nbsp; 3,665,433&nbsp;&nbsp;&nbsp; &nbsp;LASFile_3.las<br>10,691,410&nbsp;&nbsp;&nbsp; 9,940,731&nbsp;&nbsp;&nbsp; &nbsp;LASFile_4.las<br>13,522,405&nbsp;&nbsp;&nbsp; 12,672,774&nbsp;&nbsp;&nbsp; &nbsp;LDR030828_212242_0.las<br>13,058,811&nbsp;&nbsp;&nbsp; 12,157,072&nbsp;&nbsp;&nbsp; &nbsp;LDR030828_213023_0.las<br>12,244,190&nbsp;&nbsp;&nbsp; 11,502,895&nbsp;&nbsp;&nbsp; &nbsp;LDR030828_213450_0.las<br>19,820,246&nbsp;&nbsp;&nbsp; 10,414,626&nbsp;&nbsp;&nbsp; &nbsp;LDR091111_181233_1.las<br>19,424,894&nbsp;&nbsp;&nbsp; 10,193,907&nbsp;&nbsp;&nbsp; &nbsp;LDR091111_182803_1.las<br>30,451,604&nbsp;&nbsp;&nbsp; 28,680,682&nbsp;&nbsp;&nbsp; &nbsp;Lincoln.las<br>27,076,520&nbsp;&nbsp;&nbsp; 28,593,056&nbsp;&nbsp;&nbsp; &nbsp;line_27007_dd.las<br>25,753,118&nbsp;&nbsp;&nbsp; 19,594,207&nbsp;&nbsp;&nbsp; &nbsp;MARS_Sample_Filtered_LiDAR.las<br>19,820,246&nbsp;&nbsp;&nbsp; 10,414,626&nbsp;&nbsp;&nbsp; &nbsp;merrick_vertical_1.2.las<br>8,943,713&nbsp;&nbsp;&nbsp; 9,493,209&nbsp;&nbsp;&nbsp; &nbsp;MountStHelensNov202004.las<br>20,937,807&nbsp;&nbsp;&nbsp; 20,337,536&nbsp;&nbsp;&nbsp; &nbsp;MountStHelensOct42004.las<br>20,386,038&nbsp;&nbsp;&nbsp; 19,036,134&nbsp;&nbsp;&nbsp; &nbsp;ncwc000008.las<br>7,900,248&nbsp;&nbsp;&nbsp; 7,539,341&nbsp;&nbsp;&nbsp; &nbsp;PalmBeachPreHurricane.las<br>22,831,603&nbsp;&nbsp;&nbsp; 9,920,385&nbsp;&nbsp;&nbsp; &nbsp;S1C1_strip021.las<br>28,006,302&nbsp;&nbsp;&nbsp; 10,036,738&nbsp;&nbsp;&nbsp; &nbsp;SerpentMoundModelLASData.las<br>37,216,167&nbsp;&nbsp;&nbsp; 18,169,153&nbsp;&nbsp;&nbsp; &nbsp;Tetons.las<br>34,980,891&nbsp;&nbsp;&nbsp; 18,961,597&nbsp;&nbsp;&nbsp; &nbsp;USACE_Merrick_lots_of_VLRs.las<br>5,756,508&nbsp;&nbsp;&nbsp; 5,351,794&nbsp;&nbsp;&nbsp; &nbsp;xyzrgb_manuscript.las<br><br>the lasinfo details of the 1.7GB LAS file containing one swath<br><br>lasinfo Ldr100402_220229_1.laz<br>reporting all LAS header entries:<br>&nbsp; file signature:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 'LASF'<br>&nbsp; file source ID:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0<br>&nbsp; reserved (global_encoding):0<br>&nbsp; project ID GUID data 1-4:&nbsp; 0 0 0 ''<br>&nbsp; version major.minor:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.0<br>&nbsp; system identifier:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 'ALSXX'<br>&nbsp; generating software:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 'ALSXX_PP V2.69 BUILD#7 BETA'<br>&nbsp; file creation day/year:&nbsp;&nbsp;&nbsp; 92/2010<br>&nbsp; header size&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 227<br>&nbsp; offset to point data&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5697<br>&nbsp; number var. length records 4<br>&nbsp; point data format&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>&nbsp; point data record length&nbsp;&nbsp; 28<br>&nbsp; number of point records&nbsp;&nbsp;&nbsp; 66705904<br>&nbsp; number of points by return 58445315 6743224 1404140 113225 0<br>&nbsp; scale factor x y z&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.001 0.001 0.001<br>&nbsp; offset x y z&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 13000000 0 0<br>&nbsp; min x y z&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 12991192.425 588397.501 611.122<br>&nbsp; max x y z&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <a href="tel:13142242.349" target="_blank">13142242.349</a> 594146.283 3032.417<br>variable length header record 1 of 4:<br>&nbsp; reserved&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 43707<br>&nbsp; user ID&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 'LeicaGeo'<br>&nbsp; record ID&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1001<br>&nbsp; length after header&nbsp; 5120<br>&nbsp; description&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ''<br>variable length header record 2 of 4:<br>&nbsp; reserved&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 43707<br>&nbsp; user ID&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 'LeicaGeo'<br>&nbsp; record ID&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1002<br>&nbsp; length after header&nbsp; 22<br>&nbsp; description&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 'MissionInfo'<br>variable length header record 3 of 4:<br>&nbsp; reserved&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 43707<br>&nbsp; user ID&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 'LeicaGeo'<br>&nbsp; record ID&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1003<br>&nbsp; length after header&nbsp; 54<br>&nbsp; description&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 'UserInputs'<br>variable length header record 4 of 4:<br>&nbsp; reserved&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 43707<br>&nbsp; user ID&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 'LASF_Projection'<br>&nbsp; record ID&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 34735<br>&nbsp; length after header&nbsp; 56<br>&nbsp; description&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 'Projection Info'<br>&nbsp;&nbsp;&nbsp; GeoKeyDirectoryTag version 1.1.0 number of keys 6<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; key 1025 tiff_tag_location 0 count 1 value_offset 2 - GTRasterTypeGeoKey: RasterPixelIsPoint<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; key 3076 tiff_tag_location 0 count 1 value_offset 26990 - ProjLinearUnitsGeoKey: look-up for 26990 not implemented<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; key 2052 tiff_tag_location 0 count 1 value_offset 9002 - GeogLinearUnitsGeoKey: Linear_Foot<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; key 4096 tiff_tag_location 0 count 1 value_offset 5103 - VerticalCSTypeGeoKey: VertCS_North_American_Vertical_Datum_1988<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; key 4099 tiff_tag_location 0 count 1 value_offset 9002 - VerticalUnitsGeoKey: Linear_Foot<br>the header is followed by 2 user-defined bytes<br>LASzip compression (version 1.0r0 c1): POINT10 1 GPSTIME11 1<br>reporting minimum and maximum for all LAS point record entries ...<br>&nbsp; x -8807574 142242349<br>&nbsp; y 588397501 594146283<br>&nbsp; z 611122 3032417<br>&nbsp; intensity 0 255<br>&nbsp; edge_of_flight_line 0 0<br>&nbsp; scan_direction_flag 0 1<br>&nbsp; number_of_returns_of_given_pulse 1 4<br>&nbsp; return_number 1 4<br>&nbsp; classification 1 1<br>&nbsp; scan_angle_rank -26 31<br>&nbsp; user_data 161 255<br>&nbsp; point_source_ID 161 511<br>&nbsp; gps_time 511349.016753 512063.402540<br>overview over number of returns of given pulse: 51686151 10678566 3886669 454518 0 0 0<br>histogram of classification of points:<br>&nbsp;66705904 Unclassified (1)<br><br>--------------------------------------------------------------------------------------<br><br>michael's email (the graphs he mentions can be found in the archive) <br><br><a href="http://lists.osgeo.org/pipermail/liblas-devel/2011-February/001199.html" target="_blank">http://lists.osgeo.org/pipermail/liblas-devel/2011-February/001199.html</a><o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>On Fri, Feb 4, 2011 at 3:21 PM, Michael Rosen &lt;<a href="mailto:mrosen@lizardtech.com" target="_blank">mrosen@lizardtech.com</a>&gt; wrote:<o:p></o:p></p><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Here&#8217;s the summary of some LT-internal (I guess not so internal now&#8230;) benchmarking.&nbsp; Highlights:<o:p></o:p></p><p>-<span style='font-size:7.0pt'>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>I can&#8217;t really draw any conclusions about relative compression sizes:&nbsp; 2398_400 favors MG4 2:1, HGAC_Extract and AutZen favor LAZ 2:1, MtStHelens is a wash, <o:p></o:p></p><p>-<span style='font-size:7.0pt'>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>WRT extraction time, for smaller files, the MG4&#8217;s computational overhead (*) favors LAZ for all but the smallest extractions<o:p></o:p></p><p>-<span style='font-size:7.0pt'>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>For larger files, the &#8220;break even&#8221; point is much further to the right.<o:p></o:p></p><p>-<span style='font-size:7.0pt'>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>For larger files, with very small extractions, the built-in index of MG4 allows faster extractions than raw (unindexed) LAS.<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:.25in'>&nbsp;<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>The methodology here was to run &#8220;las2las&#8221; as shown before cropping out increasingly large rectangles (at full resolution)<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>I compared this with the same extraction from MG4 using a command line tool (internal) but this time, writing the output to a las file.<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>&nbsp;<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>I spot checked that the number of points written in all three cases was the same.<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>&nbsp;<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>(*) Note that the title on the graphs is not quite right.&nbsp; It&#8217;s not &#8220;Decode Time&#8221; but &#8220;Decode Time plus LAS Write Time&#8221; vs Scene Size.&nbsp; There is some speculation (based on what we were observing when omitting the output) that LT&#8217;s LAS Writer is unusually slow.&nbsp; It&#8217;s using the liblas v1.2 writer &#8230; so some here may have well-informed opinions on this.<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>&nbsp;<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Here is some raw data and some graphs:<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>&nbsp;<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 03:50 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 61,301,311 2398_400.las<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:28 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8,906,275 2398_400.laz<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:25 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4,650,992 2398_400.sid<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>&nbsp;<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:03 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 362,213,959 autzen-colorized-1.2-3.las<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:28 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 61,809,700 autzen-colorized-1.2-3.laz<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:27 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 115,857,121 autzen-colorized-1.2-3.sid<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>&nbsp;<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 03:59 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 123,876,781 Grass Lake Small.las<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:29 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 17,128,065 Grass Lake Small.laz<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:25 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 18,035,893 Grass Lake Small.sid<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>&nbsp;<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>02/02/2011&nbsp; 08:18 AM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 711,065,603 HGAC_Extract.las<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>02/02/2011&nbsp; 08:23 AM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 151,159,393 HGAC_Extract.laz<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>02/02/2011&nbsp; 08:29 AM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 269,491,108 HGAC_Extract.sid<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>&nbsp;<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 03:50 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 34,065,751 hobu.las<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:29 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 7,732,878 hobu.laz<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:24 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 9,301,431 hobu.sid<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>&nbsp;<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:00 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 185,565,975 Lincoln.las<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:29 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 28,680,682 Lincoln.laz<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:25 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 30,451,604 Lincoln.sid<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>&nbsp;<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 03:58 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 107,603,879 line_27007.las<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:30 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 22,269,252 line_27007.laz<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:25 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 24,588,596 line_27007.sid<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>&nbsp;<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 03:58 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 115,737,877 MtStHelens.las<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:30 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 9,493,209 MtStHelens.laz<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>01/28/2011&nbsp; 04:24 PM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8,943,713 MtStHelens.sid<o:p></o:p></p></div></div></div></div></div></div></div></div></blockquote></div><p class=MsoNormal><o:p>&nbsp;</o:p></p></div></div></body></html>