<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<p>To report to <a class="moz-txt-link-freetext" href="https://github.com/qgis/QGIS/issues">https://github.com/qgis/QGIS/issues</a><br>
</p>
<div class="moz-cite-prefix">Le 20/02/2023 à 18:48, Mike Stanton via
QGIS-User a écrit :<br>
</div>
<blockquote type="cite"
cite="mid:8f6005b759354111a61dccdcea05129a@geospan.com">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<meta name="Generator" content="Microsoft Word 15 (filtered
medium)">
<style>@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;}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;}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]-->
<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
<a class="moz-txt-link-rfc2396E" href="mailto:njacadieux.gitlab@gmail.com"><njacadieux.gitlab@gmail.com></a> <br>
<b>Sent:</b> Friday, February 17, 2023 7:20 AM<br>
<b>To:</b> Mike Stanton <a class="moz-txt-link-rfc2396E" href="mailto:mfstanton@geospan.com"><mfstanton@geospan.com></a><br>
<b>Cc:</b> <a class="moz-txt-link-abbreviated" href="mailto:qgis-user@lists.osgeo.org">qgis-user@lists.osgeo.org</a><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"
moz-do-not-send="true" class="moz-txt-link-freetext">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"
moz-do-not-send="true" class="moz-txt-link-freetext">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"
moz-do-not-send="true" class="moz-txt-link-freetext">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"
moz-do-not-send="true" class="moz-txt-link-freetext">QGIS-User@lists.osgeo.org</a><br>
List info: <a
href="https://lists.osgeo.org/mailman/listinfo/qgis-user"
moz-do-not-send="true" class="moz-txt-link-freetext">https://lists.osgeo.org/mailman/listinfo/qgis-user</a><br>
Unsubscribe: <a
href="https://lists.osgeo.org/mailman/listinfo/qgis-user"
moz-do-not-send="true" class="moz-txt-link-freetext">https://lists.osgeo.org/mailman/listinfo/qgis-user</a><o:p></o:p></p>
</div>
</blockquote>
</div>
</div>
</div>
<br>
<fieldset class="moz-mime-attachment-header"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
QGIS-User mailing list
<a class="moz-txt-link-abbreviated" href="mailto:QGIS-User@lists.osgeo.org">QGIS-User@lists.osgeo.org</a>
List info: <a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a>
Unsubscribe: <a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a>
</pre>
</blockquote>
<pre class="moz-signature" cols="72">--
<a class="moz-txt-link-freetext" href="http://www.spatialys.com">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
</body>
</html>