CN110113994A - 用于计算机断层摄影的数据采集 - Google Patents
用于计算机断层摄影的数据采集 Download PDFInfo
- Publication number
- CN110113994A CN110113994A CN201780048118.3A CN201780048118A CN110113994A CN 110113994 A CN110113994 A CN 110113994A CN 201780048118 A CN201780048118 A CN 201780048118A CN 110113994 A CN110113994 A CN 110113994A
- Authority
- CN
- China
- Prior art keywords
- sampling
- pixel
- data sampling
- time
- detector
- 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.)
- Pending
Links
- 238000002591 computed tomography Methods 0.000 title claims abstract description 81
- 238000005070 sampling Methods 0.000 claims abstract description 113
- 238000000034 method Methods 0.000 claims abstract description 47
- 230000005012 migration Effects 0.000 claims abstract description 33
- 238000013508 migration Methods 0.000 claims abstract description 33
- 238000003384 imaging method Methods 0.000 claims abstract description 29
- 238000005259 measurement Methods 0.000 claims description 32
- 238000004590 computer program Methods 0.000 claims description 29
- 230000004907 flux Effects 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 19
- 238000012545 processing Methods 0.000 description 18
- 238000005516 engineering process Methods 0.000 description 14
- 230000000694 effects Effects 0.000 description 10
- 238000001228 spectrum Methods 0.000 description 8
- 229910052704 radon Inorganic materials 0.000 description 7
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000013480 data collection Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010412 perfusion Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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/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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4007—Arrangements for generating radiation specially adapted for radiation diagnosis characterised by using a plurality of source units
- A61B6/4014—Arrangements for generating radiation specially adapted for radiation diagnosis characterised by using a plurality of source units arranged in multiple source-detector units
-
- 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
- A61B6/4233—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using matrix detectors
-
- 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/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
-
- 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
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/24—Measuring radiation intensity with semiconductor detectors
- G01T1/247—Detector read-out circuitry
-
- 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
-
- 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
-
- 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/54—Control of apparatus or devices for radiation diagnosis
- A61B6/545—Control of apparatus or devices for radiation diagnosis involving automatic set-up of acquisition parameters
-
- 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
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Radiology & Medical Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pulmonology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Mathematical Physics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明提供了一种用于计算机断层扫描(CT)的数据采集方法及相应的系统,基于具有多个像素(13)的探测器(12)的X射线成像系统(10),其中,数据采样是利用由探测器的相邻像素获取的测量值之间的时间偏移来执行的。
Description
技术领域
本发明所提出的技术通常涉及X射线成像,更具体地涉及用于计算机断层摄影(CT)的数据采集方法以及相应的系统,计算机程序和计算机程序产品,以及相应的计算机断层摄影(CT)系统。
背景技术
在医学应用和非破坏性测试中已经使用了诸如X射线成像等射线照相成像技术。
通常,X射线成像系统包括X射线源和X射线探测器系统。X射线源发射X射线,其穿过待成像的对象或目标,然后通过X射线探测器系统进行配准。由于一些材料比其它材料吸收更大部分的X射线,所以对象或目标的图像由此形成。
以示例性的整体X射线成像系统的简要概述开始是非常有用的,参考图10。在这个非限制性实例中,X射线成像系统10基本上包括X射线源11、X射线探测器12和相关的图像处理装置30。总的来说,X射线探测器12是被配置为用于对来自X射线源11的射线进行配准,该X射线源11可能已经被可选的X射线光学系统聚焦,并穿过物体或对象或其部分。X射线探测器12可通过合适的模拟处理和读出电子器件(其可以集成在X射线探测器12或探测器系统12中)连接到图像处理装置15,使得能通过图像处理装置15进行图像处理和/或图像重建。
通常使用的X射线成像系统的实例是计算机断层摄影(CT)系统,其可包括产生X射线的扇形或锥形束的X射线源,以及相应的X射线探测器或探测器系统,用于配准穿过患者或物体传输的X射线的部分。X射线源和探测器系统通常安装在围绕成像物体旋转的台架中
CT场的发展使得台架旋转速度变得越来越高,探测器的空间分辨率也更高;这样就增加了对足够角度采样的要求。CT中的角度采样不足导致混淆,这通过从锐边引入条纹状伪像来破坏重建图像。用于CT的过采样的几种方法,例如飞焦点(FFS)和四分之一探测器偏移(QDO)已经开发出来,并主要用作标准方法,以便降低混叠的风险。采样的数量和正弦图被采样的网格确定混叠的风险。
发明内容
期望找到一种用于管理CT系统中混叠的改进机制。
本发明的一个目的是提供一种用于计算机断层扫描(CT)的数据采集方法。
本发明的另一个目的是提供一种用于计算机断层扫描(CT)的数据采集系统。
本发明的另一目的是提供一种计算机断层摄影(CT)系统。
本发明的另一目的是提供一种用于计算机断层扫描(CT)的当由处理器执行时能进行数据采集的计算机程序。
本发明的另一目的是提供一种包括计算机可读介质的计算机程序产品,所述计算机可读介质上存储有这样的计算机程序。
通过本发明的具体实施例满足了这些和其他目的。
根据本发明所提出的技术的第一个方面,提供了一种用于计算机断层扫描(CT)的数据采集方法,基于具有多个像素的探测器的X射线成像系统,其中使用由探测器的相邻像素获取的测量值之间的时间偏移来执行数据采样。
以这种方式,诸如六边形采样图案的定制采样图案,可减少混叠或至少控制所产生的混叠的类型,而不需要依赖于标准过采样技术。
本发明提供了一种获得高效采样的鲁棒(robust)方式,能够降低CT图像中混叠伪像的风险。
根据本发明所提出的技术的第二个方面,提供了一种用于计算机断层扫描(CT)的数据采集系统,基于具有多个像素的探测器的X射线成像系统,其中,所述系统被配置为对由所述探测器的相邻像素获取的测量之间的时间偏移进行数据采样。
根据本发明所提出的技术的第三个方面,提供了一种计算机断层摄影(CT)系统,它包括根据第二个方面的系统。
根据本发明所提出的技术的第四个方面,提供了一种用于计算机断层扫描(CT)的当由处理器执行时能进行数据采集的计算机程序,采用具有多个像素的探测器的X射线成像系统,其中所述计算机程序包括指令,所述处理器在由所述处理器执行时使所述处理器执行与所述探测器的相邻像素所获取的测量值之间的时间偏移量的数据采样。
根据本发明所提出的技术的第五个方面,提供了一种计算机程序产品,该计算机程序产品包括计算机可读介质,在该介质上存储有根据第四个方面的计算机程序。
当阅读详细的说明书时,将会理解本发明的其它优点。
附图说明
具体实施例与它们的进一步的目的和优点相结合,通过参考下面的描述以及所附的附图,可以最好地理解本发明,在这些附图中:
图1A和图1B是示出正弦图和正弦图形光谱支架的示意图;分别对应于矩形采样和六边形采样的傅里叶谱的最接近混叠层。
图2是示出通过每秒像素除以半帧时间的时间偏移而获得的六边形采样网格的实例的示意图。
图3A和图3B分别示出了从具有矩形采样(图3A)和六边形采样(图3B)的高斯点采样的正弦图的2D傅立叶变换的示例的示意图。
图4是示出正弦图空间中的像素和帧之间的关系的示例的示意图。
图5是示出了用于2DCT采集的CT系统的示例的示意图,其中X射线投影测量可以例如沿着位于穿过成像物体的单个平面中的投影射线进行。
图6是示出用于3D CT采集的CT系统的实例的示意图,其中X射线投影测量可以例如沿穿过成像物体的三维子体积的投影射线进行。
图7是示出用于计算机断层摄影的数据采集和/或数据读出和/或数据处理的系统的示意图。
图8是示出根据一个实施例所述的用于数据获取的基于处理器的系统的一个示例的示意图。
图9是示出了根据一个实施例所述的计算机装置的另一个实施例的示意图。
图10是示出了整体X射线成像系统的一个实施例的示意图。
图11是示出了根据一个实施例所述的CT几何形状的一个示例的示意图。
图12是示出了每个投影线如何由角坐标θ和径向坐标r限定的一个示例的示意图。
图13是示出角度采样的一个示例的示意图,其中,探测器和辐射源在时间上是由采样间隔T隔开的两个不同的位置示出。
具体实施方式
正如所提到的,在CT场中的快速发展使得台架旋转速度变得越来越高,探测器的空间分辨率更高;这样就增加了对足够角度采样的要求。CT中的角度采样不足导致混淆,这通过从锐边引入条纹状伪像来破坏重建图像。
这里所提出的技术提出了一种独特且空前的数据获取机制,其以有效的方式来管理混叠。
根据所提出的技术的第一个方面,提供了一种用于计算机断层扫描的数据采集方法,基于具有多个像素的探测器的X射线成像系统,其中使用由探测器的相邻像素获取的测量值之间的时间偏移来执行数据采样。
以这种方式,诸如六边形采样图案的定制采样图案,可减少混叠或至少控制所产生的混叠的类型,而不需要依赖于标准过采样技术。
本发明提供了一种获得高效采样的鲁棒方式,能够降低CT图像中混叠伪像的风险。
作为示例,针对相邻像素在不同时间点处启动新帧的数据采样。
例如,数据采样包括通过由相邻像素获取的时间移位帧来连续旋转CT中的正弦图采样。
作为示例,可以执行数据采样以获得定制采样模式。
在一个特定实施例中,可以执行数据采样以获得定制采样模式,以减少混叠或至少控制所产生的混叠的类型。
例如,可以执行数据采样以获得六边形采样图案。
可选地,通过对像素施加随机的时间偏移来执行数据采样。
在一个特定实施例中,利用相邻像素获取的测量之间的半帧时间的时间偏移来执行数据采样。
例如,可以通过将每秒像素除以半帧时间的时间偏移来执行数据采样。
可选地,可以通过一部分帧时间的时间偏移每第i个像素来执行数据采样,其中i为等于或大于2的整数。
在一个特定实施例中,像素i可以被偏移(i-1)/N帧时间,其中i为2至N的整数。
通过以连续旋转方式在时间上移动由CT探测器上的相邻像素获取的帧,实现二维拉冬变换(Radon transform)的定制采样模式。在本发明的一个示例性实施方案中,可通过半帧时间将每其他像素的帧移位来实现六边形采样图案。由于CT数据的性质,六边形采样可以在理论上放松角度采样准则2倍。所获得的效果与两次角过采样的效果相同,而不增加获取的测量的数量。本发明所提出的方法简单易行,对临床高分辨率CT成像具有较大的影响。通过仿真和临床CT数据的后处理来执行评估。
在2D中的六边形采样实现了傅立叶谱的最佳平铺,从而减少了混叠的风险。这对于CT正弦图是特别有效的,因为“蝴蝶结(bowtie)”形状的光谱支持使得其适合于没有光谱重叠的平铺。理论上,正弦图的六边形采样可以将角采样需求减小2倍。
根据所提出的技术的一个方面,通过在由相邻像素获取的测量之间引入时间偏移量,可以在连续旋转CT中实现正弦图变换的六边形采样或类似的定制采样。所提出的方法将被称为时间偏移像素(TOP)法。
由于这种方法不需要对成像系统进行任何重大修改,该方法对高分辨率CT成像具有较高的实际影响,并成为当今的飞行焦斑和四分之一探测器偏移一样有用的方法。
为了更好地理解,现在将参考具体的非限制性实施例来描述本发明所提出的技术。
由于CT扫描仪以越来越高的速度工作,以满足心脏和灌注成像的要求,探测器的空间分辨率变得更高,角度采样率变得越来越关键,以便避免混叠。角混淆从对CT图像进行破坏性的锐边产生条纹。增加角度采样率需要降低探测器的帧时间,并且可以增加产生的数据量。现在将几个过采样技术应用于临床CT扫描仪中,以便减轻在径向和角度方向上的混叠。发明人已经认识到,2D拉冬变换(正弦图)为避免混叠起到很大的作用。
当2D采样时,通常在六角格栅上采样是有益的,因为这样的采样密度高于笛卡尔格栅上的采样密度。在傅里叶域中,六边形网格上的采样意味着光谱被折叠的点也位于六边形网格上。通常,这意味着傅立叶光谱可以被堆叠得更接近,而不需要与原始光谱干涉的折叠频率。对于CT,这尤其是由于正弦图的光谱内容的性质;用于正弦图的基本光谱支持具有“蝴蝶结(bowtie)”形状,特别适用于图1中所见的堆叠;本发明的六边形采样可以在理论上将角度采样需求减小2倍。2D拉冬变换的六边形采样的益处是由Paul Rattey和AllenLindgren于1981年已经描述了(采样2D拉冬变换)。
可以通过应用几种可用的过采样技术中的一种来实现正弦图中的六边形采样网格,例如在所述x,y平面和四分之一探测器中的飞行焦点(FFS),所述x,y平面和四分之一探测器在完整的旋转上具有奇数个帧偏移(否则,QDO具有与增加一半像素宽度相移的另一行像素相同的效果)。过采样技术通常出现一些并发症,例如必须跟踪焦斑运动,或者精确地控制帧的数量和台架的旋转速度之间的关系(如果你希望角样本在180度之后匹配)。
图1A和图1B是示出对应于矩形采样(图1A)和六边形采样(图1B)的正弦图的蝴蝶结(bowtie)形状的光谱支持和傅立叶光谱的最接近的无混淆平铺的示意图。
这里,除了别的之外,我们将显示在连续旋转CT中也可以实现六边形采样。本发明公开了在所述跨轴平面中的相邻像素获取的测量之间引入半帧时间的时间偏移。本发明人在现有技术中从未看到过本发明所提出的方法,现有技术中也没有进行过研究。该方法是鲁棒性的,因此对高分辨率CT成像具有较大的影响。
通常,所提出的技术提供了一种用于获取计算机断层图像的方法和相应的设备和计算机程序,具有多个像素的X射线成像系统,每个像素均获取多个帧,其中所述帧的获取在所述像素的不同子集之间具有至少一个时间偏移。换句话说,不同的像素在不同的时间点启动新帧的采样。
作为示例,具有特定的时间偏移量,可以实现六边形采样。利用其它时间偏移量,可实现其它类型的定制采样,以减少混叠或至少控制所产生的混叠的类型。
通常,像素是探测器的X射线敏感子元件。每个像素以帧的序列测量入射的X射线通量。帧是在特定时间间隔期间的总测量信号。
探测器阵列的尺寸和分段会影响CT设备的成像能力。多个探测器元件(也称为像素)在台架的旋转轴方向上,即图9的z方向使得能多片图像采集。在角度方向上(在图9中ξ)的多个探测器元件(也称为像素)能够同时在同一平面中同时测量多个投影,这被应用于扇形/锥束CT中。大多数传统的探测器是二维的(有时称为平板探测器),这意味着它们在切片方向(z)和角方向(ξ)上具有探测器元件。
对于给定的旋转位置,每个探测器元件(也称为像素)测量对于特定投影线所传输的X射线。这种测量被称为投影测量。即使探测器是二维的,许多投影线的投影测量的集合也称为正弦图,将正弦图制成三维图像。通过图像重建利用正弦图数据获得成像对象内部的图像。每个投影线(正弦图中的点)由角坐标θ和径向坐标r给出,如图10中所定义的,在由(r,θ)给出的特定坐标处的探测器元件进行测量)是正弦图的采样。正弦图中的更多采样通常导致真实正弦图的更好表示,并且因此也更精确地重建图像。
通常,台架连续旋转,并且每个探测器元件测量帧时间内的X射线通量。测量周期T在此被定义为时间间隔,在该时间间隔中,利用测量来占用特定的探测器元件。测量周期的长度可以是但不一定等于帧时间。测量周期远小于总数据采集时间,多个测量周期在整个数据采集/测量中彼此直接跟随。测量周期的长度称为时间采样间隔,采样间隔的倒数1/T称为采样频率。CT系统的角度采样间隔由台架的角速度ω=dθ/dt和时间采样间隔T给出,通过Δθ=ωT。
在图11中显示了角度采样的示意性的示例说明,其中,探测器和辐射源被示出用于由采样间隔所隔开的两个不同的位置,对应于特定探测器元件的所有投影线的径向坐标是相对于所述台架的旋转不变的。
为了对断层数据执行精确的图像重建,必须有足够量的角样本。不足的角度采样可导致图像中的伪像、混叠和较差的分辨率。
增加角采样频率的一种方法是减小时间采样间隔T。降低时间采样间隔导致产生的数据量的相应增加。
图4是示出正弦图空间中的像素和帧之间的关系的示例的示意图。
在2DCT采集中,X射线投影测量可以通过例如沿单平面内的投影射线穿过成像对象来进行,如图5所示。
在3D CT采集中,X射线投影测量可以通过例如沿投影线穿过成像物体的三维子体积来进行,如图6所示。
2D典型地意味着通过待成像的物体在一个平面内和同一平面内进行测量。在3D中,沿着与源探测器对旋转的平面正交的方向在不同的位置上执行测量。
在扇形束CT中,探测器阵列中的像素对斜线(图2中的断续线)上的正弦图进行采样,对于每个视图(角位置)。如果探测器以角速度ω旋转,帧时间为Δt,采样位置在帧之间的角度方向上移动距离ωΔt。通过对于每秒像素引入Δt/2的时间偏移量,实现正弦图的六边形采样。
图2是示出通过每秒像素除以半帧时间的时间偏移而实现的六边形采样网格的实例的示意图。
作为二维拉冬变换中六边形采样功率的演示,我们已经模拟了来自单一高斯点的高res正弦图。然后将高res正弦图向下采样2倍,到以下之上:(1)矩形网格(沿角度方向去除每秒样本)以及(2)六边形网格(在角度方向上去除不断的第二样本,但剩余的角样本移到相邻的径向位置)。然后用2D线性插值将正弦图插回到旧的高res网格上。图2示出了对应于(a)矩形采样栅和(b)六边形采样栅的上采样正弦图的傅立叶变换。由矩形采样栅(在角落中的斜线)产生的角混叠)在使用六边形采样栅时不可见。
图3A和图3B分别示出了从具有矩形采样(图3A)和六边形采样(图3B)的高斯点采样的正弦图的2D傅立叶变换的示例的示意图。图3B中的六边形采样栅去除了图3A中可见的大量角混叠。
举例来说,时间偏移像素(TOP)法的效果可以通过对实验采集的CT数据进行后处理来被证明,以便在类似于临床CT扫描仪的情况下评估混叠抑制中的改善的幅度。
同时,TOP可以与四分之一探测器偏移(QDO)和/或飞行焦点(FFS)相结合,在探测器的z方向上的时间偏移的效果可以通过仿真或实际实验而被证明。
该方法提供了鲁棒的方式来获得更有效的采样,并且可以减少CT图像中的混淆伪像的风险。理论上,该方法放宽了角度采样准则(为了避免混叠所需的视数)2倍。本发明将对临床CT数据进行评估。
通过对CT探测器上的相邻像素获取的帧进行时间移位,实现二维拉冬变换的六边形采样。由于CT数据的性质,六边形采样可以在理论上放松角度采样准则2倍。所获得的效果与两次角过采样的效果相同,然而,不需要额外的测量,这使得所提出的方法更加有效。此外,所提出的技术对高分辨率CT成像具有较大的影响。
代替每秒像素除以半帧时间的时间偏移,更一般地可能通过每第n个像素除以一部分帧时间的时间偏移。例如,相对于未偏移的像素1,像素2偏移了1/N帧时间,像素3偏移了2/N帧时间,依此类推,直到像素N偏移(N-1)/N帧时间。甚至可以对像素施加随机的时间偏移。
可替代地,也可以在所谓的z方向上施加时间偏移,即在与辐射源-探测器对旋转的平面的正交方向上。
应当理解的是,本文所述的方法和装置可以以多种方式组合和重新排列。
例如,具体实施方式可以以硬件来实现,或以通过适合的处理电路来执行的软件来实现,或硬件与软件的组合来实现。
这里描述的步骤、功能、过程、模块和/或块可以使用任何常规技术在硬件中实现,例如分立电路或集成电路技术,包括通用电子电路和专用电路。
特定示例包括一个或多个适当配置的数字信号处理器和其它已知的电子电路,例如被互连以执行专用功能的离散逻辑门或专用集成电路(ASIC)。
可选地,这里描述的至少一些步骤、功能、程序、模块和/或块可以在诸如计算机程序的软件中实现,用于由诸如一个或多个处理器或处理单元之类的适当处理电路执行。
处理电路的示例包括但不限于:一个或多个微处理器,一个或多个数字信号处理器(DSP),一个或多个中央处理单元(CPU),视频加速硬件和/或诸如一个或多个现场可编程门阵列(FPGA)的任何适当的可编程逻辑电路),或一个或多个可编程逻辑控制器(PLC)。
还应当理解的是,可以使用任何常规设备或单元的一般处理能力来重新使用所提出的技术。还可以例如通过对现有软件的重新编程或通过添加新的软件组件来重新使用现有软件。
在一个方面中,提供了一种用于在计算机断层摄影中进行采样的设备。
图7是示出用于计算机断层摄影的数据采集和/或数据读出和/或数据处理的系统的示意图。
根据一个方面,提供了一种用于计算机断层扫描(CT)的数据采集的系统20,基于具有多个像素的探测器12的X射线成像系统10,所述系统20是被配置为对于由所述探测器的相邻像素获取的测量之间的时间偏移进行数据采样。
根据另一个方面,还提供了一种计算机断层摄影(CT)系统10,其包括如这里所述的用于数据获取的系统20。
作为示例,系统20可被配置为启动采样用于相邻像素的不同时间点处的新帧
例如,系统20可被配置为通过由相邻像素获取的时间移位帧来执行连续旋转CT中的正弦图采样。
作为示例,系统20可被配置为执行数据采样以获得定制采样模式。
在一个特定实施例中,系统20被配置为执行数据采样以获得定制采样模式,以减少混叠或至少控制所产生的混叠的类型。
例如,系统20可被配置为执行数据采样以获得六边形采样图案。
或者,系统20可被配置为通过对像素施加随机化时间偏移来执行数据采样。
在一个特定实施例中,系统20被配置为执行数据采样,以在相邻像素获取的测量之间的时间偏移量为半帧时间。
例如,系统20可被配置为通过将每秒像素除以半帧时间的时间偏移来执行数据采样。
可选地,系统20可被配置为通过一部分帧时间的时间偏移每第i个像素来执行数据采样,其中i为等于或大于2的整数。
在一个特定实施例中,例如,系统20可被配置为通过(i-1)/N帧时间来偏移像素i,其中i为2至N的整数。
作为示例,每个像素13可以是探测器12的X射线敏感子元件,并且系统20可以被配置为使得每个像素13以帧的序列来测量入射的X射线通量,其中,帧是在特定时间间隔期间的总测量信号。
在一个特定实施例中,如图8所示,装置100包括处理器110和存储器120,该存储器包括处理器可执行的指令,因而处理器可操作以执行和/或控制数据获取。可选地,该装置包括用于接收输入数据并输出所产生的输出数据的输入/输出接口。
在这个特定实施例中,这里所描述的至少一些步骤、功能、过程、模块和/或块是在计算机程序中实现的,其被加载到存储器中,以供包括一个或多个处理器的处理电路执行。处理器和存储器相互连接,以实现正常的软件执行。可选的输入/输出设备也可以与处理器和/或存储器互连,以实现诸如输入参数和/或所得到的输出参数的相关数据的输入和/或输出。
术语“处理器”应在一般意义上被解释为能够在任何系统或设备执行程序代码或计算机程序指令来执行特定的处理,确定或计算任务。
因此,当执行计算机程序时,包括一个或多个处理器的处理电路被配置为执行诸如本文所述的明确定义的处理任务。
所述处理电路不必专用于仅执行上述步骤、功能、过程和/或块,也可以执行其他任务。
图9是示出了根据一个实施例的计算机装置的另一示例的示意图。
另一个方面,提供了一种用于在计算机断层摄影中改进采样的计算机程序。
更具体地,提供了一种计算机程序225、235,当由处理器210执行时,采用具有多个像素的探测器的X射线成像系统,该计算机程序用于计算机断层摄影(CT)能使数据获取,计算机程序225、235包括指令,当由处理器210执行时,导致处理器210对于由所述探测器的相邻像素获取的测量值之间的时间偏移进行数据采样。
本发明所提出的技术还提供一种包括计算机可读介质220、230的计算机程序产品,该介质上存储有这样的计算机程序。
作为示例,软件或计算机程序225、235可以被实现为计算机程序产品,其通常承载或存储在计算机可读介质220、230上,尤其是非易失性介质上。计算机可读介质可以包括一个或多个可移动或不可移动存储设备,但不限于:只读存储器(ROM)、随机存取存储器(RAM)、光盘、数字多功能光盘(DVD)、蓝光光盘、通用串行总线(USB)存储器、硬盘驱动器(HDD)存储装置、闪存、磁带或任何其他常规存储装置。因此,计算机程序可以被加载到计算机或等效处理设备的操作存储器中以供其处理电路执行。
在此呈现的方法流程可以被认为是当由一个或多个处理器执行时的计算机动作流。相应的设备、系统和/或设备可以被定义为一组功能模块,处理器对应于功能模块。在这种情况下,功能模块可被实现为在处理器上运行的计算机程序。因此,装置、系统和/或装置可替换地被定义为一组功能模块,其中所述功能模块被实现为在至少一个处理器上运行的计算机程序。
驻留在存储器中的计算机程序因此可以被组织为适当的功能模块,当由所述处理器执行时,该功能模块被配置为执行在此描述的至少部分所述步骤和/或任务。
可替代地,可以主要由硬件模块或可替换地通过硬件来实现模块。软件对硬件的程度是纯粹的实施选择。
应当理解的是,在此描述的布置可以以不同方式实现、组合和重新排列。
以上所述的实施例仅作为示例给出,应当理解的是,本发明所提出的技术不限于此。本领域技术人员将会理解,在不脱离由所附权利要求书限定的本发明范围的情况下,各种修改都可以对这些实施例进行组合和改变。具体地,在技术上可能的情况下,在不同实施例中的不同部分解决方案可以在其他配置中组合。
Claims (28)
1.一种用于计算机断层摄影(CT)的数据采集方法,基于具有多个像素(13)的探测器(12)的X射线成像系统(10),其特征在于:数据采样是利用由探测器的相邻像素获取的测量值之间的时间偏移来执行。
2.根据权利要求1所述的方法,其特征在于:新帧的数据采样在对于相邻像素的不同时间点处开始的。
3.根据权利要求1或2所述的方法,其特征在于:所述数据采样包括通过由相邻像素获取的时间移位帧以连续旋转CT方式的正弦图采样。
4.根据权利要求1至3之一所述的方法,其特征在于:所述数据采样是被执行以获得定制采样模式。
5.根据权利要求4所述的方法,其特征在于:所述数据采样是被执行以获得定制采样模式,以减少混叠或至少控制所产生的混叠的类型。
6.根据权利要求1至5之一所述的方法,其特征在于:所述数据采样是被执行以获得六边形采样图案。
7.根据权利要求1至6之一所述的方法,其特征在于:所述数据采样是通过对像素施加随机的时间偏移来执行。
8.根据权利要求1至7之一所述的方法,其特征在于:所述数据采样是利用相邻像素所获取的测量之间的半帧时间的时间偏移来执行。
9.根据权利要求8所述的方法,其特征在于:所述数据采样是通过每秒像素除以半帧时间的时间偏移来执行。
10.根据权利要求1至9之一所述的方法,其特征在于:所述数据采样是通过每第i个像素除以一部分帧时间的时间偏移来执行的,其中i是等于或大于2的整数。
11.根据权利要求10所述的方法,其特征在于:像素i是帧时间的(i-1)/N的偏移,其中i是从2至N的整数。
12.根据权利要求1至11之一所述的方法,其特征在于:每个像素(13)是所述探测器(12)的X射线敏感子元件,并且每个像素(13)将入射的X射线通量作为帧的序列进行测量,其中,帧是在特定时间间隔期间的总测量信号。
13.一种用于计算机断层扫描(CT)的数据采集的系统(20),基于具有多个像素(13)的探测器的X射线成像系统,其特征在于:所述系统(20)是被配置为以在由所述探测器的相邻像素所获取的测量值之间的时间偏移来执行数据采样。
14.根据权利要求13所述的系统,其特征在于:所述系统(20)是被配置为在对于相邻像素在不同时间点处开始对新帧的采样。
15.根据权利要求13或14所述的系统,其特征在于:所述系统(20)是被配置为通过由相邻像素所获取的时间移位帧来执行连续旋转CT中的正弦图采样。
16.根据权利要求13至15之一所述的系统,其特征在于:所述系统(20)是被配置为执行数据采样以获得定制采样模式。
17.根据权利要求16所述的系统,其特征在于:所述系统(20)是被配置为执行数据采样以获得定制采样模式,以减少混叠或至少控制所产生的混叠的类型。
18.根据权利要求13至17之一所述的系统,其特征在于:所述系统(20)是被配置为执行数据采样以获得六边形采样图案。
19.根据权利要求13至18之一所述的系统,其特征在于:所述系统(20)是被配置为通过对所述像素施加随机时间偏移来执行所述数据采样。
20.根据权利要求13至19之一所述的系统,其特征在于:所述系统(20)是被配置为以在由相邻像素获取的测量之间的半帧时间的时间偏移来执行数据采样。
21.根据权利要求20所述的系统,其特征在于:所述系统(20)是被配置为以通过每秒像素除以半帧时间的时间偏移来执行所述数据采样。
22.根据权利要求13至21之一所述的系统,其特征在于:所述系统(20)是被配置为通过每第i个像素除以一部分帧时间的时间偏移来执行数据采样,其中i是等于或大于2的整数。
23.根据权利要求22所述的系统,其特征在于:该系统被配置为将像素i除以(i-1)/n的帧时间,其中i是从2到N。
24.根据权利要求13至23之一所述的系统,其特征在于:每个像素(13)是所述探测器(12)的X射线敏感子元件,所述系统(20)是被配置为使得每个像素以帧的序列测量入射的X射线通量,其中,帧是在特定时间间隔期间的总测量信号。
25.根据权利要求13至24之一所述的系统,其特征在于:所述系统(20、100、200)包括处理器(110、210)和存储器(120,220),所述存储器包括由所述处理器执行的指令,由此处理器(110,210)可操作以执行和/或控制数据获取。
26.一种计算机断层摄影(CT)系统(10),包括根据权利要求13至25之一所述的用于数据采集的系统(20)。
27.一种计算机程序(225、235),当由处理器(210)执行时,采用具有多个像素的探测器的X射线成像系统,所述计算机程序能用于计算机断层摄影(CT)的数据采集,其中所述计算机程序(225;235)包括指令,当由处理器(210)执行时,导致处理器(210)对于由所述探测器的相邻像素获取的测量值之间的时间偏移进行数据采样。
28.一种包括计算机可读介质(220、230)的计算机程序产品,所述介质上存储有如权利要求27所述的计算机程序(225、235)。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201662373408P | 2016-08-11 | 2016-08-11 | |
US62/373,408 | 2016-08-11 | ||
PCT/SE2017/050496 WO2018030928A1 (en) | 2016-08-11 | 2017-05-15 | Data acquisition for computed tomography |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110113994A true CN110113994A (zh) | 2019-08-09 |
Family
ID=61160642
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201780048118.3A Pending CN110113994A (zh) | 2016-08-11 | 2017-05-15 | 用于计算机断层摄影的数据采集 |
Country Status (6)
Country | Link |
---|---|
US (1) | US10130313B2 (zh) |
EP (1) | EP3496616A4 (zh) |
JP (1) | JP6946417B2 (zh) |
KR (1) | KR20190033621A (zh) |
CN (1) | CN110113994A (zh) |
WO (1) | WO2018030928A1 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2579798A (en) * | 2018-12-13 | 2020-07-08 | Ucl Business Ltd | Computed tomography imaging |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5416817A (en) * | 1992-03-23 | 1995-05-16 | Siemens Aktiengesellschaft | Method of operating a fan beam computer tomography apparatus |
WO2006064403A2 (en) * | 2004-12-17 | 2006-06-22 | Koninklijke Philips Electronics, N.V. | Pulsed x-ray for continuous detector correction |
CN1969759A (zh) * | 2005-11-11 | 2007-05-30 | Ge医疗系统环球技术有限公司 | X射线ct成像方法和x射线ct设备 |
JP2007143605A (ja) * | 2005-11-24 | 2007-06-14 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
CN102778688A (zh) * | 2011-05-09 | 2012-11-14 | 西门子公司 | 用于探测x射线量子的方法和装置 |
US20140224994A1 (en) * | 2011-08-30 | 2014-08-14 | Ucl Business Plc | Radiation Detector |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5173852A (en) * | 1990-06-20 | 1992-12-22 | General Electric Company | Computed tomography system with translatable focal spot |
JP4409223B2 (ja) * | 2003-07-24 | 2010-02-03 | 東芝医用システムエンジニアリング株式会社 | X線ct装置及びx線ct用逆投影演算方法 |
US20080170658A1 (en) * | 2007-01-11 | 2008-07-17 | Jed Douglas Pack | System and method for ct imaging with increased sampling and reduced artifacts |
US9818182B2 (en) | 2012-06-20 | 2017-11-14 | Hitachi, Ltd. | X-ray CT device |
US9841514B2 (en) * | 2015-09-24 | 2017-12-12 | Prismatic Sensors Ab | X-ray detector arrangement |
JP6273241B2 (ja) | 2015-09-24 | 2018-01-31 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 放射線断層撮影方法及び装置並びにプログラム |
EP3374803B1 (en) * | 2015-11-12 | 2020-04-01 | Prismatic Sensors AB | High-resolution computed tomography using edge-on detectors with temporally offset depth-segments |
-
2017
- 2017-05-15 WO PCT/SE2017/050496 patent/WO2018030928A1/en unknown
- 2017-05-15 JP JP2019507185A patent/JP6946417B2/ja active Active
- 2017-05-15 US US15/594,800 patent/US10130313B2/en active Active
- 2017-05-15 EP EP17839916.8A patent/EP3496616A4/en not_active Withdrawn
- 2017-05-15 KR KR1020197006913A patent/KR20190033621A/ko not_active Application Discontinuation
- 2017-05-15 CN CN201780048118.3A patent/CN110113994A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5416817A (en) * | 1992-03-23 | 1995-05-16 | Siemens Aktiengesellschaft | Method of operating a fan beam computer tomography apparatus |
WO2006064403A2 (en) * | 2004-12-17 | 2006-06-22 | Koninklijke Philips Electronics, N.V. | Pulsed x-ray for continuous detector correction |
CN1969759A (zh) * | 2005-11-11 | 2007-05-30 | Ge医疗系统环球技术有限公司 | X射线ct成像方法和x射线ct设备 |
JP2007143605A (ja) * | 2005-11-24 | 2007-06-14 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
CN102778688A (zh) * | 2011-05-09 | 2012-11-14 | 西门子公司 | 用于探测x射线量子的方法和装置 |
US20140224994A1 (en) * | 2011-08-30 | 2014-08-14 | Ucl Business Plc | Radiation Detector |
Also Published As
Publication number | Publication date |
---|---|
JP2019531111A (ja) | 2019-10-31 |
WO2018030928A1 (en) | 2018-02-15 |
EP3496616A4 (en) | 2020-02-26 |
EP3496616A1 (en) | 2019-06-19 |
JP6946417B2 (ja) | 2021-10-06 |
US10130313B2 (en) | 2018-11-20 |
KR20190033621A (ko) | 2019-03-29 |
US20180042561A1 (en) | 2018-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6150940B2 (ja) | 位相コントラストctを使うことによる単色減衰コントラスト画像生成 | |
KR101632046B1 (ko) | 고해상도 컴퓨터 단층 촬영 | |
CN105682559B (zh) | 用于估计散射的方法和系统 | |
GB2462529A (en) | High-energy dual-energy CT system showing base material coefficient distribution | |
WO2003073938A1 (en) | Non-invasive stationary system | |
CN103472074A (zh) | Ct成像系统和方法 | |
JP2014232119A (ja) | 3次元陽子コンピュータ断層撮影のための高性能計算(hpc−pct) | |
JP6633201B2 (ja) | 時間オフセット深度セグメントを有するエッジオン検出器を用いた高分解能コンピュータ断層撮影 | |
US8437522B2 (en) | Motion index for medical imaging data based upon Grangeat's formula | |
US8897531B2 (en) | Method and computed tomography system for generating tomographic image datasets with integrating detector and counting detector | |
WO2017052431A1 (en) | Determining the orientation of an edge-on x-ray detector with respect to the direction of incoming x-rays | |
JP6917356B2 (ja) | 主成分材料組み合わせによる分解方法および装置 | |
JP2016087466A (ja) | X線ct装置及びその制御方法 | |
CN112763519B (zh) | 一种用拟蒙特卡罗方法模拟光子散射的方法 | |
JP2010204755A (ja) | 画像処理装置、画像再構成システム、画像処理方法およびプログラム | |
WO2014107651A1 (en) | System and method for ultra-high resolution tomorgraphic imaging | |
US20210049794A1 (en) | Method for image reconstruction of an object, in particular based on computed-tomography image reconstruction, and apparatus, system and computer program product for the same | |
CN110113994A (zh) | 用于计算机断层摄影的数据采集 | |
JP7015781B2 (ja) | 多重エネルギー(二重エネルギーを含む)コンピュータ断層撮影(ct)イメージングを行うための方法および装置 | |
US9869646B2 (en) | View triggering for radiation imaging system | |
WO2016077865A1 (en) | Single/limited angle gamma/x-ray tomography | |
JP7001176B2 (ja) | データ処理方法、プログラム、データ処理装置および陽電子放出断層撮像装置 | |
WO2011008787A1 (en) | System and method for image reconstruction by using multi-sheet surface rebinning | |
JPS6314983B2 (zh) | ||
JP6710759B2 (ja) | X線検出器用の測定回路、ならびに対応する方法およびx線撮像システム |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190809 |
|
WD01 | Invention patent application deemed withdrawn after publication |