发明名称 一种复杂凸面型物体三维重建与体积测量方法
摘要 本发明公开了一种复杂凸面型物体三维重建与体积测量方法。该方法将待测物体放置于布设有对角标志的平板上,并设置一个足够大的长方体形外包体,通过采集不同拍摄角度的待测物体图像,根据图像中提取的对角标志进行测量相机的获取图像的投影矩阵,然后利用待测物体的外轮廓和投影矩阵对初始的外包体进行切割。在不同方位进行足够多次切割之后即可获得待测物体的空间三维模型以及体积。与现有采用结构光测量技术相比,本发明的技术成本更低、操作简便灵活,能处理主动视觉无法处理的反光、透光、纹理缺乏的情况。
申请公布号 CN104154877B 申请公布日期 2016.07.06
申请号 CN201410444920.2 申请日期 2014.09.03
申请人 中国人民解放军国防科学技术大学 发明人 刘海波;姚林伸;陈圣义;张小虎;于起峰;杨玉坤;楚少臻;陈相印
分类号 G01B11/24(2006.01)I;G01B11/25(2006.01)I;G06T17/00(2006.01)I 主分类号 G01B11/24(2006.01)I
代理机构 湖南省国防科技工业局专利中心 43102 代理人 冯青
主权项 一种凸面型物体三维重建与体积测量方法,其特征在于,包括以下步骤:第一步,测量系统布设;第二步,初始的外包体数据集合<img file="dest_path_image002.GIF" wi="25" he="19" />设置;第三步,待测物体三维形貌重建;第四步,待测物体体积计算;所述第一步具体为:1.1 将待测物体放置在于黑色平板中心位置;1.2 建立测量坐标系o‑xyz坐标系o‑xyz中o点位于平板中心, o‑x轴和o‑y轴位于平板上表面所在平面;o‑z轴垂直o‑xy平面,指向平板上方;o‑x轴、o‑y轴和o‑z轴符合右手定则;1.3布设编码标志在平板上布设1个圆形对角标志和<img file="dest_path_image004.GIF" wi="19" he="19" />个方形对角标志,其中<img file="338329dest_path_image004.GIF" wi="19" he="19" />≥3,上述对角标志在o‑xyz坐标系中的三维坐标已知,且均匀分布于以待测物体为中心的圆周上;所述第二步具体为:2.1 设置一个被测物体的外包长方体设置一个被测物体的外包长方体,该长方体的顶点坐标分别为(‑X/2,‑Y/2,0)、(‑X/2,Y/2,0)、(X/2,‑Y/2,0)、(X/2,Y/2,0)、(‑X/2,‑Y/2,Z)、(‑X/2,Y/2,Z)、(X/2,‑Y/2,Z)、(X/2,Y/2,Z);2.2分割外包长方体并记录数据将外包长方体分割成<img file="dest_path_image006.GIF" wi="14" he="16" />个小的立方体单元,<img file="439491dest_path_image006.GIF" wi="14" he="16" />≥ 6000000,将立方体单元的中心坐标<img file="dest_path_image008.GIF" wi="181" he="30" />及其对应的标识码<img file="dest_path_image010.GIF" wi="26" he="24" />,记入外包体数据集合<img file="593129dest_path_image002.GIF" wi="25" he="19" />,<img file="115902dest_path_image008.GIF" wi="181" he="30" />的计算公式分别为,<img file="dest_path_image012.GIF" wi="169" he="46" />(1)<img file="dest_path_image014.GIF" wi="169" he="46" />(2)<img file="dest_path_image016.GIF" wi="129" he="46" />(3)其中,<img file="dest_path_image018.GIF" wi="75" he="48" />,<img file="dest_path_image020.GIF" wi="97" he="19" />,<img file="dest_path_image022.GIF" wi="107" he="42" />,<img file="dest_path_image024.GIF" wi="106" he="42" />,<img file="dest_path_image026.GIF" wi="107" he="42" />;2.3初始化立方体单元标识码令立方体单元对应的标识码;所述第三步具体为:3.1 待测物体图像拍摄选取拍摄位置和角度,利用测量相机拍摄一张待测物体的图像,保证所拍摄的图像包含平板上的对角标志;3.2待测物体图像轮廓提取利用已有的图像边缘检测算法,提取待测物体图像的封闭轮廓,将位于待测物体图像封闭轮廓以内及待测物体图像封闭轮廓上的图像点坐标记为集合<img file="dest_path_image028.GIF" wi="26" he="18" />,<img file="dest_path_image030.GIF" wi="113" he="35" />,其中<img file="dest_path_image032.GIF" wi="48" he="27" />为位于待测物体图像封闭轮廓以内及待测物体图像封闭轮廓上图像点像素坐标,<img file="dest_path_image034.GIF" wi="21" he="25" />为<img file="508705dest_path_image032.GIF" wi="48" he="27" />的个数;3.3 对角标志识别与定位3.3.1 对待测物体图像的封闭轮廓以外的图像区域,依次进行二值化、图像闭运算、膨胀处理,以及连通域分割,将得到的各连通域记为集合<img file="dest_path_image036.GIF" wi="50" he="30" />;其中<img file="dest_path_image038.GIF" wi="22" he="25" />为各连通域图像点像素坐标构成的子集,并按元素个数降序排列,<img file="dest_path_image040.GIF" wi="24" he="25" />为连通域的个数;3.3.2 计算集合<img file="832632dest_path_image036.GIF" wi="50" he="30" />中前<img file="dest_path_image042.GIF" wi="41" he="19" />个子集的面积<img file="dest_path_image044.GIF" wi="19" he="21" />、周长<img file="dest_path_image046.GIF" wi="17" he="25" />以及质心坐标<img file="dest_path_image048.GIF" wi="72" he="27" />,其中<img file="dest_path_image050.GIF" wi="137" he="24" />;3.3.3 在以<img file="782221dest_path_image048.GIF" wi="72" he="27" />为中心建立的一个正方形搜索区域内,采用Harris角点检测算法计算获得对角标志的亚像素图像坐标<img file="dest_path_image052.GIF" wi="74" he="27" />;3.3.4 计算<img file="dest_path_image054.GIF" wi="90" he="51" />,<img file="dest_path_image056.GIF" wi="29" he="25" />最大值所对应的对角标志,即判定为圆形对角标志,其它<img file="dest_path_image058.GIF" wi="19" he="19" />个即判定为方形对角标志,以圆形对角标志的亚像素图像坐标为起点,在图像中顺时针排序,完成所有对角标志的识别;3.4 投影矩阵计算3.4.1 利用识别出的对角标志的亚像素图像坐标<img file="544903dest_path_image052.GIF" wi="74" he="27" />,及其对应的o‑xyz坐标系中的三维坐标,利用位姿估计算法,获得测量相机在o‑xyz坐标系中的姿态矩阵<img file="dest_path_image060.GIF" wi="17" he="18" />和平移向量<img file="dest_path_image062.GIF" wi="10" he="15" />;3.4.2 利用公式(4),计算出投影矩阵<img file="dest_path_image064.GIF" wi="168" he="75" />(4);3.5外包体数据集合<img file="771704dest_path_image002.GIF" wi="25" he="19" />中立方体单元对应的标识码重置;3.5.1 如果编号<img file="dest_path_image066.GIF" wi="14" he="18" />的立方体单元对应的标识码<img file="dest_path_image068.GIF" wi="49" he="25" />,其中<img file="dest_path_image070.GIF" wi="83" he="19" />,利用公式(5)计算该立方体单元中心在测量相机成像面上的投影的像面坐标<img file="dest_path_image072.GIF" wi="51" he="26" />:<img file="dest_path_image074.GIF" wi="139" he="54" />(5)其中,<img file="dest_path_image076.GIF" wi="337" he="34" />,<img file="dest_path_image078.GIF" wi="67" he="26" />表示对<img file="dest_path_image080.GIF" wi="17" he="17" />四舍五入取整,上标<img file="dest_path_image082.GIF" wi="11" he="13" />表示向量或矩阵的转置;3.5.2 判断<img file="952019dest_path_image072.GIF" wi="51" he="26" />是否属于集合<img file="830982dest_path_image028.GIF" wi="26" he="18" />,若<img file="692628dest_path_image072.GIF" wi="51" he="26" />不属于集合<img file="340647dest_path_image028.GIF" wi="26" he="18" />,令<img file="dest_path_image084.GIF" wi="50" he="25" />;若<img file="40266dest_path_image072.GIF" wi="51" he="26" />属于集合<img file="406526dest_path_image028.GIF" wi="26" he="18" />,令<img file="806283dest_path_image068.GIF" wi="49" he="25" />;3.6 重复步骤3.1至3.5实现被测物体三维重建设置测量相机的拍摄位置,获取不同拍摄角度的待测物体图像m张,m一般取30~40,完成被测物体三维形貌的重建;所述第四步具体为:计算待测物体的体积<img file="dest_path_image086.GIF" wi="17" he="19" />,计算公式为,<img file="dest_path_image088.GIF" wi="89" he="46" />(6)。
地址 410073 湖南省长沙市砚瓦池正街47号