【点云学习PCL 】一

点云学习

说明：仅做学习使用，侵删
参考网址1

一、点云基础

0 概述

PCL（Point Cloud Library）是用于 2D/3D 图像和点云处理的大型开源跨平台的 C++ 编程库，PCL 框架实现了大量点云相关的通用算法和高效的数据结构，PCL 是 BSD 授权方式，可以免费进行商业和学术应用。相比图像数据，点云数据多了一个维度，因此能够更全面的刻画三维场景中的对象
支持多种操作系统，可以在 Windows、Linux、MacOS X、Android、部分嵌入式实时系统上运行。
应用领域：机器人、CAD/CAM、逆向工程、遥感测量、VR/AR、人机交互等
PCL 内容涉及了点云的获取、滤波、分割、配准、检索、特征提取、特征估计，表面重建、识别、模型拟合、追踪、曲面重建、可视化等。

PCL 中的所有模块和算法都是通过 Boost 共享指针来传送数据，因而避免了多次复制系统中已存在的数据

1.1 点云及其可视化

基本显示的代码逻辑：

  • 创建一个智能指针，用作指向点云文件

    pcl::PointCloud<pcl::PointXYZRGBA>::Ptr cloud(new pcl::PointCloud<pcl::PointXYZRGBA>);

  • 加载点云文件

    pcl::io::loadPCDFile("./data/pcl_logo.pcd", *cloud);

  • 创建显示的对象

    pcl::visualization::CloudViewer viewer("Cloud Viewer");

  • 对象调用显示点云的方法

    viewer.showCloud(cloud);

拓展：

  • 如果显示点云之后，判断点云是否退出，如果不退出可以处理点云

     while (!viewer.wasStopped()) {// 你可以在这里对点云做很多处理}

输出结果为：
在这里插入图片描述

1.2 点云数据

根据激光测量原理得到的点云，包含三维坐标信息和激光反射强度信息，激光反射强度与仪器的激光发射能量、波长，目标的表面材质、粗糙程度、入射角相关。根据摄影测量原理得到的点云，包括三维坐标（xyz）和颜色信息。结合两个原理的多传感器融合技术（多见于手持式三维扫描仪），能够同时得到这三种信息。

也即：三维坐标信息 ( x , y , z ) (x,y,z) (x,y,z) 强度信息 I I I（intensity）,颜色信息（RGB）

基本数据类型PointCloud

PCL的基本数据类型是PointCloud，一个PointCloud是一个C++的模板类，它包含了以下字段：

  • width(int)

    ：指定点云数据集的宽度

    • 对于无组织格式的数据集，width代表了所有点的总数
    • 对于有组织格式的数据集，width代表了一行中的总点数

  • height(int)

    ：制定点云数据集的高度

    • 对于无组织格式的数据集，值为1
    • 对于有组织格式的数据集，表示总行数

  • points(std::vector<PointT>)：包含所有PointT类型的点的数据列表

1.2.1 内置数据类型

类型主要有：

  • PointXYZ - float x, y, z
  • PointXYZI - float x, y, z, intensity
  • PointXYZRGB - float x, y, z, rgb
  • PointXYZRGBA - float x, y, z, uint32_t rgba
  • Normal - float normal[3], curvature 法线方向，对应的曲率的测量值
  • PointNormal - float x, y, z, normal[3], curvature 采样点，法线和曲率
  • Histogram - float histogram[N] 用于存储一般用途的n维直方图

1.2.2 生成PCD文件分析

以fastlio2输出的点云为例

# .PCD v0.7 - Point Cloud Data file format
VERSION 0.7
FIELDS x y z intensity normal_x normal_y normal_z curvature
SIZE 4 4 4 4 4 4 4 4
TYPE F F F F F F F F
COUNT 1 1 1 1 1 1 1 1
WIDTH 27211527
HEIGHT 1
VIEWPOINT 0 0 0 1 0 0 0
POINTS 27211527
DATA binary
̆I��E�@���     

# .PCD v.7 - Point Cloud Data file format
VERSION .7			# 版本号
FIELDS x y z rgb		#  指定一个点可以有的每一个维度和字段的名字
SIZE 4 4 4 4			# 用字节数指定每一个维度的大小。例如：
TYPE F FFF			# 用一个字符指定每一个维度的类型
COUNT 1 1 1 1			# 指定每一个维度包含的元素数目，默认情况下，如果没有COUNT，所有维度的数目被设置成1
WIDTH 640    			# 像图像一样的有序结构，有640行和480列，
HEIGHT 480    			# 这样该数据集中共有640*480=307200个点
VIEWPOINT 0 0 0 1 0 0 0		# 指定数据集中点云的获取视点 视点信息被指定为平移（tx ty tz）+ 四元数（qw qx qy qz）
POINTS 307200			# 指定点云中点的总数。从0.7版本开始，该字段就有点多余了
DATA ascii			# 指定存储点云数据的数据类型。支持两种数据类型：ascii和二进制
0.93773 0.33763 0 4.2108e+06
0.90805 0.35641 0 4.2108e+06

注解：

TYPE：用一个字符指定每一个维度的类型

  • I 表示有符号类型 int8（char）、int16（short）和 int32（int）
  • U 表示无符号类型 uint8（unsigned char）、uint16（unsigned short）和 uint32（unsigned int）
  • F 表示浮点类型

FIELDS：指定一个点可以有的每一个维度和字段的名字

  • FIELDS x y z
  • FIELDS x y z rgb
  • FIELDS x y z normal_x normal_y normal_z

PCD 文件的文件头部分必须按上面的顺序精确指定

1.2.3 点云序列化常用

  • 反序列化—将生成的点云存储到pcd文件中，并选择对应的存储类型

    将对象的状态信息转换为可以存储或传输的形式（文件、缓冲、或经由网络中发送）的过程。例如，点云在 C++ 中以类对象的形式存放，而点云存放为文件的时候是一个点 PCD、PLY 等文件，将代码中暂存的点云对象转成点 PCD 文件这一操作就叫做序列化。

    对应的代码：

    //
// Created by pf on 2023/10/25.
//
// 随机生成了 5 个点，并将之以 ASCII 形式保存（序列化）在 test_pcd.pcd 文件中：save_pcd.cpp
#include <iostream>
#include <pcl/io/pcd_io.h>
#include <pcl/point_types.h>int main(int argc, char **argv) {pcl::PointCloud<pcl::PointXYZ> cloud;cloud.width =5;cloud.height = 1;cloud.is_dense = false;cloud.points.resize(cloud.width*cloud.height);// Fill in the cloud datastd::cout << rand() << std::endl;std::cout << rand() / (RAND_MAX + 1.0f) << std::endl;std::cout << 1024 * rand() / (RAND_MAX + 1.0f) << std::endl;// 随机生成5个点for (size_t i = 0; i < cloud.points.size(); ++i) {cloud.points[i].x = 1024 * rand() / (RAND_MAX + 1.0f);cloud.points[i].y = 1024 * rand() / (RAND_MAX + 1.0f);cloud.points[i].z = 1024 * rand() / (RAND_MAX + 1.0f);}pcl::io::savePCDFileASCII("pcd_serialization.pcd", cloud);// pcl::io::savePCDFileBinary("test_pcd_binary.pcd", cloud);std::cerr << "Saved " << cloud.points.size() << " data points to test_pcd.pcd." << std::endl;for (size_t i = 0; i < cloud.points.size(); ++i)std::cerr << "    " << cloud.points[i].x << " " << cloud.points[i].y << " " << cloud.points[i].z << std::endl;return (0);
}

  • 反序列化

    和序列化相反，反序列化则是将一些已经序列化过的数据进行反序列化得到原来的对象。例如，一个点 PCD 文件，将它读取成一个点云对象这一操作就是反序列化

    对应的代码

    //
// Created by pf on 2023/10/25.
//
// 对一个点云进行反序列化操作，将之保存到 PointCloud 对象中：load_pcd.cpp
#include <iostream>
#include <pcl/io/pcd_io.h>
#include <pcl/point_types.h>int main(int argc, char **argv) {// 准备pcl::PointXYZ类型的点云pcl::PointCloud<pcl::PointXYZ>::Ptr cloud(new pcl::PointCloud<pcl::PointXYZ>);// 将pcd中的数据加载到cloud中if (pcl::io::loadPCDFile<pcl::PointXYZ>("./data/bunny.pcd", *cloud) == -1) {         //* load the filePCL_ERROR ("Couldn't read file bunny.pcd \n");return (-1);}std::cout << "Loaded "<< cloud->width * cloud->height<< " data points from test_pcd.pcd with the following fields: "<< std::endl;for (size_t i = 0; i < cloud->points.size(); ++i)std::cout << "    " << cloud->points[i].x<< " " << cloud->points[i].y<< " " << cloud->points[i].z << std::endl;return (0);
}

其他文件格式

  • PLY 是一种多边形文件格式，由 Stanford 大学的 Turk 等人设计开发
    常用的点云数据文件，ply 文件不仅可以存储点数据，而且可以存储网格数据
  • STL 是 3D Systems 公司创建的模型文件格式，主要应用于 CAD、CAM 领域
  • OBJ 是从几何学上定义的文件格式，首先由 Wavefront Technologies 开发
  • X3D 是符合 ISO 标准的基于 XML 的文件格式，表示 3D 计算机图形数据

1.3 点云的基本处理方法

在这里插入图片描述

1.3.1文件简介

common

  • pcl_common 中主要是包含了 PCL 库常用的公共数据结构和方法，比如 PointCloud 的类和许多用于表示点、曲面、法向量和特征描述等点的类型，用于计算距离、均值、协方差、角度转换以及几何变化的函数

  • angles.h 定义了标准的 C 接口的角度计算函数

  • centriod.h 定义了中心点的估算以及协方差矩阵的计算

  • commo.h 标准的 C 以及 C++ 类，是其他 common 函数的父类

  • distance.h 定义标准的 C 接口用于计算距离

  • file_io.h 定义了一些文件帮助写或者读方面的功能

  • random.h 定义一些随机点云生成的函数

  • geometry.h 定义一些基本的几何功能的函数

  • intersection.h 定义线与线相交的函数

  • norm.h 定义了标准的 C 方法计算矩阵的正则化

  • time.h 定义了时间计算的函数

  • Point_types.h 定义了所有 PCL 实现的点云的数据结构的类型

1.3.2 点云pcl和ROS之间相互转换

/*
ConstPtr：const pointer（常量指针），可理解为传递常量数据的共享指针1. 为什么要用指针而不直接传数据，因为数据太大了，值传递需要进行对数据的复制而指针引用不需要，因为涉及到右值传参和类型推断，不加引用的话类型可能发生变化2. ::constPtr 可以理解为指向常量的共享指针，::Ptr 可以理解为指向一个（可修改的）共享指针3. ::constPtr& 可以理解为对一个常量消息的共享指针的引用，::Ptr& 可以理解为对一个（可修改的）消息的共享指针的引用4. 使用共享指针，它是一种智能指针，可以自动管理他指向的内存的存储时间，不必手动释放内存
*/
void elevatedCloudCallback(const sensor_msgs::PointCloud2::ConstPtr &input_cloud) {pcl::PointCloud<PointT>::Ptr output_cloud(new pcl::PointCloud<PointT>());pcl::fromROSMsg(*input_cloud, *output_cloud);//pcl::toROSMsg(pcl::PointCloud<pcl::PointXYZ>, sensor_msgs::PointCloud2); // 逆过程
}

1.4 点云可视化工具

1.4.1 pcl_viewer工具

  • 基本使用

    进入：pcl_viewer xxxxx.pcd

    帮助：在界面中输入h，可以在控制台看到帮助信息

    退出：界面中输入q

    放大缩小：鼠标滚轮 或 Alt + [+/-]

    平移：Shift+鼠标拖拽

    旋转：Ctrl+鼠标拖拽

  • 其他技巧

    修改点颜色：数字1,2,3,4,5…9，重复按1可切换不同颜色方便观察

    放大缩小点：放大Ctrl+Shift+加号，缩小 Ctrl+减号

    保存截图：j

    显示颜色尺寸：u

    显示比例尺寸：g

    在控制列出所有几何和颜色信息：l

  • 鼠标选点打印坐标

    选点模式进入：pcl_viewer -use_point_picking bunny.pcd

    选择指定点：shift+鼠标左键

1.4.2 CloudViewer 程序的内部对象方法

//
// Created by pf on 2023/10/25.
//
#include <pcl/visualization/cloud_viewer.h>
#include <iostream>
#include <pcl/io/io.h>
#include <pcl/io/pcd_io.h>int count = 0;void viewerPsycho(pcl::visualization::PCLVisualizer &viewer);void viewerOneOff(pcl::visualization::PCLVisualizer &viewer);int main() {pcl::PointCloud<pcl::PointXYZRGBA>::Ptr cloud(new pcl::PointCloud<pcl::PointXYZRGBA>);pcl::io::loadPCDFile("./data/scans.pcd", *cloud);pcl::visualization::CloudViewer viewer("Cloud Viewer");
//    viewer.showCloud(cloud);viewer.runOnVisualizationThreadOnce(viewerOneOff);viewer.runOnVisualizationThread(viewerPsycho);while (!viewer.wasStopped()) {count++;}return 0;
}void viewerPsycho(pcl::visualization::PCLVisualizer &viewer) {static unsigned cnt = 0;std::stringstream ss;ss << "Once per viewer loop: " << cnt++;// 每次刷新时，移除text，添加新的textviewer.removeShape("text", 0);viewer.addText(ss.str(), 200, 300, "text", 0);//FIXME: possible race condition here:count++;
}void viewerOneOff(pcl::visualization::PCLVisualizer &viewer) {// 设置背景色为粉红色viewer.setBackgroundColor(1.0, 0.5, 1.0);pcl::PointXYZ o;o.x = 1.0;o.y = 0.0;o.z = 0.0;//    添加一个圆心为o，半径为0.25m的球体viewer.addSphere(o, 0.25, "sphere", 0);std::cout << "i only run once" << std::endl;
}
对应的结果

在这里插入图片描述

PCLVisualizer 雷达对象方法

方法类似，不作演示，没有点云文件。

//
// Created by pf on 2023/10/25.
//
#include <pcl/io/io.h>
#include <pcl/io/pcd_io.h>
#include <pcl/visualization/pcl_visualizer.h>int main(int argc, char **argv) {pcl::PointCloud<pcl::PointXYZ>::Ptr cloud(new pcl::PointCloud<pcl::PointXYZ>);pcl::io::loadPCDFile("./data/bunny.pcd", *cloud);pcl::PointCloud<pcl::PointXYZRGB>::Ptr cloud_milk(new pcl::PointCloud<pcl::PointXYZRGB>);pcl::io::loadPCDFile("./data/milk_color.pcd", *cloud_milk);// 创建PCLVisualizerboost::shared_ptr<pcl::visualization::PCLVisualizer> viewer(new pcl::visualization::PCLVisualizer("3D Viewer"));// 设置背景色为灰色（非必须）viewer->setBackgroundColor(0.05, 0.05, 0.05, 0);// 添加一个普通点云 (可以设置指定颜色，也可以去掉single_color参数不设置)pcl::visualization::PointCloudColorHandlerCustom<pcl::PointXYZ> single_color(cloud, 0, 255, 0);viewer->addPointCloud<pcl::PointXYZ>(cloud, single_color, "sample cloud");viewer->setPointCloudRenderingProperties(pcl::visualization::PCL_VISUALIZER_POINT_SIZE, 2, "sample cloud");// 再添加一个彩色点云及配置pcl::visualization::PointCloudColorHandlerRGBField<pcl::PointXYZRGB> rgb(cloud_milk);viewer->addPointCloud<pcl::PointXYZRGB>(cloud_milk, rgb, "sample cloud milk");viewer->setPointCloudRenderingProperties(pcl::visualization::PCL_VISUALIZER_POINT_SIZE, 3, "sample cloud milk");// 添加一个0.5倍缩放的坐标系（非必须）viewer->addCoordinateSystem(0.5);// 直到窗口关闭才结束循环while (!viewer->wasStopped()) {// 每次循环调用内部的重绘函数viewer->spinOnce();}return 0;
}

1.4.3 矩阵变换

将一个点云经过 z z z轴旋转45°，然后经过平移变换。

//
// Created by pf on 2023/10/25.
//
#include <iostream>#include <pcl/io/pcd_io.h>
#include <pcl/point_cloud.h>
#include <pcl/console/parse.h>
#include <pcl/common/transforms.h>
#include <pcl/visualization/pcl_visualizer.h>int main(int argc, char **argv) {const std::string file = "./data/bunny.pcd";pcl::PointCloud<pcl::PointXYZ>::Ptr source_cloud(new pcl::PointCloud<pcl::PointXYZ>());if (pcl::io::loadPCDFile(file, *source_cloud) < 0) {std::cout << "Error loading point cloud " << argv[1] << std::endl << std::endl;return -1;}Eigen::Matrix4f transform_1 = Eigen::Matrix4f::Identity();// 定义旋转矩阵auto m_data = float(M_PI_4);transform_1 << cos(m_data), -sin(m_data), 0, 2.5,sin(m_data), cos(m_data), 0, 0,0, 0, 0, 0,0, 0, 0, 0;// 使用仿射变换Eigen::Affine3f transform_2 = Eigen::Affine3f::Identity();transform_2.translation() << 0.8, 0.0, 0.0;transform_2.rotate(Eigen::AngleAxisf(m_data, Eigen::Vector3f::UnitZ()));pcl::PointCloud<pcl::PointXYZ>::Ptr transformed_cloud(new pcl::PointCloud<pcl::PointXYZ>());pcl::transformPointCloud(*source_cloud, *transformed_cloud, transform_2);pcl::visualization::PCLVisualizer viewer("Matrix transformation example");pcl::visualization::PointCloudColorHandlerCustom<pcl::PointXYZ>source_cloud_color_handler(source_cloud, 255, 255, 255);// whiteviewer.addPointCloud(source_cloud, source_cloud_color_handler, "original_cloud");pcl::visualization::PointCloudColorHandlerCustom<pcl::PointXYZ>transformed_cloud_color_handler(transformed_cloud, 230, 20, 20); // Red\
viewer.addPointCloud(transformed_cloud, transformed_cloud_color_handler, "transformed_cloud");// 设置坐标系系统viewer.addCoordinateSystem(0.5, "cloud", 0);// 设置背景色viewer.setBackgroundColor(0.05, 0.05, 0.05, 0); // Setting background to a dark grey// 设置渲染属性（点大小）viewer.setPointCloudRenderingProperties(pcl::visualization::PCL_VISUALIZER_POINT_SIZE, 2, "original_cloud");// 设置渲染属性（点大小）viewer.setPointCloudRenderingProperties(pcl::visualization::PCL_VISUALIZER_POINT_SIZE, 2, "transformed_cloud");//viewer.setPosition(800, 400); // Setting visualiser window positionwhile (!viewer.wasStopped()) { // Display the visualiser until 'q' key is pressedviewer.spinOnce();}return 0;
}

在这里插入图片描述

二、点云的拓扑结构

3.1 点云存储结构，K-D树

通过3D相机（雷达、激光扫描、立体相机）获取到的点云，一般数据量较大，分布不均匀，数据主要表征了目标物表面的大量点的集合，这些离散的点如果希望实现基于邻域关系的快速查找比对功能，就必须对这些离散的点之间建立拓扑关系。常见的空间索引一般是自上而下逐级划分空间的各种索引结构，包括BSP树，k-d tree、KDB tree、R tree、CELL tree、八叉树等。有了这些关系，我们就可以实现点云的降采样，计算特征向量，点云匹配，点云拆分等功能。

原理概述

k-d tree（ k-dimensional tree）是计算机科学中用于在k维空间中一些点建立关系的数据结构。它是一个包含特定约束的二叉搜索树。k-d tree对于范围搜索和最近邻居搜索非常有用。我们通常只处理三维空间的点云，因此我们所有的k-d树都是三维空间的。

k-d树的每个级别都使用垂直于相应轴的超平面沿特定维度拆分所有子级。在树的根部，所有子项都将根据第一维进行拆分（即，如果第一维坐标小于根，则它将位于左子树中，如果大于根，则显然位于右边的子树）。树中向下的每个级别都在下一个维度上划分，其他所有元素都用尽后，将返回到第一个维度。他们构建k-d树的最有效方法是使用一种分区方法，例如快速排序所用的一种方法，将中值点放置在根上，所有具有较小一维值的东西都放置在根部，而右侧则更大。然后，在左右两个子树上都重复此过程，直到要分区的最后一棵树仅由一个元素组成。

KD树的原理查资料即可。

手写kd树

//
// Created by pf on 2023/10/26.
//
#pragma once
#include <algorithm>
#include <iostream>
#include <set>
#include <stack>
#include <vector>
#include <cmath>
#include <cfloat>
struct KDNode
{int index;int axis;KDNode *left;KDNode *right;KDNode(int index, int axis, KDNode *left = nullptr, KDNode *right = nullptr){this->index = index;this->axis = axis;this->left = left;this->right = right;}
};template <class T>
class KDTree
{
private:int dimension;KDNode *root;KDNode *build(std::vector<T> &);std::set<int> visited;std::stack<KDNode *> queueNode;std::vector<T> m_data;void release(KDNode *);void printNode(KDNode *);int chooseAxis(std::vector<T> &);void dfs(KDNode *, T);// 点与点之间的距离inline double distanceT(KDNode *, T);inline double distanceT(int, T);// 点与超平面的距离inline double distanceP(KDNode *, T);// 检查父节点超平面是否在超球体中inline bool checkParent(KDNode *, T, double);public:KDTree(std::vector<T> &, int);~KDTree();void Print();int findNearestPoint(T);
};template <class T>
KDTree<T>::KDTree(std::vector<T> &data, int dim)
{dimension = dim;m_data = data;root = build(data);
}template <class T>
KDTree<T>::~KDTree()
{release(root);
}template <class T>
void KDTree<T>::Print()
{printNode(root);
}template <class T>
void KDTree<T>::release(KDNode *node)
{if (node){if (node->left)release(node->left);if (node->right)release(node->right);delete node;node = nullptr;}
}template <class T>
KDNode *KDTree<T>::build(std::vector<T> &data)
{if (data.empty())return nullptr;std::vector<T> temp = data;int mid_index = static_cast<int>(data.size() / 2);int axis = data.size() > 1 ? chooseAxis(temp) : -1;std::sort(temp.begin(), temp.end(), [axis](T a, T b){ return a[axis] < b[axis]; });std::vector<T> leftData, rightData;leftData.assign(temp.begin(), temp.begin() + mid_index);rightData.assign(temp.begin() + mid_index + 1, temp.end());KDNode *leftNode = build(leftData);KDNode *rightNode = build(rightData);KDNode *rootNode;rootNode = new KDNode(temp[mid_index].index, axis, leftNode, rightNode);return rootNode;}template <class T>
void KDTree<T>::printNode(KDNode *node)
{if (node){std::cout << "Index: " << node->index << "\tAxis: " << node->axis << std::endl;printNode(node->left);printNode(node->right);}
}template <class T>
int KDTree<T>::chooseAxis(std::vector<T> &data)
{int axis = -1;double maxVar = -1.0;size_t size = data.size();for (int i = 0; i < dimension; i++)// x数据方向的方差最大进行设置axis，也就是这个方向是最离散的{double mean = 0;double Var = 0;for (int j = 0; j < size; j++){mean += static_cast<double>(data[j][i]);}mean = mean / static_cast<double>(size);for (int j = 0; j < size; j++){Var += (static_cast<double>(data[j][i]) - mean) * (static_cast<double>(data[j][i]) - mean);}Var = Var / static_cast<double>(size);// 找到最大的这个方差对应的轴if (Var > maxVar){axis = i;// 这个就是这个轴对应最大方差的轴maxVar = Var;}}return axis;
}template <class T>
int KDTree<T>::findNearestPoint(T pt)
{visited.clear();while (!queueNode.empty())queueNode.pop();double min_dist = DBL_MAX;int resNodeIdx = -1;dfs(root, pt);while (!queueNode.empty()){KDNode *curNode = queueNode.top();queueNode.pop();double dist = distanceT(curNode, pt);if (dist < min_dist){min_dist = dist;resNodeIdx = curNode->index;}if (!queueNode.empty()){KDNode *parentNode = queueNode.top();int parentAxis = parentNode->axis;int parentIndex = parentNode->index;if (checkParent(parentNode, pt, min_dist)){if (m_data[curNode->index][parentNode->axis] < m_data[parentNode->index][parentNode->axis])dfs(parentNode->right, pt);elsedfs(parentNode->left, pt);}}}return resNodeIdx;
}template <class T>
void KDTree<T>::dfs(KDNode *node, T pt)
{if (node){if (visited.find(node->index) != visited.end())return;queueNode.push(node);visited.insert(node->index);if (pt[node->axis] <= m_data[node->index][node->axis] && node->left)dfs(node->left, pt);else if (pt[node->axis] >= m_data[node->index][node->axis] && node->right)dfs(node->right, pt);else if ((node->left == nullptr) ^ (node->right == nullptr)){dfs(node->left, pt);dfs(node->right, pt);}}
}template <class T>
double KDTree<T>::distanceT(KDNode *node, T pt)
{double dist = 0;for (int i = 0; i < dimension; i++){dist += (m_data[node->index][i] - pt[i]) * (m_data[node->index][i] - pt[i]);}return sqrt(dist);
}template <class T>
double KDTree<T>::distanceT(int index, T pt)
{double dist = 0;for (int i = 0; i < dimension; i++){dist += (m_data[index][i] - pt[i]) * (m_data[index][i] - pt[i]);}return sqrt(dist);
}template <class T>
double KDTree<T>::distanceP(KDNode *node, T pt)
{int axis = node->axis;double dist = static_cast<double>(pt[axis] - m_data[node->index][axis]);return dist >= 0 ? dist : -dist;
}template <class T>
bool KDTree<T>::checkParent(KDNode *node, T pt, double distT)
{double distP = distanceP(node, pt);return distP <= distT;
}

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.mzph.cn/news/121155.shtml

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈email:809451989@qq.com，一经查实，立即删除！

相关文章

基于XML的Web服务Java接口(JAX-WS)、Jakarta XML Web Services Eclipse 实现

基于XML的Web服务Java接口(JAX-WS)、Jakarta XML Web Services Eclipse 实现

简介 JAX-WS&#xff08;Java API for XML-Based Web Services&#xff09;&#xff0c;是创建web服务的Java编程接口&#xff0c;特别是SOAP服务。是Java XML编程接口之一&#xff0c;是Java SE 和Java EE 平台的一部分。 JAX-WS 2.0 规范是代替JAX-RPC 1.0的下一代Web服务AP…
阅读更多...
YOLOv5 onnx \tensorrt 推理

YOLOv5 onnx \tensorrt 推理

一、yolov5 pt模型转onnx code: https://github.com/ultralytics/yolov5 python export.py --weights yolov5s.pt --include onnx二、onnx 推理 import os import cv2 import numpy as np import onnxruntime import timeCLASSES [person, bicycle, car, motorcycle, airpl…
阅读更多...
stable diffusion简介和原理

stable diffusion简介和原理

Stable Diffusion中文的意思是稳定扩散&#xff0c;本质上是基于AI的图像扩散生成模型。 Stable Diffusion是一个引人注目的深度学习模型&#xff0c;它使用潜在扩散过程来生成图像&#xff0c;允许模型在生成图像时考虑到文本的描述。这个模型的出现引起了广泛的关注和讨论&am…
阅读更多...
IP 地址查询，快速查询自己的 IP 地址

IP 地址查询，快速查询自己的 IP 地址

文章目录 在线结果 在线 http://myip.top/ 结果
阅读更多...
【计算机网络笔记】Cookie技术

【计算机网络笔记】Cookie技术

系列文章目录 什么是计算机网络&#xff1f; 什么是网络协议&#xff1f; 计算机网络的结构 数据交换之电路交换 数据交换之报文交换和分组交换 分组交换 vs 电路交换 计算机网络性能&#xff08;1&#xff09;——速率、带宽、延迟 计算机网络性能&#xff08;2&#xff09;…
阅读更多...
Redis（04）| 数据结构-压缩列表

Redis（04）| 数据结构-压缩列表

压缩列表的最大特点&#xff0c;就是它被设计成一种内存紧凑型的数据结构&#xff0c;占用一块连续的内存空间&#xff0c;不仅可以利用 CPU 缓存&#xff0c;而且会针对不同长度的数据&#xff0c;进行相应编码&#xff0c;这种方法可以有效地节省内存开销。 但是&#xff0c;…
阅读更多...
腾讯云 AI 绘画：文生图、图生图、图审图 快速入门

腾讯云 AI 绘画：文生图、图生图、图审图 快速入门

腾讯云 AI 绘画是腾讯云推出的一款基于人工智能的图像生成和编辑产品&#xff0c;能够根据输入的图片或描述文本&#xff0c;智能生成与输入内容相关的图片&#xff0c;支持多样化的图片风格选择。 在本文中&#xff0c;我们将介绍如何使用腾讯云 AI 绘画的三项主要功能&#…
阅读更多...
Windows-Oracle19c 安装详解-含Navicate远程连接配置 - 同时连接Oracle11g和Oracle19c

Windows-Oracle19c 安装详解-含Navicate远程连接配置 - 同时连接Oracle11g和Oracle19c

文章目录 0 说明1 下载链接2 安装&#xff1a;一定要以管理员身份运行&#xff0c;不然后面有可能会报错。3 启动监听4. 登录Oracle4 Navicate远程连接-配置监听4.1 修改监听文件4.2 网络配置助手-配置本地监听端口4.3 Navicate连接成功 5 Navicate同时连接两个Oracle数据库 0 …
阅读更多...
鸿蒙ArkUI-X跨端应用开发，一套代码构建多平台应用

鸿蒙ArkUI-X跨端应用开发，一套代码构建多平台应用

文章目录 一、项目介绍二、技术架构三、Gitee仓库地址四、ArkUI-X开发者文档五、快速开始——环境准备1、下载DevEco Studio&#xff0c;版本V4.0 Beta2以上2、打开DevEco&#xff0c;下载相关环境配置3、配置开发环境3.1、OpenHarmony SDK3.2、安装ArkUI-X SDK3.2、Android SD…
阅读更多...
postgis ST_CoverageInvalidEdges用法

postgis ST_CoverageInvalidEdges用法

官方文档 概要 geometry ST_CoverageInvalidEdges(geometry winset geom, float8 tolerance 0); 描述 一个窗口函数&#xff0c;用于检查窗口分区中的多边形是否形成有效的多边形覆盖范围。 它返回线性指示器&#xff0c;显示每个多边形中无效边&#xff08;如果有&#x…
阅读更多...
死锁Deadlock

死锁Deadlock

定义 死锁是指两个或多个线程互相持有对方所需的资源&#xff0c;从而导致它们无法继续执行的情况。如下图所示&#xff0c;现有两个线程&#xff0c;分别是线程A及线程B&#xff0c;线程A持有锁A&#xff0c;线程B持有锁B。此时线程A想获取锁B&#xff0c;但锁B需等到线程B的结…
阅读更多...
GCC、g++、gcc的关系

GCC、g++、gcc的关系

GCC、g、gcc的关系 引言 VsCode中对编译环境进行配置的时选择编译器时发现有多种不同的编译器 GNU计划和GCC GNU的全称 GNU’s Not UNIX GNU是一个计划 Q:为什么会有这个计划 因为当时的Unix开始收费和商业闭源,有人觉得不爽→ 想要自己开发和Unix类似的→GNU计划 GUN计划目…
阅读更多...
几个常用的nosql数据库的操作方式

几个常用的nosql数据库的操作方式

dynamoDB 键 partition key&#xff1a;分区键 定义&#xff1a;分区键是用于分布数据存储的主键&#xff0c;每个项&#xff08;Item&#xff09;在表中都必须有一个唯一的分区键值。 特点&#xff1a; 唯一性&#xff1a;每个分区键值在表中必须是唯一的&#xff0c;这是因为…
阅读更多...
shell脚本变量

shell脚本变量

目录 1.变量的定义 2.shell脚本中变量的定义方法 3.变量的转译 4.Linux中命令的别名设定 5.用户环境变量的更改 6.利用命令的执行结果设定变量 7.脚本函数 1.变量的定义 1&#xff09;定义本身 变量就是内存一片区域的地址 2)变量存在的意义 命令无法操作一直变化的目…
阅读更多...
【ETL工具】Datax-ETL-SqlServerToHDFS

【ETL工具】Datax-ETL-SqlServerToHDFS

&#x1f984; 个人主页——&#x1f390;个人主页 &#x1f390;✨&#x1f341; &#x1fa81;&#x1f341;&#x1fa81;&#x1f341;&#x1fa81;&#x1f341;&#x1fa81;&#x1f341; 感谢点赞和关注 &#xff0c;每天进步一点点&#xff01;加油&#xff01;&…
阅读更多...
读高性能MySQL（第4版）笔记21_读后总结与感想兼导读

读高性能MySQL（第4版）笔记21_读后总结与感想兼导读

1. 基本信息 高性能MySQL&#xff1a;经过大规模运维验证的策略&#xff08;第4版&#xff09; High Performance MySQL, Fourth Edition [美] Silvia Botros(西尔维亚博特罗斯)&#xff1b;Jeremy Tinley(杰里米廷利) 电子工业出版社,2022年10月出版 1.1. 读薄率 书籍总字…
阅读更多...
HarmonyOS DevEso环境搭建

HarmonyOS DevEso环境搭建

DevEco Studio 3.1配套支持HarmonyOS 3.1版本及以上的应用及服务开发&#xff0c;提供了代码智能编辑、低代码开发、双向预览等功能&#xff0c;以及轻量构建工具DevEco Hvigor 、本地模拟器&#xff0c;持续提升应用及服务开发效率。 1.下载 官方网站&#xff1a; HUAWEI De…
阅读更多...
Python深度学习实战-基于tensorflow原生代码搭建BP神经网络实现分类任务（附源码和实现效果）

Python深度学习实战-基于tensorflow原生代码搭建BP神经网络实现分类任务（附源码和实现效果）

实现功能 前面两篇文章分别介绍了两种搭建神经网络模型的方法&#xff0c;一种是基于tensorflow的keras框架&#xff0c;另一种是继承父类自定义class类&#xff0c;本篇文章将编写原生代码搭建BP神经网络。 实现代码 import tensorflow as tf from sklearn.datasets import…
阅读更多...
ssh连接远程服务器，并在终端安装anaconda

ssh连接远程服务器，并在终端安装anaconda

官网下载安装&#xff1a;anaconda2023.09版本&#xff08;官网地址&#xff1a;https://www.anaconda.com/download#downloads&#xff09; wget https://repo.anaconda.com/archive/Anaconda3-2023.09-0-Linux-x86_64.sh使用阿里云镜像下载安装&#xff0c;官网下载太慢。阿…
阅读更多...
2698 求一个整数的惩罚数 （子集和，DFS）

2698 求一个整数的惩罚数 （子集和，DFS）

class Solution { public:bool dfs(int target, string s, int index, int sum) {// 只有整个字符串都被分割&#xff0c;求和&#xff0c;和看结果是不是等于targetif(index s.size()) {return sum target;}int num 0; // 在现在的子集中去依次加入余下的元素// 1 2 9 6// …
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023