CN1931098A - X-射线ct设备 - Google Patents
X-射线ct设备 Download PDFInfo
- Publication number
- CN1931098A CN1931098A CNA2006101263673A CN200610126367A CN1931098A CN 1931098 A CN1931098 A CN 1931098A CN A2006101263673 A CNA2006101263673 A CN A2006101263673A CN 200610126367 A CN200610126367 A CN 200610126367A CN 1931098 A CN1931098 A CN 1931098A
- Authority
- CN
- China
- Prior art keywords
- ray
- data
- data acquisition
- detector
- channel direction
- 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
- 239000000523 sample Substances 0.000 claims description 46
- 230000005855 radiation Effects 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000003384 imaging method Methods 0.000 abstract description 58
- 230000002093 peripheral effect Effects 0.000 abstract description 3
- 238000002591 computed tomography Methods 0.000 description 58
- 238000000034 method Methods 0.000 description 28
- 230000000007 visual effect Effects 0.000 description 28
- 230000008569 process Effects 0.000 description 13
- 230000000747 cardiac effect Effects 0.000 description 12
- 210000004072 lung Anatomy 0.000 description 11
- 230000011664 signaling Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 10
- 238000001914 filtration Methods 0.000 description 10
- 238000005286 illumination Methods 0.000 description 10
- 238000006243 chemical reaction Methods 0.000 description 9
- 238000012937 correction Methods 0.000 description 9
- 238000013480 data collection Methods 0.000 description 9
- 238000005111 flow chemistry technique Methods 0.000 description 9
- 238000010626 work up procedure Methods 0.000 description 7
- 241000446313 Lamella Species 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 6
- 238000013170 computed tomography imaging Methods 0.000 description 6
- 230000008034 disappearance Effects 0.000 description 6
- 230000035515 penetration Effects 0.000 description 6
- NCGICGYLBXGBGN-UHFFFAOYSA-N 3-morpholin-4-yl-1-oxa-3-azonia-2-azanidacyclopent-3-en-5-imine;hydrochloride Chemical compound Cl.[N-]1OC(=N)C=[N+]1N1CCOCC1 NCGICGYLBXGBGN-UHFFFAOYSA-N 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 239000012141 concentrate Substances 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 238000012805 post-processing Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 210000003050 axon Anatomy 0.000 description 2
- 210000004185 liver Anatomy 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000002685 pulmonary effect Effects 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 230000008719 thickening Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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/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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—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
-
- 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
数据采集装置(20),它具有包括在其中央部分通道细而在其外围部分通道粗或糙的多个通道宽度的X-射线探测器(24),和多个数据采集范围,它包括在通道方向上宽的数据采集范围和在通道方向上窄的数据采集范围,且能够对每个数据采集在数据采集范围内执行转换,使用X-射线数据采集装置(20)以便在通道方向上窄的数据采集范围的中央部分细的通道上执行数据采集,以此提供X-射线CT设备(100),其能够执行高分辨率成像并产生更满意的图像质量。
Description
技术领域
本发明涉及一种用于医疗用途的X-射线CT(计算机断层摄影)设备或一种用于工业用途的X-射线CT设备,具体涉及一种在诸如螺旋扫描和电影扫描或类似的所有扫描模式中在低放射曝照下能够成像或拍摄高分辨率X-射线断层摄影图像的X-射线CT设备。
背景技术
用于X-射线CT设备中的X-射线探测器,像多行X-射线探测器或以平板为代表的矩阵结构的二维X-射线区域探测器,迄今为止已经被制造成如图12(例如,参考专利文献1)中所示的固定间隔(间距)并带有固定通道宽度。
由于在通道方向上的分辨率,就是说,在传统X-射线探测器中每一通道上的每个探测器宽度/间隔即使是在采取小成像区域时仍保持不变,因而没有按比例提高实际的分辨率。因此,即使在小的成像区域中将一个像素的尺寸做的很小,也没有获得具有比根据探测器间隔所确定的基于尼奎斯特频率的分辨率更高的断层摄影图像。
【专利文献1】日本未审查的专利公开号为2000-193750
因此,即使是将成像区域做的很小试图在高分辨率下观看它,断层摄影图像也仅是模糊的,因而没有获得高分辨率的断层摄影图像。虽然通过略微加强重建函数的高频或RF区域能够实现对比度的轻微改进和分辨率的改进,但是作为相反的效果已经产生了噪音的增加和伪影的增加。
发明内容
因此,本发明的一个目的是提供一种能够成像或拍摄高分辨率X-射线断层摄影图像的X-射线CT设备。
本发明的另一个目的是提供一种当采取小成像区域时能够获得高分辨率的X-射线CT设备。
本发明进一步的目的是提供一种减少所用的X-射线量并有效地使用X-射线从而能够减少受试者对X-射线曝照的X-射线CT设备。
在第一个方面中,本发明提供一种X-射线CT设备,其包括X-射线数据采集装置,该采集装置含有X-射线发生器,和与X-射线发生器呈相对关系设置并探测X-射线的多行X-射线探测器或能围绕安置在X-射线发生器和探测器之间的旋转中心旋转的矩阵结构的二维X-射线区域探测器,从而获得穿过位于X-射线发生器和探测器之间受试者的X-射线的投影数据;图像重建装置,该图像重建装置对从X-射线数据采集装置所获得的投影数据进行图像重建;和显示装置,该显示装置显示进行了图像重建的图像,其中X-射线数据采集装置以这样的方式设定使得以通道方向看去在探测器每个外围部分的探测器通道宽度d2相对于以通道方向看去在探测器中央部分的探测器通道宽度d1变成d1<d2,或者从以通道方向看去探测器的中央部分到其外围部分所设置的多个探测器宽度(d1,d2,...di,...dn-1,dn)满足d1≤d2≤...≤di≤...≤dn-1≤dn。
在根据第一方面的X-射线CT设备中,由于通道宽度更窄的X-射线探测器通道集中在中央部分,通过使用在中央部分通道宽度窄的X-射线探测器通道执行数据采集和图像重建可进行空间高分辨率X-射线CT成像。
在第二方面,本发明提供一种X-射线CT设备,其包括X-射线数据采集装置,该采集装置含有X-射线发生器,和与X-射线发生器呈相对关系设置并探测X-射线的多行X-射线探测器或能围绕安置在X-射线发生器和探测器之间的旋转中心旋转的矩阵结构的二维X-射线区域探测器,从而获得穿过位于X-射线发生器和探测器之间受试者的X-射线的投影数据;图像重建装置,该图像重建装置对从X-射线数据采集装置所获得的投影数据进行图像重建;和显示装置,该显示装置显示进行了图像重建的图像,其中X-射线数据采集装置具有从在探测器通道方向上宽的数据采集范围l1到在通道方向上窄的数据采集范围ln的多个数据采集范围l1≥l2≥...≥li≥...ln-1≥ln,并且以这样的方式设定X-射线数据采集装置使得数据采集范围对于每个数据采集是可转换的。
在根据第二方面的X-射线CT设备中,由于在中央部分存在更窄的数据采集范围,在带有细的通道宽度和间隔的中央部分的更窄数据采集范围上实现数据的采集并执行图像的重建,从而能够进行空间高分辨率X-射线CT成像。
在第三方面,本发明提供了一种X-射线CT设备,其中当在从探测器的通道方向看去窄的数据采集范围中执行数据采集时,X-射线数据采集装置在以通道方向看去探测器中央部分的探测器通道宽度窄的部分上执行数据采集。
在根据第三方面的X-射线CT设备中,由于通道宽度更窄的探测器通道集中在中央部分并且存在更窄的数据采集范围,在带有细通道宽度并以细通道间隔的中央部分的更窄的数据采集范围上实现数据的采集并执行图像的重建,从而能够进行空间高分辨率X-射线CT成像。
在第四个方面,本发明提供了一种X-射线CT设备,其中X-射线数据采集装置具有多个执行数据采集的通道。
在根据第四方面的X-射线CT设备中,由于执行数据采集的探测器通道的数量在多种模式下是可转换的,在一种用于少量通道并以X-射线数据采集装置A/D转换器采样率为最大值的模式下高速获得关于中央部分少量探测器通道的数据,并执行图像重建,从而能够根据时间高分辨率地进行X-射线CT成像。
在第五个方面,本发明提供了一种X-射线CT设备,其中X-射线数据采集装置具有执行数据采集的多个通道以及多幅视图。
在根据第五方面的X-射线CT设备中,由于执行数据采集的探测器通道数量在多种模式下是可转换的以及执行数据采集的视图数量在多种模式下是可转换的,在X-射线数据采集装置的A/D转换器采样率的最大值和在少量视图模式下高速获得关于中央部分少量探测器通道的数据并执行图像重建,从而能够根据时间高分辨率地进行X-射线CT成像。
在第六个方面,本发明提供了一种X-射线CT设备,其中X-射线数据采集装置具有执行数据采集的多行,并且根据每个通道的位置,行的数量不同。
在根据第六方面的X-射线CT设备中,在z方向看去行数大的探测器通道集中于中央部分并且在这种情况下,在带有细通道宽度并在细通道间隔下的中央部分的窄的数据采集范围上实现数据的采集并执行图像的重建,从而能够进行空间高分辨率X-射线CT成像。
在第七个方面,本发明提供了一种X-射线CT设备,其中X-射线数据采集装置具有执行数据采集的多个数据采集采样周期。
在第八个方面,本发明提供了一种X-射线CT设备,其中X-射线数据采集装置具有执行数据采集的多个数据采集采样周期,并且根据通道位置该数据采集采样周期是不同的。
在根据第七和八方面的X-射线CT设备中,在数据采集采样周期短的模式下高速收集或采集关于中央部分少量探测器通道的数据并进行图像重建,从而根据时间高分辨率的X-射线CT成像是能够进行的。
在第九个方面,本发明提供了一种包括控制装置的X-射线CT设备,该控制装置以这样一种方式控制X-射线照射区域使得X-射线仅照射在从通道方向或其内部范围看去探测器中央部分的探测器通道宽度窄的一些范围上,在探测器通道方向或其内部范围上窄的数据采集范围,或在探测器通道方向上窄的一些数据采集范围,即,从通道方向或其内部范围看去中央部分的探测器通道宽度窄的范围。
在根据第九方面的X-射线CT设备中,由于所照射的X-射线可以最优化地更窄并且当在中央部分更窄的数据采集范围中进行数据采集时通过控制装置在通道方向上进行照射,可在低曝照放射下进行受试者的断层摄影成像。
在第十个方面,本发明提供了具有一种装置的X-射线CT设备,该装置以这样一种方式限制X-射线照射区域使得将X-射线照射到探测器通道方向上的一些范围中,该范围在从探测器的通道方向或其内部范围看去在中央部分的通道上是细的,在探测器的通道方向或其内部范围上窄的数据采集范围,或在探测器的通道方向上窄的数据采集范围,即,从通道方向或其内部范围看去中央部分的细通道的范围。
在根据第十方面的X-射线CT设备中,由于所照射的X-射线可以最优化地更窄并且当在中央部分更窄的数据采集范围中执行数据采集时通过限制装置在通道方向上进行照射,可在低曝照放射下执行受试者的断层摄影成像。
本发明可提供一种能够成像或拍摄高分辨率X-射线断层摄影图像的X-射线CT设备。而且,本发明可提供一种当采取小成像区域时能够获得高分辨率的X-射线CT设备。进一步,本发明可提供一种减少所用的X-射线量并有效地使用X-射线从而能够减少受试者对X-射线曝照的X-射线CT设备。
从所附的附图所阐明的发明优选实施例的以下描述中本发明进一步的目的和优势将更加明显。
附图说明
图1是显示了根据本发明一个实施例的X-射线CT设备的结构图。
图2是示出了X-射线发生器(X-射线管)和多行X-射线探测器的示范性视图。
图3是描述了根据本发明一个实施例的X-射线CT设备的示意性操作的流程图。
图4是显示了预处理细节的流程图。
图5是示出了三维图像重建处理细节的流程图。
图6是描述了在X-射线穿透方向上重建区域内射线投影状态的原理图。
图7是显示了X-射线探测器平面上所投影的射线的原理图。
图8是示出了在重建区域上投影数据Dr(view,x,y)的投影状态的原理图。
图9是描述了在重建区域上各个像素的反投影像素数据D2的原理图。
图10是显示了将反投影像素数据D2加到对应的所有视图的像素上以获得反投影数据D3的状态的说明图。
图11是显示了在X-射线穿透方向上投影圆形重建区域上的射线的状态原理图。
图12是示出了传统多行X-射线探测器的图。
图13是显示了中央通道形成高分辨率的多行X-射线探测器的图。
图14是示出了传统数据读取模式的图。
图15是描述了用于读取行数在内部中央部分大的模式1的图。
图16是显示了用于读取行数在内部中央部分大的模式1的图。
图17是示出了用于读取行数在内部中央部分大的模式1的图。
图18是显示了用于读取具有多种类型通道宽度和数据采集范围的多行X-射线探测器的内部中央部分上行宽总和大且细,以及读取其外围部分的行宽总和小且糙的行模式的图。
图19是示出了受试者体积大而其感兴趣区域小的方式的图。
图20是显示了在常规模式下数据传送率的图。
图21是描述了从一通道方向看去仅在中央部分获得或收集数据的数据传送率的图。
图22是显示了从该通道方向看去仅在中央部分收集数据的数据传送率的图。
图23是示出了给每个数据采集范围提供多种类型的数据采集采样周期的多行X-射线探测器的图。
图24是显示了通过通道方向准直器将X-射线照射范围与数据采集范围匹配的图。
图25是示出了由形成X-射线滤波器的线束所限定的数据采集范围的图。
图26(a)是显示了从RL方向(x方向)看去在探查图像上设置感兴趣区域的图,而图26(b)是显示了从AP方向(y方向)看去在探查图像上设置感兴趣区域的图。
图27是示出了在X-射线探测器通道宽度d和X-射线探测器通道宽度d/2之间转换的图。
图28是显示了在宽的数据采集范围和窄的数据采集范围之间转换的图。
图29是实施例2的操作的流程图。
图30(a)是显示了通道方向准直器(旋转轴偏心圆柱状系统)的图,图30(b)是显示了通道方向准直器(屏蔽板系统)的图,和图30(c)是显示了形成X-射线滤波器的线束的例子的图。
图31(a)和31(b)分别是显示了通道方向准直器控制的图。
图32是示出了增加在通道方向X-射线准直器所缺失的投影数据方式的图。
图33是显示了通道方向准直器前馈控制的图。
图34是在视角=0°时感兴趣的成像区域和照射通道范围的说明图。
图35是在视角=0°时感兴趣的成像区域、照射最小通道和照射最大通道的说明图。
图36是在视角=β时感兴趣的成像区域、照射最小通道和照射最大通道的说明图。
图37是显示了通道方向准直器反馈控制的图。
图38(a)是显示了形成X-射线滤波器32的线束的常规位置图,图38(b)是显示了在形成X-射线滤波器32的线束上位置控制(部分1)的图,而图38(c)是显示了在形成X-射线滤波器32的线束上位置控制(部分2)的图。
图39是示出了对于每个X-射线探测器通道间隔图像重建函数不同的图。
图40是显示了存在多种X-射线探测器通道宽度的重建函数卷积的流程图。
具体实施方式
通过在附图中所示的实施例将在下文中更详细地描述本发明。顺便提及的是本发明并不因此而受限。
(实施例1)
图1是根据本发明一个实施例的X-射线CT设备的构型结构图。X-射线CT设备100装备有操作控制台1、拍摄或成像工作台10和扫描台架20。
操作控制台1装备有接受操作者输入的输入装置2、执行图像重建处理或类似处理的中央处理单元3、获得或收集通过扫描台架20所得的投影数据的数据采集缓冲器5、显示从投影数据所重建的CT图像的监视器6和在其中储存程序、数据和X-射线CT图像的存储装置7。
拍摄工作台10装备有能将受试者放置其上并将受试者运入扫描台架20的空腔部分中并能从那里将受试者运出的托架12。通过内置于拍摄工作台10内的电机将托架12向上和向下移动并沿拍摄工作台10线性移动。拍摄工作台10的托架12的移动方向定义为z方向。
扫描台架20装备有X-射线管21、X-射线控制器22、片层厚度方向准直器23、多行X-射线探测器24、DAS(数据采集系统)25、控制X-射线管21或类似物围绕受试者的体轴旋转的旋转部件控制器26和与操作控制台1和拍摄工作台10交换控制信号等的掌控控制器29。通过倾斜控制器27从z方向看去扫描台架20能够向前和向后倾斜±大约30°。除了上述之外,扫描台架20具有通道方向准直器31和形成X-射线滤波器32的线束。
图2是X-射线管21和多行X-射线探测器24几何布置的说明图。
X-射线管21和多行X-射线探测器24围绕旋转中心IC旋转。当将垂直方向定为y方向、水平方向定为x方向和垂直于这两个方向的工作台行进方向定为z方向时,X-射线管21和多行X-射线探测器24的每一个旋转平面表示为xy平面。托架12的移动方向对应于z方向。
X-射线管21产生被称为“锥形线束CB”的X-射线束。当锥形线束CB的中央轴平行于y方向时,设定视角为0°。
多行X-射线探测器24具有例如对应于256行的探测器行。每个具有固定通道宽度的X-射线探测器通道在此以前已经被布置到通道方向上,并且所有通道的X-射线探测器数据总是在数据采集时读取,如图12中所示。在本实施例中,对于数据采集范围和在例如图13中所示的数据获取X-射线探测器的X-射线探测器通道宽度d这种情况下的θ角,每个探测器行具有对应于1024个通道的探测器通道。以这样的方式在多行X-射线探测器24的中央部分设置等于所有通道一半的512个通道使得即使在数据获取X-射线探测器的X-射线探测器通道宽度d/2的情况下也可读取数据。在X-射线探测器通道宽度d/2的情况下,每个探测器行具有关于数据采集范围和θ/2角的对应于1024个通道的探测器通道。
就是说,在多行X-射线探测器24中,将探测穿过受试者传输的X-射线以获得或收集X-射线探测器数据的多个通道分别布置在沿它们被旋转部件15所旋转的方向延伸的通道方向和沿关于其它们被旋转部件15所旋转的旋转轴延伸的行方向这两个方向上。如图13所示,多行X-射线探测器24具有第一区域,其中在通道方向上部署了对应于第一通道宽度d/2的多个通道,和第二区域,其中在通道方向上部署了对应于比第一通道宽度d/2大的第二通道宽度d的多个通道。在多行X-射线探测器24中形成第一区域以便对应从通道方向看去的中央部分,以及在它的外围部分形成第二区域以便从通道方向看去使第一区域置于其间。
多行X-射线探测器24和DAS 25在这种情况下具有两个数据采集模式显示如下。
(1)用于收集或获得具有从通道方向看去宽的数据采集范围(数据采集范围和角θ)内粗糙或粗大的通道宽度d的1024个通道的模式1。
(2)用于获得具有从通道方向看去窄的数据采集范围(数据采集范围和角θ/2)内细的通道宽度d/2的1024个通道的模式2。
在这种情况下,数据采集系统(DAS)25在基于模式1的数据采集时实现对具有通道宽度为d的如图13(a)所示的1到1024个通道上的所有1024个通道的所有行数据采集。在基于模式2的数据采集时,数据采集系统25执行具有通道宽度d/2的在257到768个通道上的所有1024个通道的所有行数据采集。
数据采集系统(DAS)25和多行X-射线探测器24以这样一种方式彼此电连接以便实施基于模式1的数据采集和基于模式2的数据采集。根据模式1和模式2可转换它们之间的连接。
此时在257到768个通道中,如图27中所示,在模式2中分别读取对应于每个通道宽度d/2的X-射线探测器通道的数据。在模式1中,对应于每个通道宽度d/2的X-射线探测器通道的数据分别被加在一起,之后所累加的数据作为对应于通道宽度d的X-射线探测器通道数据被读取。
这样,通过转换可使用模式1的通道方向上宽的数据采集范围和基于从模式2的通道方向看去高分辨率的X-射线探测器通道的窄的数据采集范围,例如,将模式1的通道方向上宽的数据采集范围用于肺部检查而将从模式2的通道方向看去高分辨率的窄的数据采集范围用于心脏检查,临床使用如图28中所示,从而使有效地使用各自的模式成为可能。
就是说,在本实施例中,DAS 25从多行X-射线探测器24上收集X-射线探测器数据并将X-射线探测器数据经过数据采集缓冲器5输出到中央处理单元3。如图13中所示,DAS 25将转换变到获得或收集X-射线探测器数据的区域以便从对应于多行X-射线探测器24的第一区域(数据采集范围2)和第一和第二区域(数据采集范围1)之一的通道中收集X-射线探测器数据。此处,基于从操作者输入到输入装置2中的命令,掌控控制器29传递控制信号。DAS 25将转换变到收集X-射线探测器数据的区域。当进行区域转换使得在多行X-射线探测器24的第一区域收集X-射线探测器数据时,DAS 25从第一区域内所选的布置在通道和行方向上的各个通道中获得或收集X-射线探测器数据并且从那里将同样的数据输出。另一方面,当在多行X-射线探测器24的第一和第二区域获得X-射线探测器数据时,DAS 25从第一和第二区域内所选的布置在通道和行方向上的各个通道中获得X-射线探测器数据并且从那里将同样的数据输出。关于来自第一区域的通道的各个X-射线探测器数据,DAS 25将来自第一区域内彼此邻近的多个通道的X-射线探测器数据累加以便对应第二区域内通道的通道宽度d并将累加的结果从那里输出。就是说,DAS 25将来自第一区域内彼此邻近的两个通道的X-射线探测器数据累加以便与第二区域内通道的通道宽度d相同并将结果从那里输出。将来自第二区域内通道的各个X-射线探测器数据输出而不进行累加。
以下被认为是本发明的另一个X-射线探测器。虽然如图14中所示数据在此之前已经在所有通道和行上读取,但是依据本发明一个实施例中如图15、16、17和18中所示的通道位置还可改变读取行的数量。在这种情况下,从通道方向看去在每个外围部分数量少的X-射线探测器行中读取数据,并且从通道方向看去在内部中央部分数量多的X-射线探测器行中读取数据。这样,在中央部分空间上并且高分辨率地执行数据采集。通过收集位于外围部分一些行方向数据或在行方向上的不连续数据或从行方向上看去行宽度宽和粗糙的数据,数据采集数量与传统的数据采集数量相同或相等,并且从通道方向看去可在中央部分高分辨率地实施数据采集。
从多行X-射线探测器24看去用X-射线照射并收集的X-射线探测器数据通过DAS 25进行A/D转换并且经过集电环30输入到数据采集缓冲器5。输入到数据采集缓冲器5中的数据由中央处理单元3根据存储装置7的程序进行处理,其后数据作为断层照片或断层摄影图像进行图像重建,在监视器6中将其显示。
图3是显示了根据本发明X-射线CT设备100操作要点的流程图。
在步骤S1中,X-射线管21和多行X-射线探测器24首先围绕受试者旋转。当拍摄工作台10上的托架12线性移动时,执行螺旋扫描操作。这样,工作台线性移动z方向位置Ztable(view)加到用视角view,探测器行号j和通道号i表示的X-射线探测器数据D0(view,j,i)以收集X-射线探测器数据。在传统扫描(轴向扫描)中,用放置在固定的拍摄工作台10上的托架12来收集成像数据。在本实施例中,以(2)的细通道间隔p进行数据采集。顺便提及的是,上述视角view是关于受试者从预定位置由旋转部件15在实施扫描中旋转并移动X-射线管21的角度。探测器行号j是多行X-射线探测器24中行方向上布置的每个探测器的号。通道号i是多行X-射线探测器24中通道方向上布置的每个探测器的号。X-射线探测器数据D0(view,j,i)表示当移动到预定视角view的X-射线管21对受试者应用X-射线时,通过允许以多行X-射线探测器24中探测器行号j和通道号i放置探测器收集数据以探测穿过受试者的X-射线。工作台线性移动z方向位置Ztable(view)表示在扫描实施中拍摄工作台10的托架12沿受试者的体轴方向移动的位置。
在确定了受试者的位置后,以这样的方式放置受试者使得能有效地使用位于探测器中央部分处的数据采集通道间隔p并且使受试者落在中央数据采集角θ/2内。
在步骤S2中,在X-射线探测器数据D0(view,j,i)上实施预处理并将其转换到投影数据中。如图4中所示,预处理包括步骤S21的偏移校正,步骤S22的对数变换,步骤S23的X-射线剂量校正和步骤S24的灵敏度校正。
在步骤S3中,在预处理投影数据D1(view,j,i)上施加线束硬化校正。假设在线束硬化校正S3中受到预处理S2的灵敏度校正S24的投影数据定义为D1(view,j,i)并且经线束硬化校正S3后的数据定义为D11(view,j,i),则将线束硬化校正S3表示为,例如,类似下面公式(1)的多项式形式:
D11(view,j,i)=
D1(view,j,i)·(Bo(j,i)+B1(j,i)·D1(view,j,i)+B2(j,i)·(D1(view,j,i)2
...(1)
由于独立于探测器每个j行的线束硬化校正可在此时执行,如果在拍摄或成像状态下各个数据采集系统的管电压是不同的,则可校正探测器的X-射线能量特征集合每行间的差别。
在步骤S4中,对受到线束硬化校正的投影数据D11(view,j,i)进行实施z方向(行方向)滤波的z滤波卷积处理。在步骤S4中,在对各个视角和各个数据采集系统进行预处理后,对受到在行方向上线束硬化校正的多行X-射线探测器D11(ch,row)(其中ch=1-CH,row=1-ROW)的投影数据实施行方向滤波器大小是5行的滤波,例如(w1(ch),w2(ch),w3(ch),w4(ch),w5(ch))。顺便提及的是,在此ch指代通道而row指代行。
然而,上述关系被定义为由公式(2)所给出如下形式:
校正后的探测器数据D12(ch,row)用如下所示的公式(3)表示:
顺便提及的是,当假设通道的最大值为CH并且假设行的最大值为ROW时,它们显示成下面的公式(4)和(5):
D11(ch,-1)=D11(ch,0)=D11(ch,1) ...(4)
D11(ch,ROW)=D11(ch,ROW+1)=D11(ch,ROW+2) ...(5)
当对于每个通道改变行方向滤波器系数时,根据距图像重建中心的距离可控制片层的厚度。由于通常与其重建中心相比,片层厚度在断层照片的外围部分变厚,行方向滤波器系数在中央部分和每个外围部分改变,每个行方向滤波器系数的宽度在中央部分通道的附近很大地改变,而每个行方向滤波系数的宽度在每个外围部分通道的附近较小地改变。结果,在外围部分和图像重建中央部分也还可将片层厚度变的大体平均。
通过用这种方式控制多行X-射线探测器24的中央部分通道和每个外围部分通道的行方向滤波器系数,在中央部分和外围部分处也可控制片层厚度。当依靠行方向滤波器将片层厚度变的稍微厚点时,伪影和噪声都极大地改进。这样,也可控制伪影改进的程度和噪声改进的程度。就是说,可控制三维图像重建断层照片,即,在xy平面的图像质量。作为另一个实施例,通过将行方向(z方向)滤波器系数进行反卷积滤波也可实现片层厚度薄的断层照片。
在步骤S5中,执行重建函数卷积处理。就是说,将数据进行傅立叶变换并乘以重建函数,接着将其进行反傅立叶变换。假设在重建函数卷积处理S5中,将z滤波卷积处理之后的数据定义为D12,将重建函数卷积处理之后的数据定义为D13,并且用于卷积的重建函数是Kernel(j),则重建函数卷积处理用下面公式(6)表示:
D13(view,j,i)=D12(view,j,i)*Kernel(j) ...(6)
就是说,由于重建函数Kernel(j)可对探测器的每个j行单独执行重建函数卷积处理,则可以校正每个行的噪声特征和分辨率特征上的差异。
在步骤S6中,对受到重建函数卷积处理的投影数据D13(view,j,i)可实施三维反投影处理以确定反投影数据D3(x,y)。由于在本发明中用螺旋扫描进行,在三维上将重建的图像的图像重构到垂直于z轴的平面或xy平面上。假设以下重建区域P平行于xy平面。三维反投影处理将随后参考图5进行描述。
在步骤S7中,对反投影数据D3(x,y,z)实施诸如图像滤波卷积、CT-值转化或类似的后处理以获得断层摄影图像或断层照片D31(x,y)。
假设三维反投影断层照片是D31(x,y,z),图像滤波卷积后的数据是D32(x,y,z)并且图像滤波器是Filter(z),则对应于后处理的图像滤波卷积处理可用如下的公式(7)表示:
D32(x,y,z)=D31(x,y,z)*Filter(z) ...(7)
就是说,由于对探测器的每个j行可单独进行图像滤波卷积处理,则能够校正每个行的噪声特征差异和分辨率特征差异。
这样所获得的断层照片显示在监视器6中。
图5是显示了三维反投影处理(图4的步骤S6)细节的流程图。
在本实施例中,将进行了图像重建的图像在三维上图像重建到垂直于z轴的平面或xy平面上。假设以下重建区域P平行于xy平面。
在步骤S61中,注意断层照片的图像重建所必需的所有视角(即,对应于360°的视角或对应于“180°+扇角(fan angle)”的视角)的一个,并且提取对应于重建区域P内每个像素的投影数据Dr。
如图6(a)和6(b)中所示,平行于xy平面的512×512个像素的正方形区域定义为重建区域P,且在y=0处的像素行L0,其平行于x轴,在y=63处的像素行L63,在y=127处的像素行L127,在y=191处的像素行L191,在y=255处的像素行L255,在y=319处的像素行L319,在y=383处的像素行L383,在y=447处的像素行L447,在y=511处的像素行L511分别作为多行。这样,如果提取通过将这些像素行L0到L511投影到从X-射线穿透方向看去的多行X-射线探测器24的平面上而获得的如图7中所示的在线T0到T511上的投影数据,则它们导出像素行L0到L511的投影数据Dr(view,x,y)。然而,x和y对应断层照片的每个像素(x,y)。
依据X-射线管21的X-射线焦点、各个像素和多行X-射线探测器24的几何位置确定X-射线穿透方向。然而,由于X-射线探测器数据D0(view,j,i)的z坐标z(view)被公知为伴随有X-射线探测器数据的工作台线性移动z方向位置Ztable(view),则甚至在处于加速/减速的X-射线探测器数据D0(view,j,i)的情况下也在X-射线焦点和多行X-射线探测器的数据采集几何系统中精确地确定X-射线穿透方向。
顺便提及的是,当从多行X-射线探测器探测器24的通道方向看去一些线衰落时,例如通过将像素行L0投影在从X-射线穿透方向看去的多行X-射线探测器的平面上所获得的线T0,对应的投影数据Dr(view,x,y)被假设为“0”。当它从z方向上看去衰落时,用外插地确定对应的投影数据Dr(view,x,y)。
如图8中所示,用这种方式可提取重建区域P内与各个像素有关的投影数据Dr(view,x,y)。
返回来参考图5,在步骤S62中,用锥形线束重建权重系数乘以投影数据Dr(view,x,y)以产生如图9中所示的投影数据D2(view,x,y)。
在此,锥形线束重建权重参数w(i,j)显示如下。在扇型线束图像重建的情况下,当通过连接X-射线管21的焦点和在view=βa时的重建区域P(xy平面)上的每个像素(x,y)所获得的直线与X-射线束的中央轴Bc间形成的角定义为γ,而它的相对线束通常定义为view=βb,βb导致βb=βa+180°-2γ。
假设穿过重建区域P内像素g(x,y)的X-射线束和其相对的X-射线束与重建平面P所形成的角是如下列公式(8)所指代的αa和αb,它们与取决与此的锥形线束重建权重系数ωa和ωb相乘并加在一起来确定反投影像素数据D2(0,x,y)。
D2(0,x,y)=ωa·D2(0,x,y)_a+ωb·D2(0,x,y)_b ...(8)
然而,分别将D2(0,x,y)_a定义为视角βa的投影数据,将D2(0,x,y)_b定义为视角βb的投影数据。
顺便提及的是,关于彼此相对的线束的锥形线束重建权重系数的总和导致ωa+ωb=1。
将锥形线束重建权重系数ωa和ωb乘以投影数据并加在一起使得它能够减少锥角伪影。
例如,从以下公式中确定的值能用作锥形线束重建权重系数ωa和ωb。
当假设扇型线束角的1/2为γmax时,可使用从以下公式(9)由公式(14)所确定的值。顺便提及的是,ga指代在给定方向上X-射线束的加法/乘法系数,而gb指代对应于其相对线束的X-射线束的加法/乘法系数。
ga=f(γmax,αa,βa) ...(9)
gb=f((γmax,αb,βb) ...(10)
xa=2·gaq/(gaq+gbq) ...(11)
xb=2·gbq/(gaq+gbq) ...(12)
wa=xa2·(3-2xa) ...(13)
wb=xb2·(3-2xb) ...(14)
顺便提及的是,例如,在此q=1。
假设max[]是取最大值的函数,例如,从以下公式(15)和(16)所确定的值可用作ga和gb的例子。
ga=max[0,{(π/2+γmax)-|βa|}]·|tan(αa)| ...(15)
gb=max[0,{(π/2+γmax)-|βb|}]·|tan(αb)| ...(16)
在扇型线束图像重建的情况下,重建区域P上的每个像素进一步被乘以其对应距离系数。当从X-射线管21的焦点到多行X-射线探测器24的探测器行j和通道i的距离在对应于投影数据Dr是r0并且从X-射线管21的焦点到重建区域P上每个像素的距离对应于投影数据Dr是r1时,给出距离系数为(r1/r0)2。
在平行线束图像重建的情况下,重建区域P上的每个像素可单独乘以它的对应的锥形线束重建权重系数w(i,j)。
在步骤S63中,如图10中所示,将投影数据D2(view,x,y)加到事先相关于每个像素进行清除的投影数据D3(x,y)上。
在步骤S64中,对每个断层照片的图像重建所必须的所有视角(即,对应于360°的视角和对应于“180°+扇角”的视角)重复步骤S61到S63,因而获得如图10中所示的反投影数据D3(x,y)。
顺便提及的是,重建区域P可形成为一个如图11(a)和11(b)中所示的圆形区域。
通常用如图12中所示的通道宽度d在多行X-射线探测器24的模式下对全成像视野进行成像。然而,当对观察的成像场内的小物体进行成像或拍摄时,在中央部分通道呈现如图13中所示的高分辨率的模式下执行数据采集,并且通过如上述的图像重建产生断层摄影图像。
由于在此所获的断层摄影图像是以在细通道间隔是d/2且数据采集角是θ/2时在多行X-射线探测器24的部分中所采集的投影数据为基础的,因此可获得对应小成像区域的高分辨率断层摄影图像。
就是说,对应于正常模式下通道宽度的N通道×M行,和对应于高分辨率模式的通道宽度d/2的N通道×M行可根据物体进行转换。在这种情况下,数据采集系统(DAS)25对应N通道×M行并且如果通过转换使用对应于通道宽度d的N通道和对应于通道宽度d/2的N通道,则该数据采集系统是高效的。如果物体和感兴趣区域都很小且落在对应于如图13中所示N通道×通道宽度d/2的成像区域范围内,可基于N通道×通道宽度d/2的高分辨率模式进行断层照片图像重建和断层照片图像显示。此时在图像重建中,在如图3所示的这种图像重建的流程内步骤S5的重建函数卷积处理中准备用于正常模式和高分辨率模式的图像重建函数。当以基于N通道×探测器通道宽度d/2的高分辨率模式进行数据采集时,在数据采集中增加用于采样的尼奎斯特频率。因此,当使用如图39中所示的用于高分辨率模式的图像重建函数进行图像重建时,可获得其质量适当的高分辨率图像。
顺便提及的是,作为X-射线探测器的结构,从通道方向看去的X-射线探测器的中央部分由如图13所示的具有N通道×通道宽度d/2的闪烁器和光电二极管组成。
右侧和左侧外围部分分别由N/4通道×通道宽度d的闪烁器和光电二极管组成。当中央的N通道×通道宽度d/2作为高分辨率模式是红色时,一个接一个独立读取对应于通道宽度d/2的各个通道。
然而,当用作为正常模式以N通道×通道宽度d/2读取所有通道时,用两个通道统一为一个在正常模式下读取具有位于其中央部分的对应于通道宽度d/2的各个通道的X-射线探测器。这样,FET开关被公知为用于通过转换读取X-射线探测器的闪烁器和光电二极管的输出的开关。
然而,如图19中所示,感兴趣区域小且落在具有N通道×通道宽度d/2的高分辨率模式的范围内。然而,例如,当物体的尺寸没有落在N通道×通道宽度d/2的范围内时,在具有N通道×通道宽度d/2和从行方向看去M行×行宽r的高分辨率模式下可对从通道方向看去多行X-射线探测器24内进行数据采集,反之在具有N/4通道×两点×通道宽度d和从行方向看去粗略地具有M/2行×行宽2r的正常模式下可对从通道方向看去其外部进行数据采集,如图18中所示。这种情况下在高分辨率模式下数据采集的量,即,数据采集的数量用关于正常模式中N通道×M行的以下方式表示:
N通道×M行+N/4通道×2×M/2行=1.25×N通道×M行。
通道的数量增至1.25倍。如果采取等于正常模式下的每个通道的数据采集时间,则总的数据采集时间增至1.25倍。如果希望数据采集时间符合正常模式下同样的数据采集时间,则要是每个通道的数据采集时间减少到正常模式下数据采集时间的0.8倍,就会使数据采集时间与正常模式下的数据采集时间相同。
这样,即使在从通道方向看去在多行X-射线探测器24的外围部分上粗略地实施数据采集,并且用行方向上做的很粗的投影数据的分辨率进一步进行图像重建,也不会对在高分辨率模式下感兴趣的中央区域的成像施加影响。
当在对应于如图18所示的一个视图的一行的数据中混合位于探测器通道方向外围部分的探测器通道宽度d和位于探测器通道方向中央部分的探测器通道宽度d/2时,在步骤S5的重建卷积处理中对高分辨率模式的图像重建函数进行卷积。这样,如果对应于探测器通道宽度d/2的用于高分辨率模式的图像重建函数叠加在外围部分的探测器通道宽度d的部分之上,从而产生不必要的图像噪声。对于精确来说这是不希望得到的。如下操作可将其避免。
根据目前常用的方法,将投影数据和重建函数卷积傅立叶变换到频率空间。在频率空间中乘以投影数据和重建函数的结果进行反傅立叶变换,跟着将其恢复到实空间中。然而,目前的方法不能使用仅对于探测器外围部分的投影数据不同的正常模式中的图像重建函数。
因此,在实空间中可改变在探测器外围部分和探测器中央部分所叠加或卷积的图像重建函数。以下流程处理可如图40中所示进行。
在步骤S51中,i,j,view=1。然而,假设i是1到1024的整数,假设j是1到256的整数,及假设view是1到1000的整数。
在步骤S52中,基于投影数据D12(view,j,i)确定探测器通道宽度是否是d/2并且是否被放置在高分辨率模式中。如果发现应答是是(YES),则流程处理进行到步骤S53。如果发现应答是否(NO),则流程处理进行到步骤S54。
在步骤S53中,对高分辨率模式的重建函数Kernel HR(j)进行卷积以确定重建函数卷积后的投影数据D13(view,j,i)。
在步骤S54中,对正常模式的重建函数Kernel LR(j)进行卷积以确定重建函数卷积后的投影数据D13(view,j,i)。
在步骤S55中,确定是否i=1024。如果发现应答是YES,则流程处理进行到步骤S56。
在步骤S56中,确定是否j=256。如果发现应答是YES,则流程处理进行到步骤S57。
在步骤S57中,确定是否view=1000。如果发现应答是YES,则流程处理结束。
在步骤S58中,执行i=i+1并且流程处理返回到步骤S51。
在步骤S59中,实施j=j+1并且流程处理返回到步骤S51。
在步骤S60中,实施view=view+1并且流程处理返回到步骤S51。
当用行方向上捆束(bundled)的数据抑止所采集数据的数量时,对从通道方向看去的多行X-射线探测器24的外围部分的外部进行数据采集。然而,如图17所示,可通过具有行宽度r的M/2行如此稀薄地对从通道方向看去的多行X-射线探测器24的外围部分的外部进行数据采集以至于X-射线探测器通道本身被放置在通道宽度d,M行和行宽r下而数据采集被放置在N行跳过样式中(例如,将数据采集安排在图17中N=2的两行跳过样式)。在这种情况下,即使当使用在它们邻近行的数据采集的投影数据对稀薄的且跳过行的外部外围部分的投影数据进行图像重建也不会有问题发生。为了精确,移动行方向上的投影数据。然而,如果在z方向上没有突然地改变投影数据轮廓的区域总和,则不会对高分辨率模式下断层摄影图像中央处的感兴趣区域的成像施加影响。
类似的,如作为图17另一种选择的图16中所示,在行方向将稀薄部分集中起来并且将数据采集区放置在行方向的一侧上。可选择的,数据采集区可放置在如图15中所示的行方向中央。如果在z方向上没有突然改变投影数据轮廓区域的总和,则不会对高分辨率模式下断层摄影图像的中央感兴趣区域的成像施加影响。
这样,心脏检查,作为希望细微观看感兴趣中央区域的一种临床应用被提出来了。在心脏检查中首先需要空间分辨率。图28显示了包含心脏的肺部区域成像的例子。当希望对包括肺部区域的心脏成像或拍摄时,使用X射线采集范围1在包括受试者整个肺部区域范围内执行数据采集、图像重建和图像显示。关于在此时的数据采集,在具有N通道×X-射线探测器通道宽度d的X-射线数据采集范围1上执行对应于M行的数据采集。在特别是在专门用于心脏成像时,使用X-射线数据采集范围2在以受试者心脏为中央的范围内执行数据采集、图像重建和图像显示。关于此时的数据采集,用N通道×X-射线探测器通道宽度d在X-射线数据采集范围2上执行对应于M行的数据采集。其次,需要作为心脏检查特征的时间分辨率。为此,在空间上以高分辨率模式和满意时间分辨率下高速读取感兴趣的中央区域。进一步,以高速度旋转扫描台架20的旋转部分以加快用于数据采集的扫描时间,并且用同样数量的数据采集点执行数据采集。
当行方向上数据采集范围根据如图15或16所示的通道位置进行区分时,用于控制X-射线照射区域的X-射线的方法考虑如下:
(1)基于通道方向X-射线滤波器的方法
(2)基于形成X-射线滤波器线束的方法
如图24中所示,可控制在行方向上的X-射线照射区域以便根据由通道方向准直器31在行方向上分成两个的通道位置进行区分。这样,X-射线照射区域可与数据采集范围匹配。顺便提及的是,图24是显示了作为图2中视线的从X-射线管21延伸到多行X-射线探测器24方向的图。
形成X-射线滤波器32的线束和形成X-射线滤波器32a的增加的线束如图25所示进行交迭。可控制在行方向上的X-射线照射区域以便根据通道位置进行区分。这样,X-射线照射区域能够与数据采集范围对准。在图25中,图25(a)是显示了作为图2中视线的从X-射线管21延伸到多行X-射线探测器24方向的图,图25(b)是显示了在图25(a)的A-B截面上形成X-射线滤波器32的线束和形成X-射线滤波器32a的增加的线束的剖面图,图25(c)是显示了分别在线A-B上各个部分中的布局关系的侧视图。
顺便提及的是,这种情况下存在需要将X-射线照射区域与感兴趣的诊断区域进行匹配。感兴趣的诊断区域可事先在扫描前如图26(a)和26(b)所示的探查图像中设定。
图20显示了正常模式下读取的情形。让我们假设通道数量是L,行数量是M,扫描时间是T,以及每个扫描和旋转的视图数量是L。作为一个真实扫描的例子,认为N=1024通道,M=256行,扫描时间T=0.5s,以及每个扫描和旋转的视图数量L=1000。
在正常模式下数据传送率用所给的以下公式(17)表示:
在正常模式中,在该速度下读取数据。
另一方面,图21显示了特别专用于心脏的高速读取模式的想法。
通过通道方向准直器23控制X-射线并将其放射到对应于多行X-射线探测器24的中央N/K通道的部分(例如,在N=1024和K=2的512通道)。
当对应于中央512通道的部分是n倍时,即,n=0.4,在nT=0.5秒×0.4=0.2秒下进行扫描。当在这种情况下读取对应于中央512通道的部分时,在由以下公式(18)给出的数据传送率下进行数据采集,并执行图像重建和图像显示。
顺便提及的是,由于当将图21所示的方法用于心脏成像时使用通道方向准直器23,如果使用预测缺失投影数据并对其进行校正的图像重建算法执行图像重建能获得质量更好的断层摄影图像,该算法在以后描述的实施例2中进行展示。
从通道方向看去对应于中央512通道的部分在上面的例子中并没有设定为高分辨率模式。然而,当对应于X-射线探测器通道宽度d的512通道部分在高分辨率模式下设定为具有双1024通道的X-射线探测器通道宽度d/2时,对于心脏检查它更进一步的有效。
当在此时读取位于中央部分的通道时,可由以下公式(19)所给定的数据传送率下执行数据采集,并进行图像重建和图像显示。
这样,众所周知的是具有多个数据传送模式的X-射线CT装置适于各种应用并且是高效的。
图22显示了使用形成滤波器32的X-射线线束执行在对应于从通道方向看去中央部分的通道上的高速读取的例子。
通过形成滤波器32的X-射线线束在外围部分衰减X-射线。原理上将X-射线放射到对应于多行X-射线探测器24的中央N/K通道的部分(例如,在N=1024和K=2的512通道)。
当对应于中央512通道的部分是n倍时,例如,n=0.4,在nT=0.5秒×0.4=0.2秒下进行扫描。当在这种情况下读取对应于中央512通道的部分时,在由以下公式(20)给出的数据传送率下进行数据采集,并执行图像重建和图像显示。
顺便提及的是,由于当将图22所示的方法用于心脏成像时使用形成X-射线滤波器32的线束,如果使用预测缺失投影数据并对其进行校正的图像重建算法执行图像重建能获得质量更好的断层摄影图像,该算法在以后描述的实施例2中进行展示。
从通道方向看去对应于中央512通道的部分在上面的例子中并没有设定为高分辨率模式。然而,当对应于X-射线探测器通道宽度d的512通道部分在高分辨率模式下设定为具有双1024通道的X-射线探测器通道宽度d/2时,对于心脏检查更进一步的有效。
当在此时读取位于中央部分的通道时,可由以下公式(21)所给定的数据传送率执行数据采集,并进行图像重建和图像显示。
这样,众所周知的是具有多个数据传送模式的X-射线CT装置适于各种应用并且是高效的。
图23显示了一种方法,其中由于形成X-射线滤波器32的线束不同于通道方向准直器23并且X-射线被少许放射到外围部分中,收集这种方法所获的数据并且将质量更好的断层摄影图像进行图像重建。
原理上将X-射线放射到对应于多行X-射线探测器24的中央N/K通道的部分(例如,在N=1024和K=2的512通道)上。当对应于中央512通道的部分是n倍时,例如,n=0.4,在nT=0.5秒×0.4=0.2秒下进行读取数据。顺便提及的是,此时扫描台架20的每旋转率是0.2秒且扫描速度或速率是0.2秒。例如,在外围部分512通道的情况下k=2,在kT=0.5秒×2=1.0秒下执行数据读取。
当读取对应于中央512通道的部分时,在由以下公式(22)给出的数据传送率下进行数据采集,并执行图像重建和图像显示。
在外对围部分对512通道进行数据读取时,在由以下公式(23)给出的数据传送率下进行数据采集,并执行图像重建和图像显示。
顺便提及的是,由于在将图23所示的方法用于心脏成像中使用通道方向准直器23,如果使用预测缺失投影数据并对其进行校正的图像重建算法进行图像重建能获得质量更好的断层摄影图像,该算法在实施例2中进行展示。
从通道方向看去对应于中央512通道的部分在上面的例子中并没有设定为高分辨率模式。然而,当对应于X-射线探测器通道宽度d的512通道部分在高分辨率模式下设定为具有双1024通道的X-射线探测器通道宽度d/2时,对于心脏检查它更进一步的有效。
当在此时读取位于中央部分的通道时,可由以下公式(24)所给定的数据传送率下执行数据采集,并进行图像重建和图像显示。
顺便提及的是,此时虽然有关对应于中央部分的512通道的数据在图像重建时在各个视图快速地更新,但是有关在外围部分512通道的数据却以低速度更新。然而,由于在外围部分的数据用于防止CT值变化以及伪影的发生,不会对感兴趣区域的高分辨率断层摄影图像的质量产生影响,即使它们不是稍微准确的适当的数据。
这样,存在对每个数据采集范围不同的用于数据采集采样周期的数据读取模式能够适应各种应用并产生效果。
如上所述,根据本实施例的X-射线CT装置100具有向受试者照射X-射线的X-射线管21、探测从X-射线管21照射并穿过受试者的X-射线的多行X-射线探测器24和以这样的方式移动X-射线管21和多行X-射线探测器24使得它们围绕受试者旋转的旋转部件15。基于通过实施这样的扫描使得通过旋转部件15将围绕受试者旋转的X-射线管21向受试者照射X-射线并用多行X-射线探测器24探测穿过受试者的X-射线所获得的投影数据,对受试者的断层摄影图像进行图像重建。此处,在多行X-射线探测器24中,将探测穿过受试者以产生X-射线探测器数据的多个通道分别设置在沿由旋转部件15对其进行旋转的方向延伸的通道方向上,和沿此时由旋转部件15对其进行旋转的旋转轴延伸的行方向上。多行X-射线探测器24具有在通道方向上布置对应第一通道宽度d/2的多个通道的第一区域,和布置对应比第一通道宽度d/2大的第二通道宽度d的多个通道的第二区域。在多行X-射线探测器24中,这样形成第一区域使得从通道方向看去对应中央部分,以及在外围部分这样形成第二区域以便使第一区域置入其间。因此,本实施例能够有效且容易地根据成像条件以所希望的数据采集率获得所希望的分辨率的图像。
(实施例2)
虽然实施例1显示了在高分辨率下对小成像区域进行拍摄或成像的实施例,但是实施例2显示了进一步实现低放射曝照的实施例。除了对X-射线CT装置100的操作不同于根据实施例4的操作外,实施例2与实施例1相类似。因此将不再解释重复的部分。
图29是显示了X-射线CT装置100操作要点的流程图。本实施例将解释增加通道方向准直器并因而依据人们希望重建的FOV的尺寸进行控制的实施例。
虽然屏蔽板系统(板状准直器系统)(图30(a))和屏蔽圆柱系统(旋转轴偏心柱状准直器系统)(图30(b))被认为是准直器,但是任何一个均可用于本实施例中。通过允许DAS 25读取z通道的数据已经控制Z方向(片层厚度方向)准直器控制,而在通道方向准直器31上,事先确定照射到多行X-射线探测器24的X-射线的位置,该位置依据X-射线数据采集系统的角度β(视角β)进行确定,以及事先确定将要拍摄的感兴趣区域的位置和尺寸。根据这些位置前馈控制通道方向准直器31的孔径位置和宽度。进一步,根据需要基于执行投影数据采集的DAS 25(参考图31(a)和31(b))的每个主探测器通道的数值执行在通道方向上的反馈控制。
随着用于DAS控制的CPU和用于准直器控制的CPU的性能的进步,认为用于读取对应于多行X-射线探测器24的数据的主要探测器通道和用于对通道方向准直器的孔径执行反馈控制的计算基本不会产生问题。当关于肥胖患者不能确信X-射线数据的SN时,根据基于事先视野成像的位置/尺寸所预测的通道方向准直器位置可仅执行反馈控制。
驱动系统或类似系统,例如在这种情况下控制准直器操作的脉冲马达,也被认为具有足够的响应速度。
在图29所示的整个流程中,小的成像区域可以高分辨率进行拍摄并通过后继流程具有好的图像质量。
在步骤P1中,首先执行探查图像的数据采集。
在步骤P2中,在探查图像上设置将要成像或拍摄的区域。
在步骤P3中,确定在各个z位置将要成像的轮廓区域。
在步骤P4中,在与将要成像的感兴趣区域匹配的通道方向上控制通道方向准直器。
在步骤P5中,执行扫描以进行数据采集。
在步骤P6中,执行投影数据的预处理以获得在受到探查扫描的各个z位置的所有轮廓区域信息。通道方向准直器预测从通道方向看去在每个外围部分所缺失的投影数据部分并对其进行增加。
在步骤P7中,使用增加有所缺失部分的投影数据执行图像重建处理和图像显示。
当通过如图32所示的通道方向X-射线准直器屏蔽除了在将被成像的部分的X-射线时,需要预测对应于屏蔽部分的投影数据。此处如图32(b)中所示,根据通过允许在它们被旋转360°时各个方向的Sil,Sir,Sjl和Sjr的区域进行椭圆近似而获得的轮廓之间位置关系和成像区域预测投影数据。顺便提及的是,可将投影数据以如图32(c)所示以三角形近似或圆形近似的形式加到两个侧面上。
用图33的流程图解释通道方向X-射线准直器的前馈控制。
在步骤C1中,通过根据包括X-射线管21、多行X-射线探测器24和DAS25的X-射线数据采集系统的角度β(视角β)和感兴趣的成像区域的尺寸和位置(例如中心是(xo,yo),半径是R的感兴趣的圆形区域)进行计算确定将要辐射X-射线的多行X-射线探测器24的角度范围(从最小照射通道γmin到最大照射通道γmax)或通道范围。
在步骤C2中,将通道方向准直器(其可是椭圆柱状准直器和屏蔽板类准直器)做成从最小照射通道γmin到最大照射通道γmax打开。
在步骤C3中,确认是否完成通道方向准直器控制和对应于所有视图的数据采集。
顺便提及的是,在图34中显示了最小照射通道γmin和最大照射通道γmax、包括X-射线管21、多行X-射线探测器24和DAS 25的数据采集系统以及通道方向准直器。从它的关系中可以理解,以x=FCD·sin θ,y=FCD·cosθ(此处θ:视角,FCD(焦心点距离))表示X-射线球管的位置(x,y)。
在视角=0°时的感兴趣成像区域和最小照射通道及最大照射通道之间的关系像图35中所述如下。
例如,当假设感兴趣的圆形成像区域的位置是(xo,yo)时,假设它的半径是R,假设它的视角是0°,即,X-射线焦点设置在(0,FCD),上述关系显示成用以下公式(25)、(26)和(27)表示(此处FCD:焦心点距离)。以下公式(28)和(29)源自公式(25)、(26)和(27)。如以下公式(30)和(31)所表示的,此时γ的最大值变成γmax而γ的最小值变成γmin。
在视角=β时的感兴趣成像区域和最小照射通道及最大照射通道之间的关系像图36中所述如下。
例如,当假设感兴趣的圆形成像区域的位置是(xo,yo)时,假设它的半径是R,假设它的视角是0°,就是说,X-射线焦点设置在(FCD·sinβ,FCD·cosβ),上述关系显示成用以下公式(32)、(33)和(34)表示(此处FCD:焦心点距离)。以下公式(35)源自公式(32)、(33)和(34)。从以下公式(36)和(37)中,此时γ的最大值变成γmax而γ的最小值变成γmin,如在以下公式(38)和(39)中所示。
xo′=xo·cosβ-yo·sinβ ...(36)
yo′=xo·sinβ+yo·cosβ ...(37)
下面在图37中显示通道方向X-射线准直器的反馈控制。
在步骤C1中,类似于图33的步骤C1方式,通过根据包括X-射线管21、多行X-射线探测器24和DAS 25的X-射线数据采集系统的角度β(视角β)和感兴趣的成像区域的尺寸和位置(例如中心是(xo,yo),半径是R的感兴趣的圆形区域)进行计算确定将要辐射X-射线的多行X-射线探测器24的角度范围(从最小照射通道γmin到最大照射通道γmax)或通道范围。
在步骤C2中,用类似于图34的步骤C2的方式将通道方向准直器(其可是椭圆柱状准直器和屏蔽板类准直器)做成从最小照射通道γmin到最大照射通道γmax打开。
在步骤C3中,在观察DAS 25数据的同时确定X-射线所照射数据的范围。假设从Chmin到Chmax的范围是X-射线照射的数据输入范围,则确认是否它对应步骤C1中所确定的最小照射通道γmin或最大照射通道γmax。
如果数据输入范围落在±ε的小误差范围内,则没有问题发生。然而,当它超出这一误差范围时,反馈控制进行到步骤C4。
在步骤C4中,假设为γmin-Chmin·Chang=Δγmin,γmax-Chmax·Chang=Δγmax,将校正量或值Δγmin和Δγmax加到控制量或值中。其后,反馈控制进行到步骤C5。
在步骤C5中,确认是否对应于所有视图的数据采集完成。如果没有完成,则反馈控制返回到步骤C1,此处连续地执行通道方向准直器控制和数据采集。
在这种情况下,从轮廓区域和其通道方向上的宽度做出椭圆近似。如图32中所示,通过根据受到椭圆近似的轮廓和待成像的区域之间的位置关系的在第ith个片层的各个方向上的蒙片(masked)X-射线数据认识和了解加到待成像部分或区域的左侧和右侧的投影数据Sil和Sir。通过将Sil和Sir加到投影数据的右侧和左侧并对该投影数据进行图像重建,能获得质量更好的断层摄影图像。
如上所述,提供了通道方向准直器31的本发明,所述通道方向准直器31能屏蔽照射受试者的X-射线从而调整在通道方向看去的X-射线照射范围。在扫描受试者时,掌控控制器29以这样的方式控制通道方向准直器31的孔径位置使得将来自X-射线管21的X-射线与由操作者输入到输入装置2中的受试者感兴趣区域相关地照射。基于由多行X-射线探测器24所采集的关于感兴趣区域的X-射线探测器数据,中央处理单元3允许通道方向准直器3近似关于彼此围绕感兴趣区域的区域的缺失数据并对其校正。其后,基于近似的X-射线探测器数据对受试者的断层摄影图像进行图像重建。因此有可能实现放射的低曝照。
(实施例3)
实施例3显示了将形成X-射线滤波器32的线束用于其中的例子。除了X-射线CT设备100的操作不同于实施例1外本实施例与实施例1类似。因此,将不再解释重复的部分。
虽然使用通道方向X-射线准直器31已经对实施例2进行了描述,但是即使在如图38中所示的使用形成X-射线滤波器32的线束也能产生相似效果。
图38(a)显示了通道方向X-射线准直器31正常位置的方式,即,在通道方向上其传播量是0。
图38(b)和38(c)分别显示了其中形成X-射线滤波器32的线束的传播量被认做是d1和d2。在这种情况下,可控制连接感兴趣区域中心和X-射线焦点的直线以便将其叠加到这样的直线上以便形成X-射线滤波器32的线束的X-射线穿透路径变的最短。
为了允许它们彼此重叠,使用源于公式(38)和(39)的以下公式(40)。
γmean=(γmax+γmin)/2 ...(40)
假设X-射线焦点和形成X-射线滤波器32的线束之间的距离是如图38(a)中所示的D,它在公式(41)中表示如下:
di=D·tan(γmean) ...(41)
此处di=d1或d2
根据上面X-射线CT设备100,利用本发明的X-射线CT设备可实现即使通过螺旋扫描和传统扫描(轴向扫描)就能获得高分辨率的这种CT设备。
可实现X-射线CT设备,其在仅对有限的感兴趣区域进行拍摄和图像重建时能获得高分辨率。
可实现X-射线CT设备,其当仅对有限的感兴趣区域用X-射线照射并进行拍摄时能够在低放射曝照下获得高分辨率。
具体的描述的是,提供根据本实施例的X-射线CT设备100,该设备100具有能在通道方向上调节X-射线照射分布的形成X-射线滤波器32的线束。掌控控制器29以这样的方式控制X-射线滤波器32的位置使得在对受试者进行拍摄时将来自X-射线管21的X-射线与由操作者输入到输入装置2内的受试者感兴趣区域相关地照射,因此,可实现低放射曝照。
顺便提及的是,虽然在本实施例中使用多行X-射线探测器,但是即使在使用单行X-射线探测器的X-射线CT设备情况下也能产生相似的效果。
虽然在本发明的图像重建中使用三维图像重建方法,但是也可对数据使用基于已知的Feldkamp方法的三维图像重建方法或其他三维图像重建方法或二维图像重建方法,并且在这种情况下能产生相似的效果。
虽然本发明使用螺旋扫描或传统扫描(轴向扫描),但是即使在锥形扫描的情况下也能产生相似的效果。
在不脱离本发明精神和范围内可构想处许多相差很大的不同实施例。应当理解的是除了如附属权利要求中进行的限定,本发明并非限于说明书中所描述的具体实施例。
附图说明
图1
100...X-射线CT设备
1...操作控制台
2...输入装置,3...中央处理单元,5...数据采集缓冲器,6...监视器,7...存储装置,20...扫描台架
15...旋转部件,21...X-射线管,22...X-射线控制器,23...片层厚度方向准直器,24...多行X-射线探测器,26...旋转部件控制器,29...掌控控制器,30...集电环,31通道方向准直器,32...形成X-射线滤波器的线束,
10...拍摄工作台,12...托架
图2
X-射线焦点
21...X-射线管,31...通道方向准直器
CB...X-射线束
p...重建区域
dp...探测器平面
IC...旋转通道方向的中心
24...多行X-射线探测器
图3
开始
S1...数据采集
S2...预处理
S3...线束硬化校正
S4...Z滤波卷积处理
S5...重建函数卷积处理
S6...三维反投影处理
S7...后处理
结束
图4
步骤S2
开始
S21...偏移量校正
S22...对数转换
S23...X-射线剂量校正
S24...灵敏度校正
结束
图5
步骤S6
开始三维反投影处理
S61...提取对应于重建区域P的各个像素的投影数据Dr
S62...用锥形线束重建权重系数乘以各个投影数据Dr以产生反投影数据D2
S63...将与像素有关的反投影数据D3与反投影数据D2相加
S64...对图像重建所必须的所有视角的反投影数据D2进行了相加?
结束
图6
(a)21...X-射线管
重建区域(xy平面)
原点
(b)21...X-射线管
重建区域
xz平面
图7
24...多行X-射线探测器
探测器行方向
通道方向
图8
重建区域
图9
重建区域
图11
(a)21...X-射线管
重建区域
(b)21...X-射线管
重建区域
xz平面
24...多行X-射线探测器
图12
N通道
行方向 M行
通道方向
24...多行X-射线探测器
图13
(a)N通道X通道宽度d
N通道X通道宽度d/2
通道宽度d 通道宽度d/2 通道宽度d
M行
1通道 257通道 768通道 1024通道
数据采集范围2(角θ/2)
数据采集范围1(角θ)
行方向
通道方向 24...多行X-射线探测器
(b)21...X-射线管
受试者
感兴趣区域
24...多行X-射线探测器
通道宽度d,通道宽度d/2,通道宽度d
图14
N通道,通道宽度d
通道宽度d
M行行 宽度r
行方向 数据采集范围1
通道方向
24...多行X-射线探测器
图15
N/4通道 N/4通道
通道宽度d 通道宽度d
M行,行宽度r M行,行宽度r
N通道,通道宽度d/2,
M行,行宽度r
M行 M/2行
行宽度r
行方向
通道方向
24...多行X-射线探测器
图16
N/4通道, N/4通道,
通道宽度d 通道宽度d,
M行,行宽度r M行,行宽度r
N通道,通道宽度d/2
M行,行宽度r
M行 M/2行
行宽度r
行方向
通道方向
24...多行X-射线探测器
图17
N/4通道, N/4通道,
通道宽度d 通道宽度d,
M行,行宽度r M行,行宽度r
N通道,通道宽度d/2
M行,行宽度r
M行
行宽度r
行方向
通道方向
24...多行X-射线探测器
图18
N/4通道, N/4通道,
通道宽度d 通道宽度d,
M/2行,行宽度2r M/2行,行宽度2r
N通道,通道宽度d/2
M行,行宽度r
M行 行宽度2r
行方向 行宽度r
通道方向
在行方向上读取行宽度2rX M/2行
在行方向上读取行宽度rX M行
24...多行X-射线探测器
图19
21...X-射线管
受试者
感兴趣区域
24...多行X-射线探测器
通道宽度d,通道宽度d/2,通道宽度d
图20
24...多行X-射线探测器
通道方向(N通道)
行方向
(M行)
图21
21...X-射线管,23...通道方向准直器,
通道方向(N通道)
24...多行X-射线探测器
行方向(M行),N/K通道
图22
21...X-射线管,32...形成线滤波器的X-射线束,通道方向(N通道)
24...多行X-射线探测器
行方向(M行),N/K通道
图23
21...X-射线管,32...形成滤波器的X-射线束,通道方向(N通道)24...多行X-射线探测器
位于中央部分的通道:在nT秒数据采集
位于外围部分的通道:在kT秒数据采集
行方向(M行),N/K通道
数据采集周期2T秒
数据采集周期T/2秒
图24
31...通道方向准直器
X-射线照射范围
31...通道方向准直器
23...片层厚度方向准直器
31...通道方向准直器
23...片层厚度方向准直器
31...通道方向准直器
图25
X-射线照射区域
21...X-射线管
图26
(a)头,肺部区域,肝,感兴趣区域,感兴趣区域
(b)头,肺,感兴趣区域,感兴趣区域,肝
图27
闪烁器+光电二极管
通道方向
在高分辨率模式下读取
在正常模式下读取
图28
21...X-射线管
心
肺
12...托架
24...多行X-射线探测器
X-射线数据采集范围1(N通道,M行)
X-射线数据采集范围2(N通道,M行)
图29
P1...采集探查图像
P2...设置待成像的区域
P3...在各个z位置确定待成像的轮廓区域
P4...控制通道方向上通道方向准直器与待成像的感兴趣区域匹配
P5...扫描以执行数据采集
P6...对投影数据执行预处理,并从探查扫描的所有轮廓区域信息预测投影数据缺失部分并增加它
P7...图像重建处理和显示
探查图像
获得轮廓数据
待成像区
成像
从轮廓预测投影数据并增加缺失部分
图30
(a)21...X-射线管,32...形成X-射线滤波器的线束,31...通道方向准直器,24...多行X-射线探测器
(b)21...X-射线管,32...形成X-射线滤波器的线束,31...通道方向准直器,24...多行X-射线探测器
(c)21...X-射线管,32...形成X-射线滤波器的线束,X-射线束,24...多行X-射线探测器
图31
(a)21...X-射线管,31...通道方向准直器,肺部区域,感兴趣区域(重建区域),24...多行X-射线探测器
粗的X-射线探测器,X-射线输出,细的X-射线探测器,
通道方向
(b) (b)21...X-射线管,31...通道方向准直器,肺部区域,感兴趣区域(重建区域),24...多行X-射线探测器
粗的X-射线探测器,细的X-射线探测器,X-射线输出,粗的X-射线探测器,
通道方向
图32
(a)探查图像 0°方向探查
jTH片层 特定成像区域
jTH片层 肺部区域
矩形
jTH片层
jTH片层 通道方向
(b)待成像区域
做出椭圆近似轮廓
预测是否等于Sir
预测是否等于Sir
图33
开始通道方向准直器控制
C1...计算从X-射线数据采集系统和感兴趣的成像区域(中心(xo,yo),半径R)的角度β(视角β)的每个视图的照射通道范围(从最小照射通道γmin到最大照射通道γmax)。
C2...在每个视图从γmin到γmax打开通道方向准直器的通道方向孔径
C3...所有视图完成了吗?
结束
图34
31...通道方向准直器
X-射线管位置
断层照片重建区域
待成像区域事先设为感兴趣区域
X-射线束 X-射线照射范围
N通道 1通道
照射最小通道
照射最大通道
图35
X-射线焦点
X-射线束
感兴趣区域
图36
X-射线焦点
图37
开始通道方向准直器控制
C1...计算从X-射线数据采集系统和感兴趣的成像区域(中心(xo,yo),半径R)的角度β(视角β)的每个视图的照射通道范围(从最小照射通道γmin到最大照射通道γmax)
C2...在每个视图从γmin到γmax打开通道方向准直器的通道方向孔径
C3...通过观看从DAS 25的数据输入范围Chmin到其Chmax,X-射线照射范围是否在从γmin±ε到γmax±ε的范围内?
C4...将校正值Δγmin和Δγmax加入到控制值如γmin-Chmin·Chang=Δγmin和γmax-Chmax·Chang=Δγmax
C5...所有视图完成了吗?
结束
图38
(a)21...X-射线管,32...形成X-射线滤波器的线束,24...多行X-射线探测器
(b)形成X-射线滤波器32的线束的穿透路径的短部分
肺部区域
感兴趣区域的中央
感兴趣区域(重建区域)
粗的X-射线探测器,X-射线输出,细的X-射线探测器,粗的X-射线探测器
通道方向
(c)肺部区域
感兴趣区域的中央
感兴趣区域(重建区域)
形成X-射线滤波器32的线束的穿透路径的短部分
粗的X-射线探测器,细的X-射线探测器,X-射线输出,粗的X-射线探测器
通道方向
图39
图像重建函数
在探测器通道间距d处用于正常模式的图像重建函数
在探测器通道间距d/2处用于高分辨率模式的图像重建函数
图40
开始
步骤S51...i,j,view=1
其中i=1到1024的整数
j=1到256的整数
view=1到1000的整数
步骤S52...从投影数据D12(view,j,i)判断是否是探测器通道宽度d/2并且放置在高分辨率模式下?
步骤S53...卷积用于高分辨率模式的重建函数Kernel HR(j)并确定重建函数卷积后的投影数据D13(view,j,i)。
步骤S54...卷积用于正常模式的重建函数Kernel LR(j)并确定重建函数卷积后的投影数据D13(view,j,i)。
步骤S55...i=1024?
步骤S56...j=256?
步骤S57...view=1000?
步骤S58...执行i=i+1
步骤S59...执行j=j+1
步骤S60...执行view=view+1
结束。
Claims (10)
1.一种X-射线CT设备(100)包括:
X-射线数据采集装置(20),其含有X-射线发生器(21),和设置成与X-射线发生器(21)为相对关系并探测X-射线的多行X-射线探测器(24),或者能围绕在X-射线发生器(21)和探测器(24)之间所设的旋转中心旋转的矩阵结构的二维X-射线区域探测器,从而采集穿过处于X-射线发生器(21)和探测器(24)之间的受试者的X-射线的投影数据;
图像重建装置(3),其对从X-射线数据采集装置(20)所采集的投影数据进行图像重建;以及
显示装置(6),其显示经过图像重建的图像,
其中以这样的方式配置X-射线数据采集装置(20)使得从通道方向看去在探测器(24)的每个外围部分的探测器通道宽度d2相对于从通道方向看去在探测器中央部分的探测器通道宽度d1变成d1<d2,或者以通道方向看去从探测器(24)的中央部分到其外围部分所设置的多个探测器宽度(d1,d2,...di,...dn-1,dn)满足d1≤d2≤...≤di≤...≤dn-1≤dn。
2.一种X-射线CT设备(100)包括:
X-射线数据采集装置(20),其含有X-射线发生器(21)和设置成与X-射线发生器(21)为相对关系并探测X-射线的多行X-射线探测器(24),或者能围绕在X-射线发生器(21)和探测器(24)之间所设的旋转中心旋转的矩阵结构的二维X-射线区域探测器,从而采集穿过处于X-射线发生器(21)和探测器(24)之间的受试者的X-射线的投影数据;
图像重建装置(3),其对从X-射线数据采集装置(20)所采集的投影数据进行图像重建;以及
显示装置(6),其显示经过图像重建的图像,
其中X-射线数据采集装置(20)具有从在探测器(24)通道方向上宽的数据采集范围l1到在通道方向上窄的数据采集范围ln的多个数据采集范围l1≥l2≥...≥li≥...ln-1≥ln,并且X-射线数据采集装置以这样的方式配置使得数据采集范围对每个数据采集是可转换的。
3.根据权利要求1或2的X-射线CT设备(100),其中当从探测器(24)的通道方向看去在窄的数据采集范围执行数据采集时,X-射线数据采集装置(20)在从通道方向看去在探测器(24)中央部分的探测器通道宽度窄的部分处执行数据采集。
4.根据权利要求1到3任何一项的X-射线CT设备(100),其中X-射线数据采集装置(20)具有执行数据采集的多个通道。
5.根据权利要求1到3任何一项的X-射线CT设备(100),其中X-射线数据采集装置(20)具有执行数据采集的多个通道和多个视图。
6.根据权利要求1到5任何一项的X-射线CT设备(100),其中X-射线数据采集装置(20)具有执行数据采集的多个行,并且行的数量根据每个通道位置是不同的。
7.根据权利要求1到6任何一项的X-射线CT设备(100),其中X-射线数据采集装置(20)具有执行数据采集的多个数据采集采样周期。
8.根据权利要求1到7任何一项的X-射线CT设备(100),其中X-射线数据采集装置(20)具有执行数据采集的多个数据采集采样周期,并且数据采集采样周期根据通道位置是不同的。
9.根据根据权利要求1到8任何一项的X-射线CT设备(100),进一步包括控制装置(26),其以这样的方式控制X-射线照射区域使得X-射线仅照射到从通道方向或它的内部范围看去探测器(24)的中央部分的探测器通道宽度窄的一些范围、在探测器(24)的通道方向或它的内部范围上窄的数据采集范围、或在探测器(24)的通道方向上窄的一些数据采集范围,即,从通道方向或它的内部范围看去中央部分的探测器通道宽度窄的范围。
10.根据权利要求1到8任何一项的X-射线CT设备(100),进一步包括装置(23,31,32),其以这样的方式限制X-射线照射区域使得X-射线照射到探测器(24)的通道方向上的某些范围,该范围在从探测器(24)的通道方向或它的内部范围看去在中央部分通道上是细的,在探测器(24)的通道方向或它的内部范围上窄的数据采集范围,或在探测器(24)的通道方向上窄的数据采集范围,即,从通道方向或它的内部范围看去中央部分的细通道的范围。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005208235A JP5011482B2 (ja) | 2005-07-19 | 2005-07-19 | X線ct装置 |
JP2005-208235 | 2005-07-19 | ||
JP2005208235 | 2005-07-19 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1931098A true CN1931098A (zh) | 2007-03-21 |
CN1931098B CN1931098B (zh) | 2010-09-29 |
Family
ID=37307291
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2006101263673A Expired - Fee Related CN1931098B (zh) | 2005-07-19 | 2006-07-19 | X-射线ct设备 |
Country Status (5)
Country | Link |
---|---|
US (1) | US7522695B2 (zh) |
EP (1) | EP1745744A3 (zh) |
JP (1) | JP5011482B2 (zh) |
KR (1) | KR20070011176A (zh) |
CN (1) | CN1931098B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102971644A (zh) * | 2010-09-14 | 2013-03-13 | 株式会社东芝 | Mo准直器及使用该Mo准直器的X射线检测器、X射线检查装置以及CT装置 |
CN104856714A (zh) * | 2014-02-21 | 2015-08-26 | 上海西门子医疗器械有限公司 | Ct扫描参数指示方法、装置及ct机 |
Families Citing this family (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8804899B2 (en) | 2003-04-25 | 2014-08-12 | Rapiscan Systems, Inc. | Imaging, data acquisition, data transmission, and data distribution methods and systems for high data rate tomographic X-ray scanners |
WO2006129282A2 (en) * | 2005-05-31 | 2006-12-07 | Arineta Ltd. | Graded resolution field of view ct scanner |
JP5179007B2 (ja) * | 2005-12-01 | 2013-04-10 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置およびそのx線ct画像再構成方法 |
CN1989908B (zh) * | 2005-12-28 | 2010-05-19 | 株式会社东芝 | X射线ct装置及其控制方法 |
DE102007020879A1 (de) * | 2006-05-10 | 2009-04-02 | Gachon University Of Medicine & Science Industry-Academic Cooperation Foundation | Verfahren und Vorrichtung für die äußerst schnelle Symmetrie- und SIMD- gestützte Projektion/Rückprojektion für die 3D-PET-Bildrekonstruktion |
JP5383005B2 (ja) * | 2007-05-08 | 2014-01-08 | キヤノン株式会社 | X線ct撮影装置 |
JP2009082250A (ja) * | 2007-09-28 | 2009-04-23 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
US8553959B2 (en) * | 2008-03-21 | 2013-10-08 | General Electric Company | Method and apparatus for correcting multi-modality imaging data |
DE102008021639A1 (de) * | 2008-04-30 | 2009-12-03 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Vorrichtung und Verfahren zum Erzeugen einer CT-Rekonstruktion eines Objekts mit einem hochaufgelösten interessierenden Objektbereich |
US8300765B2 (en) * | 2008-08-04 | 2012-10-30 | Koninklijke Philips Electronics N.V. | Interventional imaging and data processing |
JP5159543B2 (ja) * | 2008-09-30 | 2013-03-06 | 株式会社東芝 | X線ct装置 |
JP5461106B2 (ja) * | 2009-08-25 | 2014-04-02 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
US20110211667A1 (en) * | 2010-02-26 | 2011-09-01 | Abdelaziz Ikhlef | De-populated detector for computed tomography and method of making same |
JP5725981B2 (ja) * | 2010-06-16 | 2015-05-27 | 株式会社東芝 | 医用画像表示装置及びx線コンピュータ断層撮影装置 |
JP5657614B2 (ja) | 2011-08-26 | 2015-01-21 | 富士フイルム株式会社 | 放射線検出器および放射線画像撮影装置 |
US8942341B2 (en) * | 2011-09-01 | 2015-01-27 | General Electric Company | Method of dose reduction for CT imaging and apparatus for implementing same |
JP6109560B2 (ja) * | 2012-12-21 | 2017-04-05 | 東芝メディカルシステムズ株式会社 | X線コンピュータ断層撮影装置 |
US9076563B2 (en) * | 2013-06-03 | 2015-07-07 | Zhengrong Ying | Anti-scatter collimators for detector systems of multi-slice X-ray computed tomography systems |
JP5677543B2 (ja) * | 2013-10-03 | 2015-02-25 | キヤノン株式会社 | X線ct撮影装置、x線ctの撮影方法、x線センサ、及び再構成装置 |
JP6342175B2 (ja) * | 2014-02-10 | 2018-06-13 | キヤノンメディカルシステムズ株式会社 | X線ct装置 |
US10126254B2 (en) * | 2014-12-18 | 2018-11-13 | Toshiba Medical Systems Corporation | Non-uniform photon-counting detector array on a fourth-generation ring to achieve uniform noise and spectral performance in Z-direction |
US10695011B2 (en) * | 2018-08-08 | 2020-06-30 | General Electric Company | X-ray collimator for imaging system |
CN111652951B (zh) * | 2020-05-07 | 2023-06-06 | 中国工程物理研究院材料研究所 | 一种稀疏角度快中子ct成像方法 |
JP7466401B2 (ja) | 2020-07-31 | 2024-04-12 | キヤノンメディカルシステムズ株式会社 | 医用画像診断装置 |
JP7141441B2 (ja) * | 2020-12-18 | 2022-09-22 | キヤノンメディカルシステムズ株式会社 | X線コンピュータ断層撮影装置 |
Family Cites Families (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SE347859B (zh) * | 1970-11-30 | 1972-08-14 | Medinova Ab | |
GB1577172A (en) * | 1976-07-15 | 1980-10-22 | Tokyo Shibaura Electric Co | Tomographing device |
DE3112891A1 (de) * | 1980-04-02 | 1982-01-28 | General Electric Co., New York, N.Y. | "rechnergesteuertes tomographiegeraet und -verfahren" |
JPS59181136A (ja) * | 1983-03-17 | 1984-10-15 | シーメンス・アクチェンゲゼルシャフト | 断層撮影システム |
US4638499A (en) * | 1984-08-06 | 1987-01-20 | General Electric Company | High resolution collimator system for X-ray detector |
US5033075A (en) * | 1988-05-18 | 1991-07-16 | Rad/Red Laboratories Inc. | Radiation reduction filter for use in medical diagnosis |
FR2654917A1 (fr) * | 1989-11-24 | 1991-05-31 | Gen Electric Cgr | Appareil de radiologie avec filtre d'homogeneisation. |
US5054041A (en) * | 1990-03-19 | 1991-10-01 | General Electric Company | High precision x-ray collimator |
JPH04179100A (ja) * | 1990-11-13 | 1992-06-25 | Toshiba Corp | X線絞り装置 |
JP3468372B2 (ja) * | 1992-09-07 | 2003-11-17 | 株式会社日立メディコ | 定位的放射線治療装置 |
JP3763611B2 (ja) * | 1996-07-12 | 2006-04-05 | 株式会社東芝 | X線ctスキャナ |
IL119033A0 (en) * | 1996-08-07 | 1996-11-14 | Elscint Ltd | Multi-slice detector array |
JP3828967B2 (ja) * | 1996-10-30 | 2006-10-04 | 株式会社東芝 | X線ctスキャナ |
US6023494A (en) | 1996-12-19 | 2000-02-08 | General Electric Company | Methods and apparatus for modifying slice thickness during a helical scan |
US5974109A (en) * | 1997-11-07 | 1999-10-26 | General Electric Company | Methods and apparatus for cell ganging in a multislice computed tomography system |
US6115448A (en) | 1997-11-26 | 2000-09-05 | General Electric Company | Photodiode array for a scalable multislice scanning computed tomography system |
WO1999030616A1 (en) * | 1997-12-16 | 1999-06-24 | Koninklijke Philips Electronics N.V. | Computer tomography device |
US5982846A (en) * | 1998-04-13 | 1999-11-09 | General Electric Company | Methods and apparatus for dose reduction in a computed tomograph |
US6275562B1 (en) | 1998-04-28 | 2001-08-14 | General Electric Company | Apparatus and methods for performing scalable multislice computed tomography scan |
US6198791B1 (en) | 1998-08-25 | 2001-03-06 | General Electric Company | Scalable multislice imaging system |
JP3977558B2 (ja) | 1998-10-20 | 2007-09-19 | 株式会社東芝 | 2次元アレイ型x線検出器、x線検出器の製造方法及び2次元アレイ型x線検出器を用いたctスキャナ装置 |
DE69936769T2 (de) | 1998-12-30 | 2008-04-30 | General Electric Co. | Bilddickeselektion für mehrschichtbildgerät |
DE19925395B4 (de) * | 1999-06-02 | 2004-11-25 | Siemens Ag | Verfahren zum Betrieb eines Computertomographie(CT)-Gerätes |
JP2001120534A (ja) * | 1999-10-28 | 2001-05-08 | Hitachi Medical Corp | マルチスライス型x線ct装置 |
JP2001187046A (ja) | 1999-12-27 | 2001-07-10 | Ge Medical Systems Global Technology Co Llc | マルチスライスx線ct装置及びその制御方法 |
DE10015191A1 (de) | 2000-03-27 | 2001-10-25 | Siemens Ag | Detektor für Röntgen-Computertomograph |
US6385278B1 (en) * | 2000-04-28 | 2002-05-07 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for region of interest multislice CT scan |
US6700948B2 (en) * | 2000-12-12 | 2004-03-02 | Ge Medical Systems Global Technology Company, Llc | Low-cost, multislice CT detector with multiple operating modes |
JP3942142B2 (ja) * | 2000-12-15 | 2007-07-11 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 放射線断層撮影装置およびその方法 |
JP4654551B2 (ja) * | 2001-08-27 | 2011-03-23 | 株式会社島津製作所 | Ct装置 |
US6895077B2 (en) * | 2001-11-21 | 2005-05-17 | University Of Massachusetts Medical Center | System and method for x-ray fluoroscopic imaging |
JP4041025B2 (ja) | 2003-07-15 | 2008-01-30 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線分布調整フィルタ装置およびそれを用いたx線ct装置 |
JP3942178B2 (ja) * | 2003-07-29 | 2007-07-11 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ctシステム |
JP4041040B2 (ja) * | 2003-09-08 | 2008-01-30 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 放射線断層撮影装置 |
JP4342274B2 (ja) * | 2003-11-04 | 2009-10-14 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
-
2005
- 2005-07-19 JP JP2005208235A patent/JP5011482B2/ja not_active Expired - Fee Related
-
2006
- 2006-07-19 US US11/489,065 patent/US7522695B2/en active Active
- 2006-07-19 EP EP20060253787 patent/EP1745744A3/en not_active Ceased
- 2006-07-19 KR KR1020060067632A patent/KR20070011176A/ko not_active Application Discontinuation
- 2006-07-19 CN CN2006101263673A patent/CN1931098B/zh not_active Expired - Fee Related
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102971644A (zh) * | 2010-09-14 | 2013-03-13 | 株式会社东芝 | Mo准直器及使用该Mo准直器的X射线检测器、X射线检查装置以及CT装置 |
CN104856714A (zh) * | 2014-02-21 | 2015-08-26 | 上海西门子医疗器械有限公司 | Ct扫描参数指示方法、装置及ct机 |
Also Published As
Publication number | Publication date |
---|---|
US20070019779A1 (en) | 2007-01-25 |
JP5011482B2 (ja) | 2012-08-29 |
US7522695B2 (en) | 2009-04-21 |
CN1931098B (zh) | 2010-09-29 |
EP1745744A3 (en) | 2008-01-30 |
EP1745744A2 (en) | 2007-01-24 |
JP2007020906A (ja) | 2007-02-01 |
KR20070011176A (ko) | 2007-01-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1931098A (zh) | X-射线ct设备 | |
CN1968654A (zh) | 放射线断层像摄影装置 | |
CN1945630A (zh) | 图像重构方法和计算x线断层成像设备 | |
CN1917811A (zh) | 断层摄影图象的再构成方法及断层摄影装置 | |
CN101032408A (zh) | X线ct设备 | |
CN1955725A (zh) | X射线ct系统 | |
CN101049243A (zh) | X射线ct装置 | |
CN1324327C (zh) | 用于获得与环境内产生的辐射有关的信息的方法 | |
CN1249987C (zh) | 辐射图像数据处理设备和辐射图像数据处理方法 | |
CN1305011C (zh) | 图像处理方法和图像处理设备 | |
CN1965762A (zh) | Ct系统束硬化后处理方法及ct系统 | |
CN1688254A (zh) | X线断层摄影装置 | |
CN1268035A (zh) | 利用螺旋扫描ct再现立体图像 | |
CN1309548A (zh) | 用于螺旋扫描立体计算的x断层照相系统的改进的探测器阵列的几何形状 | |
CN1307943C (zh) | 逆投影方法和x射线计算机化断层摄影装置 | |
CN1493258A (zh) | 图像处理设备与超声波诊断设备 | |
CN1400806A (zh) | 采用自适应二值化的图象处理方法和设备 | |
CN1121608A (zh) | 位置检测装置及其所用的位置指示器 | |
CN1658241A (zh) | 信号处理方法与设备以及记录介质 | |
CN1460452A (zh) | 三维反投影方法和一种x射线计算机层析成像装置 | |
CN1721815A (zh) | 光学对象识别装置 | |
CN101065782A (zh) | 在网格表面和立体对象上产生和测量表面线的系统和方法以及网格分割技术(“曲线测量”) | |
CN1115873C (zh) | 高清晰度电视系统 | |
CN1254102C (zh) | 信息及图像信号处理装置、系数种类数据生成装置及方法 | |
CN1533551A (zh) | 信号处理设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100929 Termination date: 20210719 |