CN111476856A - 一种多谱ct成像方法 - Google Patents
一种多谱ct成像方法 Download PDFInfo
- Publication number
- CN111476856A CN111476856A CN202010271896.2A CN202010271896A CN111476856A CN 111476856 A CN111476856 A CN 111476856A CN 202010271896 A CN202010271896 A CN 202010271896A CN 111476856 A CN111476856 A CN 111476856A
- Authority
- CN
- China
- Prior art keywords
- reference component
- correction parameter
- correction
- hardening artifact
- different energies
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Radiology & Medical Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种多谱CT成像方法,该方法包括:基于参考组分对成像物体进行硬化伪影重建校正得到不同能量下成像物体硬化伪影校正后的重建图像衰减系数;基于参考组分得到成像物体在不同能量下各组分的衰减系数;根据成像物体在不同能量下各组分的衰减系数和不同能量下成像物体硬化伪影校正后的重建图像衰减系数,通过DCM算法进行组分表征,得到成像物体各组分体积表征图像。采用本发明能够减少所需成像物体信息和噪声,获得较好的组分表征结果。
Description
技术领域
本发明涉及CT应用技术领域,特别涉及一种多谱CT成像方法。
背景技术
现有的多谱电子计算机断层扫描(CT)成像方法,有投影域分解、迭代过程分解、图像域分解等方法,投影域分解的方法使用非线性模型将多能谱投影分解为各基材料的投影,之后将各投影使用常规CT重建方法进行重建就可得到各基材料的重建图像。迭代过程分解为基于多能谱下X射线的非线性衰减模型,对多能谱投影进行重建过程时通过迭代修正策略求解与能谱无关或近似无关的图像信息,投影域分解的方法要考虑能谱信息,能谱估计的仪器和时间损耗相对较大,并且能谱波动无法计算在内,迭代过程分解的不足是计算效率问题。图像域分解方法为先对多能谱投影进行CT重建,然后对重建图像进行线性组合可得到各基材料的分解图像。现有技术常规CT仪器发出的X射线为连续谱,所以在对投影数据进行重建时会产生硬化伪影,这也是本文要解决的问题。
使用图像域方法进行材料组分表征过程中,硬化伪影校正效果决定了最终组分表征的正确性。硬化伪影校正方法可分为硬件校正和软件校正。在硬件校正方面,有多种方法。如使用滤波片减小能谱宽度,具有良好单色性的同步辐射光源,可以通过设置多个电子学系统阈值滤除低能量信息的光子计数探测器。软件校正方法,有多项式拟合法、蒙特卡罗校正法、双能校正法等。多项式拟合法是利用多项式拟合的方法得到材质厚度和多色投影之间的关系,并将其校正为线性关系,进而校正硬化伪影;蒙特卡罗校正(Monte Carlo,MC)法是借助MC模拟软件模拟来得到不同物质在宽能谱下的衰减曲线,进而利用非线性函数法将多能投影校正为单能投影数据;双能校正法的思路是对衰减系数按照基材料进行分解,分解成空间依赖分量与能量依赖分量,然后重建基材料图像。
1)利用滤波片虽然可以减小能谱宽度,但是会降低投影信噪比,导致重建图像噪声大,分解误差增大。同步辐射光源具有良好的单色性,但是其造价昂贵,目前难以普及。光子计数探测器可以通过设置多个电子学系统阈值滤除低能量信息,但是因为在采集过程中探测的窄谱通道光子数目减少,信噪比较低,此外,光子计数探测器收集透射光子过程中易受脉冲堆积、电荷共享和非相干散射等影响,加大了窄谱投影数据的噪声,成像效果不佳,实际应用受到限制。
2)多项式拟合法需要与成像物体相同材质的体模,蒙特卡罗校正法的模拟所需时间较长,使用双能校正法需要两个射线源和两个探测器或两次扫描,对硬件要求较高,这制约了其实际应用。
由上述可以确定,现有技术存在的问题是:硬化伪影校正方式中的硬件校正因为引入噪声或成本问题导致实际应用受限,软件校正受限于需要获取额外信息或仿真时间较长。
发明内容
有鉴于此,本发明提供一种多谱CT成像方法,能够减少所需成像物体信息和噪声,获得较好的组分表征结果。
本发明提供了一种多谱CT成像方法,该方法包括:
基于参考组分对成像物体进行硬化伪影重建校正得到不同能量下成像物体硬化伪影校正后的重建图像衰减系数;
基于参考组分得到成像物体在不同能量下各组分的衰减系数;
根据成像物体在不同能量下各组分的衰减系数和不同能量下成像物体硬化伪影校正后的重建图像衰减系数,通过数据约束模型(DCM)算法进行组分表征,得到成像物体各组分体积表征图像。
由上述的技术方案可见,本发明提供了一种实现多谱CT成像方法,基于参考组分对成像物体进行硬化伪影重建校正,并获取各组分衰减系数,最后使用DCM算法获得各组分体积结构表征图像。相比于硬件校正和软件校正硬化伪影,引入参考组分可以减少所需成像物体信息和噪声,并且参考组分信息可以获取成像物体各组分的衰减系数,获得较好的组分表征结果。
附图说明
图1为本发明提供的一种多谱CT成像方法的流程示意图。
图2为硬化伪影校正参数α,β划分示意图。
图3为重建图像模体直径上的衰减系数曲线示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面结合附图和具体实施例对本发明进行详细描述。
本发明的目的在于提供一种实现多谱CT成像方法,关键在于引入参考组分,成像物体是混合物体,参考组分是成像物体其中的一种组分。本发明基于参考组分求解硬化伪影最优校正参数,并利用硬化伪影最优校正参数对成像物体进行硬化伪影重建校正得到不同能量下成像物体硬化伪影校正后的重建图像衰减系数。同时,基于参考组分得到成像物体在不同能量下各组分的衰减系数。最后根据成像物体在不同能量下各组分的衰减系数和不同能量下成像物体硬化伪影校正后的重建图像衰减系数,通过DCM算法进行组分表征,得到成像物体各组分体积表征图像,完成本发明的多谱CT成像方法。
图1为本发明提供的一种多谱CT成像方法的流程示意图,该方法包括:
步骤11、基于参考组分对成像物体进行硬化伪影重建校正得到不同能量下成像物体硬化伪影校正后的重建图像衰减系数。
其中,基于参考组分对成像物体进行硬化伪影重建校正,具体包括:
S111、从成像物体中选择参考组分和参考组分模体;
所述成像物体为混合物体,根据混合物体组分选择参考组分,根据几何结构选择参考组分模体;
S112、对参考组分模体和成像物体进行投影数据采集;
S113、根据所采集的投影数据对成像物体进行硬化伪影重建校正。
本步骤具体包括:
基于不同能量下的参考组分模体投影数据进行CT重建和硬化伪影校正,得到不同能量下的参考组分模体硬化伪影校正后的重建图像衰减系数,将所述不同能量下的参考组分模体硬化伪影校正后的重建图像衰减系数带入硬化伪影校正指标求解相应能量下的参考组分硬化伪影最优校正参数;
利用相应能量下的参考组分硬化伪影最优校正参数,并基于相应能量下的成像物体投影数据进行CT重建和硬化伪影校正,得到不同能量下成像物体硬化伪影校正后的重建图像衰减系数。
步骤12、基于参考组分得到成像物体在不同能量下各组分的衰减系数。
本步骤具体包括:
S121、基于参考组分得到参考组分在不同能量下的衰减系数;
本步骤具体实现为:根据不同能量下的参考组分模体硬化伪影校正后的重建图像衰减系数得到参考组分在不同能量下的衰减系数。
S122、根据参考组分在不同能量下的衰减系数,得到成像物体各组分在不同能量下的衰减系数。
步骤13、根据成像物体在不同能量下各组分的衰减系数和不同能量下成像物体硬化伪影校正后的重建图像衰减系数,通过DCM算法进行组分表征,得到成像物体各组分体积表征图像。
至此,完成了本发明的多谱CT成像方法。
其中,在步骤S113中,本发明基于参考组分求解硬化伪影最优校正参数,这个最优校正参数是确定的值,并非一个范围,因此,可以达到更好的组分表征结果。
参考组分硬化伪影校正参数包括第一校正参数和第二校正参数,所述第一校正参数在第一预定范围内,所述第二校正参数在第二预定范围内,从第一预定范围内确定第一最优校正参数,从第二预定范围内确定第二最优校正参数,第一最优校正参数和第二最优校正参数组成参考组分硬化伪影最优校正参数,不同能量对应不同参考组分硬化伪影最优校正参数;
同一能量下,将参考组分模体硬化伪影校正后的重建图像衰减系数带入硬化伪影校正指标求解参考组分硬化伪影最优校正参数,具体包括:
A、在第一预定范围内等分取任意一个第一校正参数值,作为第一初始值;在第二预定范围内等分取任意一个第二校正参数值,作为第二初始值;
B、固定第二初始值,求解第一最优校正参数;固定第一最优校正参数,求解第二最优校正参数;
如果求解的第一最优校正参数与第一初始值相比有变化,或第二最优校正参数与第二初始值相比有变化,则循环执行步骤B,在第二预定范围内再次选择一个第一校正参数值进行固定,求解第一最优校正参数;固定第一最优校正参数,求解第二最优校正参数;
如果求解的第一最优校正参数与第一初始值相比无变化,且第二最优校正参数与第二初始值相比无变化,则结束循环,确定第一初始值为第一最优校正参数,第二初始值为第二最优校正参数。
其中,固定第二初始值,求解第一最优校正参数,具体包括:
将在第一预定范围内等分选取的r个第一校正参数值,分别与第二初始值相结合,并基于参考组分投影数据进行CT重建和硬化伪影校正,得到r个参考组分模体硬化伪影校正后的重建图像;
针对每个重建图像,将重建图像上参考组分模体柱体直径上等间隔选择的q个点将其衰减系数带入硬化伪影校正指标,将硬化伪影校正指标最小对应的值,作为第一最优校正参数。
同理,固定第一最优校正参数,求解第二最优校正参数的原理与上述方法相同,不再赘述。
为清楚说明本发明,下面在另一实施例中进行详细阐述。
1、根据混合物体组分选择参考组分,根据几何结构选择参考组分模体。
例如,成像物体是由镁(Mg)和铝(Al)混合而成,参考组分为其中的一种物质Al。所选择的参考组分及模体需要满足以下(1)(2)(3)这三个条件。具体地,
在选择参考组分时,需注意参考组分与混合物体各m组分在最低电压时的衰减系数应满足:
其中fm(E)为m组分衰减系数随能量变化的函数,Emax,l为最低能量谱时对应的最大电压,式(1)中的d的经验值为0.8。满足式(1)可防止根据参考组分获取其余组分等效衰减系数时出现较大误差。
为便于评价硬化伪影校正效果,参考组分模体可选择任意规则形状柱体(如圆柱体、正方形柱体等),柱体切面最大长度L1与成像物体最大长度L2满足:
以保证硬化伪影校正效果。
2、对参考组分模体和成像物体进行投影数据采集。
采集参考组分模体的投影数据,即,采集参考组分模体的多张投影,形成一个投影序列。不同能量对应不同的参考组分模体投影序列。
采集成像物体的投影数据,即,采集成像物体的多张投影,形成一个投影序列。不同能量对应不同的成像物体投影序列。
3、基于不同能量下的参考组分模体投影数据进行CT重建和硬化伪影校正,得到不同能量下的参考组分模体硬化伪影校正后的重建图像衰减系数。
1)硬化伪影校正方法
与能量相关的衰减系数可以表示为光电效应μpe和康普顿散射μcs:
μ=μpe+μcs (4)
在5-200keV能量范围内,以下公式近似有效:
其中τ为光电吸收的相互作用效率,ρ为物体的质量密度,Z为原子序数,E为入射原子的能量。康普顿散射在低能量时不明显,光电效应在使用的能量区域为主要效应,因此在考虑校正硬化效应时,仅考虑光电效应就够了。
有效能量Eeff与衰减系数之间可近似为线性关系,这种关系可以表示为:
式(6)中E0为未硬化射线的效能量,Eeff,n为第n个像素的射线有效能量,μi为第i个像素的衰减系数,μ表示有效能量随衰减系数的变化,根据式(5)和式(6),可得:
进而可得:
进而可得在有效能量Eeff下像素n的衰减系数为:
其中α=η/E0。
在重建算法中,SART(simultaneous algebraic reconstruction technique联合代数重建算法)迭代公式为:
将式(9)带入式(10)可得到迭代公式:
其中α∈[0,1],β∈[2.5,3.5]为硬化伪影校正参数。α∈[0,1]为第一校正参数,由成像物体的材料性质和几何结构决定,β∈[2.5,3.5]为第二校正参数,也称为能量依赖参数。
2)将不同能量下的参考组分模体投影数据分别结合公式(11),对参考组分模体进行CT重建和硬化伪影校正,则得到不同能量下的参考组分模体硬化伪影校正后的重建图像。重建图像区域上每一个点都对应一个衰减系数值,重建图像区域对应一个平均后的衰减系数值。
4、将所述不同能量下的参考组分模体硬化伪影校正后的重建图像衰减系数带入硬化伪影校正指标求解相应能量下的参考组分硬化伪影最优校正参数。
已知α∈[0,1],β∈[2.5,3.5],参考组分硬化伪影最优校正参数(α,β)可通过求解下式获得:
其中,a点和b点为参考组分模体重建图像边缘的衰减系数值,在直径上等间隔取q点,cq为各等间隔点的衰减系数值。
式(11)中为降低参数搜索求解难度,在搜索过程中在α的取值范围[0,1]中进行r等分,在β的取值范围[2.5,3.5]中进行s等分,其中,r和s都为自然数。如图2所示,图2为硬化伪影校正参数α,β划分示意图。
图3为重建图像模体直径上的衰减系数曲线示意图。其中a点和b点为参考组分模体重建图像边缘的衰减系数值,在直径上等间隔取q点,cq为各等间隔点的衰减系数值,其中,q为自然数。
在求解最优点(α,β)时,进行如下循环:
1).先在(α,β)中r*s任选一点作为初始值
2).固定β,对α求解:
3).固定α,对β求解:
在上述求解最优点(α,β)过程中,具体方法可以是:
针对步骤1)
由于在第一预定范围[0,1]内,有r个α值;在第二预定范围[2.5,3.5]内,有s个β值。从r个α值中任意选择一个作为第一初始值,从s个β值中任意选择一个作为第二初始值。
针对步骤2)
固定第二初始值βm,则第二初始值对应r个α值(α1…αr),则出现
(α1,βm),(α2,βm)…(αr,βm),m∈[1,s]。
将(α1,βm),(α2,βm)…(αr,βm)分别代入公式(11),基于参考组分投影数据,对参考组分模体进行CT重建和硬化伪影校正,得到r个参考组分模体硬化伪影校正后的重建图像。
针对r中每个重建图像,将重建图像上参考组分模体柱体直径上等间隔选择的q个衰减系数带入硬化伪影校正指标,将硬化伪影校正指标最小对应的值,作为第一最优校正参数。
针对步骤3)
然后固定第一最优校正参数αn,n∈[1,r],则第一最优校正参数对应s个β值(β1…βs),则出现(αn,β1),(αn,β2)…(αn,βs)。
将(αn,β1),(αn,β2)…(αn,βs)分别代入公式(11),基于参考组分投影数据,对参考组分模体进行CT重建和硬化伪影校正,得到s个参考组分模体硬化伪影校正后的重建图像。
针对s中每个重建图像,将重建图像上参考组分模体柱体直径上等间隔选择的q个衰减系数带入硬化伪影校正指标,将硬化伪影校正指标最小对应的值,作为第二最优校正参数。
针对步骤4)
此时,判断
如果求解的第一最优校正参数与第一初始值相比有变化,或第二最优校正参数与第二初始值相比有变化,则循环执行步骤2)和3),在第二预定范围内再次选择一个第一校正参数值进行固定,求解第一最优校正参数;固定第一最优校正参数,求解第二最优校正参数;
如果求解的第一最优校正参数与第一初始值相比无变化,且第二最优校正参数与第二初始值相比无变化,说明此时已经是最优结果,则结束循环,确定第一初始值为第一最优校正参数,第二初始值为第二最优校正参数。
5、对成像物体进行硬化伪影重建校正。
具体为:将不同能量下的成像物体投影数据分别结合代入(α,β)的公式(11),对成像物体进行CT重建和硬化伪影校正,即可对成像物体进行硬化伪影校正,从而得到不同能量下成像物体硬化伪影校正后的重建图像衰减系数。
6、基于参考组分得到成像物体在不同能量下各组分的衰减系数。
1)根据不同能量下的参考组分模体硬化伪影校正后的重建图像衰减系数得到参考组分在不同能量下的衰减系数。
2)根据参考组分在不同能量下的衰减系数,得到成像物体各组分在不同能量下的衰减系数。
μ(1,1)、μ(1,2)···μ(1,k)带入g1(μ)可得:
Ek,eff=g1(μ(1,k)) (15)
式中Ek,eff为第k个能谱的有效能量。将Ek,eff带入fm(E)中:
μ(m,k)=fm(Ek,eff) (16)
其中μ(m,k)为m组分在第k个能谱中的衰减系数。
7、使用DCM进行组分表征。
具体为:根据成像物体在不同能量下各组分的衰减系数和不同能量下成像物体硬化伪影校正后的重建图像衰减系数,通过DCM算法进行组分表征,得到成像物体各组分体积表征图像。
DCM分解的基本假设是:
1、样品的任意一个体元由包含孔隙在内的M种组分构成。
Ikv、某一能量下的CT实验中,任一体元对X射线的吸收等于各组分的吸收之和。
上述假设可表示为:
式(17)、(18)、(19)可写作矩阵形式为:
式中μ(m,k)代表k能量下m种组分的衰减系数,fm表示m种组分的体积分数,μk表示k能量下成像物体硬化伪影校正后的重建图像衰减系数,对上式进行求解,可得:
求解上式即可获得的m种组分的体积结构表征图像。
综上,本发明提供的基于参考组分的多谱CT成像方法,基于参考组分求解相应能量下的参考组分硬化伪影最优校正参数,并利用该参考组分硬化伪影最优校正参数对成像物体进行硬化伪影重建校正,得到不同能量下成像物体硬化伪影校正后的重建图像衰减系数。同时,基于参考组分得到成像物体在不同能量下各组分的衰减系数。最后,根据上述两方面,通过DCM算法进行组分表征,得到成像物体各组分体积表征图像。
以上所述仅为本发明的较佳实施例而已,并不用以限定本发明的包含范围,凡在本发明技术方案的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种多谱电子计算机断层扫描CT成像方法,其特征在于,该方法包括:
基于参考组分对成像物体进行硬化伪影重建校正得到不同能量下成像物体硬化伪影校正后的重建图像衰减系数;
基于参考组分得到成像物体在不同能量下各组分的衰减系数;
根据成像物体在不同能量下各组分的衰减系数和不同能量下成像物体硬化伪影校正后的重建图像衰减系数,通过数据约束模型DCM算法进行组分表征,得到成像物体各组分体积表征图像。
2.根据权利要求1所述的方法,其特征在于,基于参考组分对成像物体进行硬化伪影重建校正,具体包括:
从成像物体中选择参考组分和参考组分模体;
对参考组分模体和成像物体进行投影数据采集;
根据所采集的投影数据对成像物体进行硬化伪影重建校正。
4.根据权利要求2所述的方法,其特征在于,根据所采集的投影数据对成像物体进行硬化伪影重建校正,具体包括:
基于不同能量下的参考组分模体投影数据进行CT重建和硬化伪影校正,得到不同能量下的参考组分模体硬化伪影校正后的重建图像衰减系数,将所述不同能量下的参考组分模体硬化伪影校正后的重建图像衰减系数带入硬化伪影校正指标求解相应能量下的参考组分硬化伪影最优校正参数;
利用相应能量下的参考组分硬化伪影最优校正参数,并基于相应能量下的成像物体投影数据进行CT重建和硬化伪影校正,得到不同能量下成像物体硬化伪影校正后的重建图像衰减系数。
5.根据权利要求1所述的方法,其特征在于,基于参考组分得到成像物体在不同能量下各组分的衰减系数,具体包括:
基于参考组分得到参考组分在不同能量下的衰减系数;
根据参考组分在不同能量下的衰减系数,得到成像物体各组分在不同能量下的衰减系数。
6.根据权利要求5所述的方法,其特征在于,基于参考组分得到参考组分在不同能量下的衰减系数,具体包括:
根据不同能量下的参考组分模体硬化伪影校正后的重建图像衰减系数得到参考组分在不同能量下的衰减系数。
7.根据权利要求4所述的方法,其特征在于,参考组分硬化伪影校正参数包括第一校正参数和第二校正参数,所述第一校正参数在第一预定范围内,所述第二校正参数在第二预定范围内,从第一预定范围内确定第一最优校正参数,从第二预定范围内确定第二最优校正参数,第一最优校正参数和第二最优校正参数组成参考组分硬化伪影最优校正参数,不同能量对应不同参考组分硬化伪影最优校正参数;
同一能量下,将参考组分模体硬化伪影校正后的重建图像衰减系数带入硬化伪影校正指标求解参考组分硬化伪影最优校正参数,具体包括:
A、在第一预定范围内等分取任意一个第一校正参数值,作为第一初始值;在第二预定范围内等分取任意一个第二校正参数值,作为第二初始值;
B、固定第二初始值,求解第一最优校正参数;固定第一最优校正参数,求解第二最优校正参数;
如果求解的第一最优校正参数与第一初始值相比有变化,或第二最优校正参数与第二初始值相比有变化,则循环执行步骤B,在第二预定范围内再次选择一个第一校正参数值进行固定,求解第一最优校正参数;固定第一最优校正参数,求解第二最优校正参数;
如果求解的第一最优校正参数与第一初始值相比无变化,且第二最优校正参数与第二初始值相比无变化,则结束循环,确定第一初始值为第一最优校正参数,第二初始值为第二最优校正参数。
8.根据权利要求7所述的方法,其特征在于,固定第二初始值,求解第一最优校正参数,具体包括:
将在第一预定范围内等分选取的r个第一校正参数值,分别与第二初始值相结合,并基于参考组分投影数据进行CT重建和硬化伪影校正,得到r个参考组分模体硬化伪影校正后的重建图像;
针对每个重建图像,将重建图像上参考组分模体柱体直径上等间隔选择的q个衰减系数带入硬化伪影校正指标,将硬化伪影校正指标最小对应的值,作为第一最优校正参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010271896.2A CN111476856B (zh) | 2020-04-08 | 2020-04-08 | 一种多谱ct成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010271896.2A CN111476856B (zh) | 2020-04-08 | 2020-04-08 | 一种多谱ct成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111476856A true CN111476856A (zh) | 2020-07-31 |
CN111476856B CN111476856B (zh) | 2023-06-06 |
Family
ID=71751207
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010271896.2A Active CN111476856B (zh) | 2020-04-08 | 2020-04-08 | 一种多谱ct成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111476856B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112581556A (zh) * | 2020-12-25 | 2021-03-30 | 上海联影医疗科技股份有限公司 | 多能ct图像硬化校正方法、装置、计算机设备和存储介质 |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1235811A (zh) * | 1998-05-14 | 1999-11-24 | 通用电气公司 | 对于高衰减物体产生的虚影的校正 |
US6249595B1 (en) * | 1998-01-22 | 2001-06-19 | General Electric Company | Iterative reconstruction for EPI |
US20040208290A1 (en) * | 2003-04-17 | 2004-10-21 | Masatake Nukui | X-ray CT system and beam-hardening post-processing method |
CN103876772A (zh) * | 2014-03-20 | 2014-06-25 | 中北大学 | 一种多谱成像方法和装置 |
US20150164456A1 (en) * | 2012-08-30 | 2015-06-18 | Kabushiki Kaisha Toshiba | Medical image processing apparatus and x-ray computed tomography apparatus |
CN105608720A (zh) * | 2016-01-30 | 2016-05-25 | 上海联影医疗科技有限公司 | 计算机断层成像伪影校正方法及装置 |
US20160278733A1 (en) * | 2013-04-08 | 2016-09-29 | Hitachi Medical Corporation | X-ray ct apparatus, reconstruction arithmetic device, and reconstruction arithmetic method |
WO2016158138A1 (ja) * | 2015-04-01 | 2016-10-06 | 株式会社日立製作所 | X線ct装置、再構成演算装置およびx線ct画像生成方法 |
CN106204673A (zh) * | 2016-06-29 | 2016-12-07 | 中国人民解放军信息工程大学 | 基于能谱滤波和图像残差重投影的ct图像金属伪影校正方法 |
KR101698033B1 (ko) * | 2016-02-26 | 2017-01-19 | 연세대학교 산학협력단 | Ct 영상에서의 콘빔 아티팩트 보정 장치 및 방법 |
US20170061654A1 (en) * | 2015-08-24 | 2017-03-02 | Analogic Corporation | Projection data correction and computed tomography value computation |
KR20170026984A (ko) * | 2015-09-01 | 2017-03-09 | 연세대학교 산학협력단 | 엑스레이 컴퓨터 단층촬영 환경에서 빔 경화현상에 의한 인공물 보정방법 |
CN107895390A (zh) * | 2017-08-28 | 2018-04-10 | 深圳联影医疗科技有限公司 | 电子计算机断层扫描图像数据处理系统及方法 |
GB201810788D0 (en) * | 2018-06-29 | 2018-08-15 | Biocompatibles Uk Ltd | Biodegradable polymer |
CN110811660A (zh) * | 2019-10-25 | 2020-02-21 | 赛诺威盛科技(北京)有限公司 | 一种校正ct射线束硬化伪影的方法 |
CN110827370A (zh) * | 2019-11-09 | 2020-02-21 | 中北大学 | 一种非等厚构件的多能ct循环迭代重建方法 |
-
2020
- 2020-04-08 CN CN202010271896.2A patent/CN111476856B/zh active Active
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6249595B1 (en) * | 1998-01-22 | 2001-06-19 | General Electric Company | Iterative reconstruction for EPI |
CN1235811A (zh) * | 1998-05-14 | 1999-11-24 | 通用电气公司 | 对于高衰减物体产生的虚影的校正 |
US20040208290A1 (en) * | 2003-04-17 | 2004-10-21 | Masatake Nukui | X-ray CT system and beam-hardening post-processing method |
US20150164456A1 (en) * | 2012-08-30 | 2015-06-18 | Kabushiki Kaisha Toshiba | Medical image processing apparatus and x-ray computed tomography apparatus |
US20160278733A1 (en) * | 2013-04-08 | 2016-09-29 | Hitachi Medical Corporation | X-ray ct apparatus, reconstruction arithmetic device, and reconstruction arithmetic method |
CN103876772A (zh) * | 2014-03-20 | 2014-06-25 | 中北大学 | 一种多谱成像方法和装置 |
WO2016158138A1 (ja) * | 2015-04-01 | 2016-10-06 | 株式会社日立製作所 | X線ct装置、再構成演算装置およびx線ct画像生成方法 |
US20170061654A1 (en) * | 2015-08-24 | 2017-03-02 | Analogic Corporation | Projection data correction and computed tomography value computation |
KR20170026984A (ko) * | 2015-09-01 | 2017-03-09 | 연세대학교 산학협력단 | 엑스레이 컴퓨터 단층촬영 환경에서 빔 경화현상에 의한 인공물 보정방법 |
CN105608720A (zh) * | 2016-01-30 | 2016-05-25 | 上海联影医疗科技有限公司 | 计算机断层成像伪影校正方法及装置 |
KR101698033B1 (ko) * | 2016-02-26 | 2017-01-19 | 연세대학교 산학협력단 | Ct 영상에서의 콘빔 아티팩트 보정 장치 및 방법 |
CN106204673A (zh) * | 2016-06-29 | 2016-12-07 | 中国人民解放军信息工程大学 | 基于能谱滤波和图像残差重投影的ct图像金属伪影校正方法 |
CN107895390A (zh) * | 2017-08-28 | 2018-04-10 | 深圳联影医疗科技有限公司 | 电子计算机断层扫描图像数据处理系统及方法 |
GB201810788D0 (en) * | 2018-06-29 | 2018-08-15 | Biocompatibles Uk Ltd | Biodegradable polymer |
CN110811660A (zh) * | 2019-10-25 | 2020-02-21 | 赛诺威盛科技(北京)有限公司 | 一种校正ct射线束硬化伪影的方法 |
CN110827370A (zh) * | 2019-11-09 | 2020-02-21 | 中北大学 | 一种非等厚构件的多能ct循环迭代重建方法 |
Non-Patent Citations (8)
Title |
---|
DENG, JUNJUN ET AL.: "Beam Hardening Correction Using an Attenuation Coefficient Decomposition Approach", 《IEEE NUCLEAR SCIENCE SYMPOSIUM AND MEDICAL IMAGING CONFERENCE》 * |
HEUSSER T ET AL.: "Prior-based artifact correction (PBAC) in computed tomography", 《MEDICAL PHYSICS》 * |
PENGCHONG JIN ET AL.: "A method for simultaneous image reconstruction and beam hardening correction", 《2013 IEEE NUCLEAR SCIENCE SYMPOSIUM AND MEDICAL IMAGING CONFERENCE (2013 NSS/MIC)》 * |
李磊等: "基于加权TV正则化的X射线CT系统能谱估计方法", 《光谱学与光谱分析》 * |
王文杰 等: "基于参考组分的双能CT成像方法", 《CT理论与应用研究》 * |
谢德华: "基于参考组分的多能CT组分表征算法研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》 * |
赵云松 等: "双能谱CT的迭代重建模型及重建方法", 《电子学报》 * |
陈思宇: "工业锥束CT图像环状及硬化伪影校正技术研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112581556A (zh) * | 2020-12-25 | 2021-03-30 | 上海联影医疗科技股份有限公司 | 多能ct图像硬化校正方法、装置、计算机设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN111476856B (zh) | 2023-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Szczykutowicz et al. | Dual energy CT using slow kVp switching acquisition and prior image constrained compressed sensing | |
CN108010098B (zh) | 一种双能谱ct基材料图像迭代重建方法 | |
Ruchala et al. | Megavoltage CT on a tomotherapy system | |
Elbakri et al. | Segmentation-free statistical image reconstruction for polyenergetic x-ray computed tomography with experimental validation | |
US10064591B2 (en) | System, method and computer readable medium for preview of low-dose x-ray projection and tomographic images | |
US9600866B2 (en) | Projection data de-noising | |
Zhang et al. | A practical reconstruction method for dual energy computed tomography | |
CN110175957B (zh) | 多能量ct基材料物质分解方法 | |
Ding et al. | Image‐domain multimaterial decomposition for dual‐energy CT based on prior information of material images | |
Rakvongthai et al. | Spectral CT using multiple balanced K-edge filters | |
Evans et al. | Experimental implementation of a polyenergetic statistical reconstruction algorithm for a commercial fan-beam CT scanner | |
KR20200059712A (ko) | 엑스선 영상 산란 보정 시스템 및 방법 | |
CN111476856A (zh) | 一种多谱ct成像方法 | |
Chen et al. | Line integral alternating minimization algorithm for dual-energy x-ray CT image reconstruction | |
CN113298905B (zh) | 基于基础材料分解的多能谱ct重建方法及计算机可读介质 | |
Chung et al. | Numerical algorithms for polyenergetic digital breast tomosynthesis reconstruction | |
Zhao et al. | K-edge eliminated material decomposition method for dual-energy X-ray CT | |
CN116183647A (zh) | 一种物质识别方法 | |
Kim et al. | Accuracy improvement of interior region of interest reconstruction using field modulation CT acquisition | |
Zhang | Basis vector model method for proton stopping power estimation using dual-energy computed tomography | |
Altunbas | 10 Image corrections for scattered radiation | |
Abella et al. | A new statistical image reconstruction algorithm for polyenergetic X-ray CT | |
Huh et al. | Fast kVp-switching dual energy CT for PET attenuation correction | |
Jo et al. | Scatter correction using a primary modulator for dual energy digital radiography: A Monte Carlo simulation study | |
Humphries et al. | Comparison of image and data domain methods for three-material decomposition in dual-energy CT |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |