一种基于双视角点云配准的类圆形水果纵横径测量方法

1.本发明涉及一种类圆形水果纵横径测量方法，尤其是涉及了一种基于双视角点云配准的类圆形水果纵横径测量方法。


背景技术：

2.我国是世界上的水果产量大国，近年水果总产量达到25688.35万吨，其中常见类圆形水果苹果、柑橘及甜瓜总产量分别达到了4242.54万吨、4584.54万吨及1355.70万吨。对类圆形水果进行品质分级有利于提升其商品化品质，促进类圆形水果在流通中按质论价、优质优价。类圆形水果的纵径及横径是果形的重要参数之一(gb/t 12947
‑
2008鲜柑橘、gb/t 10651
‑
2008鲜苹果)，同时也是品质分级的重要依据。与传统人工测量耗时耗力相比较，采用机器视觉技术通过双视角点云配准的方法测量类圆形水果的纵横径，具有非接触且快速无损等优势，且有重要的实际应用价值。
3.精确的类圆形水果纵横径尺寸的获得，需要对类圆形水果进行纵剖与横剖，但是采用剖切的方法一方面会对类圆形水果造成严重损伤，同时也无法用于类圆形水果的分级生产过程中。
4.利用二维图像测量水果纵横径外形的常见方法主要有最小外接法及边界点测量法：章文英等(2001)(章文英,应义斌.苹果图象的低层处理及尺寸检测[j].浙江农业学报,2001(04):38
‑
41.)利用链符法搜索并细化轮廓线后，利用最小外接矩形求苹果的纵径与横径，结果显示测量尺寸与实际尺寸相关系数达到了0.955；魏新华等(2011)(沈宝国,魏新华,尹建军.基于最小外接圆法的苹果直径检测技术[j].农机化研究,2011,33(12):131
‑
134.)对12个样本苹果采集了144幅图像，利用最小外接圆法测量苹果直径，拟合方程的相关系数达到了0.988；龚昌来等(2013)(龚昌来,罗聪,黄杰贤,杨冬涛.基于机器视觉的柚子尺寸检测系统的研究[j].农机化研究,2013,35(11):22
‑
25.)利用机器视觉获得目标图像的上、下、左、右4个边界点，通过边界点来估算柚子的纵径与横径。基于二维图像测量水果纵横径的方法能够快速实现水果纵横径的计算，但是其对于表面具有较强反光特性的类圆形水果光照环境设置等提出了较高的要求。
[0005]
近年来，随着三维重建技术以及消费级图像采集设备的发展，有学者利用生成水果全表面点云的方式进行三维重建测量水果尺寸或体积：yawe等(2020)(yawei,w,yifeic.fruitmorphologicalmeasurement based on three
‑
dimensional reconstruction.agronomy，2020,10,455)将梨置于旋转台上，利用9张图片得到了梨的三维点云图，并计算了其三维尺寸。但是该方法存在匹配点计算量大，点云数量多等困难，耗时长，不适于生产现场。此外，蜜柚等表面较为光滑的类圆形水果无明显特征点可用来进行信息匹配。


技术实现要素：

[0006]
为了解决背景技术中存在的问题和需求，本发明提供了一种基于双视角点云配准
的类圆形水果纵横径测量方法。
[0007]
本发明的技术方案如下：
[0008]
本发明包括以下步骤：
[0009]
1)搭建双视角采集系统：所述双视角采集系统包括置物平台、类圆形水果、两个rgbd相机、同步线、数据传输线、标定板和主机；
[0010]
两个rgbd相机和主机均放置于置物平台上，两个rgbd相机之间用同步线连接，两个rgbd相机分别通过数据传输线与主机相连，类圆形水果或者标定板放置于置物平台上，类圆形水果或者标定板放置在两个rgbd相机和主机的一侧；两个rgbd相机的光轴均朝向类圆形水果或者标定板，并且两个rgbd相机的光轴成夹角布置；
[0011]
2)建立世界坐标系：记其中一个rgbd相机为主相机c1，另一个rgbd相机为副相机c2，以主相机c1的光心为世界坐标系的原点，以主相机c1的光轴向外为世界坐标系的z轴，世界坐标系的x轴与y轴构成的平面平行于主相机c1的成像平面，x轴与y轴垂直；
[0012]
3)标定主相机c1、副相机c2位姿：将标定板置于置物平台上，利用双视角采集系统对标定板进行采集，获取主相机c1与副相机c2各自采集的标定板彩色图像和对应的标定板深度图像，通过双目标定方法对主相机c1与副相机c2各自采集的标定板彩色图像进行处理后，获得将副相机c2坐标转换至主相机c1坐标的初始转换矩阵t
12
，由主相机c1与副相机c2各自采集的标定板彩色图像和对应的标定板深度图像分别组成主相机标定板点云和副相机标定板点云；
[0013]
4)主相机、副相机标定板点云粗配准并配对，获得粗配准标定板点云的所有正确匹配点对；
[0014]
5)优化初始转换矩阵t
12
；
[0015]
6)重复步骤5)迭代优化转换矩阵t
12
，获得优化转换矩阵
[0016]
7)获取类圆形水果的主、副点云：将类圆形水果置于置物平台上，控制两个rgbd相机同步采集类圆形水果的主点云与副点云；
[0017]
8)类圆形水果的主、副点云精配准：利用步骤6)得到的优化转换矩阵将类圆形水果的主点云与副点云配准得到类圆形水果的精配准水果点云；
[0018]
9)精配准水果点云下采样：使用voxelgrid滤波器对精配准水果点云进行下采样，获得精配准下采样水果点云；
[0019]
10)精配准下采样水果点云滤波：设置点云世界坐标系的x轴的取值范围为(x
min
，x
max
)，y轴的取值范围为(y
min
，y
max
)，z轴的取值范围为(z
min
，z
max
)，根据x轴、y轴和z轴的取值范围利用conditionand算法除去精配准下采样水果点云的外部噪声点云，得到精配准滤波水果点云；
[0020]
11)精配准滤波水果点云离群点去除：取精配准滤波水果点云中的任意一空间点k，获得当前空间点k相邻的50个空间点，计算出当前空间点k与相邻的50个空间点的平均距离，利用高斯分布选取1.5倍平均距离范围外的空间点作为当前空间点k的离群点并进行去除，遍历精配准滤波水果点云中的所有空间点，得到去离群点精配准水果点云；
[0021]
12)纵横径测量：采用有向最小包围盒obb算法对去离群点精配准水果点云进行处理后绘制出类圆形水果的最小obb包围盒，提取最小obb包围盒轴向与类圆形水果纵径方向及横径方向一致的纵径及横径，得到类圆形水果的纵径h及横径w。
[0022]
所述步骤4)具体为：
[0023]
4.1)利用初始转换矩阵t
12
将副相机标定板点云与主相机标定板点云进行配准得到粗配准标定板点云；
[0024]
4.2)在粗配准点云中，任取主相机标定板点云中的一个空间点，在副相机标定板点云中寻找与当前空间点的间距小于误差阈值d的空间点，如果存在，则两个空间点作为一个正确匹配点对；遍历副相机标定板点云，寻找当前空间点的对应的所有正确匹配点对；如果不存在，则当前空间点没有正确匹配点对；
[0025]
4.3)遍历主相机标定板点云中的所有空间点，获得粗配准标定板点云的所有正确匹配点对。
[0026]
所述步骤5)具体为：
[0027]
5.1)计算主相机、副相机点集合的质心：将粗配准标定板点云的所有正确匹配点中的主相机标定板点云和副相机标定板点云分别记为主相机点集合p和副相机点集合q，分别满足p＝(p1，p2，
…
，p
i
，
…
，p
n
，)和q＝(q1，q2，
…
，q
i
，
…
，q
n
)，其中p
i
表示主相机点集合p中第i个空间点，q
i
表示副相机点集合q中第i个空间点，n表示粗配准标定板点云中正确匹配点对的对数，主相机点集合p和副相机点集合q的质心q分别为：
[0028][0029][0030]
5.2)计算主相机、副相机点集合中以质心为原点的空间点坐标：主相机点集合p中以质心p为原点的空间点坐标p
′
i
以及副相机点集合q中以质心q为原点的空间点坐标q
′
i
分别为：
[0031]
p
′
i
＝p
i
‑
p
[0032]
q
′
i
＝q
i
‑
q
[0033]
5.3)计算第i对正确匹配点对的距离误差e
i
，满足e
i
＝|p
i
‑
(r
·
q
i
+t)|，其中r表示当前次迭代的旋转矩阵，t表示当前次迭代的平移矩阵，则所有正确匹配点对的距离误差平方和e(r，t)为：
[0034][0035]
5.4)对所有正确匹配点对的距离误差平方和进行svd分解后获得优化转换矩阵t
′
12
。
[0036]
所述步骤6)具体为：
[0037]
不断重复步骤5)进行迭代，直至达到以下迭代终止条件之一时停止迭代；迭代终止条件1：达到最大迭代次数，迭代终止条件2：所有正确匹配点对的距离误差平方和e(r，t)小于给定的终止误差值；得到优化转换矩阵
[0038]
所述初始转换矩阵t
12
由初始旋转矩阵r
12
与初始平移矩阵t
12
组成。
[0039]
本发明的有益效果为：
[0040]
本发明对表面具有较强反光特性的类圆形水果光照环境要求低，能够避免对于类圆形水果的损伤；同时通过迭代得到的优化转换矩阵进行点云精配准，有利于解决多视角点云配准过程中匹配点计算量大，点云数量多、耗时长等困难，能够做到省时省力，实现了类圆形水果的纵横径测量，有利于类圆形水果的在线生产。
附图说明
[0041]
图1是本发明的整体流程图。
[0042]
图2是本发明的双视角采集系统布置图。
[0043]
图3是本发明标定过程的标定板粗配准点云图。
[0044]
图4是本发明标定过程的标定板精配准点云图。
[0045]
图5是本发明类圆形水果的主点云图。
[0046]
图6是本发明类圆形水果的副点云图。
[0047]
图7是本发明类圆形水果的精配准滤波水果点云图。
[0048]
图8是本发明类圆形水果的去离群点精配准水果点云图。
[0049]
图9是本发明类圆形水果的最小obb包围盒图。
[0050]
图中：1、置物平台，2、类圆形水果，3、rgbd相机，4、同步线，5、数据传输线，6、主机。
具体实施方式
[0051]
下面结合附图和实施例对本发明作进一步说明。
[0052]
本发明选取类圆形水果蜜柚实施例如下：
[0053]
如图1所示，本发明包括以下步骤：
[0054]
1)如图2所示，搭建双视角采集系统：双视角采集系统包括置物平台(1)、类圆形水果(2)、两个rgbd相机(3)、同步线(4)、数据传输线(5)、标定板和主机(6)；
[0055]
两个rgbd相机(3)和主机(6)均放置于置物平台(1)上，两个rgbd相机(3)之间用同步线连接，用于控制两个rgbd相机(3)同步采集点云，两个rgbd相机(3)分别通过数据传输线(5)与主机(6)相连，类圆形水果(2)或者标定板放置于置物平台(1)上，类圆形水果(2)或者标定板放置在两个rgbd相机(3)和主机(6)的一侧；两个rgbd相机(3)的光轴均朝向类圆形水果(2)或者标定板，并且两个rgbd相机(3)的光轴成夹角布置；
[0056]
2)建立世界坐标系：记其中一个rgbd相机(3)为主相机c1，另一个rgbd相机(3)为副相机c2，以主相机c1的光心为世界坐标系的原点，以主相机c1的光轴向外为世界坐标系的z轴，世界坐标系的x轴与y轴构成的平面平行于主相机c1的成像平面，水平向右为x轴，竖直向下为y轴；
[0057]
3)标定主相机c1、副相机c2位姿：将标定板置于置物平台(1)上，利用双视角采集系统对标定板进行采集，获取主相机c1与副相机c2各自采集的标定板彩色图像和对应的标定板深度图像，通过双目标定方法对主相机c1与副相机c2各自采集的标定板彩色图像进行处理后，获得将副相机c2坐标转换至主相机c1坐标的初始转换矩阵t
12
，初始转换矩阵t
12
由初始旋转矩阵r
12
与初始平移矩阵t
12
组成。由主相机c1与副相机c2各自采集的标定板彩色图像
和对应的标定板深度图像分别组成主相机标定板点云和副相机标定板点云；
[0058]
初始转换矩阵t
12
、初始旋转矩阵r
12
与初始平移矩阵t
12
分别为：
[0059][0060][0061]
t
12
＝[
‑
0.5935 0.00787 0.05856]
[0062]
4)主相机、副相机标定板点云粗配准并配对，获得粗配准标定板点云的所有正确匹配点对；
[0063]
步骤4)具体为：
[0064]
4.1)利用初始转换矩阵t
12
将副相机标定板点云与主相机标定板点云进行配准得到粗配准标定板点云，如图3所示；
[0065]
4.2)在粗配准点云中，任取主相机标定板点云中的一个空间点，在副相机标定板点云中寻找与当前空间点的间距小于误差阈值d的空间点，如果存在，则两个空间点作为一个正确匹配点对；遍历副相机标定板点云，寻找当前空间点的对应的所有正确匹配点对；如果不存在，则当前空间点没有正确匹配点对；
[0066]
4.3)遍历主相机标定板点云中的所有空间点，获得粗配准标定板点云的所有正确匹配点对。
[0067]
5)优化初始转换矩阵t
12
；
[0068]
步骤5)具体为：
[0069]
5.1)计算主相机、副相机点集合的质心：将粗配准标定板点云的所有正确匹配点中的主相机标定板点云和副相机标定板点云分别记为主相机点集合p和副相机点集合q，分别满足p＝(p1，p2，
…
，p
i
，
…
，p
n
，)和q＝(q1，q2，
…
，q
i
，
…
，q
n
)，其中p
i
表示主相机点集合p中第i个空间点，q
i
表示副相机点集合q中第i个空间点，n表示粗配准标定板点云中正确匹配点对的对数，主相机点集合p和副相机点集合q的质心q分别为：
[0070][0071][0072]
5.2)计算主相机、副相机点集合中以质心为原点的空间点坐标：主相机点集合p中以质心p为原点的空间点坐标p
′
i
以及副相机点集合q中以质心q为原点的空间点坐标q
′
i
分别为：
[0073]
p
′
i
＝p
i
‑
p
[0074]
q
′
i
＝q
i
‑
q
[0075]
5.3)计算第i对正确匹配点对的距离误差e
i
，满足e
i
＝|p
i
‑
(r
·
q
i
+t)|，其中r表示当前次迭代的旋转矩阵，t表示当前次迭代的平移矩阵，则所有正确匹配点对的距离误差平方和e(r，t)为：
[0076][0077]
5.4)对所有正确匹配点对的距离误差平方和进行svd分解后获得优化转换矩阵t
′
12
。
[0078]
6)重复步骤5)迭代优化转换矩阵t
12
，获得优化转换矩阵
[0079]
步骤6)具体为：
[0080]
不断重复步骤5)进行迭代，直至达到以下迭代终止条件之一时停止迭代；迭代终止条件1：达到最大迭代次数，迭代终止条件2：所有正确匹配点对的距离误差平方和e(r，t)小于给定的终止误差值；得到优化转换矩阵优化转换矩阵由优化后的旋转矩阵r
*
和优化后的平移矩阵t
*
组成；将副相机标定板点云与主相机标定板点云进行配准得到的精配准标定板点云，如图4所示；
[0081]
优化转换矩阵优化后的旋转矩阵r
*
和优化后的平移矩阵t
*
分别为：
[0082][0083][0084]
t
*
＝[
‑
0.73497 0.016825
ꢀ‑
0.03818]
[0085]
7)获取类圆形水果的主、副点云：将类圆形水果(2)置于置物平台(1)上，控制两个rgbd相机(3)同步采集类圆形水果(2)的主点云与副点云，如图5和6所示；
[0086]
8)类圆形水果的主、副点云精配准：利用步骤6)得到的优化转换矩阵将类圆形水果(2)的主点云与副点云配准得到类圆形水果(2)的精配准水果点云；
[0087]
9)精配准水果点云下采样：使用voxelgrid滤波器对精配准水果点云进行下采样，以减少点云数量，获得精配准下采样水果点云；
[0088]
10)精配准下采样水果点云滤波：设置点云世界坐标系的x轴的取值范围为(x
min
,x
max
)，y轴的取值范围为(y
min
,y
max
)，z轴的取值范围为(z
min
,z
max
)，根据x轴、y轴和z轴的取值范围利用conditionand算法除去精配准下采样水果点云的外部噪声点云，得到精配准滤波水果点云，如图7所示；
[0089]
11)精配准滤波水果点云离群点去除：取精配准滤波水果点云中的任意一空间点k，获得当前空间点k相邻的50个空间点，计算出当前空间点k与相邻的50个空间点的平均距
离，利用高斯分布选取1.5倍平均距离范围外的空间点作为当前空间点k的离群点并进行去除，遍历精配准滤波水果点云中的所有空间点，得到去离群点精配准水果点云，如图8所示；
[0090]
12)纵横径测量：采用有向最小包围盒obb算法对去离群点精配准水果点云进行处理后绘制出类圆形水果的最小obb包围盒，如图9所示，提取最小obb包围盒轴向与类圆形水果纵径方向及横径方向一致的纵径及横径，得到类圆形水果的纵径h及横径w。
[0091]
表1列出的是根据本实施例选取直径65mm圆球连续进行10次测量获得的直径数据，其中圆球真值取10次游标卡尺测量均值为65.24mm。
[0092]
表1 65mm圆球直径测量数据表
[0093][0094]