Python解析Las(点云)格式文件

最新推荐文章于 2024-05-09 22:47:04 发布
最新推荐文章于 2024-05-09 22:47:04 发布
阅读量2.2w 收藏 125
点赞数 34
分类专栏： python
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_22034353/article/details/102506188
版权

Python解析Las(点云)格式文件

tips:

  • 本文代码基于python3编写
  • 代码仓库
  • 推荐使用laspy或者plcpy包，本文做基础研究。

一、点云

  • 点云：

    在逆向工程中通过测量仪器得到的产品外观表面的点数据集合也称之为点云，通常使用三维坐标测量机所得到的点数量比较少，点与点的间距也比较大，叫稀疏点云；而使用三维激光扫描仪或照相式扫描仪得到的点云，点数量比较大并且比较密集，叫密集点云。

  • 点云格式：

    *.las、*.pcd、*.txt

二、Las格式

1. 简介

las文件是一个二进制文件，其中定义的数据类型与C语言中数据类型一致。到目前为止，las共有6版分别是：

2. Las规范数据类型

数据类型字节
char1 byte
unsigned char1 byte
short2 bytes
unsigned short2 bytes
long4 bytes
unsigned long4 bytes
double8 bytes

3. las文件整体构成

las1.0~las1.2las1.3~las1.4
公共头公共头
变长记录域变长记录域
点数据域点数据域
扩展变长记录域

4. 公共头不同版本构成

las1.0~las1.2las1.3~las1.4
文件标识（“LASF”）文件标识（“LASF”）
保留字节文件ID
GUID数据全球编码
版本信息GUID数据
系统标识版本信息
飞行时间系统标识
头部大小文件创建时间
数据偏移头部大小
变长记录域数目数据偏移
点数据格式、长度、数目、不同回波点数目变长记录域数目
X、Y、Z刻度因子点数据格式、长度
X、Y、Z偏移值老格式点数目、不同回波点数目
X、Y、Z最大最小值X、Y、Z刻度因子
X、Y、Z偏移值
X、Y、Z最大最小值
波形数据包记录起始位置
扩展变长记录起始、数目
点数目、不同回波点数目

5. 点记录格式说明

在las1.0版本中定义了点数据格式0，其一共20字节数据，在las1.0~las1.4的版本中点数据格式1到5都是在点数据格式0基础上增添字段。具体字段说明可参见不同版本的官方文档。

在las1.4版本中增加了点格式6，其一共30字节数据，在las1.4版本中点格式7到10都是在点数据格式6基础上增添字段。具体字段说明可参见不同版本的官方文档。

备注：

在点格式0~5中Return Number、Number of Returns (given pulse)、Scan Direction Flag、Edge of Flight Line等字段共占1字节数据。

在点格式6~10中Return Number、Number of Returns、 Classification Flags、 Scan Direction Flag、 Edge of Flight Line等字段共占2字节数据。

这些字段的值都是通过位计算获得。

三、Python解析二进制

Python解析二进制文件使用内置struct包

  1. Python数据类型与C语言数据类型对应关系，参考链接。

    FormatC TypePython字节数
    xpad byteno value1
    ccharbytes of length 11
    bsigned charinteger1
    Bunsigned charinteger1
    ?_Boolbool1
    hshortinteger2
    Hunsigned shortinteger2
    iintinteger4
    Iunsigned intinteger4
    llonginteger4
    Lunsigned longinteger4
    qlong longinteger8
    Qunsigned long longinteger8
    ffloatfloat4
    ddoublefloat8
    schar[]bytes1
    pchar[]bytes1
    Pvoid *integer

  2. 简单使用struct包

    import struct

# 创建float型二进制数据
a = struct.pack('f', 12.34)
print(f"转换的二进制：{a}", )
# 创建char[]类型二进制数据
b = struct.pack('4s', 'test'.encode('utf-8'))
print(f"转换的二进制：{b}", )

# 写入到二进制文件
f = open('test.dat', 'wb+')
f.write(a)
f.write(b)
f.close()

# 读取二进制文件
f = open('test.dat', 'rb')
a = f.read(4)
print(f"读取的二进制：{a}", a)
a = struct.unpack('f', a)
print(f"转换为Python数据类型：{a}", )
b = f.read(4)
print(f"读取的二进制：{b}", )
b = struct.unpack('4s', b)
print(f"转换为Python数据类型：{b}", )
f.close()

四、读取las文件并提取XYZ坐标值

1. 读取公共头

  • 新建head.py，定义las1.0~las1.4版本不同的公共头的类,具体实现参考代码

    class OneZeroHeader(object):
  """
  1.0版本
  """
  def __init__():
      pass
  
  def read_header(f):
      pass

class OneOneHeader(object):
  """
  1.1版本
  """
  def __init__():
      pass
  
  def read_header(f):
      pass

class OneTwoHeader(object):
  """
  1.2版本
  """
  def __init__():
      pass
  
  def read_header(f):
      pass

class OneThreeHeader(object):
  """
  1.3版本
  """
  def __init__():
      pass
  
  def read_header(f):
      pass

class OneFourHeader(object):
  """
  1.4版本
  """
  def __init__():
      pass
  
  def read_header(f):
      pass

def get_header(f, version):
  """
  根据不同版本读取不同头
  """
  if version == (1, 0):
      new_header = OneZeroHeader()
  elif version == (1, 1):
      new_header = OneOneHeader()
  elif version == (1, 2):
      new_header = OneTwoHeader()
  elif version == (1, 3):
      new_header = OneThreeHeader()
  elif version == (1, 4):
      new_header = OneFourHeader()
  else:
      raise Exception("未找到对应文件版本")

  new_header.read_header(f)
  return new_header

  • 新建test_read_las.py，读取公共头数据

    
import struct
from head import get_header


def get_version(f):
    f.read(24)
    version_major, version_minor = struct.unpack('2B', f.read(2))
    print(f"las版本:{version_major}.{version_minor}")
    return version_major, version_minor


if __name__ == '__main__':
    f = open('test.las', 'rb')
    version = get_version(f)
    header = get_header(f, version)
    print(header.__dict__)

2. 读取点数据

  • 新建point.py，定义不同点格式。由本文只需要提取xyz值，所有只需要xyz值，其他数据解析跳过

    """
定义点数据0~10格式解析
"""
import struct


class Point(object):

    def __init__(self, x_s_f, y_s_f, z_s_f, x_o, y_o, z_o):
        self.x_scale_factor = x_s_f
        self.y_scale_factor = y_s_f
        self.z_scale_factor = z_s_f
        self.x_offset = x_o
        self.y_offset = y_o
        self.z_offset = z_o
    
    def get_offset_bytes(self, point_data_record_format):
        """
        根据不同的点格式跳过的字节数
        :param point_data_record_format:点格式0~10
        :return:
        """
        # x,y,z共占12字节
        data_format = {
            0: 8,  # 点格式0 共20字节
            1: 16,  # 点格式1 共28字节
            2: 14,  # 点格式2 共26字节
            3: 22,  # 点格式3 共34字节
            4: 45,  # 点格式4 共57字节
            5: 51,  # 点格式5 共63字节
            6: 18,  # 点格式6 共30字节
            7: 24,  # 点格式7 共36字节
            8: 26,  # 点格式8 共38字节
            9: 47,  # 点格式9 共59字节
            10: 55,  # 点格式10 共67字节
        }
    
        offset_bytes = data_format.get(point_data_record_format, None)
        if offset_bytes is None:
            raise Exception(f"不存在当前的点格式{point_data_record_format}")
        return offset_bytes
    
    def read_point(self, f, offset_to_point_data, point_data_record_format, num):
        """
        读取当前文件中的点数据
        :param f:  数据文件
        :param offset_to_point_data: 点数据开始读取的地方
        :param point_data_record_format:  点数据格式0~10
        :param num: 读取多少个点
        :return: 读取的数据点
        """
        offset_bytes = self.get_offset_bytes(point_data_record_format)
        f.seek(offset_to_point_data)
        points = list()
    
        i = 0
        while i < num:
            point_bytes = f.read(12)
            x_record, y_record, z_record = struct.unpack_from('3l', point_bytes)
            x = x_record * self.x_scale_factor + self.x_offset
            y = y_record * self.y_scale_factor + self.y_offset
            z = z_record * self.z_scale_factor + self.z_offset
            i += 1
            f.read(offset_bytes)
            points.append((x, y, z))
    
        return points

  • test_read_las.py中读取点数据

    import struct
from head import get_header
from point import Point


def get_version(f):
    f.read(24)
    version_major, version_minor = struct.unpack('2B', f.read(2))
    print(f"las版本:{version_major}.{version_minor}")
    return version_major, version_minor


if __name__ == '__main__':
    f = open('test.las', 'rb')
    version = get_version(f)
    header = get_header(f, version)
    print(header.__dict__)

    points = Point(header.x_scale_factor,
                   header.y_scale_factor,
                   header.z_scale_factor,
                   header.x_offset,
                   header.y_offset,
                   header.z_offset,
                   )

    data = points.read_point(f, header.offset_to_point_data,
                             header.point_data_record_format,
                             header.number_of_point_records)

    print(data)

3. 输出XYZ坐标值到文件

在test_read_las.py将解析出来的点数据按照一定的格式（x,y,z）输出到txt文件。

import struct
from head import get_header
from point import Point


def get_version(f):
    f.read(24)
    version_major, version_minor = struct.unpack('2B', f.read(2))
    print(f"las版本:{version_major}.{version_minor}")
    return version_major, version_minor


if __name__ == '__main__':
    f = open('test.las', 'rb')
    version = get_version(f)
    header = get_header(f, version)
    print(header.__dict__)

    points = Point(header.x_scale_factor,
                   header.y_scale_factor,
                   header.z_scale_factor,
                   header.x_offset,
                   header.y_offset,
                   header.z_offset,
                   )

    data = points.read_point(f, header.offset_to_point_data,
                             header.point_data_record_format,
                             header.number_of_point_records)

    f.close()

    with open("point.txt", 'a+') as f:
        for item in data:
            f.write(f'{item[0]},{item[1]},{item[2]}\n')

总结

  • 本文完成了对las文件的读取并输出到文件
  • 下一篇将讲解多进程分块加速读取las文件并输出到对应文件
那些年踩过的坑丶
关注
  • 34
    点赞
  • 踩
  • 125
    收藏
    觉得还不错? 一键收藏
  • 32
    评论
专栏目录
Python点云处理（一）点云数据读取与写入
weixin_43894929的博客
07-23 5873
利用python语言，实现不同格式点云数据的读取与写入
Python Las与Pcd文件的读取与保存、格式转换
weixin_43896283的博客
10-26 556
Las与Pcd文件的读取与保存、格式转换
点云las格式转换为pcd
02-13
点云PCL库支持的pcd格式包含点云强度等信息，但一般扫描仪得到的多为las格式，这里提高一个转换工具，亲测速度很快
用 Python 读取 .las 文件并修改 RGB 信息
lalala6983的博客
07-15 566
.las 文件的读写及 RGB 信息的修改
使用Python读取las点云，写入las点云，无损坐标精度
最新发布
qq_27816785的博客
05-09 750
在 LAS 文件中，您实际上将存储一个整数：123456789123，该整数将在读取时使用比例因子转换为 123456789.123。因此存储时，将该值偏移 123450000，实际存的是6789123。las 格式原生支持的属性字段，与las版本密切相关。* las 存储数据时，需要设置 scales 和 offsets，否则会出现精度问题。读写点云文件时，着实被坐标精度问题难到了，顺便就仔细学习了下las格式到底是什么来头。* 实测版本 laspy v2.3，2.x版本应该都可用，但不适合1.x版本。
点云格式解读LAS/PCD
热门推荐
lhm_19960601的博客
07-27 1万+
1.LAS点云格式 从本质上来说,LAS格式是一种二进制文件格式。其目的是提供一种开放的格式标准,允许不同的硬件和软件提供商输出可互操作的统一格式。现在LAS格式文件已成为LiDAR数据的工业标准格式   LAS文件按每条扫描线排列方式存放数据,包括激光点的三维坐标、多次回波信息、强度信息、扫描角度、分类信息、飞行航带信息、飞行姿态信息、项目信息、GPS信息、数据点颜色信息等。LAS格式定义...
LAS格式解析
天凉好丶秋的博客
04-18 3507
LAS格式解析1. LAS的起源2. LAS的组成部分3. 各部分所包含的内容3.1 公共头新的改变功能快捷键合理的创建标题，有助于目录的生成如何改变文本的样式插入链接与图片如何插入一段漂亮的代码片生成一个适合你的列表创建一个表格设定内容居中、居左、居右SmartyPants创建一个自定义列表如何创建一个注脚注释也是必不可少的KaTeX数学公式新的甘特图功能，丰富你的文章UML 图表FLowchart流程图导出与导入导出导入 1. LAS的起源   为了方便各大硬件、软件厂商和用户之间进行数据的交互，ASP
利用laspy读取写入las文件
xinxiangwangzhi_的博客
03-03 5240
首先要了解las数据格式的基本构成，其次利用开源库计算，laspy提供了两种方式 读取las文件 import laspy import numpy as np #参考链接： # 用Matlab处理LAS点云（1）——LAS文件概述--知乎 # https://pythonhosted.org/laspy/tut_part_1.html # https://laspy.readthedocs.io/en/latest/examples.html # 按间距中的绿色按钮以运行脚本。 """ Laspy can
【Lidar】Laspy库介绍+基础函数（读取、可视化、保存、旋转、筛选、创建点云数据）
Remote Sensing迷途小书童
12-11 2157
laspy是一个Python库，主要用于读取、修改和创建LAS点云文件。该库兼容Python 2.6+和3.5+，并且可以处理LAS版本1.0-1.3的文件。
LAStool读取和保存.las格式点云
03-25
VS2013编译LAStool，读取和保存las格式数据，使用pcl点云库进行显示和数据处理
.las格式点云文件
12-02
las格式点云文件，测试数据。我就是想上传一个测试数据，实在没啥描述的，非让我写够50个字，这不是扯淡嘛。
.las文件读写实验数据
06-12
用于.las文件读写的实验数据.las格式的地面激光雷达数据.las格式的激光雷达数据，读取.las文件提供源数据。
las格式转pcd和txt
05-23
激光点云las格式转pcd(point cloud libray格式)和txt格式
livox-SDK lvx格式转pcd、las、txt
06-07
Livox SDK提供了C++和Python两种语言的API，可以方便地读取LVX文件并按需导出到其他格式。转换过程一般包括以下步骤： 1. **初始化Livox SDK**：在代码中设置必要的配置，如设备连接参数，然后启动SDK。 2. **加载...
点云转换成3dtiles源码
12-11
在处理点云数据时，我们需要理解如何解析不同格式的点云文件，如LAS或XYZ等。 其次，3DTiles是Cesium开源JavaScript库推出的一种高效、流式加载的三维地理空间数据格式。它将大规模三维模型分割成小块（tiles），每...
libLAS-1.8.1.tar.gz
10-19
《Python环境下libLAS库解析LAS点云数据详解》 在当今的三维空间数据处理领域，点云数据的应用越来越广泛，特别是在地理信息系统（GIS）、遥感、建筑建模、自动驾驶等领域。其中，.las（LAS）文件格式作为一种常用...
点云体素下采样(抽稀)（附open3d python代码）
三维点云技术探索
11-29 1561
参数voxel_size为体素大小，该参数越小，降采样得到的结果的点越多，采样后的点排列看上去会很整齐均匀，但是位置和之前相比有轻微的移动。2. 每个被占用的体素通过平均内部所有点的坐标而生成一个点。体素下采样使用常规体素网格从输入点云创建均匀下采样的点云。它通常用作多点云处理任务的预处理步骤。因为下采样之后的点是通过平均计算得到的。该方法的缺陷是会改变原有点的坐标，1. 点被储存到体素中。
基于python将激光点云.las文件投影为深度图/强度图/反强度图
moringacui的博客
03-19 7038
基于python实现激光点云投影为深度图/强度图/反强度图，提供了代码与数据集供测试。
python读取las点云数据
05-02
要读取LAS点云数据，可以使用Python中的Laspy库。 首先需要安装Laspy库。可以使用以下命令进行安装： ``` pip install laspy ``` 然后可以使用以下代码读取LAS文件： ```python import laspy # 打开LAS文件 inFile = laspy.file.File("example.las", mode="r") # 读取点云数据 points = inFile.points # 输出点云数据 for point in points: print(point) ``` 在上面的代码中，我们使用`laspy.file.File()`方法打开LAS文件，使用`mode="r"`参数以只读模式打开文件。然后我们可以使用`inFile.points`属性读取点云数据。最后，我们可以使用循环来输出点云数据。 除了读取点云数据外，Laspy库还提供了许多其他的功能，例如读取点云属性、写入点云数据等。这些功能可以通过查看Laspy文档来学习。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 32
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值