CN113017668A - 一种基于非柱面探测器的ct图像重建方法及ct扫描仪 - Google Patents

一种基于非柱面探测器的ct图像重建方法及ct扫描仪 Download PDF

Info

Publication number
CN113017668A
CN113017668A CN202110199513.XA CN202110199513A CN113017668A CN 113017668 A CN113017668 A CN 113017668A CN 202110199513 A CN202110199513 A CN 202110199513A CN 113017668 A CN113017668 A CN 113017668A
Authority
CN
China
Prior art keywords
cylindrical
detector
channel
cone angle
central
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
Application number
CN202110199513.XA
Other languages
English (en)
Other versions
CN113017668B (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.)
FMI Technologies Inc
Original Assignee
FMI Technologies Inc
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 FMI Technologies Inc filed Critical FMI Technologies Inc
Priority to CN202110199513.XA priority Critical patent/CN113017668B/zh
Publication of CN113017668A publication Critical patent/CN113017668A/zh
Application granted granted Critical
Publication of CN113017668B publication Critical patent/CN113017668B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pulmonology (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明属于医学影像技术领域，具体涉及一种基于非柱面探测器的CT图像重建方法及CT扫描仪。其中，CT图像重建方法，包括以下步骤：S1、采用数据插值方法将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面数据；S2、根据对应于柱面坐标系的柱面数据进行CT图像重建。本发明的CT图像重建方法，针对非柱面探测器，采用数据插值方法将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面数据，能够保证得到的采样点与非柱面探测器的采样点尽量接近，减少了空间分辨率的损失和混叠伪影的增加。

Description

一种基于非柱面探测器的CT图像重建方法及CT扫描仪
技术领域
本发明属于医学影像技术领域，具体涉及一种基于非柱面探测器的CT图像重建方法及CT扫描仪。
背景技术
第三代CT系统的主要组成部分包括Tube(球管)、Collimator(限束器)和Detector(探测器)。球管发出X光，经过限束器限制形成锥形的光束，锥形光束照射到探测器上，经过探测器转化为电信号并经过数据采集和转换单元转换为数字信息存储在图像处理系统中，图像处理系统经过一系列的校正算法和图像重建算法生成图像，显示在显示器上。从成本和技术成熟程度来说，主流探测器都采用了由很多个探测器模块排列在一个弧形或者多边形的面上来形成整个探测器，每个探测器模块中规则排列一个探测器单元形成的矩阵。
现有探测器像素设计一般都为在扫描平面内(X方向)探测器像素较大，而Z方向探测器像素较小。在扫描平面内像素较大，是为了兼顾几何探测效率和空间分辨率而设计。扫描平面内可以使用探测器1/4偏移做到主射线和相隔180度机架旋转角的共轭射线结合做到2倍采样率；或者，进一步和X射线焦点平面内摆动技术结合做到4倍采样率；这些技术可以做到扫描平面内的奈奎斯特采样频率大幅提升，增加了平面内的空间分辨率，并减少了混叠伪影。
随着探测器技术的发展，出现了排数达到256排，甚至达到320排，Z方向扫描范围达到160mm的探测器。这些探测器形状复杂，都属于非柱面探测器，其形状多接近球面探测器或者其它扩展形状，如公开号为US10492746的美国专利文献公开的多个球面组成的探测器，以及存在中心多排是柱面探测器而外围是球面或接近球面的探测器；它们一个共同的特点是每个像素在通道方向以及Z方向的角度都是不均匀的；分布的复杂性造成了重建算法的复杂程度。
然而，现有标准的3D重建算法都是基于柱面探测器进行的；非柱面探测器的几何分布比较复杂，直接从其进行重建，计算难度很高。因此，如何针对非柱面探测器进行CT图像重建是本领域亟需解决的难题。
发明内容
基于现有技术中存在的上述缺点和不足，本发明的目的之一是至少解决现有技术中存在的上述问题之一或多个，换言之，本发明的目的之一是提供满足前述需求之一或多个的一种基于非柱面探测器的CT图像重建方法及CT扫描仪。
为了达到上述发明目的，本发明采用以下技术方案：
一种基于非柱面探测器的CT图像重建方法，包括以下步骤：
S1、采用数据插值方法将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面数据；
S2、根据对应于柱面坐标系的柱面数据进行CT图像重建。
作为优选方案，对于非柱面探测器，假设探测器的排数为M，每一排的通道数为N；
各个通道的空间坐标标记为X方向的扇角坐标γ(m，n)和Z方向的锥角坐标α(m，n)，γ(m，n)表示X方向各通道与中心通道ISOch之间的夹角，α(m，n)表示Z方向各排与探测器Z方向中心线之间的夹角，m为排的标号，取值为1～M之间的整数，n为排内的通道标号，取值为1～N之间的整数。
作为优选方案，所述步骤S1，包括：
计算中心排的张角为γ(middle，N)-γ(middle，1)，计算通道的平均张角
Figure BDA0002948066120000021
其中，middle表示中心排；
以通道的平均张角Δγ和中心通道ISOch为基准，计算每个通道的柱面扇角坐标为γcyl(m，n)＝(n-ISOch)×Δγ；其中，中心通道
Figure BDA0002948066120000031
或ISOch＝i，i为离中心通道最近的通道标号；
Figure BDA0002948066120000032
表示1/4探测器偏移；
根据每个通道的扇角坐标γ(m，n)至柱面扇角坐标γcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面扇角数据。
作为优选方案，所述步骤S1，包括：
以探测器的中心排的各通道的扇角作为每排对应通道的柱面扇角坐标，即γcyl(m，n)＝γ(middle，n)；
根据扇角坐标γ(m，n)至柱面扇角坐标γcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面扇角数据。
作为优选方案，所述步骤S1，包括：
以中心通道ISOch的总锥角除以排数，计算通道的平均锥角
Figure BDA0002948066120000033
其中，中心通道
Figure BDA0002948066120000034
或ISOch＝i，i为离中心通道最近的通道标号；
Figure BDA0002948066120000035
表示1/4探测器偏移；
以通道的平均锥角Δα和探测器Z方向中心线为基准，计算每排的柱面锥角坐标
Figure BDA0002948066120000036
根据锥角坐标α(m，n)至柱面锥角坐标αcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面锥角数据。
作为优选方案，所述步骤S1，包括：
以探测器的各排中心通道的锥角作为每排的柱面锥角坐标，即αcyl(m，n)＝α(m，ISOch)；
根据锥角坐标α(m，n)至柱面锥角坐标αcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面锥角数据。
作为优选方案，所述步骤S1，包括：
以探测器中心范围的部分通道对应的平均锥角作为每个通道的柱面锥角坐标：
Figure BDA0002948066120000041
其中，中心通道
Figure BDA0002948066120000042
或ISOch＝i，i为离中心通道最近的通道标号；
Figure BDA0002948066120000043
表示1/4探测器偏移；n1和n2为探测器中心范围的部分通道对应的起始通道标号和终止通道标号；
根据锥角坐标α(m，n)至柱面锥角坐标αcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面锥角数据。
作为优选方案，所述数据插值方法为线性插值或拉格朗日插值。
作为优选方案，所述步骤S2中，采用滤波反投影算法进行CT图像重建。
本发明还提供一种CT扫描仪，采用如上任一方案所述的CT图像重建方法。
本发明与现有技术相比，有益效果是：
本发明的CT图像重建方法，针对非柱面探测器，采用数据插值方法将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面数据，能够保证得到的采样点与非柱面探测器的采样点尽量接近，减少了空间分辨率的损失和混叠伪影的增加。
附图说明
图1是本发明实施例的CT图像重建方法得到的MTF50％对比图；
图2是本发明实施例的CT图像重建方法得到的MTF10％对比图；
图3是本发明实施例的CT图像重建方法得到的CT图像的对比图。
具体实施方式
以下通过具体实施例对本发明的技术方案作进一步解释说明。
本发明实施例的基于非柱面探测器的CT图像重建方法，包括以下步骤：
S1、采用数据插值方法将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面数据；
S2、根据对应于柱面坐标系的柱面数据进行CT图像重建。
其中，在从非柱面数据转换为柱面数据的过程中，因为两者的三维坐标不可能完全相同，因此需要进行数据插值，数据插值一般会采用线性插值或者多点的拉格朗日插值，这种插值相当于低通滤波器，会造成空间分辨率的损失和混叠的增加。
对于非柱面探测器，假设探测器的排数为M，每一排的通道数为N；
各个通道的空间坐标采用X方向的扇角坐标γ(m，n)和Z方向的锥角坐标α(m，n)，m为排的标号，取值为1～M之间的整数，n为排内的通道标号，取值为1～N之间的整数。
对于X方向的扇角坐标，上述步骤S1可以采用以下两种方案：
第一种方案：
计算探测器的中心排张开的扇角(即张角)再除以通道数得到Δγ，然后在中心通道两边设置柱面探测器的通道，使得每两个相邻通道之间的扇角为Δγ；在锥角方向可以采用类似的思想来进行；如此，能够保证得到的采样点与非柱面探测器的采样点尽量接近，因而能够减少插值造成的损失，减少了空间分辨率的损失和混叠伪影的增加。
具体地，以探测器的中心排为基准，计算中心排的张角为γ(middle，N)-γ(middle，1)，计算通道的平均张角
Figure BDA0002948066120000061
其中，middle表示中心排；
以通道的平均张角Δγ和中心通道ISOch为基准，计算每个通道的柱面扇角坐标为γcyl(m，n)＝(n-ISOch)×Δγ；其中，中心通道
Figure BDA0002948066120000062
或ISOch＝i，i为离中心通道最近的通道标号；
Figure BDA0002948066120000063
表示1/4探测器偏移，其目的是使得互相共轭的射线能够交叉增加采样率；其它特殊目的如飞焦点的ISOCh(中心通道)也有不同的选择方式。
根据每个通道的扇角坐标γ(m，n)至柱面扇角坐标γcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面扇角数据。
第二种方案：
以探测器的中心排的各通道的扇角作为每排对应通道的柱面扇角坐标，即γcyl(m，n)＝γ(middle，n)；
根据扇角坐标γ(m，n)至柱面扇角坐标γcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面扇角数据。
对于Z方向的锥角坐标，上述步骤S1可以如下三种方案：
第三种方案：
以中心通道ISOch的总锥角除以排数，计算得到通道的平均锥角：
Figure BDA0002948066120000064
其中，中心通道
Figure BDA0002948066120000065
或ISOch＝i，i为离中心通道最近的通道标号；
Figure BDA0002948066120000066
表示1/4探测器偏移；
以通道的平均锥角Δα和探测器Z方向中心线为基准，计算每排的柱面锥角坐标
Figure BDA0002948066120000071
根据锥角坐标α(m，n)至柱面锥角坐标αcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面锥角数据。
第四种方案：
以探测器的各排中心通道的锥角作为每排的柱面锥角坐标，即αcyl(m，n)＝α(m，ISOch)；
根据锥角坐标α(m，n)至柱面锥角坐标αcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面锥角数据。
第五种方案：
以探测器中心范围的部分通道对应的平均锥角作为每个通道的柱面锥角坐标：
Figure BDA0002948066120000072
其中，中心通道
Figure BDA0002948066120000073
或ISOch＝i，i为离中心通道最近的通道标号；
Figure BDA0002948066120000074
表示1/4探测器偏移；n1和n2为探测器中心范围的部分通道对应的起始通道标号和终止通道标号；
根据锥角坐标α(m，n)至柱面锥角坐标αcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面锥角数据。
针对非柱面探测器，扇角的坐标转换采用上述第一种方案和第二种方案中的任一种，锥角的坐标转换采用上述第三种方案至第五种方案中的任一种；按照上述方案计算的新的采样点和非柱面探测器的采样点足够接近，都能够达到减少插值损失以减少空间分辨率损失增加混叠伪影的风险。
具体地，本发明实施例的图像重建过程，按照以下步骤来进行操作：
(1)采集到的原始数据进行预处理，得到物体沿每个射线路径上的总吸收量P(m，n)；
(2)计算每个射线路径对应的非球面探测器的各个像素的扇角γ(m，n)和锥角α(m，n)；
(3)根据本发明实施例上述的数种方案选取一种，例如选取最中间一排的探测器像素所对应γ,计算出一个平均的Δγ；然后在中心通道两边依次计算柱面探测器像素对应的扇角γcyl(m，n)；对于锥角方向的处理，如果原始的探测器其排列使得各排像素在Z方向是排列整齐的，可以对锥角方向不用插值；如果原始的探测器其排列在Z方向是不对齐的，则可以按照本发明实施例上述方案所述的方法计算αcyl(m，n)，进行插值；
(4)对步骤(1)得到的数据P(m，n)按照从非柱面坐标系(γ(m，n)，α(m，n))到柱面坐标系(γcyl(m，n)，αcyl(m，n))的线性插值或双线性插值得到Pcyl(m，n)。
(5)按照现有标准的柱面探测器的重排以及滤波反投影算法进行重建，得到CT图像。
以通道方向插值为例。通道方向的插值影响平面内的分辨率，平面内的分辨率通常使用MTF(Modulation Transfer Function，调制传递函数)来量化。如图1和2所示，显示了分别使用本发明提到的四种插值方法重建得到的MTF对比结果：图中的a采用中心排的扇角作为柱面的扇角(第二种方案)；b采用中心排所有通道的平均扇角为Δγ(第一种方案)；c采用中心模块的中心像素的张角为Δγ；d采用中心通道ISOch所在通道的张角Δγ；图中横坐标d表示距图像上ISO中心的距离。由图1和2可知，a和b的表现较为接近，c和d的表现较为接近，四种方法在ISO中心的表现较为相同，但随着距ISO中心的距离变大，a和b的效果要远远优于c和d。相应地，重建得到的CT图像如图3所示，重建图像DFOV(Display Field OfView，重建视野)为150mm，从图像可以看出，a和b表现较为接近，在整个DFOV范围内图像均很清晰；c和d表现较为接近，随着距ISO中心的距离增大，图像逐渐变模糊。
本发明实施例还提供一种CT扫描仪，采用如上任一方案所述的CT图像重建方法。实现配置有非柱面探测器的CT扫描仪的图像重建更加高效。
以上所述仅是对本发明的优选实施例及原理进行了详细说明，对本领域的普通技术人员而言，依据本发明提供的思想，在具体实施方式上会有改变之处，而这些改变也应视为本发明的保护范围。

Claims (10)

1.一种基于非柱面探测器的CT图像重建方法，其特征在于，包括以下步骤：
S1、采用数据插值方法将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面数据；
S2、根据对应于柱面坐标系的柱面数据进行CT图像重建。
2.根据权利要求1所述的CT图像重建方法，其特征在于，对于非柱面探测器，假设探测器的排数为M，每一排的通道数为N；
各个通道的空间坐标标记为X方向的扇角坐标γ(m，n)和Z方向的锥角坐标α(m，n)，γ(m，n)表示X方向各通道与中心通道ISOch之间的夹角，α(m，n)表示Z方向各排与探测器Z方向中心线之间的夹角，m为排的标号，取值为1～M之间的整数，n为排内的通道标号，取值为1～N之间的整数。
3.根据权利要求2所述的CT图像重建方法，其特征在于，所述步骤S1，包括：
计算中心排的张角为γ(middle，N)-γ(middle，1)，计算通道的平均张角
Figure FDA0002948066110000011
其中，middle表示中心排；
以通道的平均张角Δγ和中心通道ISOch为基准，计算每个通道的柱面扇角坐标为γcyl(m，n)＝(n-ISOch)×Δγ；其中，中心通道
Figure FDA0002948066110000012
或ISOch＝i，i为离中心通道最近的通道标号；
Figure FDA0002948066110000013
表示1/4探测器偏移；
根据每个通道的扇角坐标γ(m，n)至柱面扇角坐标γcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面扇角数据。
4.根据权利要求2所述的CT图像重建方法，其特征在于，所述步骤S1，包括：
以探测器的中心排的各通道的扇角作为每排对应通道的柱面扇角坐标，即γcyl(m，n)＝γ(middle，n)；
根据扇角坐标γ(m，n)至柱面扇角坐标γcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面扇角数据。
5.根据权利要求3或4所述的CT图像重建方法，其特征在于，所述步骤S1，包括：
以中心通道ISOch的总锥角除以排数，计算通道的平均锥角
Figure FDA0002948066110000021
其中，中心通道
Figure FDA0002948066110000022
或ISOch＝i，i为离中心通道最近的通道标号；
Figure FDA0002948066110000023
表示1/4探测器偏移；
以通道的平均锥角Δα和探测器Z方向中心线为基准，计算每排的柱面锥角坐标
Figure FDA0002948066110000024
根据锥角坐标α(m，n)至柱面锥角坐标αcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面锥角数据。
6.根据权利要求3或4所述的CT图像重建方法，其特征在于，所述步骤S1，包括：
以探测器的各排中心通道的锥角作为每排的柱面锥角坐标，即αcyl(m，n)＝α(m，ISOch)；
根据锥角坐标α(m，n)至柱面锥角坐标αcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面锥角数据。
7.根据权利要求3或4所述的CT图像重建方法，其特征在于，所述步骤S1，包括：
以探测器中心范围的部分通道对应的平均锥角作为每个通道的柱面锥角坐标：
Figure FDA0002948066110000031
其中，中心通道
Figure FDA0002948066110000032
或ISOch＝i，i为离中心通道最近的通道标号；
Figure FDA0002948066110000033
表示1/4探测器偏移；n1和n2为探测器中心范围的部分通道对应的起始通道标号和终止通道标号；
根据锥角坐标α(m，n)至柱面锥角坐标αcyl(m，n)的数据插值，以将非柱面探测器采集的非柱面数据转换为对应于柱面探测器的柱面锥角数据。
8.根据权利要求1所述的CT图像重建方法，其特征在于，所述数据插值方法为线性插值或拉格朗日插值。
9.根据权利要求1所述的CT图像重建方法，其特征在于，所述步骤S2中，采用滤波反投影算法进行CT图像重建。
10.一种CT扫描仪，其特征在于，采用如权利要求1-9任一项所述的CT图像重建方法。
CN202110199513.XA 2021-02-23 2021-02-23 一种基于非柱面探测器的ct图像重建方法及ct扫描仪 Active CN113017668B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110199513.XA CN113017668B (zh) 2021-02-23 2021-02-23 一种基于非柱面探测器的ct图像重建方法及ct扫描仪

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110199513.XA CN113017668B (zh) 2021-02-23 2021-02-23 一种基于非柱面探测器的ct图像重建方法及ct扫描仪

Publications (2)

Publication Number Publication Date
CN113017668A true CN113017668A (zh) 2021-06-25
CN113017668B CN113017668B (zh) 2022-06-14

Family

ID=76460993

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110199513.XA Active CN113017668B (zh) 2021-02-23 2021-02-23 一种基于非柱面探测器的ct图像重建方法及ct扫描仪

Country Status (1)

Country Link
CN (1) CN113017668B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1028684A (ja) * 1996-04-10 1998-02-03 Siemens Ag スパイラル動作におけるマルチライン検出器コンピュータ断層撮影装置に対する画像再構成方法
JP2005334230A (ja) * 2004-05-26 2005-12-08 Ge Medical Systems Global Technology Co Llc 放射線ｃｔ装置およびそれを用いた画像再構成方法
US20080056432A1 (en) * 2006-08-30 2008-03-06 General Electric Company Reconstruction of CT projection data
CN103300874A (zh) * 2012-03-14 2013-09-18 Ge医疗系统环球技术有限公司 图像生成方法和装置、放射线断层摄影装置及程序
CN103961122A (zh) * 2013-01-31 2014-08-06 通用电气公司 用于非等γ角CT系统中数据转换的方法和装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1028684A (ja) * 1996-04-10 1998-02-03 Siemens Ag スパイラル動作におけるマルチライン検出器コンピュータ断層撮影装置に対する画像再構成方法
JP2005334230A (ja) * 2004-05-26 2005-12-08 Ge Medical Systems Global Technology Co Llc 放射線ｃｔ装置およびそれを用いた画像再構成方法
US20080056432A1 (en) * 2006-08-30 2008-03-06 General Electric Company Reconstruction of CT projection data
CN103300874A (zh) * 2012-03-14 2013-09-18 Ge医疗系统环球技术有限公司 图像生成方法和装置、放射线断层摄影装置及程序
US20130243152A1 (en) * 2012-03-14 2013-09-19 Ge Medical Systems Global Technology Company, Llc Image generating method, image generating apparatus and radiation tomographic imaging apparatus, and program therefor
CN103961122A (zh) * 2013-01-31 2014-08-06 通用电气公司 用于非等γ角CT系统中数据转换的方法和装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨昆等: "数据插值对正电子发射断层成像设备的图像重建影响的研究", 《物理学报》 *
杨芳等: "一种基于测量空间样条插值扩展的扩散光学层析图像改善方法", 《光学学报》 *

Also Published As

Publication number Publication date
CN113017668B (zh) 2022-06-14

Similar Documents

Publication Publication Date Title
CN102413767B (zh) 用于微分相位对比扇形线束ct、锥形线束ct以及混合锥形线束ct的方法和装置
US8290233B2 (en) X-ray computed tomography apparatus and image processing method
US8116426B2 (en) Computed tomography device and method using circular-pixel position-adaptive interpolation
US7656991B2 (en) Image reconstruction method and X-ray CT apparatus
US7864916B2 (en) X-ray computer tomography apparatus
CN101036164A (zh) 采用与体素相关插值的图像重建
Hasan et al. Hybrid-collaborative Noise2Noise denoiser for low-dose CT images
EP1031943A2 (en) Efficient cone-beam reconstruction system using circle-and-line orbit data
IL100312A (en) Method and device for recreating a three-dimensional computerized tomography image of a complete data set on a conical beam projection
JP2009160408A (ja) 磁気共鳴イメージングにおける歪み補正方法及び装置
CN103065279A (zh) 圆轨道锥束计算机断层摄影（ct）中用于大幅度地减轻伪影的方法以及系统
CN103679706B (zh) 一种基于图像各向异性边缘检测的ct稀疏角度重建方法
CN110363825B (zh) 一种z向高分辨率的CT扫描方式和图像重建方法
CN103300874B (zh) 图像生成方法和装置、放射线断层摄影装置及程序
CN114886444B (zh) 一种cbct成像重建方法
JP2022520977A (ja) 医用画像処理システムとその画像再構成方法のためのコリメータ
CN108283503B (zh) 一种ct机、扫描方法及图像重建方法
JP2004166923A (ja) Ｘ線コンピューター断層撮影装置
CN113017668B (zh) 一种基于非柱面探测器的ct图像重建方法及ct扫描仪
US7242004B2 (en) Image correction method, image correction apparatus, and image correction program
CN114359126A (zh) 三维图像的几何校正方法、装置、电子设备及存储介质
CN102599887B (zh) 一种基于螺旋扫描轨道的光学投影断层成像方法
US20120093279A1 (en) Medical image processing apparatus, x-ray computed tomography apparatus, and medical image processing method
CN110176045A (zh) 一种由单能ct图像生成双能ct图像的方法
JP3583567B2 (ja) Ｘ線断層撮像装置

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