Think about whether there are cheap and non-intrusive ways to guard against overflow. about laspy HOT 10 CLOSED

import os, struct, sys, glob
import laspy

if __name__ == "__main__":
    wkt = 'PROJCS["NAD_1983_UTM_Zone_15N",GEOGCS["GCS_NAD_1983_2011",DATUM["NAD_1983_2011",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433],AUTHORITY["EPSG","4269"]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-93.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0],VERTCS["NAVD_1988",VDATUM["North_American_Vertical_Datum_1988"],PARAMETER["Vertical_Shift",0.0],PARAMETER["Direction",1.0],UNIT["Meter",1.0]],AUTHORITY["EPSG","26915"]]'

    # Read las files
    for las_file in sorted(glob.glob(sys.argv[1])):
        print("Attempting %s..." % os.path.basename(las_file))

        f = laspy.file.File(las_file, mode='rw')

        # WARNING: Removing all VLRs
        inVLRs = []

        # OGC Coordinate System WKT
        new_vlr = laspy.header.VLR(user_id = "LASF_Projection",
                                   record_id = 2112,
                                   VLR_body = wkt,
                                   description = "OGC Coordinate System WKT")

        inVLRs.append(new_vlr)
        f.header.vlrs = inVLRs

        # Set WKT Global Encoding Bit
        f.header.wkt = 1

        f.close()

Traceback (most recent call last):
File "D:\path\to\above\code.py" line...
f.header.vlrs = inVLRs
File "C:\path\to\laspy\header.py", line 1275, in set_vlrs
    self.reader.set_vlrs(value)
File "C:\path\to\laspy\base.py", line 948, in set_vlrs
    dat_part_2 = self.data_provider.fileref.read(current_size - old_offset)
OverflowError: Python int too large to convert to C long

for las_file in sorted(glob.glob(sys.argv[1])):
    print("Attempting %s..." % os.path.basename(las_file))

    f = laspy.file.File(las_file, mode='r')
    h = f.header
    inVLRs = []

    # OGC Coordinate System WKT
    new_vlr = laspy.header.VLR(user_id = "LASF_Projection",
                               record_id = 2112,
                               VLR_body = wkt,
                               description = "OGC Coordinate System WKT")

    # Initialize output file for writing data
    out_file = os.path.join(os.path.dirname(las_file), os.path.splitext(las_file)[0] + "_HEADER.las")
    out_f = laspy.file.File(out_file, mode='w', header = h)

    inVLRs.append(new_vlr)
    out_f.header.vlrs = inVLRs

    # Set WKT Global Encoding Bit
    out_f.header.wkt = 1

    # Set GPS Global Encoding Bit
    out_f.header.gps_time_type = 1

    out_f.points = f.points
    f.close()
    out_f.close()

out_f.points = f.points