comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: Open Source LiDAR LAS file manipulator for IDL and ENVI
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: Open Source LiDAR LAS file manipulator for IDL and ENVI [message #92787] Sun, 28 February 2016 16:07
Dick Jackson is currently offline  Dick Jackson
Messages: 347
Registered: August 1998
Senior Member
On Thursday, 25 September 2014 01:45:51 UTC-7, Antoine Cottin wrote:
> Dick,
>
> Glad you enjoyed it.
>
> Thank you for the couple of lines, very nice and useful. I'll integrate them to the next release.
>
> We actually had in the past something like that implemented already. It works very well with "small" LAS file. But when you load a full data set, several Go of data, then the loading time is just a pain and softwares like LASTools/lasview (Open Source), CloudCompare (Open Source) or QT Viewer are much more efficient for quick visualisation.
>
> Cheers.
> Antoine

Hi again,

I learned that the Fleurdelas library (and the required Sazerac library) were updated, and some of the calling sequences have changed to be more convenient. Here's an update to my earlier instructions (only changed in the first four lines, which became three):

[...] for those who are curious about what's actually *in* that MiltonKeynes.las file (or any other LAS file they may come across), I offer these eight lines of code:

root = File_dirname(File_dirname(Routine_filepath('fleurdelas__defi ne')))
inputFile = FILEPATH('MiltonKeynes.las', ROOT_DIR=root, SUBDIRECTORY=['data'])
lasobj = obj_new('fleurdelas', INPUTFILE=inputFile)
pts = lasobj.getData(/ALL)
xyz = Transpose([[pts.east],[pts.north],[pts.elev]])
rgb = Transpose([[pts.r],[pts.g],[pts.b]])/256U
o = IDLgrPolygon(xyz, Style=0, Vert_Colors=rgb)
XObjView, o

(remember to click the Reset button at left after resizing the window)

BTW, for thicker points, do something like this:
o.SetProperty,thick=3
... then click the display window again to refresh it.

For those who want to try this, the simple way is:
1) use Git to clone (IDL Workbench: "Import...") these two repositories ("projects"):
https://github.com/carbomap/sazerac.git
https://github.com/carbomap/fleurdelas.git
2) build both projects (compile all files... sometimes requires compiling a few times, without resetting... I think there's a minor compiling bug in IDL)
3) paste those eight lines into the IDL console

Thanks yet again, Antoine!
-Dick

Dick Jackson Software Consulting Inc.
Victoria, BC, Canada
www.d-jackson.com
[Message index]
 
Read Message
Previous Topic: problem with idl program
Next Topic: FILE_LINES change in 8.4

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Wed Oct 08 14:21:25 PDT 2025

Total time taken to generate the page: 0.00228 seconds