Julia wrapper for LibLAS, a library for reading and writing the LAS LiDAR format.
Note: Like LibLAS itself, this project is currently in hibernation or maintenance mode. For reading and writing LAS files in pure Julia, have a look at LasIO.jl.
pkg> add https://github.com/visr/LibLAS.jlYou need to have the liblas_c shared library in your path, this is currently not added to the package. See the LibLAS website for how to obtain it.
using LibLAS
# testfile downloaded from http:https://www.liblas.org/samples/
filename = joinpath(dirname(pathof(LibLAS)), "..", "test", "libLAS_1.2.las")
reader = LibLAS.create(LibLAS.LASReader, filename)
header = LibLAS.lasheader(reader)
n = LibLAS.pointrecordscount(header)
bounds = [LibLAS.min_x(header), LibLAS.min_y(header), LibLAS.min_z(header),
LibLAS.max_x(header), LibLAS.max_y(header), LibLAS.max_z(header)]
function mean_coordinates(header, reader)
n = LibLAS.pointrecordscount(header)
x_sum = 0.0
y_sum = 0.0
z_sum = 0.0
# looping over all points
for i = 1:n
p = LibLAS.nextpoint(reader)
x = LibLAS.xcoord(p)
y = LibLAS.ycoord(p)
z = LibLAS.zcoord(p)
x_sum += x
y_sum += y
z_sum += z
end
x_avg = x_sum / n
y_avg = y_sum / n
z_avg = z_sum / n
x_avg, y_avg, z_avg
end
x_avg, y_avg, z_avg = mean_coordinates(header, reader)
# jumping to a specific point
p = LibLAS.point_at(reader, 101)
intens = LibLAS.intensity(p)For more documentation on specific functions see the LibLAS C API reference. Function names are adapted to follow Julian naming conventions. They have not all been tested. Please open an issue/PR if you find any mistakes.