<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=utf-8">
<meta name="Generator" content="Microsoft Word 15 (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:DengXian;
panose-1:3 0 5 9 0 0 0 0 0 0;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"\@DengXian";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
span.EmailStyle21
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@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" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal">Hi Nicolas, thank you for the well researched and detailed answer. I was running on the latest build and it had a the new feature for Cloud Optimization COPC. That feature doubled the number of open files for the process and I dropped
back to 3.22 to get my processing done. It appears that during the indexing process of 3.28 the COPC files are left open even after moving on to another LAZ file. The 3.22 level did not create the COPC files and everything worked out for the images.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b>From:</b> Nicolas Cadieux <njacadieux.gitlab@gmail.com> <br>
<b>Sent:</b> Friday, February 17, 2023 7:20 AM<br>
<b>To:</b> Mike Stanton <mfstanton@geospan.com><br>
<b>Cc:</b> qgis-user@lists.osgeo.org<br>
<b>Subject:</b> Re: [Qgis-user] Two LAZ file possible issues iwth some new data we received to processs<o:p></o:p></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Hi,<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">I did run into this layer number limite in the past. If you check the user list or if you Google “QGIS, maximum amount of layers… “ you will find multiple post. The issue, at the time was caused by limites in the OS. (Windows 7 at the
time), OSX had and easy work around. Same thing for the various flavours of Linux.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">At the time, I either combined files when using QGIS or opened them in Saga that did not have this problem. I believe this was because Saga, at least back then, loaded all the file into memory and did not keep a handle on the files once
is was open.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">This is the type of thing that got me into learning Python so that I could create loops and speed up operations on LiDAR data.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">It’s normal (or at least expected) for large batches of file to have bugs in them like missing crs in headers. File transfers can be difficult to handle. Talk with the provider to see if they have the same issue. If you know a bit of
Python, I could send you a script that reads the LiDAR header directly. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">LasPy is a good little library for this <a href="https://laspy.readthedocs.io/en/latest/intro.html#header">https://laspy.readthedocs.io/en/latest/intro.html#header</a> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">You could find tools in WhiteboxTools also. <o:p></o:p></p>
<div>
<p class="MsoNormal">Nicolas Cadieux<o:p></o:p></p>
<div>
<p class="MsoNormal"><a href="https://gitlab.com/njacadieux">https://gitlab.com/njacadieux</a><o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><br>
<br>
<o:p></o:p></p>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<p class="MsoNormal" style="margin-bottom:12.0pt">Le 16 févr. 2023 à 17:50, Mike Stanton via QGIS-User <<a href="mailto:qgis-user@lists.osgeo.org">qgis-user@lists.osgeo.org</a>> a écrit :<o:p></o:p></p>
</blockquote>
</div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<p class="MsoNormal"> <o:p></o:p></p>
<p class="MsoNormal">I have LAZ files for a square mile of the local area where my office is and when it is loaded under QGIS I hit two problems. It is a bit over 2000 files and I seem to have hit a capacity limit where it stops processing the images at about
200 files. The message being output for the files is “Failed to Calculate Statistics of the Point Cloud of XXXXXXXXXXXX.LAZ”. Breaking it into smaller pieces gets the project around the scale limits of the product because the doing only the rejected files
gets me more files processed, but it is not convenient considering I expect to be processing a lot more of these files.<o:p></o:p></p>
<p class="MsoNormal"> <o:p></o:p></p>
<p class="MsoNormal">The other problem is that QGIS tells me I don’t have a CRS for the images and the image has a field that says “No CRS Found” in the log file and what I have is data in the header that says “COMPD_CS [“Projected”, PROJCS “UTM 15N”, GEOGCS
[WGS84 …. “ but the strings have unpaired square brackets round the relevant text. Is this what may be causing the error to occur? <o:p></o:p></p>
<p class="MsoNormal">_______________________________________________<br>
QGIS-User mailing list<br>
<a href="mailto:QGIS-User@lists.osgeo.org">QGIS-User@lists.osgeo.org</a><br>
List info: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a><br>
Unsubscribe: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a><o:p></o:p></p>
</div>
</blockquote>
</div>
</div>
</div>
</body>
</html>