如何通过PDAL Python将LAS文件多维度写入Zarr文件变量？

问题分析 #

1. gdalopts不生效：PDAL Python API中gdalopts参数需要传入键值对字符串，而非字典，原代码的字典形式无法被正确解析。
2. 无法批量写入多维度：GDAL Writer默认单次仅处理一个维度，需为每个目标维度创建独立的Writer阶段。

解决方案 #

1. 修正gdalopts传递格式 #

将原字典形式的gdalopts改为字符串（比如"ARRAY_NAME=Intensity"），确保参数被PDAL正确识别。

2. 批量创建多维度Writer #

为每个需要导出的维度（如Z、Intensity、ReturnNumber）创建单独的GDAL Writer，每个Writer指定对应的dimension和自定义的ARRAY_NAME，最后将所有阶段加入Pipeline执行。

修改后的代码示例 #

%%time

import os
import pdal

# 定义需要导出的维度列表
target_dimensions = ['Intensity', 'Z', 'ReturnNumber']

# Loop through each LAS file and process it
for las_file in las_files:
    base_name = os.path.basename(las_file)
    base_no_ext = os.path.splitext(base_name)[0]
    output_file = f'zarr/{base_no_ext}.zarr'
    
    # 基础处理阶段
    stage_1 = pdal.Reader.las(filename=las_file)
    stage_2 = pdal.Filter.sort(dimension="X")
    
    # 创建多维度Writer阶段
    writer_stages = []
    for dim in target_dimensions:
        # 每个维度对应一个Writer，设置独立的ARRAY_NAME
        writer = pdal.Writer.gdal(
            gdaldriver='zarr',
            output_type='mean',
            resolution=1.0,
            filename=output_file,
            data_type='Float32',
            gdalopts=f'ARRAY_NAME={dim}',
            dimension=dim
        )
        writer_stages.append(writer)
    
    # 构建并执行Pipeline，包含所有阶段
    pipeline = pdal.Pipeline([stage_1, stage_2] + writer_stages)
    pipeline.execute()

说明 #

• 每个Writer阶段对应一个维度，通过ARRAY_NAME指定Zarr中数组的名称，避免覆盖。
• 所有Writer共享同一个Zarr输出路径，PDAL会自动将多个维度数组写入同一个Zarr存储中。
• 可根据需求修改target_dimensions列表，添加或删除需要导出的LAS维度。

内容的提问来源于stack exchange，提问作者Thomas Cobban QG