一种快速实现X辐射图像分离的方法

1、    一种快速实现X辐射图像分离的方法        摘要通过改变入射的X射线能谱，从而得到物体内部不同密度成分的分布像的方法，是X辐射成像领域的一个重要课题。本文提出了一种新的像分离方法，即通过影像灰度进行密度分布像重建，从本质上对这种算法进行了探讨，证明了这种算法的可行性，并且进行了两种成分的像分离实验，取得了较好的结果。这种方法的突出特点是减少了运算量和运算时间，可以快速实现像分离，有一定的实用性。关键词X射线能谱密度解析成像快速分离算法非线性校正A Fast Sep

2、arating Method for X-ray ImagingWang Chunyan(Dept of Biomedical Engineering, Capital University of Medical Sciences, Beijing, 100054)Li Zheng, Kang Kejun, Gao Wenhuan, Liu Yongkang(Dept of Engineering Physics,Tsinghua University, Beijing 100084)AbstractIt is an important issue to get the density dis

3、tributed image inside an object by changing energy spectrum of incident X-ray. In this paper, the authors brought up a new method for image separation, i.e. by reconstruction of density distributed imaging with gray levels of images. We discussed the nature of the algorithm, and proved the feasibili

4、ty of the algorithm. Finally. Some experimental to separate images of two components results are given. From the results we can couclude:the features of this method are: the time of calculation can be reduced a lot, and the calculation is simplified; the images can be separated faster. Hopefully, th

5、is method can be, applied in practice.Key words:X-ray; Energy spectrum; Density analytic imaging; Fast separation algorithm; Non-Linear correction 从50年代提出利用变化X射线能谱实现物体密度解析成像的想法至今，双能乃至多能X射线像分离的成像方法仍是一个较为热门的话题。通过分析不同的像从而得到不同密度的分布仍然具有诱人的前景1，2。本文提出了密度解析成像的快速分离算法，从本质上对这种算法进行了探讨，证明了这种算法的可行性，并且进行了两种成分的像分离实

6、验。从实验结果可以看出，这种算法是完全可行性的。1像分离算法为了研究方便，选取包含两种物质的物体作为研究对象，我们姑且认为这两种物质为医学上常见的骨组织和软组织，厚度分别是tb和ts，骨组织对高、低两种能量的X射线的吸收系数分别是bh和bl，软组织对高、低两种能量的X射线的吸收系数分别是sh和sl，入射的高、低两种能量X射线的强度分别是Ih0和Il0，探测器上得到的高、低两种能量X射线的强度分别是Ih和Il，则有3(1)两边取对数得(2)从(2)式解出：(3)(4)(5)(6)式中bltb为骨组织低能像中的吸收值，在已知两种成分的吸收系数的情况下，对像中的每一象素都用(3)式处理，即可得到一幅

7、骨组织的像。同理，可得到一幅软组织的像。两种成分的厚度分布亦可得到。但是，在实际应用中，往往不知道待测材料的吸收系数，这时，一种只需从像上的灰度分布就可以进行像分离的方法就显得很有意义了。下面给出这种算法的原理。当tb=0时，(2)式变为则(7)同理，当ts=0时，(2)式变为则  (8)软组织像(9)骨组织像(10)(7)、(8)两式说明，软组织对高低两种能量的X射线的吸收系数之比等于tb=0处高、低能量像灰度值之比(HL)。同理，骨组织对高低两种能量的X射线的吸收系数之比等于ts=0处高、低能量像灰度值之比(HL)。因此，在不知道材料吸收系数的情况下，只要物体中存在ts=0(只有

8、骨组织)和tb=0(只有软组织)的点，就可以按(9)、(10)两式分别得到两种物质的分布像。当物体中不存在ts=0(只有骨组织)和tb=0(只有软组织)的点时，可以选ts0(与骨组织相比软组织含量很少)和tb0(与软组织相比骨组织含量很少)的点，用这些点的灰度值对整幅像进行处理，(7)、(8)、(9)、(10)式近似成立，得到近似结果。根据成像系统分析，可以近似导出灰度Gln(I0I)，则H、L分别代表灰度值。在软组织的影像中，任一象素值(灰度)sl.ts等于该点低能量影像灰度L减去该点高能量影像灰度H与(LH)的乘积H.(LH)，其中(LH)为骨组织吸收值远大于软组织吸收值的点，可认为是影像

9、中最亮点低、高能量影像灰度之比，值是对整个影像各象素都相同的常数。对骨组织影像有类似的结果。2能量限制条件由方程组(2)可看出，若方程组可求解的话，必须有(11)因此，为了从物体的两次成像中分离出两种成分的像，我们必须根据两种成分的吸收系数曲线，正确地选择两束X射线的能量范围，所以(11)式给出了像分离的能量限制条件。3算法的物理解释像分离快速算法的内容是对灰度进行加权处理，然后相减。下面，我们从物理原理上对这一算法进行探讨。1中，横轴表示射线透射的物体厚度，纵轴表示影像灰度，直线h、l分别表示高、低两种能量的射线对同一物体产生的灰度。1灰度与透射厚度的关系从X射线物理学我们知道，对于某种确定

10、的物体，其X射线的衰减系数与射线能量有关，所以用不同能量的X射线透射物体，经探测器得到的灰度是不同的，以因子(lh)遍乘H，相当于把直线h以原点为轴转一角度(-)而与直线l重合，则该种物质的象元灰度值在高、低能量影像中调整为相同。当透射路径上有两种物质时，由于这两种物质在高、低能量下的吸收系数的比值不同，以一种物质在不同能量下的吸收系数的比值对灰度作上述处理，实际只是将一种物质的象元灰度值调整为相同，减影后，则只剩另一种物质造成的灰度。下面以骨组织和软组织的减影为例进行说明。我们用Ls、Lb分别表示低能射线同一透射路径上软组织和骨组织的灰度值，用Hs、Hb分别表示高能射线同一透射路径上软组织和

11、骨组织的灰度值。(12)由(2)式，仅对骨组织(13)(14)由此可看出，对于既含软组织，又含骨组织的某一象素，经过上述处理，骨组织引起的灰度完全消除，只剩下软组织引起的灰度。对整幅像的各个象素均作同样处理，就可实现两种成分的像分离。但是我们也看到，在软组织的吸收值表达式中，包含(blbh)，所以，软组织灰度受到了骨组织的影响。这种分离算法的优点是：(1)象元灰度值直接由物质的吸收值给出，符合影像的形成机理和影像检查目的；(2)与其他解方程组的算法相比，该算法仅涉及加、减、乘法运算，减少了运算量和运算时间，可以快速实现像分离，有一定的实用性。4实验验证实验装置为清华大学粒子技术与辐射成像国家专

12、业实验室的FXS60型微焦点X光机，通过改变管电压的方法产生高、低能X射线，实验用管电压为120kV和60kV，管电流为60A，实验材料为铝和有机玻璃，在射线的透射路径上，铝的厚度为1.5mm，有机玻璃的厚度为060mm。表1有机玻璃和铝的分离结果(mm)有机真实值铝真实值有机测算值铝测算值有机玻璃误差铝误差01.51.21.5-1.2051.57.51.7-2.5-0.2101.511.81.6-1.8-0.1151.515.81.5-0.80201.519.61.40.40.1251.522.11.62.9-0.1301.527.81.42.20.1351.530.71.64.3-0.14

13、01.534.51.75.5-0.2451.538.31.46.70.1501.542.61.37.40.2551.545.71.29.30.3601.549.81.310.20.22是有机玻璃的真实值和分离实测值的对比曲线，横轴为透射路径上的有机玻璃厚度，可以看出：随着路径上有机玻璃厚度的增加，其实测值也增加，基本是线性变化。3是铝的真实值和分离实测值的对比曲线，横轴为透射路径上的有机玻璃厚度，可以看出实测值始终在真实值附近，且不受有机玻璃厚度的影响。这说明采用快速像分离算法进行像分离，其像的输入输出基本是线性的，因而是完全可行的。2有机玻璃的测算值与真实值比较曲线真实值测算值3铝的测算值与

14、真实值比较 真实值测算值46给出了运用快速像分离算法的分离效果。实验装置和实验条件同前，实验样品上，脸是用六层5mm厚的有机玻璃做成，眼睛、鼻子、嘴都是用1.5mm厚的铝做成，这个样品虽然比较简单，但大致可以说明快速像分离的效果。4是像分离前的样品的X光照片，5是有机玻璃被减去后的铝影像，6是减去铝影像后的有机玻璃的影像。可以看出这种像分离算法是完全可行的。 4铝和有机玻璃原5铝的减影6有机玻璃减影5结束语本文讨论了利用变化能谱实现X辐射像分离的一般原理和方法，提出了一种新的像分离方法，即通过影像灰度进行密度分布像重建。并且用实测像进行了检验，取得了较好的结果。这种方法的突出特点是减少了运算量

15、和运算时间，可以快速实现像分离，有一定的实用性。基金项目：国家自然科学基金资助项目(19705006).作者简介：王春燕(1958)，女，副教授。作者单位：王春燕首都医科大学生物医学工程系(北京100054)李政康克军高文焕刘永康清华大学工程物理系粒子技术与辐射成像国家专业实验室(北京100084)参考文献1Hiroshi Tsutsui, Tetsuro Ohtsuchi, Koichi Ohmori, et al. X-ray energy separation method using a cdTe semiconductor X-ray imaging sensor and photon counting method. IEEE Tran Nucl Sci,1993，40(1)：402Shoji Kido, Morgan DR, Sones R