CN116664707B - 一种ct图像迭代重建方法、装置及计算机设备 - Google Patents

一种ct图像迭代重建方法、装置及计算机设备 Download PDF

Info

Publication number
CN116664707B
CN116664707B CN202210146430.9A CN202210146430A CN116664707B CN 116664707 B CN116664707 B CN 116664707B CN 202210146430 A CN202210146430 A CN 202210146430A CN 116664707 B CN116664707 B CN 116664707B
Authority
CN
China
Prior art keywords
value
image
angle
projection
image value
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.)
Active
Application number
CN202210146430.9A
Other languages
English (en)
Other versions
CN116664707A (zh
Inventor
何艾静
杨书
张康平
孙宇
张文宇
吴宏新
王亚杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BEIJING LANGSHI INSTRUMENT CO LTD
Original Assignee
BEIJING LANGSHI INSTRUMENT CO LTD
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by BEIJING LANGSHI INSTRUMENT CO LTD filed Critical BEIJING LANGSHI INSTRUMENT CO LTD
Priority to CN202210146430.9A priority Critical patent/CN116664707B/zh
Publication of CN116664707A publication Critical patent/CN116664707A/zh
Application granted granted Critical
Publication of CN116664707B publication Critical patent/CN116664707B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/60—Rotation of whole images or parts thereof
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10081—Computed x-ray tomography [CT]

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明公开了一种CT图像迭代重建方法、装置及计算机设备，该方法包括，根据目标对象在初始投影角度下的先验图像值和径迹矩阵，得到在对应角度的投影理论值，基于投影理论值和相应测量值确定偏差值，基于偏差值更新先验图像值，以单位旋转角度和更新后的先验图像值旋转，得到旋转单位投影角度后的图像值；将旋转单位投影角度后的图像值作为旋转单位投影角度的先验图像，并将单位投影角度作为初始投影角度，返回执行基于先验图像值和径迹矩阵得到在初始投影角度下的投影理论值的步骤至基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值的步骤，直至旋转角度达到预设角度值，基于各单位旋转角度对应图像值得到CT图像重建结果。

Description

一种CT图像迭代重建方法、装置及计算机设备
技术领域
本发明涉及图像重建技术领域，具体涉及一种CT图像迭代重建方法、装置及计算机设备。
背景技术
随着锥形束CT(Computed Tomography，电子计算机断层扫描，简称CT)技术的产生及发展，使得口腔颌面临床诊断及治疗水平得到了很大的提高，然而辐射剂量对被照患者的潜在危害性也受到广泛关注。在保证图像质量和诊疗效果的同时，尽可能减少受检者的辐射剂量、合理使用低剂量已成为当前CT重建领域重要研究方向。
CT重建算法可分为迭代类重建算法和解析类重建算法，解析重建算法具有分辨率高和成像速度快等优点，但易受统计波动的影响，投影数量若不足，重建图像的质量会明显下降，因此为了保证重建效果需保证投影数据的完备性，这对扫描过程的辐射剂量要求较高。迭代重建算法可以有效克服解析类重建算法的固有问题，能在投影数量不足的条件下重建出较高质量的图像，可有效减少投影数据的获取、减少扫描时间、从而减少患者受照剂量，作为目前低剂量成像的实现手段之一被广泛研究。
但迭代类重建算法存在占用较大的计算资源以及计算实时性较差的缺点。传统迭代算法中需要计算每个旋转角度的径迹矩阵，即在每个旋转角度下，每条射线穿过各个重建体素内的长度，对于探测器偏置的锥形束CT扫描系统，存在物理模型复杂、待求量较大的问题，导致迭代重建过程计算量大、重建时间长。
发明内容
因此，本发明要解决的技术问题在于克服现有迭代重建算法计算量大、重建时间长的缺陷，从而提供一种CT图像迭代重建方法、装置及计算机设备。
根据第一方面，本发明实施例公开了一种CT图像迭代重建方法，包括：获取目标对象在初始投影角度下的先验图像值及径迹矩阵；基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值；基于所述投影理论值和在所述初始投影角度下的测量值确定偏差值；基于所述偏差值更新所述先验图像值，得到更新后的先验图像值；基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值；将旋转单位投影角度后的图像值作为旋转单位投影角度的先验图像，并将所述单位投影角度作为初始投影角度，返回执行基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值的步骤至基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值的步骤，直至当前旋转角度达到预设角度值；基于所述当前旋转角度中各单位旋转角度对应的旋转单位投影角度后的图像值得到CT图像重建结果。
可选地，在基于所述当前旋转角度中各单位旋转角度对应的旋转单位投影角度后的图像值得到CT图像重建结果之前，所述方法还包括：基于相邻两次单位旋转角度对应的旋转单位投影角度后的图像值判断迭代过程是否满足预设收敛条件；若迭代过程满足预设收敛条件，则执行基于所述当前旋转角度中各单位旋转角度对应的旋转单位投影角度后的图像值得到CT图像重建结果的步骤。
可选地，所述方法还包括：若迭代过程不满足预设收敛条件，则返回执行基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值的步骤至基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值的步骤。
可选地，所述获取目标对象在初始投影角度下的先验图像值及径迹矩阵，包括：获取所述目标对象的投影数据和初始投影角度；基于所述投影数据和预设的重建算法得到所述先验图像值；基于所述先验图像值得到初始投影角度下的径迹矩阵。
可选地，所述基于所述投影数据和预设的重建算法得到所述先验图像值，包括：对所述投影数据进行加权和滤波处理，得到第一图像值；对所述第一图像值进行加权反投影得到所述先验图像值。
可选地，所述基于所述先验图像值得到初始投影角度下的径迹矩阵，包括：基于CT射线源和像素点得到穿过所述先验图像集的射线长度；基于各所述射线长度得到所述径迹矩阵。
可选地，所述基于所述偏差值更新所述先验图像值，得到更新后的先验图像值，包括：基于所述偏差值和所述径迹矩阵中的每一条射线得到各所述先验图像值的校正值；基于各所述校正值得到所述更新后的先验图像值。
根据第二方面，本发明实施例还公开了一种CT图像迭代重建装置，包括：获取模块，用于获取目标对象在初始投影角度下的先验图像值及径迹矩阵；投影理论值模块，用于基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值；偏差模块，用于基于所述投影理论值和在所述初始投影角度下的测量值确定偏差值；更新模块，用于基于所述偏差值更新所述先验图像值，得到更新后的先验图像值；第一旋转模块，用于基于所述旋转角度和更新后的先验图像值，得到旋转角度后的图像值；第二旋转模块，用于将旋转单位投影角度后的图像值作为旋转单位投影角度的先验图像，并将所述单位投影角度作为初始投影角度，返回执行基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值的步骤至基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值的步骤，直至当前旋转角度达到预设角度值；输出模块，用于基于所述当前旋转角度中各单位旋转角度对应的旋转单位投影角度后的图像值得到CT图像重建结果。
根据第三方面，本发明实施例还公开了一种计算机设备，包括：至少一个处理器；以及与所述至少一个处理器通信连接的存储器；其中，所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器执行如第一方面或第一方面任一可选实施方式所述的CT图像迭代重建方法的步骤。
根据第四方面，本发明实施方式还公开了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现如第一方面或第一方面任一可选实施方式所述的CT图像迭代重建方法的步骤。
本发明技术方案，具有如下优点：
本发明提供的CT图像迭代重建方法、装置及计算机设备，该方法包括：获取目标对象在初始投影角度下的先验图像值及径迹矩阵；基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值；基于所述投影理论值和在所述初始投影角度下的测量值确定偏差值；基于所述偏差值更新所述先验图像值，得到更新后的先验图像值；基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值；将旋转单位投影角度后的图像值作为旋转单位投影角度的先验图像，并将所述单位投影角度作为初始投影角度，返回执行基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值的步骤至基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值的步骤，直至当前旋转角度达到预设角度值；基于所述当前旋转角度中各单位旋转角度对应的旋转单位投影角度后的图像值得到CT图像重建结果。根据重建图像值作为先验图像值，并在初始投影角度下计算径迹矩阵，在其他角度下根据前一角度的图像值作为先验矩阵，只需要进行一次径迹矩阵的计算，减少了计算量，降低了物理模型的复杂度。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中CT图像迭代重建方法的一个具体示例的流程图；
图2为本发明实施例中CT图像迭代重建装置的一个具体示例的原理框图；
图3为本发明实施例中计算机设备的一个具体示例图；
图4为本发明实施例中CT图像迭代重建方法的一个具体示例的示意图；
图5为本发明实施例中CT图像迭代重建方法的一个具体示例的示意图；
图6为本发明实施例中CT图像迭代重建方法的一个具体示例的示意图；
图7为本发明实施例中CT图像迭代重建方法的一个具体示例的示意图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，还可以是两个元件内部的连通，可以是无线连接，也可以是有线连接。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
此外，下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
本发明实施例公开了一种CT图像迭代重建方法，如图1所示，该方法包括如下步骤：
步骤101，获取目标对象在初始投影角度下的先验图像值及径迹矩阵。
示例性地，目标对象是指需要进行CT成像的对象，例如可以是人体的某一部位，例如是脑部、胸部、肺部、心脏以及腹部等。先验图像值为在初始投影角度下每一个点的像素值的集合，对应在CT图像中也可以称为体素，先验图像值也相当于在初始投影角度下的各个像素点的初始值，可以根据实际情况进行赋值。径迹矩阵为在对应的投影角度下，射线源发出的每条射线穿过各个先验图像值的长度，在初始投影角度下，根据先验图像值可以计算得到穿过每一个先验图像值的射线长度，所有的射线长度共同构成初始投影角度下的径迹矩阵。本发明实施例对目标对象的类型和先验图像值的大小形式不作限定，本领域技术人员可以根据实际需要确定。
步骤102，基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值。
示例性地，根据预设的先验图像值和计算得到的径迹矩阵，可以计算得到径迹矩阵中每一条射线对应的投影理论值，其中投影理论值可以通过以下公式计算得到，
其中为射线aij所对应的投影理论值，aij为径迹矩阵中的射线，/>为先验图像值，i为径迹矩阵中的射线编号，j为先验图像值中的像素编号，J为总的先验图像值数量，i、j为正整数。
步骤103，基于所述投影理论值和在所述初始投影角度下的测量值确定偏差值。
示例性地，根据扫描目标对象所得的测量值，计算测量值与投影理论值的差值，得到该角度下所有射线的投影值偏差Δi，其中pi为测量值，/>为投影理论值。
步骤104，基于所述偏差值更新所述先验图像值，得到更新后的先验图像值。示例性地，在进行CT图像扫描时，在初始投影角度下可以获取实际的测量值，根据测量值与对应投影理论值之间的差值，可以根据投影理论值和差值计算得到先验图像值中所有的像素点的真实测量值，且不需要再一次进行测量，只需测量一次即可，减少了目标对应的CT扫描时间。
步骤105，基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值。
示例性地，单位旋转角度为对应CT图像扫描过程中，需要进行总旋转角度的扫描精度，例如总旋转角度可以是扫描旋转180度，其中单位旋转角度可以是每0.02度进行一次扫描。根据单位旋转角度将更新后的先验图像值进行旋转，根据更新后的先验图像值和单位旋转角度，由公式(1)计算得到旋转单位投影角度后的理论图像值，根据步骤103中计算得到偏差值，更新计算得到旋转单位投影角度后的图像值。本发明实施例对单位旋转角度和总旋转角度的数值大小不作限定，本领域技术人员可以根据实际需要确定。
如图5所示为先验图像旋转的示意图，先验图像中的每一个像素点跟随投影角度同步旋转，每次放射源-探测器系统在旋转单位投影角度时，更新的先验图像矩阵跟着旋转相应角度，利用当前重建的图像值X(0)，进行坐标旋转单位投影角度后的插值，得到下一投影角度下的重建值插值后的各射线穿过各重建体素的长度与关系，与初始投影角度下的一致。
本发明实施例以FDK算法和CBCT系统(Cone beasm锥形束，简称CBCT)进行介绍，在CBCT系统中通常旋转一圈采集的帧数都大于365帧(很多大于600帧)，因此单次的旋转角度(单位旋转角度)会小于1度，在这样的小旋转角度下，对每一重建切片层的像素来说，旋转之前的坐标为(x1,y1,z1,)，旋转角度θ之后的坐标为(x2,y2,z2,)，其中x2＝x1sinθ+y1cosθ，y2＝x1cosθ-y1sinθ，z1＝z2，因此只需考虑X、Y轴方向的插值变换，利用旋转前的坐标和图像值，进行旋转后坐标的插值计算，得到旋转后的图像值。插值方法有很多，此处边界像素可以采用最近邻插值法(变换后的像素的图像值等于距它最近的变换前像素的图像值)，中间像素可以选用双线性插值法，其原理是要想得到旋转后坐标为(x,y)的P点的图像值f(x,y)，先找到旋转前离其最近的四个像素点Q11(x1,y1)、Q12(x1,y2)、Q21(x2,y1)、Q22(x2,y2)及其图像值，计算公式为
将每一层切面的重建图像值旋转插值之后的结果整合到一起，得到旋转之后的重建三维像素矩阵如图6为插值计算示意图，图7为三维像素矩阵重建示意图。
步骤106，将旋转单位投影角度后的图像值作为旋转单位投影角度的旋转单位投影角度后的先验图像，并将所述单位投影角度作为初始投影角度，返回执行步骤102至步骤105，直至当前旋转角度达到预设角度值。
示例性地，在旋转单位投影角度后得到图像值作为下一次旋转单位投影角度的先验图像值，即前一次得到的图像值，作为下一次投影的先验图像值，重复执行上述步骤102至步骤105，直至总旋转角度达到预设值，如上述步骤105中的180度。
步骤107，基于所述当前旋转角度中各单位旋转角度对应的旋转单位投影角度后的图像值得到CT图像重建结果。示例性地，在达到总旋转角度后，根据各个单位旋转角度下的图像值，共同构成CT图像的重建结果。
本发明提供的CT图像迭代重建方法，包括：获取目标对象在初始投影角度下的先验图像值及径迹矩阵；基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值；基于所述投影理论值和在所述初始投影角度下的测量值确定偏差值；基于所述偏差值更新所述先验图像值，得到更新后的先验图像值；基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值；将旋转单位投影角度后的图像值作为旋转单位投影角度的先验图像，并将所述单位投影角度作为初始投影角度，返回执行基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值的步骤至基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值的步骤，直至当前旋转角度达到预设角度值；基于所述当前旋转角度中各单位旋转角度对应的旋转单位投影角度后的图像值得到CT图像重建结果。根据重建图像值作为先验图像值，并在初始投影角度下计算径迹矩阵，在其他角度下根据前一角度的图像值作为先验矩阵，只需要进行一次初始投影角度下径迹矩阵的计算，减少了计算量，降低了物理模型的复杂度。
作为本发明一个可选实施方式，在步骤107之前，所述方法还包括：基于相邻两次单位旋转角度对应的旋转单位投影角度后的图像值判断迭代过程是否满足预设收敛条件；若迭代过程满足预设收敛条件，则执行步骤107。
示例性地，在旋转角度达到总旋转角度后，需要判断得到的图像重建值是否满足预设收敛条件，其中预设收敛条件由相邻两次的图像值计算得到，通过以下公式计算得到，
其中X(n)为第n次旋转单位投影角度后的图像值，X(n-1),为第n-1次旋转单位投影角度后的图像值，当ρ的值满足预设阈值时，执行步骤107，其中预设阈值可以是ρ≤0.02，本发明实施例对预设阈值的数据大小不作限定，本领域技术人员可以根据实际需要确定。
作为本发明一个可选实施方式，在步骤107之前，所述方法还包括：若迭代过程不满足预设收敛条件，则返回执行步骤102至步骤105。示例性地，若不满足将此重建结果作为下一次迭代的初值X(k+1)＝X(k)并重复步骤102至步骤105，直到满足设定的收敛条件，输出对应的重建结果。
作为本发明一个可选实施方式，所述步骤101，包括：获取所述目标对象的投影数据和初始投影角度；基于所述投影数据和预设的重建算法得到所述先验图像值；基于所述先验图像值得到初始投影角度下的径迹矩阵。
作为本发明一个可选实施方式，步骤101中，所述基于所述投影数据和预设的重建算法得到所述先验图像值，包括：对所述投影数据进行加权和滤波处理，得到第一图像值；对所述第一图像值进行加权反投影得到所述先验图像值。
示例性地，在获取先验图像值和相应的径迹矩阵时，需要根据目标对象的实际需要投影的范围等构造先验图像值。预设的重建算法可以是FDK算法，先验图像值的计算过程如下：
对投影数据预加权处理；
进行滤波；
进行加权反投影计算；
其中β为初始投影角度，p(β,a,b)，为投影数据，γ为CBCT系统的扇形角度，k为CBCT系统的锥角，R为圆轨道半径，g(a)为相应的滤波函数，(a,b)为探测器坐标，(x,y,z)为重建点坐标。f(x,y,z)为最终重建的三维图像数据，将f(x,y,z)作为初始投影角度下的先验图像值X(0)。
作为本发明一个可选实施方式，步骤101中，所述基于所述先验图像值得到初始投影角度下的径迹矩阵，包括：基于CT射线源和像素点得到穿过所述先验图像集的射线长度；基于各所述射线长度得到所述径迹矩阵。
示例性地，径迹矩阵可以采用点-点模型，将探测器上的像素和射线源都分别视为一个点，两点的连线穿过中间重建图像值的长度可建立坐标系求得，如图4所示，为第i条射线穿过4*4个体素，其中这条射线穿过每一个体素的长度即为径迹矩阵，例如，如图4所示，ai1＝射线i穿过像素1内的长度＝L5，ai2＝射线i穿过体素1内的长度＝0。
作为本发明一个可选实施方式，所述步骤104包括：基于所述偏差值和所述径迹矩阵中的每一条射线得到各所述先验图像值的校正值；基于各所述校正值得到所述更新后的先验图像值。
示例性地，沿每条径迹矩阵的射线反投影于其经过的像素，在迭代初始值的基础上依次校正理论值，
其中，为校正后的像素值，/>为校正前的像素值，aij为径迹矩阵中的射线，Δi测量值与投影理论值的差值，i为径迹矩阵中的射线编号，j为先验图像值中的像素编号，J为总的先验图像值数量，像素值为图像值中的各像素值。
本发明实施例还公开了一种CT图像迭代重建装置，如图2所示，该装置包括：
获取模块201，用于获取目标对象在初始投影角度下的先验图像值及径迹矩阵。示例性地，详细内容见上述方法实施例中步骤101的内容，此处不再赘述。
投影理论值模块202，用于基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值。示例性地，详细内容见上述方法实施例中步骤102的内容，此处不再赘述。
偏差模块203，用于基于所述投影理论值和在所述初始投影角度下的测量值确定偏差值。示例性地，详细内容见上述方法实施例中步骤103的内容，此处不再赘述。
更新模块204，用于基于所述偏差值更新所述先验图像值，得到更新后的先验图像值。示例性地，详细内容见上述方法实施例中步骤104的内容，此处不再赘述。
第一旋转模块205，用于基于所述旋转角度和更新后的先验图像值，得到旋转角度后的图像值。示例性地，详细内容见上述方法实施例中步骤105的内容，此处不再赘述。
第二旋转模块206，用于将旋转单位投影角度后的图像值作为旋转单位投影角度的先验图像，并将所述单位投影角度作为初始投影角度，返回执行基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值的步骤至基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值的步骤，直至当前旋转角度达到预设角度值。示例性地，详细内容见上述方法实施例中步骤106的内容，此处不再赘述。
输出模块207，用于基于所述当前旋转角度中各单位旋转角度对应的旋转单位投影角度后的图像值得到CT图像重建结果。示例性地，详细内容见上述方法实施例中步骤107的内容，此处不再赘述。
本发明提供的CT图像迭代重建装置，包括，获取模块201，用于获取目标对象在初始投影角度下的先验图像值及径迹矩阵；投影理论值模块202，用于基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值；偏差模块203，用于基于所述投影理论值和在所述初始投影角度下的测量值确定偏差值；更新模块204，用于基于所述偏差值更新所述先验图像值，得到更新后的先验图像值；第一旋转模块205，用于基于所述旋转角度和更新后的先验图像值，得到旋转角度后的图像值；第二旋转模块206，用于将旋转单位投影角度后的图像值作为旋转单位投影角度的先验图像，并将所述单位投影角度作为初始投影角度，返回执行基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值的步骤至基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值的步骤，直至当前旋转角度达到预设角度值；输出模块207，用于基于所述当前旋转角度中各单位旋转角度对应的旋转单位投影角度后的图像值得到CT图像重建结果。根据重建图像值作为先验图像值，并在初始投影角度下计算径迹矩阵，在其他角度下根据前一角度的图像值作为先验矩阵，只需要进行一次计算，减少了计算量，降低了物理模型的复杂度。
作为本发明一个可选实施方式，所述装置还包括：第一判断模块，用于基于相邻两次单位旋转角度对应的旋转单位投影角度后的图像值判断迭代过程是否满足预设收敛条件，若迭代过程满足预设收敛条件，则触发输出模块207，执行基于所述当前旋转角度中各单位旋转角度对应的旋转单位投影角度后的图像值得到CT图像重建结果过程。示例性地，详细内容见上述方法实施例中步骤107的内容，此处不再赘述。
作为本发明一个可选实施方式，所述装置还包括：第二判断模块，用于判断若迭代过程不满足预设收敛条件，则返回触发投影理论值模块202、偏差模块203、更新模块204和第一旋转模块205，并执行对应的步骤。示例性地，详细内容见上述方法实施例中步骤107的内容，此处不再赘述。
作为本发明一个可选实施方式，所述获取模块201，包括：获取子模块，用于获取所述目标对象的投影数据和初始投影角度；第一计算模块，用于基于所述投影数据和预设的重建算法得到所述先验图像值；第二计算模块，用于基于所述先验图像值得到初始投影角度下的径迹矩阵。示例性地，详细内容见上述方法实施例中步骤101的内容，此处不再赘述。
作为本发明一个可选实施方式，步骤获取模块201，所述基于所述投影数据和预设的重建算法得到所述先验图像值，包括：第一图像值模块，用于对所述投影数据进行加权和滤波处理，得到第一图像值；第二图像值模块，用于对所述第一图像值进行加权反投影得到所述先验图像值。示例性地，详细内容见上述方法实施例中步骤101的内容，此处不再赘述。
作为本发明一个可选实施方式，所述获取模块201中，所述基于所述先验图像值得到初始投影角度下的径迹矩阵，包括：第一径迹矩阵模块，用于基于CT射线源和像素点得到穿过所述先验图像集的射线长度；第二径迹矩阵模块，用于基于各所述射线长度得到所述径迹矩阵。示例性地，详细内容见上述方法实施例中步骤101的内容，此处不再赘述。
作为本发明一个可选实施方式，所述更新模块204包括：第一校正子模块，用于基于所述偏差值和所述径迹矩阵中的每一条射线得到各所述先验图像值的校正值；第二校正子模块，用于基于各所述校正值得到所述更新后的先验图像值。示例性地，详细内容见上述方法实施例中步骤104的内容，此处不再赘述。
本发明实施例还提供了一种计算机设备，如图3所示，该计算机设备可以包括处理器301和存储器302，其中处理器301和存储器302可以通过总线或者其他方式连接，图3中以通过总线连接为例。
处理器301可以为中央处理器(Central Processing Unit，CPU)。处理器301还可以为其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现场可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等芯片，或者上述各类芯片的组合。
存储器302作为一种非暂态计算机可读存储介质，可用于存储非暂态软件程序、非暂态计算机可执行程序以及模块，如本发明实施例中的CT图像迭代重建方法对应的程序指令/模块。处理器301通过运行存储在存储器302中的非暂态软件程序、指令以及模块，从而执行处理器的各种功能应用以及数据处理，即实现上述方法实施例中的CT图像迭代重建方法。
存储器302可以包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需要的应用程序；存储数据区可存储处理器301所创建的数据等。此外，存储器302可以包括高速随机存取存储器，还可以包括非暂态存储器，例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施例中，存储器302可选包括相对于处理器301远程设置的存储器，这些远程存储器可以通过网络连接至处理器301。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
所述一个或者多个模块存储在所述存储器302中，当被所述处理器301执行时，执行如图1所示实施例中的CT图像迭代重建方法。
上述计算机设备具体细节可以对应参阅图1所示的实施例中对应的相关描述和效果进行理解，此处不再赘述。
本领域技术人员可以理解，实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)、随机存储记忆体(Random AccessMemory，RAM)、快闪存储器(Flash Memory)、硬盘(Hard Disk Drive，缩写：HDD)或固态硬盘(Solid-State Drive，SSD)等；所述存储介质还可以包括上述种类的存储器的组合。
虽然结合附图描述了本发明的实施例，但是本领域技术人员可以在不脱离本发明的精神和范围的情况下作出各种修改和变型，这样的修改和变型均落入由所附权利要求所限定的范围之内。

Claims (10)

1.一种CT图像迭代重建方法，其特征在于，包括：
获取目标对象在初始投影角度下的先验图像值及径迹矩阵，所述先验图像值为所述目标对象在初始投影角度下每一个点的像素值的集合，所述径迹矩阵为在所述初始投影角度下射线源发出的每条射线穿过各个先验图像值的长度；
基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值；
基于所述投影理论值和在所述初始投影角度下的测量值确定偏差值；
基于所述偏差值更新所述先验图像值，得到更新后的先验图像值；
基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值；
将旋转单位投影角度后的图像值作为旋转单位投影角度的先验图像，并将所述单位投影角度作为初始投影角度，返回执行基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值的步骤至基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值的步骤，直至当前旋转角度达到预设角度值；
基于所述当前旋转角度中各单位旋转角度对应的旋转单位投影角度后的图像值得到CT图像重建结果。
2.根据权利要求1所述的方法，其特征在于，在基于所述当前旋转角度中各单位旋转角度对应的旋转单位投影角度后的图像值得到CT图像重建结果之前，所述方法还包括：
基于相邻两次单位旋转角度对应的旋转单位投影角度后的图像值判断迭代过程是否满足预设收敛条件；
若迭代过程满足预设收敛条件，则执行基于所述当前旋转角度中各单位旋转角度对应的旋转单位投影角度后的图像值得到CT图像重建结果的步骤。
3.根据权利要求2所述的方法，其特征在于，所述方法还包括：
若迭代过程不满足预设收敛条件，则返回执行基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值的步骤至基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值的步骤。
4.根据权利要求1所述的方法，其特征在于，所述获取目标对象在初始投影角度下的先验图像值及径迹矩阵，包括：
获取所述目标对象的投影数据和初始投影角度；
基于所述投影数据和预设的重建算法得到所述先验图像值；
基于所述先验图像值得到初始投影角度下的径迹矩阵。
5.根据权利要求4所述的方法，其特征在于，所述基于所述投影数据和预设的重建算法得到所述先验图像值，包括：
对所述投影数据进行加权和滤波处理，得到第一图像值；
对所述第一图像值进行加权反投影得到所述先验图像值。
6.根据权利要求4所述的方法，其特征在于，所述基于所述先验图像值得到初始投影角度下的径迹矩阵，包括：
基于CT射线源和像素点得到穿过所述先验图像值的射线长度；
基于各所述射线长度得到所述径迹矩阵。
7.根据权利要求1所述的方法，其特征在于，所述基于所述偏差值更新所述先验图像值，得到更新后的先验图像值，包括：
基于所述偏差值和所述径迹矩阵中的每一条射线得到各所述先验图像值的校正值；
基于各所述校正值得到所述更新后的先验图像值。
8.一种CT图像迭代重建装置，其特征在于，包括：
获取模块，用于获取目标对象在初始投影角度下的先验图像值及径迹矩阵，所述先验图像值为所述目标对象在初始投影角度下每一个点的像素值的集合，所述径迹矩阵为在所述初始投影角度下射线源发出的每条射线穿过各个先验图像值的长度；
投影理论值模块，用于基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值；
偏差模块，用于基于所述投影理论值和在所述初始投影角度下的测量值确定偏差值；
更新模块，用于基于所述偏差值更新所述先验图像值，得到更新后的先验图像值；
第一旋转模块，用于基于单位旋转角度和更新后的先验图像值，得到旋转角度后的图像值；
第二旋转模块，用于将旋转单位投影角度后的图像值作为旋转单位投影角度的先验图像，并将所述单位投影角度作为初始投影角度，返回执行基于所述先验图像值和径迹矩阵计算得到在所述初始投影角度下的投影理论值的步骤至基于单位旋转角度和更新后的先验图像值，得到旋转单位投影角度后的图像值的步骤，直至当前旋转角度达到预设角度值；
输出模块，用于基于所述当前旋转角度中各单位旋转角度对应的旋转单位投影角度后的图像值得到CT图像重建结果。
9.一种计算机设备，其特征在于，包括：至少一个处理器；以及与所述至少一个处理器通信连接的存储器；其中，所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器执行如权利要求1-7任一所述的CT图像迭代重建方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1-7中任一项所述的CT图像迭代重建方法的步骤。
CN202210146430.9A 2022-02-17 2022-02-17 一种ct图像迭代重建方法、装置及计算机设备 Active CN116664707B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210146430.9A CN116664707B (zh) 2022-02-17 2022-02-17 一种ct图像迭代重建方法、装置及计算机设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210146430.9A CN116664707B (zh) 2022-02-17 2022-02-17 一种ct图像迭代重建方法、装置及计算机设备

Publications (2)

Publication Number Publication Date
CN116664707A CN116664707A (zh) 2023-08-29
CN116664707B true CN116664707B (zh) 2024-02-13

Family

ID=87726549

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210146430.9A Active CN116664707B (zh) 2022-02-17 2022-02-17 一种ct图像迭代重建方法、装置及计算机设备

Country Status (1)

Country Link
CN (1) CN116664707B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102243766A (zh) * 2011-06-08 2011-11-16 无锡引速得信息科技有限公司 一种ct图像重建加速的新方法
CN103810733A (zh) * 2014-02-28 2014-05-21 南方医科大学 一种稀疏角度x射线ct图像的统计迭代重建方法
CN106355624A (zh) * 2016-08-26 2017-01-25 中国科学院武汉物理与数学研究所 一种基于变采样率的肺部超极化气体动态成像方法
CN110702707A (zh) * 2019-10-16 2020-01-17 四川轻化工大学 一种获取核废物桶层析γ扫描图像的方法
CN111476854A (zh) * 2020-04-03 2020-07-31 上海联影医疗科技有限公司 图像重建方法、装置、终端及存储介质
WO2022000192A1 (zh) * 2020-06-29 2022-01-06 深圳先进技术研究院 一种ct图像的构建方法、ct设备以及存储介质

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7042974B2 (en) * 2004-05-20 2006-05-09 Eastman Kodak Company Method for reconstructing tomographic images
US9261467B2 (en) * 2013-06-14 2016-02-16 General Electric Company System and method of iterative image reconstruction for computed tomography
CN105593905B (zh) * 2013-09-30 2019-07-23 皇家飞利浦有限公司 针对完全3d迭代ct重建中的图像质量优化用于对正则化参数的局部调节的方法
US9848844B2 (en) * 2015-10-09 2017-12-26 Carestream Health, Inc. Iterative reconstruction process

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102243766A (zh) * 2011-06-08 2011-11-16 无锡引速得信息科技有限公司 一种ct图像重建加速的新方法
CN103810733A (zh) * 2014-02-28 2014-05-21 南方医科大学 一种稀疏角度x射线ct图像的统计迭代重建方法
CN106355624A (zh) * 2016-08-26 2017-01-25 中国科学院武汉物理与数学研究所 一种基于变采样率的肺部超极化气体动态成像方法
CN110702707A (zh) * 2019-10-16 2020-01-17 四川轻化工大学 一种获取核废物桶层析γ扫描图像的方法
CN111476854A (zh) * 2020-04-03 2020-07-31 上海联影医疗科技有限公司 图像重建方法、装置、终端及存储介质
WO2022000192A1 (zh) * 2020-06-29 2022-01-06 深圳先进技术研究院 一种ct图像的构建方法、ct设备以及存储介质

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
An improved algebraic reconstruction technique for reconstructing tomographic gamma scanning image;Honglong Zheng;《ElsevierNuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment》;第906卷;77-82 *
An Improved OSEM Iterative Reconstruction Algorithm for Transmission Tomographic Gamma Scanning;Aijing He;《Applied Radiation and Isotopes》;第142卷(第01期);51-55 *
层析γ扫描透射图像迭代重建算法的改进研究;何艾静;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑 》(第02期);C040-63 *

Also Published As

Publication number Publication date
CN116664707A (zh) 2023-08-29

Similar Documents

Publication Publication Date Title
US8705828B2 (en) Methods and apparatus for super resolution scanning for CBCT system and cone-beam image reconstruction
US8218715B2 (en) Methods and systems to facilitate correcting gain fluctuations in iterative image reconstruction
US5251128A (en) Motion artifact reduction in projection imaging
JP5348855B2 (ja) 対象の画像再構成方法およびその方法を実施するための装置
CN104318536B (zh) Ct图像的校正方法及装置
JP6066596B2 (ja) Ｘ線撮像における散乱補正の方法及びシステム
CN106232007B (zh) X射线ct装置和处理装置
CN104644200B (zh) 减少计算机断层扫描图像重构中伪像的方法和装置
US8885903B2 (en) Method and apparatus for statistical iterative reconstruction
JP6050029B2 (ja) 断層写真法画像の再構成の方法及びシステム
WO2017128892A1 (zh) 计算机断层成像伪影校正方法及系统
US20240029207A1 (en) Systems and methods for adaptive blending in computed tomography imaging
CN113344876A (zh) 一种ct和cbct间可变形配准方法
CN116664707B (zh) 一种ct图像迭代重建方法、装置及计算机设备
CN110063739B (zh) Spect和pet的图像校正方法、装置以及电子设备
US6307908B1 (en) System and method for data interpolation in a multislice x-ray computed tomography system
US7379525B2 (en) Method and system for efficient helical cone-beam reconstruction
JP7403585B2 (ja) 計算機式断層写真法画像再構成のためのシステム及び方法
US20230145920A1 (en) Systems and methods for motion detection in medical images
CN111839574B (zh) Ct超低剂量自动三维定位扫描方法及系统
CN106651985A (zh) Ct图像的重建方法和装置
CN107341836B (zh) 一种ct螺旋扫描图像重建方法及装置
CN114903510B (zh) 一种双能ct成像方法、装置及计算机设备
CN110599563B (zh) 一种自适应nlm修正的ct重建方法及装置
TW201822717A (zh) 斷層合成影像邊緣假影抑制方法

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