PDAL管线导入LAS点云至PostgreSQL/PostGIS时RGB等数据缺失问题

我来帮你搞定这个问题，咱们先理清楚核心原因，再一步步解决：

问题根源 #

你现在遇到两个关键问题：

1. 用ST_AsEWKT(PC_EXPLODE(pa)::geometry)查询时，::geometry类型转换只会提取点的空间坐标，完全不会保留点云的其他属性——这是geometry类型的特性，它只存位置信息。
2. 虽然你加了output_dims，但可能没包含所有需要的属性，而且查询方式不对，导致看不到这些属性。

分步解决方案 #

1. 补全PDAL管线的output_dims参数 #

从你用pdal info查到的结果来看，LAS文件里有很多属性（Red、Green、Blue、Intensity、Classification等等），你需要把这些全部加到output_dims里，而且要注意字段名的大小写（LAS的字段名是大小写敏感的，比如要写Red而不是RED）：

{
  "pipeline":[
    {
      "type":"readers.las",
      "filename":"F:\\Drive\\KAP\\KAP\\99.526\\pointclouds\\201406071257\\201406071257.las",
      "spatialreference":"EPSG:32635"
    },
    {
      "type":"filters.chipper",
      "capacity":600
    },
    {
      "type":"writers.pgpointcloud",
      "connection":"host='localhost' dbname='dbname' user='user' password='pw'",
      "schema":"pointcloud_99_526",
      "table":"pc_201406071257",
      "compression":"dimensional",
      "srid":"32635",
      "output_dims":"X,Y,Z,Red,Green,Blue,Intensity,Classification,ReturnNumber,NumberOfReturns,ScanDirectionFlag,EdgeOfFlightLine,ScanAngleRank,UserData,PointSourceId"
    }
  ]
}

建议你直接把pdal info返回的所有属性字段都加进去，避免遗漏。

2. 用正确的SQL查询所有属性 #

别再把PC_EXPLODE的结果转成geometry了，改用PC_Get函数来提取每个属性值，这样就能拿到所有信息：

基础查询方式

SELECT
  ST_AsEWKT(PC_Explode(pa)) AS point_geometry,
  PC_Get(PC_Explode(pa), 'Red') AS red_channel,
  PC_Get(PC_Explode(pa), 'Green') AS green_channel,
  PC_Get(PC_Explode(pa), 'Blue') AS blue_channel,
  PC_Get(PC_Explode(pa), 'Intensity') AS intensity,
  PC_Get(PC_Explode(pa), 'Classification') AS classification
FROM pointcloud_99_526.pc_201406071257;

更高效的查询方式（避免重复调用PC_Explode）

WITH exploded_points AS (
  SELECT PC_Explode(pa) AS single_pc_point
  FROM pointcloud_99_526.pc_201406071257
)
SELECT
  ST_AsEWKT(single_pc_point) AS point_geometry,
  PC_Get(single_pc_point, 'Red') AS red_channel,
  PC_Get(single_pc_point, 'Green') AS green_channel,
  PC_Get(single_pc_point, 'Blue') AS blue_channel,
  PC_Get(single_pc_point, 'Intensity') AS intensity,
  PC_Get(single_pc_point, 'Classification') AS classification
FROM exploded_points;

3. 验证数据库中的点云维度 #

先确认所有属性都已经成功写入数据库，你可以用这个查询查看点云的维度信息：

SELECT PC_Summary(pa) FROM pointcloud_99_526.pc_201406071257 LIMIT 1;

这个结果会列出点云包含的所有维度，检查是否有你需要的Red、Green等字段。

额外小提示 #

• 如果你用的是PostGIS 3.0及以上版本，可以用更简洁的ST_PointCloudExtract函数来提取属性：

  SELECT
    ST_AsEWKT(ST_PointCloudExtract(pa, 'geometry')) AS point_geometry,
    ST_PointCloudExtract(pa, 'Red') AS red_channel,
    ST_PointCloudExtract(pa, 'Green') AS green_channel
  FROM pointcloud_99_526.pc_201406071257;
  

• 注意PDAL和PGPointCloud的版本兼容性，你用的PDAL 1.8.0建议搭配PGPointCloud 1.2.0及以上版本，避免出现奇怪的兼容性问题。

内容的提问来源于stack exchange，提问作者Spatial Digger