Lasclip crash at polygon 507

[bobbysmob...@gmail.com]

Jochen,

I need to clip out classes 4 and 5 that reside within a polygon SHP.  I have a SHP containing ~120K polygons to work through.

My command is:

C:\LAStools\bin\lasclip64 -v -i *.laz -merged -poly G:\DATA_FOR_PROCESSING\C6_AOI.shp -split RAND -keep_class 4 5 -odir %DELIVERY_DIRECTORY%\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG -olaz

Version of LAStools was downloaded on January 24, 2023.


The clip operation works fine until clip polygon 507 is encountered, then the script crashes, see below.


.
.
.
creating clip polygon 501 took 0.001 sec.
** create output file [X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\7482020.laz]
removed points 56970911667 (445). outputting 57 points took 0.177 sec.
creating clip polygon 502 took 0.001 sec.
** create output file [X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\139277.laz]
removed points 56970911573 (273). outputting 151 points took 0.17 sec.
creating clip polygon 503 took 0.001 sec.
** create output file [X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\4527452.laz]
removed points 56970911724 (41). outputting 0 points took 0.388 sec.
creating clip polygon 504 took 0.002 sec.
** create output file [X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\4963274.laz]
removed points 56970911690 (138). outputting 34 points took 0.163 sec.
creating clip polygon 505 took 0.001 sec.
** create output file [X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\8412267.laz]
removed points 56970911722 (219). outputting 2 points took 0.176 sec.
creating clip polygon 506 took 0.002 sec.
** create output file [X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\9059566.laz]
ERROR: cannot open file '755255.laz' for read
ERROR: cannot open file '755255.laz' for read
ERROR: could not open lasreaderlas for file '755255.laz'
removed points 56970911724 (0). outputting 0 points took 0.062 sec.
creating clip polygon 507 took 0.002 sec.
** create output file [X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\8671774.laz]
ERROR: cannot open file 'X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\8671774.laz' for write
ERROR: cannot open file 'X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\8671774.laz' for write
ERROR: cannot open laswriterlas with file name 'X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\8671774.laz'
ERROR: could not open laswriter


The issue is not tile 755255.laz.  I have rearranged the order of the polygons processed by manipulating the order of the polygons encountered via the SHP attribute (i.e. RAND).

The script fails once it reaches clip polygon 507 every time, regardless of the tile being processed.

Interested in our thoughts regarding this issue.

Thank you,
Bobby 

[Jochen Rapidlasso]

Hi Bobby,

just as a fast reply: I remember, SHP files have problems if getting too large - I need to do some research on that.

Can you share your data so I can test with the original data?

Best regards,

Jochen @rapidlasso

[bobbysmob...@gmail.com]

Jochen,

I can send you the SHP via WeTransfer as it is large.  Shall I send it to sup...@rapidlasso.de?

In the mean time I will create a smaller version of SHP and see if it works.

Thanks,

Bobby

[Jochen Rapidlasso]

Hi Bobby,

yes, please send the link to support ät rapidlasso.de 

Thanks,

Jochen @rapidlasso

[bobbysmob...@gmail.com]

Jochen,

I ran the task again after excluding 50% of the original SHP polys and it crashed at polygon 507.

creating clip polygon 505 took 0.001 sec.

** create output file [X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\2379729.laz]
removed points 56970911724 (0). outputting 0 points took 0.058 sec.
creating clip polygon 506 took 0.001 sec.
** create output file [X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\196279.laz]
ERROR: cannot open file U670795.laz' for read
ERROR: cannot open file '670795.laz' for read
ERROR: could not open lasreaderlas for file '670795.laz'
removed points 56970911724 (0). outputting 0 points took 0.056 sec.
creating clip polygon 507 took 0.001 sec.
** create output file [X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\1262045.laz]
ERROR: cannot open file 'X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\1262045.laz' for write
ERROR: cannot open file 'X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\1262045.laz' for write
ERROR: cannot open laswriterlas with file name 'X:\Projects_2022\AREA_GG\20230310_DEL1\STEP-6_SPECIFIC_CLASSES\VEG_INSIDE_BLDG\1262045.laz'

ERROR: could not open laswriter

Thanks,

Bobby

[bobbysmob...@gmail.com]

RESOLVED, sort of.

On a lark I referenced a 2021 version of lasclip and it works.  So there is something not right in the lasclip code subsequent to version 210330.

C:\LAStools_2021\bin>lasclip -version

WARNING: your license has expired
Please note that LAStools is not "free" (see http://lastools.org/LICENSE.txt)
contact 'martin....@rapidlasso.com' to clarify licensing terms if needed.
LAStools (by mar...@rapidlasso.com) version 210330 (unlicensed)

C:\LAStools_2021\bin>

Thanks,

Bobby

[Jochen Rapidlasso]

Hi Bobby,

this issue is fixed now with LAStools version 230330.

We support now > 500 shapes even in "-split" mode.

Best regards,

Jochen @rapidlasso

[Bobby Riley]

Wonderful!  Thank you!

Bobby

--
Download LAStools at
https://rapidlasso.de
Manage your settings at
https://groups.google.com/g/lastools/membership
---
You received this message because you are subscribed to the Google Groups "LAStools - efficient tools for LiDAR processing" group.
To unsubscribe from this group and stop receiving emails from it, send an email to lastools+u...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/lastools/ce4d847a-1b5a-42d6-871e-49e76e9ca51bn%40googlegroups.com.

--

Sent from my iPhone. Please excuse brevity and typos.