三维扫描的精度验证方法二——拟合球体

已于 2022-03-09 13:53:03 修改
阅读量1k 收藏 6
点赞数
分类专栏： 三维扫描 文章标签： 深度学习 cnn 人工智能
于 2021-12-18 16:21:20 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/sunnyrainflower/article/details/122013660
版权
三维扫描 专栏收录该内容
18 篇文章 200 订阅 ¥19.90 ¥99.00
订阅专栏

借助于CloudCompare软件

1、用三维扫描仪扫描一个标准球体，已知球体的直径，用CloudCompare打开

2、找到CloudCompare软件中的拟合球体选项，计算球体的直径，然后与实际的尺寸做比对，下图为得到的拟合的球体

下图为操作步骤

大胡子大叔
关注
  • 0
    点赞
  • 踩
  • 6
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 2
    评论
专栏目录
订阅专栏
点云拟合球体
qq_40041064的博客
07-26 702
在很多情况下，需要根据点云来拟合球体，本博文主要介绍各种方法的拟合情况及优缺点，希望对各位小伙伴有所帮助！
三维扫描的精度验证方法之一——平面拟合（平面度）
大胡子大叔的专栏
12-18 2375
需要借助halcon软件 1、三维扫描仪扫描一个特别平的平面，如标定板、量块……等 2、将生成的三维点云数据对齐三维的轴线（三维点云的基本操作，不做详解，网上很多资料） 3、对齐后，查看点云在另一个方向上的波动情况，是否收敛于某个区域 下图为点云的平面对齐三维的轴的效果图 4、如何查看：在halcon的变量窗口，找到轴的数据，在变量上右击，点击“像函数一样绘图”，即可查看波浪图像，效果图如下图 注：需要注意点云的单位，别搞错了 本文为原创文章，未经同意请勿转载，谢谢配合 ..
Halcon3D学习笔记
十层楼
03-12 4822
gen_object_model_3d_from_points 可以将离散的坐标数据转换为3D模型，数据来源可以是CSV等文件 write_object_model_3d 保存3D模型到我本地文件，格式多样 read_object_model_3d 可以将保存到本地的3D模型读取出来 select_points_object_model_3d 可以选择离散点数据范围，用于分割高度或分...
Halcon 用点来拟合平面
weixin_39354845的博客
06-12 5734
项目要求用多个点来拟合一个平面，然后再用其他平面上的点来计算这个点到平面的距离，halcon 有现成的拟合函数。
python实现三维拟合的方法
09-19
今天小编就为大家分享一篇python实现三维拟合的方法，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
基于三维激光扫描仪测距精度的分析
06-14
三维激光扫描仪是传统测量的技术革新,研究三维激光扫描仪基本...在20 m范围内,点位中误差低于6 mm,满足三维激光扫描仪的标称精度。扫描仪距离与中误差之间存在二次函数关系,随着扫描距离增加中误差增大,测距精度降低。
task.zip_ployfit_三维_三维插值_二维插值_拟合公式
07-15
本案例中，"task.zip_ployfit_三维_三维插值_二维插值_拟合公式"是一个项目，它涉及到使用ployfit进行三维数据的拟合，并进行三维和二维的插值操作。以下是关于这些知识点的详细解释： 1. **ployfit函数**：ployfit...
python-三维拟合的方法.zip
11-21
本主题将深入探讨如何利用Python实现三维拟合的方法，包括理论基础、使用库以及实际应用示例。 首先，我们需要了解什么是三维拟合。在统计学和数学中，三维拟合是指寻找一个三维函数来最佳地描述给定的三维数据点集...
Python实现二维曲线拟合的方法
01-01
如下所示： from numpy import * import numpy as np import matplotlib.pyplot as plt plt.close() fig=plt.figure() plt.grid(True) plt.axis([0,10,0,8]) #列出数据 point=[[1,2],[2,3],[3,6],[4,7],[6,5],[7,3]...
直线、平面、球体的拟合算法
06-29
本文是将若干离散点，通过最小二乘法，常用曲面的数学算法。
空间离散点球体最小二乘拟合SphereFit
04-06
空间离散点球体最小二乘拟合MFC源代码
快速球体拟合 [非迭代]：使用精确解法使用几何拟合来拟合 3D 数据集。 算法快速且非迭代-matlab开发
05-30
使用几何拟合的精确解拟合3D数据。 有关算法的详细信息，请参阅并引用“Sumith YD，“使用精确解的球体快速几何拟合算法”，2015，arXiv:1506.02776”。 该算法适用于球体、半球和球体的任何部分。
球面拟合-基于最小二乘法
04-12
在三维空间中对当前数据集的散点进行球体拟合获得球体描述，球体中心坐标，球体半径。
球体拟合（最小二乘法）：将球体拟合到一组噪声数据。 不需要宽弧或许多点。-matlab开发
05-30
给定一组数据点，该函数以最小二乘法计算数据的中心和半径。 最小二乘方程用于将矩阵求逆为 3x3，而不是直接在数据集上进行。 不需要大弧或许多数据点。 假设点不是奇异的（共面的）和真实的...... 在 R2010b 上创建，但应该适用于所有版本。
最小二乘法拟合球及其相关代码实现
woniu199166的博客
01-19 6453
当我们手中握有大量的数据时，对于二维的数据，我们会对他们进行直线拟合、对数拟合，圆曲线的拟合等等。这些拟合的方法都是运用的了非常古老而又非常有效的方法，即最小二乘法。 今天给大家介绍一种三维球面数据的拟合方法，该方法也是运用的最小二乘的方法。旨在使拟合的半径在均方意义下误差达到最小。公式推导设拟合后的球面的球心为(x_0,y_0,z_0)及半径r。 对于每一点拟合后估计的值与实际值的差值为:  则
CloudCompare——实现点云由粗到精的配准
热门推荐
点云侠的博客
08-28 2万+
CloudCompare——实现点云由粗到精的配准
PCL RANSAC拟合球体（C++详细过程版）
点云侠的博客
03-20 1033
C++详细过程版实现的RANSAC拟合球体。
球面数据拟合算法简介
Inspire
12-05 1万+
当我们手中握有大量的数据时，对于二维的数据，我们会对他们进行直线拟合、对数拟合，圆曲线的拟合等等。这些拟合的方法都是运用的了非常古老而又非常有效的方法，即最小二乘法。 今天给大家介绍一种三维球面数据的拟合方法，该方法也是运用的最小二乘的方法。旨在使拟合的半径在均方意义下误差达到最小。
ransac球体拟合python
z504727099的博客
04-11 706
下面展示一些 内联代码片。 // A code block var foo = 'bar'; camera_cx = 637.166 camera_cy = 369.465 camera_fx = 611.223; camera_fy = 610.982; import random import numpy as np def estimate_sphere(data): x1=data[0][0] y1=data[0][1] z1=data[0][2] x2 = d
c++最小二乘球体拟合
最新发布
01-02
最小二乘球体拟合是一种用于拟合三维数据点到球体模型的算法。它通过最小化数据点到拟合球体的距离平方和来确定最佳的球体参数。以下是一个使用C++实现最小二乘球体拟合的示例代码： ```cpp #include <iostream> #include <vector> #include <cmath> struct Point3D { double x; double y; double z; }; struct Sphere { double centerX; double centerY; double centerZ; double radius; }; Sphere fitSphere(const std::vector<Point3D>& points) { double sumX = 0.0; double sumY = 0.0; double sumZ = 0.0; double sumX2 = 0.0; double sumY2 = 0.0; double sumZ2 = 0.0; double sumXY = 0.0; double sumXZ = 0.0; double sumYZ = 0.0; double sumX3 = 0.0; double sumY3 = 0.0; double sumZ3 = 0.0; double sumXY2 = 0.0; double sumXZ2 = 0.0; double sumYZ2 = 0.0; for (const auto& point : points) { double x = point.x; double y = point.y; double z = point.z; sumX += x; sumY += y; sumZ += z; sumX2 += x * x; sumY2 += y * y; sumZ2 += z * z; sumXY += x * y; sumXZ += x * z; sumYZ += y * z; sumX3 += x * x * x; sumY3 += y * y * y; sumZ3 += z * z * z; sumXY2 += x * y * y; sumXZ2 += x * z * z; sumYZ2 += y * z * z; } double N = points.size(); double D11 = sumX2 + sumY2 + sumZ2; double D12 = sumX * sumY + sumX * sumZ + sumY * sumZ; double D13 = sumX * sumY2 + sumX * sumZ2 + sumY * sumZ2; double D14 = sumX2 * sumY + sumX2 * sumZ + sumY2 * sumZ; double D21 = D12; double D22 = sumX2 + sumY2 + sumZ2; double D23 = sumX2 * sumY + sumX2 * sumZ + sumY2 * sumZ; double D24 = sumX * sumY2 + sumX * sumZ2 + sumY * sumZ2; double D31 = D13; double D32 = D23; double D33 = sumX2 + sumY2 + sumZ2; double D34 = sumX * sumY + sumX * sumZ + sumY * sumZ; double D41 = D14; double D42 = D24; double D43 = D34; double D44 = N; double D = D11 * (D22 * D33 * D44 + D23 * D34 * D42 + D24 * D32 * D43 - D24 * D33 * D42 - D23 * D32 * D44 - D22 * D34 * D43) - D12 * (D21 * D33 * D44 + D23 * D34 * D41 + D24 * D31 * D43 - D24 * D33 * D41 - D23 * D31 * D44 - D21 * D34 * D43) + D13 * (D21 * D32 * D44 + D22 * D34 * D41 + D24 * D31 * D42 - D24 * D32 * D41 - D22 * D31 * D44 - D21 * D34 * D42) - D14 * (D21 * D32 * D43 + D22 * D33 * D41 + D23 * D31 * D42 - D23 * D32 * D41 - D22 * D31 * D43 - D21 * D33 * D42); double Dx = D11 * (D22 * D33 * sumX + D23 * D34 * sumY + D24 * D32 * sumZ - D24 * D33 * sumY - D23 * D32 * sumZ - D22 * D34 * sumX) - D12 * (D21 * D33 * sumX + D23 * D34 * sumZ + D24 * D31 * sumY - D24 * D33 * sumZ - D23 * D31 * sumY - D21 * D34 * sumX) + D13 * (D21 * D32 * sumX + D22 * D34 * sumZ + D24 * D31 * sumY - D24 * D32 * sumZ - D22 * D31 * sumY - D21 * D34 * sumX) - D14 * (D21 * D32 * sumY + D22 * D33 * sumX + D23 * D31 * sumZ - D23 * D32 * sumX - D22 * D31 * sumZ - D21 * D33 * sumY); double Dy = D11 * (D22 * D33 * sumY + D23 * D34 * sumZ + D24 * D32 * sumX - D24 * D33 * sumZ - D23 * D32 * sumX - D22 * D34 * sumY) - D12 * (D21 * D33 * sumY + D23 * D34 * sumX + D24 * D31 * sumZ - D24 * D33 * sumX - D23 * D31 * sumZ - D21 * D34 * sumY) + D13 * (D21 * D32 * sumY + D22 * D34 * sumX + D24 * D31 * sumZ - D24 * D32 * sumX - D22 * D31 * sumZ - D21 * D34 * sumY) - D14 * (D21 * D32 * sumZ + D22 * D33 * sumY + D23 * D31 * sumX - D23 * D32 * sumY - D22 * D31 * sumX - D21 * D33 * sumZ); double Dz = D11 * (D22 * D33 * sumZ + D23 * D34 * sumX + D24 * D32 * sumY - D24 * D33 * sumX - D23 * D32 * sumY - D22 * D34 * sumZ) - D12 * (D21 * D33 * sumZ + D23 * D34 * sumY + D24 * D31 * sumX - D24 * D33 * sumY - D23 * D31 * sumX - D21 * D34 * sumZ) + D13 * (D21 * D32 * sumZ + D22 * D34 * sumY + D24 * D31 * sumX - D24 * D32 * sumY - D22 * D31 * sumX - D21 * D34 * sumZ) - D14 * (D21 * D32 * sumY + D22 * D33 * sumZ + D23 * D31 * sumX - D23 * D32 * sumZ - D22 * D31 * sumX - D21 * D33 * sumY); double Dxx = D11 * (D22 * D33 * sumX2 + D23 * D34 * sumXY + D24 * D32 * sumXZ - D24 * D33 * sumXY - D23 * D32 * sumXZ - D22 * D34 * sumX2) - D12 * (D21 * D33 * sumX2 + D23 * D34 * sumXZ + D24 * D31 * sumXY - D24 * D33 * sumXZ - D23 * D31 * sumXY - D21 * D34 * sumX2) + D13 * (D21 * D32 * sumX2 + D22 * D34 * sumXZ + D24 * D31 * sumXY - D24 * D32 * sumXZ - D22 * D31 * sumXY - D21 * D34 * sumX2) - D14 * (D21 * D32 * sumXY + D22 * D33 * sumX2 + D23 * D31 * sumXZ - D23 * D32 * sumX2 - D22 * D31 * sumXZ - D21 * D33 * sumXY); double Dyy = D11 * (D22 * D33 * sumY2 + D23 * D34 * sumYZ + D24 * D32 * sumY2 - D24 * D33 * sumYZ - D23 * D32 * sumY2 - D22 * D34 * sumY2) - D12 * (D21 * D33 * sumY2 + D23 * D34 * sumY2 + D24 * D31 * sumYZ - D24 * D33 * sumY2 - D23 * D31 * sumYZ - D21 * D34 * sumY2) + D13 * (D21 * D32 * sumY2 + D22 * D34 * sumY2 + D24 * D31 * sumYZ - D24 * D32 * sumY2 - D22 * D31 * sumYZ - D21 * D34 * sumY2) - D14 * (D21 * D32 * sumYZ + D22 * D33 * sumY2 + D23 * D31 * sumY2 - D23 * D32 * sumY2 - D22 * D31 * sumY2 - D21 * D33 * sumYZ); double Dzz = D11 * (D22 * D33 * sumZ2 + D23 * D34 * sumZ2 + D24 * D32 * sumYZ - D24 * D33 * sumZ2 - D23 * D32 * sumYZ - D22 * D34 * sumZ2) - D12 * (D21 * D33 * sumZ2 + D23 * D34 * sumYZ + D24 * D31 * sumZ2 - D24 * D33 * sumYZ - D23 * D31 * sumZ2 - D21 * D34 * sumZ2) + D13 * (D21 * D32 * sumZ2 + D22 * D34 * sumYZ + D24 * D31 * sumZ2 - D24 * D32 * sumYZ - D22 * D31 * sumZ2 - D21 * D34 * sumZ2) - D14 * (D21 * D32 * sumYZ + D22 * D33 * sumZ2 + D23 * D31 * sumZ2 - D23 * D32 * sumZ2 - D22 * D31 * sumZ2 - D21 * D33 * sumYZ); double Dxy = D11 * (D22 * D33 * sumXY + D23 * D34 * sumY2 + D24 * D32 * sumYZ - D24 * D33 * sumY2 - D23 * D32 * sumYZ - D22 * D34 * sumXY) - D12 * (D21 * D33 * sumXY + D23 * D34 * sumYZ + D24 * D31 * sumY2 - D24 * D33 * sumYZ - D23 * D31 * sumY2 - D21 * D34 * sumXY) + D13 * (D21 * D32 * sumXY + D22 * D34 * sumYZ + D24 * D31 * sumY2 - D24 * D32 * sumYZ - D22 * D31 * sumY2 - D21 * D34 * sumXY) - D14 * (D21 * D32 * sumY2 + D22 * D33 * sumXY + D23 * D31 * sumYZ - D23 * D32 * sumXY - D22 * D31 * sumYZ - D21 * D33 * sumY2); double Dxz = D11 * (D22 * D33 * sumXZ + D23 * D34 * sumYZ + D24 * D32 * sumZ2 - D24 * D33 * sumYZ - D23 * D32 * sumZ2 - D22 * D34 * sumXZ) - D12 * (D21 * D33 * sumXZ + D23 * D34 * sumZ2 + D24 * D31 * sumYZ - D24 * D33 * sumZ2 - D23 * D31 * sumYZ - D21 * D34 * sumXZ) + D13 * (D21 * D32 * sumXZ + D22 * D34 * sumZ2 + D24 * D31 * sumYZ - D24 * D32 * sumZ2 - D22 * D31 * sumYZ - D21 * D34 * sumXZ) - D14 * (D21 * D32 * sumYZ + D22 * D33 * sumXZ + D23 * D31 * sumZ2 - D23 * D32 * sumXZ - D22 * D31 * sumZ2 - D21 * D33 * sumYZ); double Dyz = D11 * (D22 * D33 * sumYZ + D23 * D34 * sumZ2 + D24 * D32 * sumY2 - D24 * D33 * sumZ2 - D23 * D32 * sumY2 - D22 * D34 * sumYZ) - D12 * (D21 * D33 * sumYZ + D23 * D34 * sum

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

大胡子大叔

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值