C++中从LibLAS迁移至PDAL时，从头编写LAS文件的文件头配置问题

我之前从LibLAS迁移到PDAL的时候也卡在了LAS文件头配置这一步，PDAL的思路和LibLAS不太一样，但其实用LasHeader和LasWriter就能轻松实现自定义头的从头写入，给你详细说下：

核心思路 #

PDAL中没有直接的SetHeader方法，而是通过创建自定义LasHeader对象，再将其作为选项传递给LasWriter来实现的。整个流程分为：构建头信息 → 关联点表结构 → 写入点数据。

完整代码示例 #

#include <pdal/LasHeader.hpp>
#include <pdal/LasWriter.hpp>
#include <pdal/PointTable.hpp>
#include <pdal/PointView.hpp>
#include <pdal/Options.hpp>

int main() {
    // 1. 构建自定义LAS文件头
    pdal::LasHeader lasHeader;
    
    // 设置LAS版本（这里以1.2为例）
    lasHeader.setVersion(pdal::LasHeader::Version::e1_2);
    
    // 设置点格式（LAS 1.2对应格式0/1/2，这里选格式0）
    lasHeader.setPointFormatId(pdal::LasHeader::PointFormat::e0);
    
    // 可选：设置坐标范围（提前设置能优化写入性能，也可以后续自动计算）
    lasHeader.setMinX(0.0);
    lasHeader.setMaxX(1000.0);
    lasHeader.setMinY(0.0);
    lasHeader.setMaxY(1000.0);
    lasHeader.setMinZ(0.0);
    lasHeader.setMaxZ(50.0);
    
    // 可选：预估点数量（帮助PDAL预分配资源）
    lasHeader.setPointCount(10000);
    
    // 可选：设置坐标参考系统（比如EPSG:4326）
    lasHeader.setSpatialRef("EPSG:4326");

    // 2. 创建PointTable并关联头的点格式
    pdal::PointTable table;
    lasHeader.setup(table); // 这一步会根据头的点格式自动添加对应的字段
    
    // 创建PointView来存储点数据
    pdal::PointViewPtr view = pdal::PointView::create(table);
    
    // 3. 向PointView添加点数据（示例：添加一个点）
    pdal::PointRef point(*view, 0);
    point.setField(pdal::Dimension::Id::X, 100.0);
    point.setField(pdal::Dimension::Id::Y, 200.0);
    point.setField(pdal::Dimension::Id::Z, 30.0);
    point.setField(pdal::Dimension::Id::Intensity, 255);
    // 如果是格式1/2，还可以添加GPS时间、RGB等字段
    // point.setField(pdal::Dimension::Id::GpsTime, 123456.789);

    // 4. 配置LasWriter并写入文件
    pdal::Options writerOptions;
    writerOptions.add("filename", "output.las");
    writerOptions.add("header", lasHeader); // 关键：将自定义头传入Writer
    
    pdal::LasWriter writer;
    writer.setOptions(writerOptions);
    writer.prepare(table);
    writer.write(view);
    
    return 0;
}

关键细节说明 #

• 头与点表的关联：lasHeader.setup(table)是必须的，它会根据你设置的点格式，自动在PointTable中添加对应的维度字段（比如X/Y/Z/Intensity等），确保点数据和头格式匹配。
• 版本与点格式兼容性：不同LAS版本支持的点格式不同，比如LAS 1.2支持格式0-2，LAS 1.4支持格式0-6、10等，设置时要注意对应关系。
• 其他头信息配置：如果需要设置文件源ID、项目ID、生成软件等信息，可以调用lasHeader.setFileSourceId()、lasHeader.setProjectId()、lasHeader.setSoftwareId("YourAppName")等方法。
• 自动计算范围：如果预先不知道坐标范围，可以跳过setMinX/MaxX等设置，PDAL会在写入点数据时自动计算并更新头信息。

内容的提问来源于stack exchange，提问作者Mathieu Poulin