如何用Python读取仅含单强度值的PLY点云文件

方法一：直接读取并访问强度属性 #

如果你的PLY文件头部正确定义了intensity字段（比如ASCII格式的PLY头部包含property float intensity），Open3D其实已经读取了该数据，只是默认可视化不会显示。你可以通过以下方式访问并可视化：

import open3d as o3d
import numpy as np

path = './example.ply'
pcd = o3d.io.read_point_cloud(path)

# 检查并读取强度属性
if 'intensity' in pcd.point:
    intensities = pcd.point.intensity
    print(f"成功读取{len(intensities)}个强度值")
    
    # 将强度值映射为灰度颜色用于可视化
    norm_intensities = (intensities - np.min(intensities)) / (np.max(intensities) - np.min(intensities))
    pcd.colors = o3d.utility.Vector3dVector(np.stack([norm_intensities, norm_intensities, norm_intensities], axis=1))

o3d.visualization.draw_geometries([pcd])

方法二：手动解析PLY文件（适用于非标准格式） #

如果PLY文件的强度字段定义不规范，Open3D无法自动识别，可以手动解析文件内容：

import open3d as o3d
import numpy as np

path = './example.ply'

# 读取并解析PLY文件
with open(path, 'r') as f:
    lines = f.readlines()

# 定位数据起始行
data_start_idx = lines.index('end_header\n') + 1

# 提取点坐标和强度值
points = []
intensities = []
for line in lines[data_start_idx:]:
    x, y, z, intensity = map(float, line.strip().split())
    points.append([x, y, z])
    intensities.append(intensity)

# 创建Open3D点云对象并赋值
pcd = o3d.geometry.PointCloud()
pcd.points = o3d.utility.Vector3dVector(np.array(points))
pcd.point.intensity = o3d.utility.VectorXdVector(np.array(intensities))

# 映射强度为可视化颜色
norm_intensities = (np.array(intensities) - np.min(intensities)) / (np.max(intensities) - np.min(intensities))
pcd.colors = o3d.utility.Vector3dVector(np.stack([norm_intensities]*3, axis=1))

o3d.visualization.draw_geometries([pcd])

注意事项 #

• 确保PLY文件头部包含property float intensity字段（ASCII格式），或二进制格式的字段定义符合规范，Open3D才能自动识别强度属性。
• 强度值范围较大时，可视化前需归一化到0-1区间，否则颜色显示会异常。

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