Convert PointCloud to Geotiff raster file

Question

Here is my problem, I have point cloud data in .las format. I need to find a way to create a Geotiff file from this cloud. The pixel size of the tiff file is 2cm. The cloud has been resampled in the same way.

I was thinking of using the GDAL library to create the tiff file but saw that it needed ordered data, so my las file is scattered.

The problem is how to convert the scattered data of the las into a 2D array containing the elevations (Z) in an ordered way. All without needing several for loops that slow down the process too much.

I started from the solution found in this topic: writing array to raster

The numpy array looks like :

array([[6.48913764e+05, 6.83548129e+06, 8.25270000e+01],
       [6.48914073e+05, 6.83548840e+06, 8.20520000e+01],
       [6.48917150e+05, 6.83549002e+06, 8.22260000e+01],
       ...,
       [6.49070029e+05, 6.83557728e+06, 8.40550000e+01],
       [6.49070563e+05, 6.83557851e+06, 8.44870000e+01],
       [6.49070563e+05, 6.83557851e+06, 8.44870000e+01]])

Where the first colomn is X coordinates, the second is Y coordinates and the last for Z. I can't transfer files, but a cloud contains about 500k observations.

Thanks for your answer

Answer 1

you could consider using opalsDTM developed by TU Vienna: https://opals.geo.tuwien.ac.at/html/stable/ModuleDTM.html