Merge las files using laspy (basic example)

merge_las_files.py

#!/usr/bin/env python3

import os
import laspy
import shutil

file_ext = ".las"
output_file = "output.las"

if os.path.exists(output_file):
    os.remove(output_file)
    print(output_file + " was removed")

# Merges las files with identical headers, VLRs and point formats

for f in os.listdir('.'):
    if f.endswith(file_ext) and f != output_file:
        with laspy.open(f) as las:
            if not os.path.exists(output_file):
                shutil.copy(f, output_file)
                print(output_file + " was created as a copy of " + f)
                continue
            with laspy.open(output_file, mode="a") as writer:
                for points in las.chunk_iterator(2_000_000):
                    writer.append_points(points)
            print(f + " points have been added to " + output_file)