复杂点阵镂空结构CT尺寸测量方法，首先通过复杂点阵镂空结构CT成像质量影响因素分析、CT检测策略制定和确定检测参数、CT图像处理技术研究，获得最佳的图像处理方法，然后通过图像边界提取技术、图像边界提取理论计算方法分析和边界提取理论计算方法的研究以及理论算法的验证，获得精确提取复杂点阵镂空结构内部轮廓边界的理论计算方法，最后通过设计尺寸测量误差校准标准试件，应用复杂点阵镂空结构内部轮廓边界的提取算法，对复杂点阵镂空结构CT检测尺寸测量误差进行校准，建立一种复杂点阵镂空结构CT尺寸测量方法，解决复杂点阵镂空结构检测难题。
1.复杂点阵镂空结构CT尺寸测量方法，其特征在于包括如下步骤：步骤1：对复杂点阵镂空结构进行CT扫描，获得复杂点阵镂空结构CT图像；
步骤2：将CT图像进行滤波降噪，逐级统计复杂点阵镂空结构CT图像灰度值的数值，拟合一条开口向上的曲线，取曲线波谷处的灰度值，并记为阈值ε；
步骤3：将复杂点阵镂空结构CT图像采用重叠的方式拆分成多个图像块，并将图像块按从上到下、从左到右的顺序排成行，每一图像块内部像素按从上到下、从左到右的顺序排成列，构成新的图像矩阵φ；
步骤4：新的图像矩阵φ中的每一个元素减去该元素所在列的平均值，实现新的图像矩阵φ的归一化，得到矩阵blocks；
步骤5：根据得到的阈值，利用K-SVD算法更新字典和稀疏模型解算原理，对矩阵blocks进行稀疏编码和影像重构，再将重构的blocks矩阵的每个列向量转换为图像块，从而得到滤波后的CT图像；
步骤6：利用形态学理论开运算制作背景图像，然后将滤波后的CT图像减去背景图像；
步骤7：将CT图像进行边界提取，计算每个连通区域的中心坐标；
步骤8：将各个连通区域以各个连通区域中心为原点，以长半轴、短半轴、长轴旋转角度、中心点X坐标、中心点Y坐标，确定画椭圆所需的6个参数，并绘制椭圆，记录椭圆边线的(x,y)坐标，进而计算焦点坐标；
步骤9：判断图像各像素是否在椭圆内，如果像素在椭圆区域外，则舍去原始图像中的像素，如果在椭圆区域内，则保留原始图像中的像素，遍历图像所有像素；
步骤10：制作复杂点阵镂空结构CT图像边界提取标准试件，采用三坐标测量机对复杂点阵镂空结构CT图像边界提取标准试件进行尺寸测量，获得标定值；
步骤11：对复杂点阵镂空结构CT图像边界提取标准试件进行CT扫描，得到CT图像，然后重复步骤2-步骤10，获得CT图像提取边界的尺寸测量值，与标定值作差，获得尺寸测量误差；
步骤12：对尺寸测量误差进行校准，然后采用校准后的CT图像边界提取方法对点阵镂空结构内部轮廓尺寸进行测量，获得复杂点阵镂空结构内部轮廓尺寸。
2.根据权利要求1所述的复杂点阵镂空结构CT尺寸测量方法，其特征在于：通过采用微焦点CT设备或常规焦点CT设备，对所述的步骤1的复杂点阵镂空结构进行CT扫描。
3.根据权利要求1或2所述的复杂点阵镂空结构CT尺寸测量方法，其特征在于：所述的步骤2中将CT图像进行滤波降噪的方法为采用稀疏表示的图像滤波方法。
4.根据权利要求1或2所述的复杂点阵镂空结构CT尺寸测量方法，其特征在于：所述的步骤3将复杂点阵镂空结构CT图像采用重叠的方式拆分成的多个图像块为8×8大小。
5.根据权利要求1或2所述的复杂点阵镂空结构CT尺寸测量方法，其特征在于：所述的焦点坐标的计算方法为：(81)对于沿顺时针旋转角度较小的椭圆，首先计算椭圆的焦距c式中，a代表椭圆的长半轴，b代表椭圆的短半轴；
(83)对于沿顺时针旋转角度较大的椭圆，计算椭圆的焦距c
6.根据权利要求1或2所述的复杂点阵镂空结构CT尺寸测量方法，其特征在于：所述的判断图像各像素是否在椭圆内的方法为：若图像中像素到椭圆焦点的距离大于长轴长，则在椭圆外，若等于长轴长，则在椭圆上，若小于长轴长，则在椭圆内。
7.根据权利要求1或2所述的复杂点阵镂空结构CT尺寸测量方法，其特征在于：所述的对尺寸测量误差进行校准为采用双校准试块尺寸校准方法。
8.根据权利要求1或2所述的复杂点阵镂空结构CT尺寸测量方法，其特征在于：所述的步骤3中采用重叠的方式拆分成的图像块均代表一个64×1的列向量。
9.根据权利要求1或2所述的复杂点阵镂空结构CT尺寸测量方法，其特征在于：所述的步骤5中将重构的blocks矩阵的每个列向量转换为图像块，从而得到滤波后的CT图像的方法为：将重构的blocks矩阵的每个64×1的列向量转换为8×8的图像块，从而得到滤波后的CT图像。
复杂点阵镂空结构CT尺寸测量方法
技术领域
[0001] 本发明涉及一种复杂点阵镂空结构CT尺寸测量方法，该方法属于CT检测技术尺寸测量技术领域，适用于测量复杂点阵镂空结构节点连接杆直径范围在φ0.5～2.0mm，主要用于测量复杂点阵镂空结构内部轮廓尺寸。
背景技术
[0002] 复杂点阵镂空结构是通过在零件表面、零件内部设计具有一定分布规律的三维点阵曲面替代零件原始的实体区域、支撑结构，实现零部件的形状优化具有高度的空间对称性，可将外部载荷均匀分解，在实现建中的同时保证了承载能力。
[0003] 针对复杂点阵镂空结构检测而言，CT检测技术与其他无损检测技术相比其优势在于：CT检测技术不受试件材料种类、形状结构及表面状况等限制；能给出与复杂点阵结构的几何结构、组分及密度特性相对应的三维图像，且成像直观、目标特征不受周围细节的遮挡；可方便地测量内部轮廓尺寸等。所以基于原理CT检测技术可有效检测复杂点阵镂空结构内部轮廓尺寸，现阶段应用CT检测复杂镂空点阵结构时，对CT图像采用等值面阈值方法(ISO50)提取CT检测图像上关键结构内部轮廓边界，但是CT检测图像边界提取不准确，提取边界混淆、连通，影响尺寸测量精度。其原因在于两个方面：一是由于受点阵结构复杂性的影响，CT检测成像质量比较差，图像存在伪影、信噪比低、图像边界模糊等原因，二是由于受点阵结构表面不线性、不规则形状的影响，造成现用边界分割方法提取边界不准确。提高CT检测复杂点阵镂空结构内部轮廓尺寸测量精度，需解决的两个问题：
[0004] (1)由于复杂点整镂空结构特点造成CT检测成像质量差。(2)对于增材制造结构成形时表面不线性、不规则的特性和CT图像质量不高，不能精确提取边界，造成尺寸测量误差。
[0005] CT检测提高尺寸测量精度的主要途径有两个：一是可通过提高CT检测系统设备硬件水平，二是可改进CT检测工艺方法和理论计算方法提高精度。对硬件的提高是有限的，但通过检测方法的研究和理论计算方法提高测量精度的方法是简单有效的。
发明内容
[0006] 本发明解决的技术问题是：克服现有技术的不足，提供了复杂点阵镂空结构CT尺寸测量方法，解决了复杂点阵镂空结构CT检测图像成像质量差和边界提取精度差的问题，提高复杂点阵镂空结构内部轮廓尺寸测量精度，核心技术包括复杂点阵镂空结构CT检测图像的成像质量控制技术、复杂点阵镂空结构边界提取技术。
[0007] 本发明的技术解决方案是：复杂点阵镂空结构CT尺寸测量方法，包括如下步骤：
[0008] 步骤1：对复杂点阵镂空结构进行CT扫描，获得复杂点阵镂空结构CT图像；
[0009] 步骤2：将CT图像进行滤波降噪，逐级统计复杂点阵镂空结构CT图像灰度值的数值，拟合一条开口向上的曲线，取曲线波谷处的灰度值，并记为阈值ε；
[0010] 步骤3：将复杂点阵镂空结构CT图像采用重叠的方式拆分成多个图像块，并将图像块按从上到下、从左到右的顺序排成行，每一图像块内部像素按从上到下、从左到右的顺序排成列，构成新的图像矩阵φ，其中，每一个64×1的列向量代表一个图像块；
[0011] 步骤4：新的图像矩阵φ中的每一个元素减去该元素所在列的平均值，实现新的图像矩阵φ的归一化，得到矩阵blocks；
[0012] 步骤5：根据得到的阈值，利用K-SVD算法更新字典和稀疏模型解算原理，对矩阵blocks进行稀疏编码和影像重构，再将重构的blocks矩阵的每个64×1的列向量转换为8×8的图像块，从而得到滤波后的CT图像；
[0013] 步骤6：利用形态学理论开运算制作背景图像，然后将滤波后的CT图像减去背景图像；
[0014] 步骤7：将CT图像进行边界提取，计算每个连通区域的中心坐标；
[0015] 步骤8：将各个连通区域以各个连通区域中心为原点，以形状先验信息为参数画椭圆，记录椭圆边线的(x,y)坐标，进而计算焦点坐标；
[0016] 步骤9：判断图像各像素是否在椭圆内，如果像素在椭圆区域外，则舍去原始图像中的像素，如果在椭圆区域内，则保留原始图像中的像素，遍历图像所有像素；
[0017] 步骤10：制作复杂点阵镂空结构CT图像边界提取标准试件，采用三坐标测量机对复杂点阵镂空结构CT图像边界提取标准试件进行尺寸测量，获得标定值；
[0018] 步骤11：对复杂点阵镂空结构CT图像边界提取标准试件进行CT扫描，得到CT图像，然后重复步骤2-步骤10，获得CT图像提取边界的尺寸测量值，与标定值作差，获得尺寸测量误差；
[0019] 步骤12：对尺寸测量误差进行校准，然后采用校准后的CT图像边界提取方法对点阵镂空结构内部轮廓尺寸进行测量，获得复杂点阵镂空结构内部轮廓尺寸。
[0020] 所述的对步骤1的复杂点阵镂空结构进行CT扫描为通过采用微焦点CT或常规焦点CT设备。
[0021] 所述的步骤2中将CT图像进行滤波降噪的方法为采用稀疏表示的图像滤波方法。
[0022] 所述的步骤3将复杂点阵镂空结构CT图像采用重叠的方式拆分成的多个图像块为8×8大小。
[0023] 所述的形状先验信息为包括长半轴、短半轴、长轴旋转角度。
[0024] 所述的焦点坐标的计算方法为：
[0025] (81)对于沿顺时针旋转角度较小的椭圆，首先计算椭圆的焦距c
[0026]
[0027] 式中，a代表椭圆的长半轴，b代表椭圆的短半轴；
[0028] (82)求解焦点坐标包括
[0029]
[0030] 式中，α代表椭圆长轴与x轴的交角；
[0031] (83)对于沿顺时针旋转角度较大的椭圆，计算椭圆的焦距c
[0032]
[0033] 式中，a代表椭圆的长半轴，b代表椭圆的短半轴；
[0034] (84)求解焦点坐标包括
[0035] 所述的判断图像各像素是否在椭圆内的方法为：
[0036] 若图像中像素到椭圆焦点的距离大于长轴长，则在椭圆外，若等于长轴长，则在椭圆上，若小于长轴长，则在椭圆内。
[0037] 所述的对尺寸测量误差进行校准为采用双校准试块尺寸校准方法。
[0038] 本发明与现有技术相比的优点在于：
[0039] (1)本发明采用CT测量复杂点阵镂空结构内部轮廓尺寸，进行了误差校准，内部轮廓尺寸测量精度高；
[0040] (2)本发明提供了复杂点阵镂空结构CT检测图像的处理方法，可获得最佳的复杂点阵镂空结构CT检测图像，保证了尺寸测量精度；
[0041] (3)本发明提供了复杂点阵镂空结构CT检测图像边界提取模型，可精确提取复杂点阵镂空结构的边界，保证了尺寸测量精度；
[0042] (4)本发明可以检测复杂点阵镂空结构内部缺陷，如：节点连接杆断裂变形、内部困粉等问题，可以保证复杂点阵镂空结构内部质量。
附图说明
[0043] 图1为复杂点阵镂空结构节点连接杆单元设计图；
[0044] 图2为复杂点阵镂空结构CT图像边界提取校准试件。
[0045] 图3为稀疏表示的图像滤波方法进行滤波后的CT图像；
[0046] 图4为基于稀疏表示的复杂点阵镂空结构CT图像滤波流程图；
[0047] 图5为形态学理论开运算制作的背景图像；
[0048] 图6为伪影去除后的CT图像；
[0049] 图7为求取连通区域中心示意图；
[0050] 图8为画椭圆并求出焦点坐标示意图；
[0051] 图9为求解沿顺时针旋转角度较小的椭圆焦点坐标示意图；
[0052] 图10为求解沿顺时针旋转角度较大的椭圆焦点坐标示意图；
[0053] 图11为复杂点阵镂空结构边界提取实验结果示意图；
具体实施方式
[0054] 本发明针对现有技术的不足，提出一种复杂点阵镂空结构CT尺寸测量方法，1)通过复杂点阵镂空结构CT成像质量影响因素分析、CT检测策略制定和确定检测参数、CT图像处理技术研究，获得最佳的图像处理方法。2)通过图像边界提取技术、图像边界提取理论计算方法分析和边界提取理论计算方法的研究以及理论算法的验证，获得精确提取复杂点阵镂空结构内部轮廓边界的理论计算方法。3)通过设计尺寸测量误差校准标准试件，应用复杂点阵镂空结构内部轮廓边界的提取算法，对复杂点阵镂空结构CT检测尺寸测量误差进行校准。4)建立一种复杂点阵镂空结构CT尺寸测量方法，解决复杂点阵镂空结构检测难题。下面结合附图对本发明进行详细的解释和说明，本发明包括如下步骤。
[0055] 步骤1：采用复杂点阵镂空结构作为CT尺寸测量载体。图1所示的点阵节点连接杆单元以阵列的方式设计的结构。
[0056] 步骤2：采用采用复杂点阵镂空结构边界提取标准试件作为校准的标准物质，试件圆柱尺寸为φ2mm、φ3mm、φ4mm、φ5mm，制作试件时应注意选择与被检测复杂点阵镂空结构一致的材料，以保证射线衰减系数保持一致，规避在CT检测图像上出现由于材料不同带来的伪像。图2所示的复杂点阵镂空结构CT图像边界提取校准试件设计图。然后采用三坐标测量机对复杂点阵镂空结构CT图像边界提取标准试件进行尺寸测量，获得准确可靠的标定值。
[0057] 步骤3：采用微焦点CT或常规焦点CT设备对步骤1的复杂点阵结构进行CT扫描，获得复杂点阵镂空结构CT图像。
[0058] 步骤4：将步骤3获得的CT图像采用稀疏表示的图像滤波方法进行滤波降噪。逐级统计复杂点阵镂空结构CT图像灰度值的数值，拟合一条开口向上的曲线，取曲线的波谷处的灰度值为阈值ε。
[0059] 步骤5：将复杂点阵镂空结构CT图像采用重叠的方式拆分成8×8大小的图像块，并将图像块按从上到下，从左到右的顺序排成行，每一图像块内部像素按从上到下，从左到右的顺序排成列，构成新的图像矩阵φ，其中每一个64×1的列向量代表一个图像块。
[0060] 步骤6：将步骤5获得的新图像矩阵φ中的每一个元素减去该元素所在列的平均值，实现矩阵φ的归一化，从而得到新的矩阵blocks；
[0061] 步骤7：根据第4步骤得到的阈值，利用K-SVD算法更新字典和稀疏模型解算原理，对矩阵blocks进行稀疏编码和影像重构，再将重构的blocks矩阵的每个64×1的列向量转换为8×8的图像块，从而得到滤波后的CT图像，如图3所示的采用稀疏滤波方法滤波后的图像。如图4所示基于稀疏表示的复杂点阵镂空结构CT图像滤波流程图
[0062] 步骤8：利用形态学理论开运算制作背景图像，如图5所示的形态学理论开运算制作的背景图像。然后用步骤7获得的原始图像减去背景图像，去除复杂点阵镂空结构CT图像伪影，如图6所示的伪影去除后的CT图像。
[0063] 步骤9：将步骤8获得的CT图像进行形状先验的复杂点阵镂空结构CT图像边界提取。求连通区域并标记，计算每个连通区域的中心坐标，如图7所示求取连通区域中心示意图。
[0064] 步骤10：将步骤9所获得的连通区域以各连通区域中心为原点，以形状先验信息为参数画椭圆。以长半轴、短半轴、长轴旋转角度(顺时针)、中心点X坐标、中心点Y坐标、画椭圆所需的点数6个参数绘制椭圆(其中前三个参数是形状先验信息)，并记录椭圆边线的(x,y)坐标，如图8所示画椭圆并求出焦点坐标示意图。由于图像中蝶状区域是有旋转角度的椭圆，固需要分两种情况计算焦点坐标。对于沿顺时针旋转角度较小的椭圆为第一种情况，对于沿顺时针旋转角度较大的椭圆为第二种情况。
[0065] 第一种情况(见图9):
[0066] 1)对于沿顺时针旋转角度较小的椭圆，首先计算椭圆的焦距c
[0067]
[0068] 式中，a代表椭圆的长半轴，b代表椭圆的短半轴,c代表椭圆的焦距。
[0069] 2)求解1号焦点坐标：
[0070]
[0071] 式中a代表椭圆长轴与x轴的交角，x轴、y轴可以为CT图像中的任一直角坐标系。
[0072] 3)求解2号焦点坐标：
[0073]
[0074] 式中，α代表椭圆长轴与x轴的交角；
[0075] 第二种情况(见图10)：
[0076] 1)对于沿顺时针旋转角度较大的椭圆，首先计算椭圆的焦距c
[0077]
[0078] 式中，a代表椭圆的长半轴，b代表椭圆的短半轴。
[0079] 2)求解1号焦点坐标：
[0080]
[0081] 式中a代表椭圆长轴与x轴的交角。
[0082] 3)求解2号焦点坐标：
[0083]
[0084] 步骤11：判断图像各像素是否在椭圆内，求解原图像上的椭圆区域。椭圆有一个重要特征即是椭圆上任一点到两个焦点的距离等于长轴长。可根据这个特性，遍历图像中的像素，图像中任一像素到椭圆焦点的距离大于长轴长，则在椭圆外；等于长轴长，则在椭圆上；小于长轴长，则在椭圆内。如果像素在椭圆区域外，则舍去原始图像中的像素，继续往下遍历，如果在椭圆区域内，则保留原始图像中的像素，直至所有像素都遍历完为止。复杂点阵镂空结构边界提取实验结果见图11。
[0085] 步骤12：将步骤2的复杂点阵镂空结构边界提取标准试件进行CT扫描。获得CT图像后，执行步骤4-11的流程进行CT图像边界提取。对提取的CT图像进行尺寸测量，获得CT图像尺寸测量值。
[0086] 步骤13：将步骤12获得的CT图像测量值与步骤2获得的三坐标机测量的标定值作差，获得尺寸测量误差。
[0087] 步骤14：如果尺寸测量误差不大于±0.1mm则复杂点阵镂空结构尺寸测量完毕。如果尺寸测量误差大于±0.1mm则进行14步骤，进行尺寸误差校准。
[0088] 步骤14：采用双校准试块尺寸校准方法对步骤13获得的尺寸测量误差进行校准后，获得复杂点阵镂空结构尺寸测量值，则测量完毕。
[0089] 本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
