如何使用laspy在Python中正确读取LiDAR .las文件的GPS_time

这种问题本质是LAS文件中GPS时间的存储基准与你预期的Unix时间戳不一致，你的数据集里的GPS_time实际是截断了前导1（即减去了10^9）的Unix时间戳，以下是更高效的原生处理方案：

1. 高效修正：使用numpy向量化运算替代列表推导 #

laspy返回的gps_time是numpy数组，直接用向量运算修正比列表推导效率高得多（尤其适合大体积LAS数据）：

import laspy
import numpy as np

path = r"\\xxxx\las_data.las"
las_data = laspy.read(path)

# 直接对numpy数组做批量运算，无需循环
correct_gps_time = las_data.gps_time + 1000000000
# 转为整数类型（按需）
correct_gps_time = correct_gps_time.astype(np.int64)

2. 严谨溯源：从LAS头文件确认时间基准 #

如果想彻底避免手动修正的不确定性，可以先查看LAS文件的头信息，确认GPS时间的存储规则：

# 查看GPS时间的存储类型
print(las_data.header.gps_time_type)
# 查看GPS周相关参数（若有）
print(f"GPS周: {las_data.header.gps_week}, 周内偏移: {las_data.header.gps_week_offset}")

如果头文件中包含标准GPS周信息，可以通过GPS基准时间（1980年1月6日对应的Unix时间戳315964800）计算完整时间：

gps_epoch = 315964800
seconds_per_week = 604800
full_gps_time = gps_epoch + las_data.header.gps_week * seconds_per_week + las_data.gps_time

但如果你的数据集确实是人为截断了Unix时间戳的前导1，第一种向量化修正的方法最直接高效。

内容的提问来源于stack exchange，提问作者Valerii Tarasova