基于分段硬化曲线的ＣＴ射束硬化校正方法技术

本发明专利技术公开了一种基于分段硬化曲线的ＣＴ射束硬化校正方法，获取ＣＴ系统对被检测物体扫描的硬化数据，建立拟合坐标系，选择不含常数项的多项式ｆ（ｘ）对该数据进行最小二乘法拟合，得射束硬化曲线；求ｆ（ｘ）＝０、ｆ′（ｘ）＝０和ｆ″（ｘ）＝０的实数解；将多项式的拟合区间Ｌ均分，分别对区间（０，Ｌ／３）和（２Ｌ／３，Ｌ）中的多项式曲线进行振荡性判断，若有振荡，将其振荡部分替换成幂函数曲线；若无振荡，满足要求；确定过原点的校正直线斜率；对各幅多色投影图像的每个像素，将其灰度代入硬化曲线方程计算得到对应的穿越长度；将该穿越长度代入校正直线方程计算，得到该像素射束硬化校正后的近似单色投影灰度。本发明专利技术增强了硬化曲线的稳定性，方法简单，便于实现，可有效提高射束硬化校正的可靠性。

【技术实现步骤摘要】

【技术保护点】
基于分段硬化曲线的ＣＴ射束硬化校正方法，其特征在于包括下述步骤：　（ａ）获取所使用的ＣＴ系统对被检测物体进行扫描的一组硬化数据，其中每一射线穿越长度对应一个多色投影灰度；　（ｂ）建立以多色投影灰度为Ｘ轴、射线穿越长度为Ｙ轴的拟合坐标系，选择４～７次不含常数项的多项式ｆ（ｘ）＝ａ↓［１］ｘ＋ａ↓［２］ｘ↑［２］＋…＋ａ↓［ｎ］ｘ↑［ｎ］对这组硬化数据进行最小二乘法拟合，所得的多项式为一条过原点的曲线，即射束硬化曲线；其中，ａ↓［１］～ａ↓［ｎ］为多项式系数，ｆ（ｘ）为射线穿越长度，ｘ为多色投影灰度；　（ｃ）根据射束硬化曲线位于第一象限、呈凹状且单调递增的特点，对上一步所得的多项式ｆ（ｘ），分别求出ｆ（ｘ）＝０、ｆ′（ｘ）＝０和ｆ″（ｘ）＝０的实数解；　（ｄ）将多项式的拟合区间Ｌ均分为三段，从左到右分别设为（０，Ｌ／３）、［Ｌ／３，２Ｌ／３］和（２Ｌ／３，Ｌ），分别对区间（０，Ｌ／３）和（２Ｌ／３，Ｌ）中的多项式曲线进行振荡性判断，若区间内有上一步所得的实数解，则认为该区间的多项式曲线有振荡，需将其振荡部分替换成曲线性质较稳定的幂函数曲线ｇ（ｘ）＝ａ.ｘ↑［ｂ］；若区间没有上一步所得的实数解，则认为该区间的多项式曲线无振荡，满足要求；最终的硬化曲线由未被替换的多项式曲线与幂函数曲线替换段联合组成；　（ｅ）确定过原点的校正直线斜率：若多项式曲线在区间（０，Ｌ／３）无振荡，则校正直线的斜率直接取为多项式的一次项系数；若多项式曲线在区间（０，Ｌ／３）有振荡而被替换，由于幂函数在原点的切线斜率为０，此时对位于区间（０，Ｌ／３）左端１／４～１／２的硬化数据进行不含常数项的最小二乘法直线拟合，所得直线的斜率即为所需斜率；　（ｆ）对各幅多色投影图像的每个像素，将其灰度代入硬化曲线方程计算得到对应的穿越长度；　（ｇ）将该穿越长度代入校正直线方程计算，得到该像素射束硬化校正后的近似单色投影灰度。...

【技术特征摘要】

【专利技术属性】
技术研发人员：黄魁东，张定华，卜昆，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：87[中国|西安]