PDAL库读写多种格式的点云

最新推荐文章于 2024-03-14 11:35:38 发布
最新推荐文章于 2024-03-14 11:35:38 发布
阅读量2.7k 收藏 16
点赞数 6
分类专栏： 点云读写 C++ PDAL 文章标签： PDAL 点云读写
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/jacken123456/article/details/109780699
版权
C++ 同时被 3 个专栏收录
34 篇文章 0 订阅
订阅专栏
点云读写
7 篇文章 2 订阅
订阅专栏
PDAL
2 篇文章 1 订阅
订阅专栏

 

// PDALTest.cpp : 定义控制台应用程序的入口点。
//
#pragma execution_character_set("utf-8")  //解决中文

#include "stdafx.h"

#include <iostream>

#include <memory>
#include <pdal/PointTable.hpp>
#include <pdal/PointView.hpp>
#include <pdal/Options.hpp>
#include <pdal/io/BufferReader.hpp>
#include <pdal/StageFactory.hpp>
#include <pdal/Stage.hpp>

#include <pdal/Filter.hpp>
#include <pdal/filters/DividerFilter.hpp>



#include <pdal/io/LasReader.hpp>
#include <pdal/io/LasHeader.hpp>

#include <pdal/io/PlyReader.hpp>
#include <pdal/io/PlyWriter.hpp>

#include <pdal/io/TextReader.hpp>
#include <pdal/io/TextWriter.hpp>

#include <pdal/PluginDirectory.hpp>
#include <pdal/PluginManager.hpp>

using namespace std;

struct las_Point
{
	double x;
	double y;
	double z;
	int intensity;
};

vector<las_Point> pointCloud;

void fillView(pdal::PointViewPtr view)
{
	for (int i = 0; i < pointCloud.size(); ++i)
	{
		view->setField(pdal::Dimension::Id::X, i, pointCloud[i].x);
		view->setField(pdal::Dimension::Id::Y, i, pointCloud[i].y);
		view->setField(pdal::Dimension::Id::Z, i, pointCloud[i].z);
		view->setField(pdal::Dimension::Id::Intensity, i, pointCloud[i].intensity);
	}
}

//pdal库读写ply
void WriteText()
{
	using namespace pdal;
	using namespace pdal::Dimension;

	std::string fileName = "d:/pdal_test.txt";

	Options xjOptions;
	xjOptions.add("filename", fileName);

	PointTable table;
	table.layout()->registerDim(Dimension::Id::X);
	table.layout()->registerDim(Dimension::Id::Y);
	table.layout()->registerDim(Dimension::Id::Z);
	table.layout()->registerDim(Dimension::Id::Intensity);
	
	PointViewPtr view(new PointView(table));
	fillView(view);
	BufferReader xjBufferReader;
	xjBufferReader.addView(view);

	StageFactory factory;
	std::string r_drivername = factory.inferReaderDriver(fileName);
	std::string w_drivername = factory.inferWriterDriver(fileName);
	Stage *writer = factory.createStage(w_drivername);
	//Stage *writer = factory.createStage("writers.e57");
	writer->addOptions(xjOptions);

	writer->setInput(xjBufferReader);
	writer->setOptions(xjOptions);
	writer->prepare(table);

	writer->execute(table);
}

//pdal库读写ply
void readText()
{
	using namespace pdal;
	using namespace pdal::Dimension;

	pdal::TextReader text_reader;
	string xjFileName = "d:/pdal_test.txt";
	pdal::Option las_opt("filename", xjFileName);//参数1："filename"（键）
	pdal::Options las_opts;
	las_opts.add(las_opt);
	las_opts.add("skip", 1);   //忽略第一行 ，一定要设置
	las_opts.add("header", "X, Y, Z, Intensity");  //读取的文件头，一定要设置
	pdal::PointTable table;

	text_reader.setOptions(las_opts);
	text_reader.prepare(table);

	pdal::PointViewSet point_view_set = text_reader.execute(table);
	pdal::PointViewPtr point_view = *point_view_set.begin();
	pdal::Dimension::IdList dims = point_view->dims();

	for (pdal::PointId idx = 0; idx < point_view->size(); ++idx)
	{
		las_Point cur_point;
		cur_point.x = point_view->getFieldAs<double>(pdal::Dimension::Id::X, idx);
		cur_point.y = point_view->getFieldAs<double>(pdal::Dimension::Id::Y, idx);
		cur_point.z = point_view->getFieldAs<double>(pdal::Dimension::Id::Z, idx);
		int Pintensity = point_view->getFieldAs<int>(pdal::Dimension::Id::Intensity, idx);
		int PointSourceID = point_view->getFieldAs<int>(pdal::Dimension::Id::PointSourceId, idx);
		double GPStime = point_view->getFieldAs<double>(pdal::Dimension::Id::GpsTime, idx);
		int red = point_view->getFieldAs<int>(pdal::Dimension::Id::Red, idx);
		int green = point_view->getFieldAs<int>(pdal::Dimension::Id::Green, idx);
		int blue = point_view->getFieldAs<int>(pdal::Dimension::Id::Blue, idx);

		pointCloud.push_back(cur_point);
	}
}

//pdal库读写ply
void readPly()
{
	using namespace pdal;
	using namespace pdal::Dimension;

	pdal::PlyReader ply_reader;
	string xjFileName = "d:/pdal_test.ply";
	pdal::Option las_opt("filename", xjFileName);//参数1："filename"（键）
	pdal::Options las_opts;
	las_opts.add(las_opt);

	pdal::PointTable table;

	ply_reader.setOptions(las_opts);
	ply_reader.prepare(table);

	pdal::PointViewSet point_view_set = ply_reader.execute(table);
	pdal::PointViewPtr point_view = *point_view_set.begin();
	pdal::Dimension::IdList dims = point_view->dims();

	for (pdal::PointId idx = 0; idx < point_view->size(); ++idx)
	{
		las_Point cur_point;
		cur_point.x = point_view->getFieldAs<double>(pdal::Dimension::Id::X, idx);
		cur_point.y = point_view->getFieldAs<double>(pdal::Dimension::Id::Y, idx);
		cur_point.z = point_view->getFieldAs<double>(pdal::Dimension::Id::Z, idx);
		int Pintensity = point_view->getFieldAs<int>(pdal::Dimension::Id::Intensity, idx);
		int PointSourceID = point_view->getFieldAs<int>(pdal::Dimension::Id::PointSourceId, idx);
		double GPStime = point_view->getFieldAs<double>(pdal::Dimension::Id::GpsTime, idx);
		int red = point_view->getFieldAs<int>(pdal::Dimension::Id::Red, idx);
		int green = point_view->getFieldAs<int>(pdal::Dimension::Id::Green, idx);
		int blue = point_view->getFieldAs<int>(pdal::Dimension::Id::Blue, idx);

		pointCloud.push_back(cur_point);
	}
}

//pdal库读写ply
void WritePly()
{
	using namespace pdal;
	using namespace pdal::Dimension;

	pdal::PlyReader ply_reader;
	
	std::string fileName = "d:/pdal_test.ply";

	Options xjOptions;
	xjOptions.add("filename", fileName);
	xjOptions.add("storage_mode", "ascii");   //Valid values are ‘ascii’, ‘little endian’, ‘big endian’. [Default: “ascii”]
	//xjOptions.add("dims", "Y=int32_t, X,Red=char");  //默认所有维度为存储类型

	//xjOptions.add("faces", true);   //Write a mesh as faces in addition to writing points as vertices. [Default: false]
	//xjOptions.add("sized_types", true);
	xjOptions.add("precision", "g-style");
	
	PointTable table;
	table.layout()->registerDim(Dimension::Id::X);
	table.layout()->registerDim(Dimension::Id::Y);
	table.layout()->registerDim(Dimension::Id::Z);
	table.layout()->registerDim(Dimension::Id::Intensity);
	PointViewPtr view(new PointView(table));
	fillView(view);
	BufferReader xjBufferReader;
	xjBufferReader.addView(view);

	StageFactory factory;
	std::string r_drivername = factory.inferReaderDriver(fileName);
	std::string w_drivername = factory.inferWriterDriver(fileName);
	Stage *writer = factory.createStage(w_drivername);
	writer->addOptions(xjOptions);

	writer->setInput(xjBufferReader);
	writer->setOptions(xjOptions);
	writer->prepare(table);
	writer->execute(table);
}

//pdal库读写las
void WriteLASPC()
{
	using namespace pdal;
	using namespace pdal::Dimension;

	std::string fileName = "d:/pdal_test.las";
	Options xjOptions;
	xjOptions.add("filename", fileName);
	
	double xoffset = 0.0, yoffset = 0.0, zoffset = 0.0;

	xjOptions.add("offset_x", xoffset);
	xjOptions.add("offset_y", yoffset);
	xjOptions.add("offset_z", zoffset);
	xjOptions.add("scale_x", 0.001);
	xjOptions.add("scale_y", 0.001);
	xjOptions.add("scale_z", 0.001);

	//xjOptions.add("type", "filters.splitter");
	//xjOptions.add("scale_z", "100");
	//xjOptions.add("scale_z", "638900.0");
	//xjOptions.add("scale_z", "835500.0");

	PointTable table;
	table.layout()->registerDim(Dimension::Id::X);
	table.layout()->registerDim(Dimension::Id::Y);
	table.layout()->registerDim(Dimension::Id::Z);
	//table.layout()->registerDim(Dimension::Id::Intensity);

	PointViewPtr view(new PointView(table));
	fillView(view);

	BufferReader xjBufferReader;
	xjBufferReader.addView(view);
	
	StageFactory factory;
	std::string r_drivername = factory.inferReaderDriver(fileName);
	std::string w_drivername = factory.inferWriterDriver(fileName);
	//Stage *writer = factory.createStage(w_drivername);
	Stage *writer = factory.createStage("writers.las");
	
	writer->setInput(xjBufferReader);
	writer->setOptions(xjOptions);
	writer->prepare(table);
	writer->execute(table);
}

//pdal库读写las
void readlasPC()
{
	string xjFileName = "F:\\pointCloud_data\\zhangjiang8\\pointcloud\\181013_030701-11-16-47-028.las";//中文路径可能报错
	pdal::Option las_opt("filename", xjFileName);//参数1："filename"（键）
	pdal::Options las_opts;
	las_opts.add(las_opt);

	pdal::PointTable table;
	pdal::LasReader las_reader;



	las_reader.setOptions(las_opts);
	las_reader.prepare(table);

	pdal::PointViewSet point_view_set = las_reader.execute(table);

	pdal::PointViewPtr point_view = *point_view_set.begin();
	pdal::Dimension::IdList dims = point_view->dims();

	pdal::LasHeader las_header = las_reader.header();

	unsigned int PointCount = las_header.pointCount();
	//char* projstr = table.spatialRef().getWKT(pdal::SpatialReference::eCompoundOK).c_str();
	double scale_x = las_header.scaleX();
	double scale_y = las_header.scaleY();
	double scale_z = las_header.scaleZ();
	double offset_x = las_header.offsetX();
	double offset_y = las_header.offsetY();
	double offset_z = las_header.offsetZ();
	double xmin = las_header.minX();
	double xmax = las_header.maxX();
	double ymin = las_header.minY();
	double ymax = las_header.maxY();
	double zmin = las_header.minZ();
	double zmax = las_header.maxZ();

	unsigned int n_features = las_header.pointCount();
	int count = 0;
	for (pdal::PointId idx = 0; idx < point_view->size(); ++idx)
	{
		las_Point cur_point;
		cur_point.x = point_view->getFieldAs<double>(pdal::Dimension::Id::X, idx);
		cur_point.y = point_view->getFieldAs<double>(pdal::Dimension::Id::Y, idx);
		cur_point.z = point_view->getFieldAs<double>(pdal::Dimension::Id::Z, idx);
		cur_point.intensity = point_view->getFieldAs<int>(pdal::Dimension::Id::Intensity, idx);
		int PointSourceID = point_view->getFieldAs<int>(pdal::Dimension::Id::PointSourceId, idx);
		double GPStime = point_view->getFieldAs<double>(pdal::Dimension::Id::GpsTime, idx);
		int red = point_view->getFieldAs<int>(pdal::Dimension::Id::Red, idx);
		int green = point_view->getFieldAs<int>(pdal::Dimension::Id::Green, idx);
		int blue = point_view->getFieldAs<int>(pdal::Dimension::Id::Blue, idx);

		
		pointCloud.push_back(cur_point);
	}


	/*double deltaX = xmax - xmin;
	cout << "deltaX=" << fixed << setprecision(4) << deltaX << endl;

	double ymin = las_header.minY();
	double ymax = las_header.maxY();
	double deltaY = ymax - ymin;
	cout << "deltaY=" << fixed << setprecision(4) << deltaY << endl;

	double zmin = las_header.minZ();
	double zmax = las_header.maxZ();
	double deltaZ = zmax - zmin;
	cout << "deltaZ=" << fixed << setprecision(4) << deltaZ << endl;

	unsigned int n_features = las_header.pointCount();
	cout << "PointCount=" << n_features << endl;*/
}

int main()
{
	readlasPC();
	//readPly();
	//readText();

	//WriteLASPC();
	//WritePly();
	//WriteText();
	dividerLAS();

	system("pause");
    return 0;
}

 

chen_jared
关注
专栏目录
pdal读取.las文件(将las数据读取成XYZRGB格式)//las转txt
swjtuJW的博客
08-05 3202
学习笔记：pdal读取.las文件(将las数据读取成XYZRGB格式)//las转txt 本着对c++里一些库的学习态度，整理了一下在实际工程项目中写下的代码分享给用得到的人；有不对的地方希望指正，勿喷。 pdal可以读取.las数据中的详细信息；由于本人只用到点云的xyzrgb信息，故将.las直接转换成: pcl::PointCloud<PointXYZRGB> 读取数据用到了pdal库,写入数据用到了pcl库! 一、pdal是什么？ PDAL是点云数据处理的库。这是一个C/C++开源库
LiDAR点云处理软件
weixin_30640291的博客
02-18 1888
转自https://www.cnblogs.com/l-kong/ 1、主要介绍当前主流的激光雷达数据处理库、开源软件、商业软件等。 英文主页： http://opentopo.sdsc.edu/tools/listTools2、日期 工具名称 工具类型 关键词 评分2017年1月12日 planlauf /地形说明： planlauf / TERRAIN是Windows 7/8/10易于...
PDAL点云库官方文档
11-19
该文档是直接从PDAL官网上下载下来的。包括了PDAL库的相关介绍、代码例子等，可以快速的入门pdal库。
LASTools、PDAL、CGAL点云读取效率
March的专栏
11-22 1268
在本文之前，让我们先悼念一位大师：Martin Isenburg，德国人，是rapidlasso的所有者，也是LASTools、LASzip的创造者，于2021年9月9日在家中自杀，天才总是孤独的，也是具有社会担当的，给我们留下宝贵的财富，愿大师在另一世界再无痛苦。 本文测试的是读取las以及压缩后laz文件的效率问题，主要测试LASTools、PDAL以及CGAL，PCL没有读取las文件的接口，这里就不进行测试了。 测试数据分别是包含60万、2千...
osg+pdal大数据点云的LOD可视化
Diasepan的博客
05-21 1298
大数据点云可视化
点云：PDAL读写点云——VS2017，Qt5.9.8
累了就要打游戏
08-18 4729
1 配置环境 1.1 属性——C/C++——常规——附加包含目录，添加头文件路径C:\OSGeo4W64\include 1.2 属性——链接器——输入——附加依赖项，添加lib文件C:\OSGeo4W64\lib\pdalcpp.lib和C:\OSGeo4W64\lib\pdal_util.lib 2 分享给有需要的人，代码质量勿喷 2.1 头文件PDALtest.h #pr...
读取PCD文件（一）
m0_62246041的博客
09-25 1325
读取PCD是为了使用数据集中的数据而不是在实车上使用。 一般来说PCD文件时用bin格式二进制存储，有利于读取速度，但是必须用编程方法进行读取。 根据pcl官网的介绍分为header和正文，header里有如下一些项目 VERSION版本 FIELDS数据各个维度的详细内容 SIZE每个维度内容的存储占用字节数 TYPE带符号整型，无符号整型，浮点型等 COUNT每个维度有多少元素 WIDTH对于格式化的...
点云读取文件
09-10
例如，对于PCL库，可以使用`pcl::io::loadPCDFile`函数读取PCD文件，而对于LAS/LAZ格式，可以使用libLAS或PDAL（Point Data Abstraction Library）库。 在处理点云数据时，通常会涉及以下步骤： 1. **数据预处理**...
cloudcompare pdal
最新发布
05-24
而PDAL则是一款基于开源的点云数据处理库，支持多种点云数据格式的读写和转换，并提供了众多的点云数据处理工具和算法。 两者在功能上有所重叠，但也有一些区别。相对来说，CloudCompare更加注重于点云数据的可视化...
LAS文件读取库
11-23
详情配置请参考博客https://blog.csdn.net/wokaowokaowokao12345/article/details/50765989
LAS 文件标题阅读器：读取 LAS 文件的公共标题块和可变记录长度-matlab开发
05-30
读取 LAS 文件的公共标题块和可变记录长度 可用于机载和地面 LiDAR/激光扫描仪文件格式 激光雷达数据文件 - *.LAS 运行 LASreader.m 作者 - Sandeep Sasidharan 联系方式 - sandeep40@gmail.com； sasidhar.sandeep@gmail.com 主页 - http://sandeepsasidharan.webs.com 更多关于 Las 格式 - http://www.asprs.org/a/society/committees/standards/asprs_las_format_v12.pdf
LASTools读取操作
02-28 1935
LASTools可以读快速las数据，对于1.2版本包括点的r、g、b颜色属性，但是其取值在0-65535范围内（2的16）。需要将其转换成正常的0-255范围内，其示例如下： #include "lasreader.hpp" #include "laswriter.hpp" #include<fstream> #include<iostream> #include<iomanip> using namespace std; void main() { LASr
PDAL点云处理库介绍
小白学视觉
05-26 1106
点击上方“小白学视觉”，选择加"星标"或“置顶”重磅干货，第一时间送达PDAL是点云数据处理的库。这是一个C/C++开源库，用于点云数据的转换和处理。尽管该库中许多工具的重点和发展都起源于激光雷达点云数据的处理，但它也不限于激光雷达数据。什么是PDAL？PDAL是点云数据处理的库。这是一个C/C++开源库，用于点云数据的转换和处理。尽管该库中许多工具的重点和发展都起源于激光雷达点云数据的处理，但它...
LAS点云的读写以及投影信息的读写（LASlib）
March的专栏
10-27 5043
本文使用使用LASlib库完成点云及投影信息的读写。具体安装方式可自行网上搜索，推荐使用vcpkg安装，非常方便。 对于库的使用下面直接上代码： 头文件： #include <LASlib/lasreader.hpp> #include <LASlib/laswriter.hpp> 1、las文件的读 先是LASreadOpener对象构造，设置文件名，然后进行点云的打开。 LASreadOpener lasreadopener; lasreadope
PDAL库点云读写环境配置及基本操作
zou_ys88的博客
05-13 1959
点云处理库 pdal 环境配置 读写操作
PDAL: Point Data Abstraction Library
gitblog_00034的博客
03-14 528
PDAL: Point Data Abstraction Library PDAL是一个开源的C++库，用于处理点云数据。它提供了一个灵活的数据抽象层，可以读取、转换、过滤、聚合和可视化大量的点云数据。此外，PDAL还支持多种输入和输出格式，包括LAS、LAZ、XYZ、BAG等。 PDAL的功能和用途 数据读取和写入：PDAL可以从各种源读取点云数据，并将其导出到不同的格式。 数据转换和过滤：P...
点云：点云格网化 （师弟师妹）简单教学版
累了就要打游戏
11-18 5051
一、原理与方法 点云格网化是点云数据处理的常用方法，是一种降维或压缩的方法。（论文里烂大街了） 基本原理是将点放入不同的格网，进而以格网为基本操作单元，对点云进行分块、抽稀、语义分割等。二维格网化是投影，类似像素；三维格网是体素。 基本方法是根据范围，划分成规则或不规则的格网，将不同的点分别放入不同的格网。下面是规则格网的划分。 （1） 式中，Xmax、Xmin、Ymax、Ymin、Zmax、Zmin分别代表点云区域X、Y、Z坐标.
PDAL中filters.divider过滤器对大点云文件进行分割
Jared_Chen的博客
11-19 593
//pdal库将las文件分割为多个文件 void dividerLAS() { using namespace pdal; using namespace pdal::Dimension; std::string in_filename = "F:\\pointCloud_data\\zhangjiang8\\pointcloud\\181013_030701-11-16-47-028.las";//中文路径可能报错 std::string fileName = "d:/out_#.la...
las数据的读写及注意事项总结
06-28 8730
1.C#使用laslib或liblas读写las数据 2.Python使用laspy或pylas读写las数据 3.C++使用laslib读写las数据
pdal LasField
06-09
PDAL LasField是一个类，用于描述LAS格式文件中的一种数据字段（field）。LAS格式是一种常见的激光雷达数据格式，其中包含了点云数据中的各种属性信息，例如点的位置、强度、回波时间等。LasField类用于描述LAS格式文件中的一个字段，包括该字段的名称、数据类型、数据长度等信息。LasField类通常用于PDAL中读取和写入LAS格式文件时，用于解析和生成LAS格式文件中的各个数据字段。例如，可以使用LasField类来读取LAS格式文件中的X、Y、Z坐标、强度和分类等信息。
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值