一种校正CT射线束硬化伪影的方法技术

本发明专利技术公开了一种校正CT射线束硬化伪影的方法，该方法基于X射线能谱分析对原图像进行补偿，消除X射线束的硬化伪影；即，根据已知的扫描条件，获得X射线的能谱分布，计算等效能量；然后，进行扫描，根据建像条件重建原始图像，将原始图像分割成软组织、骨组织和空气部分，分别对不同部分进行能谱投影和等效能量投影；计算等效能量投影与能谱投影的差值；对投影差值进行反投影得到补偿图像，将补偿图像与原始图像相加得到射线束硬化伪影校正后的目标图像。

A method of correcting CT beam hardening artifacts

【技术实现步骤摘要】
一种校正CT射线束硬化伪影的方法
本专利技术涉及一种校正CT射线束硬化伪影的方法，本专利技术属于图像处理

技术介绍
CT是1969年由英国电子工程师亨斯菲尔德.G.N.专利技术，1972年问世的。CT不同于普通X线成像，它是用X射线束对人体层面进行扫描，取得信息，经计算机处理而获得的重建图像，是数字成像而不是模拟成像。CT所显示的断层解剖图像，其密度分辨力明显优于X线成像，使X线成像不能显示的解剖结构及其病变得以显影，从而显著扩大了人体的检查范围，提高了病变检出率和诊断的准确率。CT成像基本原理是用X射线束对人体检查部位一定厚度的层面进行扫描，由探测器接收透过该层面的X射线束，转变为可见光后，由光电转换器转变为电信号，再经模拟/数字转换器转为数字信号，输入计算机进行图像成像处理。在X射线束穿透人体器官或组织时，由于人体器官或组织是由多种不同成分、不同密度的物质构成的，所以，实际上各点对X射线束的吸收系数是不同的，同一种物质对不同能量的光子的衰减系数也不同，而X射线束的能谱是多色的，在X射线穿过人体器官或组织的过程中，低能量的光子通常容易被吸收衰减的快，高能量的光子吸收的慢衰减的慢，因此在射线穿过物质过程中射线的能谱发生变化，射线随着穿过的路径变长，越来越不易吸收，出现了X射线束硬化现象。由于X射线的硬化现象，导致重建的CT图像中出现伪影，这种伪影称为射线束硬化伪影。通常，射线束硬化伪影分两类，一类是：在扫描均匀的组织时，产生的呈“杯状”或“盆状”的伪影；另一类是：扫描多种不同密度物体时，在致密组织之间形成的暗色的条带。这些伪影都极易引起医生误诊，影响医生诊断。通常，为了消除这两类伪影需要进行两步处理，首先，通过采集水模数据模拟人体的软组织，通过将水模图像校正成均匀图像从而确定多项式的系数，利用多项式拟合的方法校正射线束穿过均匀软组织产生的射线束硬化伪影；第二步，将图像中的骨头分割出来，对骨头部分的图像通过再次多项式拟合的方式消除骨头引起的射线束硬化伪影。为消除射线束硬化引起的伪影，需要进行两步校正操作。另外，通过多项式拟合的方式确定多项式系数来源于经验，缺少理论依据，因此，这种校正伪影的方法的鲁棒性比较差。
技术实现思路
为了更好地校正CT射线束硬化伪影，本专利技术的目的是提供一种基于X射线能谱分析补偿的方式校正CT射线束硬化伪影的方法。为实现上述目的，本专利技术采用以下技术方案：一种校正CT射线束硬化伪影的方法，具体步骤如下：S1：获取X光源能谱分布，计算等效能量；S2：根据预设的扫描部位、扫描强度对人体进行X射线束扫描，图像重建，得到待校正的带有射线束硬化伪影的原始图像I0；S3：将原始图像I0分割出软组织部分、骨组织部分和空气部分，得到分割图像I1；S4：计算分割图像I1中X射线穿过软组织和骨组织的路径长度；S5：根据S4计算出的X射线穿过软组织和骨组织的路径长度，进行等效能量投影得到等效能量投影P1；等效能量投影P1计算公式如下：P1＝μb(e)Lb+μs(e)Ls其中：P1为等效能量对应的投影；μb(e)为等效能量值对应的骨头衰减系数，该衰减系数通过公开资料可以查得；Lb为射线穿过骨头的路径长度；μs(e)为等效能量对应的软组织的衰减系数，该衰减系数通过公开资料可以查得；Ls为射线穿过软组织的路径长度；S6：根据S4计算出的射线穿过软组织和骨组织的路径长度，进行X射线束投影得到能谱投影P2；能谱投影P2计算公式如下：其中：P2为射线能谱投影的结果；S(i)e为第i个能量值对应的谱分布；μf(i)为第i个能量值对应的滤过材料的衰减系数，该衰减系数通过公开资料可以查得；Lf为射线穿过滤过材料的路径长度；μb(i)为第i个能量值对应的骨头衰减系数，该衰减系数通过公开资料可以查得；Lb为射线穿过骨头的路径长度；μs(i)为第i个能量值对应的软组织的衰减系数，该衰减系数通过公开资料可以查得；Ls为射线穿过软组织的路径长度；S7：计算等效能量投影P1与能谱投影P2的投影差值：P＝P1-P2；S8：将投影差值进行滤波反投影得到补偿图像I2；S9：将原始图像I0与补偿图像I2相加得到校正后的消除射线束硬化伪影的目标图像I。附图说明图1是本专利技术校正CT射线束硬化伪影的方法流程图；图2是120KV管电压对应的X光源能谱分布图；图3是本专利技术具体实施例中的原始图像；图4是本专利技术具体实施例中的分割图像；图5是本专利技术具体实施例中的补偿图像；图6是本专利技术具体实施例中校正后的消除X射线束硬化伪影后的图像。具体实施方式下面结合附图和实施例对本专利技术的结构及特征进行详细说明。需要说明的是，可以对此处公开的实施例做出各种修改，因此，说明书中公开的实施例不应该视为对本专利技术的限制，而仅是作为实施例的范例，其目的是使本专利技术的特征显而易见。本专利技术校正CT射线束硬化伪影的方法是基于X射线能谱分析对原图像进行补偿，消除X射线束的硬化伪影；即，首先，根据已知的扫描条件，获得X射线的能谱分布，计算等效能量；然后，进行扫描，根据建像条件重建原始图像，将原始图像分割成软组织、骨组织和空气部分，分别对不同部分进行能谱投影和等效能量投影；计算等效能量投影与能谱投影的投影差值；对投影差值进行反投影得到补偿图像，将补偿图像与原始图像相加得到射线束硬化伪影校正后的目标图像。如图1所示，本专利技术校正CT射线束硬化伪影的方法具体步骤如下：S1：获取X光源能谱分布，计算等效能量对于任何一台CT机，其出厂后即可确定其发射X射线的球管的管电压，根据球管管电压即可获取该CT机X光源能谱分布。如图2所示，为120KV管电压对应的X光源能谱分布图。获取CT机X光源能谱后，以摆放在扫描视野中心的直径为20cm的水模作为模拟扫描对象，进行投影计算水模的能谱投影值，将该投影值作为目标投影；再以X光能谱中最高能量的40％-60％为能量，以0.1为步长计算单能量投影，计算X光谱中每个能量对应的水模投影结果；与目标投影均方差最接近的单能量投影对应的能量值即为等效能量。S2：根据预设的扫描部位、扫描强度对人体进行X射线束扫描，图像重建，得到待校正的带有射线束硬化伪影的原始图像I0。图3为对人体头部进行扫描，然后进行图像重建得到的待校正的带有射线束硬化伪影的原始图像。S3：将原始图像I0分割出软组织部分、骨组织部分和空气部分，得到分割图像I1。通过CT值阈值将原始图像I0分割出软组织部分、骨组织部分和空气部分。在本专利技术具体实施例中，如图4所示，将CT值在[-200，100]之间的图像定义为软组织；将CT值大于100的图像定义为骨组织；将CT值小于-200的图像定义为空气。S4：计算分割图像I1中X射线穿过软组织和骨组织的路径长度。根据CT机的几何结构，模拟计算从球管发射出的X光射线到探测器检测单元穿过本文档来自技高网...

【技术保护点】
1.一种校正CT射线束硬化伪影的方法，具体步骤如下：/nS1：获取X光源能谱分布，计算等效能量；/nS2：根据预设的扫描部位、扫描强度对人体进行X射线束扫描，图像重建，得到待校正的带有射线束硬化伪影的原始图像I

【技术特征摘要】
1.一种校正CT射线束硬化伪影的方法，具体步骤如下：
S1：获取X光源能谱分布，计算等效能量；
S2：根据预设的扫描部位、扫描强度对人体进行X射线束扫描，图像重建，得到待校正的带有射线束硬化伪影的原始图像I0；
S3：将原始图像I0分割出软组织部分、骨组织部分和空气部分，得到分割图像I1；
S4：计算分割图像I1中X射线穿过软组织和骨组织的路径长度；
S5：根据S4计算出的X射线穿过软组织和骨组织的路径长度，进行等效能量投影得到等效能量投影P1；
等效能量投影P1计算公式如下：
P1＝μb(e)Lb+μs(e)Ls
其中：P1为等效能量对应的投影；μb(e)为等效能量值对应的骨头衰减系数，该衰减系数通过公开资料可以查得；Lb为射线穿过骨头的路径长度；μs(e)为等效能量对应的软组织的衰减系数，该衰减系数通过公开资料可以查得；Ls为射线穿过软组织的路径长度；
S6：根据S4计算出的射线穿过软组织和骨组织的路径长度，进行X射线束投影得到能谱投影P2；
能谱投影P2计算公式如下：



其中：P2为射线能谱投影的结果；S(i)e为第i个能量值对应的谱分布；μf(i)为第i个能量值对应的滤过材料的衰减系数，该衰减系数通过公开资料可以查得；Lf为射线穿过滤过材料的路径长度；...

【专利技术属性】
技术研发人员：王秀清，任彦，
申请(专利权)人：赛诺威盛科技北京有限公司，
类型：发明
国别省市：北京;11