基于CT序列图像的旋转体积测量方法.pdfVIP

第25卷第9期 强 激 光 与 粒 子 束 V01．25，No．9 2013年9月 HIGHPOWERLASERANDPARTICLEBEAMS Sep．，2013 文章编号：1001—4322(2013)09—2423—06 基于CT序列图像的旋转体积测量方法” 邹永宁1， 伍立芬2， 王 珏1， 段伟伟2 (1．重庆大学光电技术与系统教育部重点实验室ICT研究中心，重庆400044； 2．重庆大学自动化学院，重庆400044) 摘要： 针对工程应用中常见的轴向变化剧烈工件，研究一种Zernike矩亚像素层间拟合轮廓旋转积分 的高精度体积测量算法。首先用Zernike矩方法提取亚像素级轮廓，其次采用多项式拟合切片轮廓，并进行等 相角间隔采样；然后对等相位的轮廓点采用多项式拟合并通过外插法自动预测工件顶端，将各拟合轮廓旋转一 个相角，对旋转块积分并累加得到待测工件的体积。分别对体积相等的3个试块进行CT扫描测量实验，结果 表明该算法的测量误差低于0．7％；对实际含内腔工件进行测量实验，结果表明该算法测量精度是Canny台体 法的1．5倍。 关键词：CT；体积测量；亚像素；旋转积分；内腔体积 中圈分类号：TP391 文献标志码： A doi：10。3788／HPLP2423 实现工业CT图像的高精度定量化分析和测量是工业CT技术发展的方向之一。由于受工业CT系统层 厚的影响，断层图像在垂直方向分辨率较低；针对轴向轮廓非线性变化的工件，现有的三维体积测量方法测量 精度不高，无法满足实际需要，因此，研究高精度CT图像体积测量方法十分必要。CT图像体积测量方法有体 数据直接测量和基于切片图像测量两种；体数据直接测量的方法有3DFacet模型[21法：用一个三元三次多项 式来拟合当前体素的三维邻域灰度，把二阶方向导数的零交叉点作为物体的三位的三维边缘。基于二维切片 图像的测量方法包括灰度插值法、形状插值法[3-83、数学形态学方法[gqz]等；基于形状插值的方法运算繁冗复 杂，抗噪性较差；基于数学形态学的方法边缘定位精度只有像素级别的精度。ChenPing[1朝等提出了一种CT 图像精确尺寸的测量方法，虽然精度较高，但只能测量物体的长度不能测量体积。卿杰[61等提出的基于层间轮 廓拟合的形状插值方法只对简单柱形工件有效。王珏[11等提出的台体累加方法，解决了工件缺顶问题，当待测 工件轴向非线性变化很大时精度不高。本文针对大多数工件为旋转体这一特点，研究了轴向轮廓强烈非线性 工件的三维高精度测量算法，首先获得各层切片图像亚像素轮廓点并进行多项式拟合，计算二维轮廓点的曲率 半径，对轮廓点按照等相位准则进行层间匹配与拟合，将拟合曲线旋转一定角度，最终得到工件的精确体积。 1算法原理 1．1 Zernike矩边缘提取 Zernike矩是正交矩，它具有平移、缩放和旋转不变性，常用于数字图像的特征检测、分类和识别n4。18|。 置作为实际边缘的位置。首先定义图像f(p，口)的挖阶m次Zernike矩为 1厂(1D，O)R√ID)eMp@dO’ (1) k=掣r 7c J0 J0 m 式中：72，m为正整数，是正交矩的阶数，并且满足胛≥0，竹一，7z为偶数，且I1≤竹；V。。(JD，目)一R。。e，枷是正交 复数多项式，R—cP，孝示径向多项式，定义为Rnm cp一‘,,-蚤ImlⅢi_五王娄主}≥亨产墨兰蠡。 在单位圆中建立理想的阶跃边缘模型，圆内边缘两侧的灰度值分别为h和^+愚，点为灰度差，z是圆心到边缘 的垂直距离，9是边缘与z轴的夹角。将边缘顺时针旋转妒角后，使其与y轴平行，正交复数多项式V。。，V。。和V。。 *收稿日期：2012-12—24；修订日期：2013-04—22 基金项目：重庆市科技攻关项目(CSTC2009AC3