CN1056506C - 用于双束计算机ｘ射线断层扫描术的稳定的可间断滤波器 - Google Patents
用于双束计算机ｘ射线断层扫描术的稳定的可间断滤波器 Download PDFInfo
- Publication number
- CN1056506C CN1056506C CN91105322A CN91105322A CN1056506C CN 1056506 C CN1056506 C CN 1056506C CN 91105322 A CN91105322 A CN 91105322A CN 91105322 A CN91105322 A CN 91105322A CN 1056506 C CN1056506 C CN 1056506C
- Authority
- CN
- China
- Prior art keywords
- filter
- wave filter
- signal
- interrupted
- output
- 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.)
- Expired - Fee Related
Links
- 230000009977 dual effect Effects 0.000 title abstract 2
- 238000002591 computed tomography Methods 0.000 title description 16
- 238000001914 filtration Methods 0.000 claims abstract description 6
- 238000005070 sampling Methods 0.000 claims description 26
- 239000000523 sample Substances 0.000 claims description 21
- 238000013170 computed tomography imaging Methods 0.000 claims description 11
- 230000003595 spectral effect Effects 0.000 abstract description 2
- 230000010354 integration Effects 0.000 description 19
- 230000008859 change Effects 0.000 description 12
- 238000000034 method Methods 0.000 description 12
- 238000003384 imaging method Methods 0.000 description 11
- 238000005516 engineering process Methods 0.000 description 9
- 239000003990 capacitor Substances 0.000 description 8
- 238000006243 chemical reaction Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 238000002156 mixing Methods 0.000 description 3
- 230000001172 regenerating effect Effects 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000010009 beating Methods 0.000 description 2
- 238000007600 charging Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000012071 phase Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000008929 regeneration Effects 0.000 description 2
- 238000011069 regeneration method Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 241000446313 Lamella Species 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000012937 correction Methods 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
- 238000001514 detection method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000008384 inner phase Substances 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000009418 renovation Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000013519 translation Methods 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/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/405—Source units specially adapted to modify characteristics of the beam during the data acquisition process
-
- 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/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- 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/4021—Arrangements for generating radiation specially adapted for radiation diagnosis involving movement of the focal spot
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
一种用于双束CT系统的数据采集滤波器用两个频谱滤波器根据X射线束的状态交替接收来自每一探测器的信号。当一个滤波器对探测信号滤波时，另一个滤波器处于“保持”状态，其输出和内部量值都被冻结。每一滤波器只对一种波束状态下产生的信号进行有效的滤波，而不受在另一波束状态下产生的信号的影响，也不受在另一波束状态下经过的时间的影响。处于保持状态的滤波器的输出是一常数，可以在该周期的任何时刻对它取样。
Description
本发明涉及计算机X射线断层扫描(CT)系统，特别涉及用于对CT系统采集的数据进行处理的抗混杂滤波器。
在计算机X射线断层扫描系统中，X射线源经过校正形成具有一定的扇形束角的扇形束。扇形束指向被称作“成象平面”的直角坐标系统的X-Y平面内，穿过被成象的物体射向排列在成象平面内的X射线探测器阵列。
探测器阵列是由许多探测单元组成的，单元之间的间距近似等于单元的宽度。每一探测单元测量由X射线源投射至特定接收单元的射线束的透射辐射强度。透射辐射强度取决于X射线束沿途因被成象物体阻挡而造成的衰减。
X射线源和探测器阵列可以在成象平面内围绕被成象物体在机架上旋转，使得扇形束与被成象物体的交角不断变化。当机架旋转时，就得到了构成一个投影群的许多投影，每个投影由若干探测单元的强度信号组成，这些强度信号是在探测单元扫过以投影角为中心的很小的机架旋转角时接收到的。
所得到的断层扫描投影群典型地以数值形式储存起来，以便进行计算机处理，根据熟知的重建算法，“重建”出片层图象。利用扇形束重建方法，可以将扇形束投影的投影群直接重建成图象，或者将投影的强度数据储存在平行束中，并根据平行束再生方法加以再生。再生的断层造影图象可以在传统的阴极射线管CRT上显示，也可利用计算机控制的摄像机转换到胶片载体上。
机架的连续旋转在每一探测器上产生一个恒定变化的信号，这些变化的信号与同这些探测器相关的X射线束在其束角改变时产生的衰减变化相对应。将这一信号相对和每一投影角相关的机架旋转增量求积分，即得到该投影角的探测信号。然后将达一积分值加以保持，利用数据采集系统(DAS)进行取样，并将其转换成数字化的探测值，以便于计算机的储存和再生。
探测信号的积分，使探测器灵敏度增加，同时也给探测器信号提供固有的带宽限制，从而阻止了探测器信号取样时由DAS引起的“混杂”。众所周知，混杂信号是信号赝象，它是由存在于取样信号中的具有高于取样速率之半频率的频率分量产生的。
为了提供合适的时间供DAS对探测器信号取样，每一探测器单元中需用两个这样的积分器。当其中一个积分器从探测单元积分新的现行数据时，另一积分器正保持着先前积分数据的值以供取样。这种两积分器系统被称作“积分和转储电路”。它具有如下的优点：即既提供一个精确限定的积分时间，又为DAS提供一个足够大的取样窗口。
尽管如此，积分和转储电路对于由积分电路中积分电容数值的改变而产生的增益的变化是十分敏感的。更有甚者，专门用于联接两个积分器的分立元件转接开关漏电十分严重，这就需要首先对探测信号进行前置放大。这一独立前置放大器增益的变化会对积分和转储电路遭受的增益变化施加影响。增益变化会在断层扫描的再生图像中产生不能接受的拖影、0环、斑点或其他赝象。
像低通滤波器这样的连续波滤波器可以用来代替积分器。在连续波滤波器中，滤波器的输出将以加权滚动的方式反映出先前的探测信号。采用具有适当频率截止点的低通滤波器，能够消除混杂。
在连续波滤波器工作状态下，只要对机架的位置作适当的修正，取样可以在任何时间进行，因此一个探测器只需要一个单个滤波器。这就无需采用与积分和转储电路的分立元件转接开关相联的独立前置放大器。
采用连续波滤波器改善了积分和转储电路的增益稳定性。(与积分电路不同)在连续波滤波器中可以建立起直流反馈。因此，滤波器的增益可以由一个单个电阻来确定，而不是由积分和转储电路的电容及与之相联的前置放大器的电阻两者来确定。一般说来，电阻比电容更稳定，单个元件比多个元件更稳定。
在一些CT成象技术中，当机架旋转时，X射线束的两个束强度或两个束位置间迅速变换。例如，在双能量扫描技术中，加到X射线管上的电压可以改变，以获得两束具有不同频谱的X射线，从而生成两种图象，将两种图象加以比较，有助于区分各种不同的组织。另外，在“点跳动”扫描技术中，X射线发射点可以相对于机架“跳动”，以产生具有稍许不同角度的两束射线，从而增加X射线成象的分辨率。这后一种技术在1990年6月20日提交的美国申请序号为07/540,995的专利“具有可变换焦点的计算机断层扫描系统”中有详尽的描述，上述专利与本发明具有相同的受让人，这里列举出来以作参考。
在这些双束技术中，随着机架的旋转，束的状态迅速变换，以减少由于病人移动对采集数据一致性的影响。假如束的状态只是在机架转动一周之后才变换的话，那么病人的移动通常会在相应的探测信号之间引起更多的变动。
随着X射线束迅速变换。来自探测器的信号也将变化，必须相对于束的变换将变化的信号同步分开，才能采集到分别与每种束的状态相关的两组不同的数据。将探测器信号分离成与每一束的状态相关的一种方法是采用前面所述的积分和转储电路的双积分器。第一积分器调节到只有在第一束状态下才积分(随即采集数据)，第二积分器调节到只有在第二束状态下才积分。两个积分器的输出产生了两组不同的分别与各个束状态相关的探测器数据。
遗憾的是，这一方法仍存在如上所说的积分和转储电路与增益灵敏度相关的缺点，亦即电路的增益是由一个与每一积分电路相关的相对不稳定的电容和一个与独立的前置放大器相关的电阻确定的。
前面所述的连续波滤波器不适用于这些双束技术，因为滤波器产生一个连续输出，该输出是先前探测器信号的函数，而与束状态无关。连续波滤波器的连续属性过去被采用是基于它们能够灵活采样的优点，但却阻碍了双束信号的彻底分离。
本发明提供一种探测器采集滤波器，它与双束CT技术结合使用，既具有连续波滤波器的稳定性，又具有积分和转储滤波器同步取样的能力。尤其是数据采集滤波器含有第一和第二可间断滤波器，它们或者工作于滤波状态，或者工作于保持状态。当滤波器处于滤波状态时，它们的输出取决于仅在当前和上一个滤波状态下探测信号的频率分量。当处于保持状态时，滤波器的输出是上一个滤波状态期间滤波器的最后输出值。当CT机的X射线束处于第一状态时，第一可间断滤波器处于滤波状态，第二可间断滤波器处于保持状态。反之，当X射线束处于第二状态时，情形正好相反。在一种实施例中，转接开关提供一个单一取样输出，它被交替联接至处于保持状态的滤波器上。
因此，本发明的一个目的就是提供一种装置，它能对由两个迅速交谈的X射线束状态产生的间断探测器信号进行滤波。
优选的是，两个可间断滤波器的每一个共用一个前置放大器以接收探测器信号，并含有一个可转接的反馈电路单元，用于接收来自处于滤波状态的可间断滤波器的反馈信号。
因此，本发明的另一个目的就是提供一种增益稳定的抗混杂的滤波器。单独享用的反馈电路单元为第一和第二可间断滤波器提供精确的增益匹配，可间断滤波器可以利用单一的电阻反馈元件加以控制，使其比诸如在积分和转储滤波器电路中采用的电容和电阻反馈元件具有更高的稳定性。
本发明的再一个目的是提供一种采用了上述抗混杂滤波器的Cr成象系统。
根据本发明的一个方面，一种用于双束CT系统的数据采集滤波器，所述的双束CT系统具有在第一和第二周期被激活的第一和第二X射线束，该数据采集滤波器接收探测器信号并产生取样输出，包括一个与第一和第二周期同步用以产生时钟信号的时钟；其特征在于包括：
用于接收时钟信号(90、92)的第一和第二可间断滤波器(19)，所述可间断滤波器具有产生取样输出(82)的滤波器输出端，每个滤波器具有滤波状态和保持状态，在滤波状态下每一滤波器的滤波输出只取决于当前和前一个滤波状态下信号(62)的频率分量，而在保持状态下，滤波输出取决于前一个滤波状态的最后输出值；
其中第一可间断滤波器在第一周期处于滤波状态，在第二周期处于保持状态，而第二可间断滤波器在第一周期处于保持状态在第二周期处于滤波状态。
根据本发明的另一个方面，一种CT成象系统，具有在第一和第二周期激发的第一和第二X射线束，一个探测器(18)接收该第一和第二X射线束，用以产生探测器信号(62)，其特征在于，采用了一个上述数据采集滤波器。
本发明的上述目的以及其他目的和优点将会在随后的描述中看到。在描述中将会涉及一些附图，这些附图是本发明的组成部分，并以图示方式给出了本发明的最佳实施例。这样的实施例不一定代表本发明的全部内容，但是必定涉及解释本发明范畴的权利要求书。
图1是适用于本发明的CT系统示意图；
图2是由图1系统产生的X射线扇形束的详细图，图中示出了与之有关的相应角度和坐标轴；
图3是用图1CT系统采用焦点“跳动”方法获得的投影数据的射线角T与半径r的曲线图，并示出了与两束状态有关的周期
图4(a)是运用图3的跳动技术由探测单元产生的探测器信号的代表性曲线，图中示出了两种束状态数据的交织；
图4(b)是类似于图4(a)的曲线，它表示出由本发明的一个可间断滤波器接收图4(a)的信号时看到的有效输入；
图4(c)是本发明的可间断滤波器接收图4(a)的输入信号时的代表性输出曲线；
图5是本发明的可间断滤波器的简化框图；
图6是用于使图5的可间断滤波器与CT机双束状态同步的时钟波形图；
图7是本发明的可间断滤波器的详细方框图。
参看图1，CT机架16是″第三代″CT扫描装置的代表，它包括一个X射线源10，X射线源能够从焦点11投射出一个扇形X射线束24，经过被成象物体12到达探测器阵列18上。探测器阵列18由许多探测单元26组成，这些探测单元共同探测由X射线穿过被成象物体12所得到的投影图象。与每一探测器相关联的是可间断滤波器19，下面将对它作更全面的介绍。机架16围绕位于被成象物体12内的旋转中心14旋转。
适用于本发明的CT扫描装置的控制系统具有同机架相关联的控制模块28，控制模块28包括X射线控制器30、机架马达控制器32以及数据采集系统(″DAS″)34三部分。X射线控制器30为X射线源10提供电源和定时信号，并控制X射线管焦点11的位置。机架马达控制器32控制机架16的转速和位置。数据采集系统34经过可间断滤波器19，在探测器阵列18的探测器单元26上对投影数据取样。并将已经滤了波的数据转换成为以后计算机处理所需的数字字。
X射线控制器30和机架马达控制器32被联接到计算机36上。计算机36是一台通用小型计算机，如Data General Eclipse MV/7800C。
DAS34与图象再生装置38相联接。图象再生装置38接收来自DAS34的取样的和数字化的信号，并根据熟知的方法进行高速图象再生处理。图象再生装置38可以是一个阵列处理机，如弗吉尼亚明星技术公司(Star Technologies of Virginia)制造的阵列处理机。
计算机36通过操作控制台40接收指令和扫描参数。操作控制台40通常为一阴极射线管(CRT)显示屏和键盘，操作人员可以输入扫描参数，显示再生图象和其他来自计算机36的其它信息。大容量存储装置42用于存储CT成象系统的操作程序以及图象数据，以供操作人员备用。
参看图2，与特定探测单元26相关联的扇形束24的那一部分可以用通过X射线焦点11的中心和特定探测单元26的中心的射线20加以标记。射线20用一条径线来加以描述，该径线距旋转中心14的垂直距离长为″r″，与固定在空间的任意参考轴22构成的旋转角为T。
如图3所示每一射线20的r和T可以绘制成r-T关系图。其中T为水平坐标轴，r为垂直坐标轴。参看图3，对于某一投影，在数据采集的开始，给定探测器26的射线20处于r-T关系图上的位置50。为了简单起见，图3只示出了一条射线20的起点位置50。但是，在现有技术中应当知道，一个投影通常包含有近千条射线20和相应的强度测量数据。
当机架16旋转时，在第一个周期I内，射线20的位置由起点50沿机架旋转角dT的水平线52水平地在r-T关系图上移动。水平线52对应于由于机架16的旋转而引起的T的增加。在水平线52上沿射线20方向X射线辐射的变化强度由探测器单元26接收。在机架16转过dT角后，X射线管10的焦点11相对于机架16的位置发生偏移，将射线20的位置移至起点位置54，此时r增加了，T减少了。在机架旋转的第二个周期II内，射线20由该起点位置54沿着水平线56再次水平移动距离dT。在这一起始位置末尾，焦点11相对于机架移回到原来的位置，由于机架交错旋转，射线20在r-T关系图上处于新的位置58，它与起点位置50具有相同的r值，但是具有位移T。随着X射线焦点11在第一周期I的第一状态和第二周期II的第二状态之间运动，这一过程重复进行。这两种状态由X射线焦点11相对于机架16的相对位置所确定。
参看图4(a)，从探测器26取出的信号62的振幅在时间t＝nΔt时发生突变，其中n为整数，Δt为机架转动dT所需的时间。这些不连续点60的产生是由于在两种状态之间X射线焦点11不连续运动的结果。亦即在第一周期I和第二周期II之间，X射线束24相对于被成象物体12的方位发生变化的结果。如上所述，当连续波滤波器接收到这一不连续信号时，它将试图把第1周期和第II周期的数据混合起来，以期达到使第I周期数据与第II周期数据“模糊”的效果。
现在参看图5，本发明的可间断滤波器19接收在探测器输入端64的探测器信号62，探测器信号62来自探测器阵列18中的探测单元26。还有另外一些独立的可间断滤波器分别与其余各个探测器单元26相关(未示出)。探测器信号62由差分前置放大器66放大，该差分前置放大器66在其反相输入端接收探测器信号62，并将探测器信号62的电平提高以减少联接至前置放大器66输出端的分立元件转接开关68的漏电电流和注入电流的影响。前置放大器66可以是任何现成的具有低输入电流、低偏压和低偏压漂移的商品运算放大器。前置放大器66的非反相输入端接地。
电阻72为该前置放大器66提供局部反馈，以提高其稳定性，它将前置放大器66输出端联接至其反向输入端，这是具有普通技术的人员都能理解的。
前置放大器66的输出信号被单刀单掷(“SPST”)分立元件转接开关68与70的“动片”端接收，转接开关68与70的“定片”端分别接至低通滤波器子单元74和76的输入端。
滤波器子单元74和76的输出端随后接至SPST转接开关84和86的“动片”端，而它们的“定片”端联在一起，并通过反馈元件88与前置放大器66的反相输入端相联接。因此，当转接开关84接通时，输向前置放大器66的反馈信号由滤波器子单元74的输入端提供，而当转接开关86接通时，输向前置放大器66的反馈信号由滤波器子单元76的输出端提供。
因此，应当理解，当转接开关68和84接通时，滤波器子单元74呈反馈状与前置放大器66联接，形成第一滤波器；而当转接开关70和86接通时，滤波器子单元76呈反馈状与前置放大器66联接，形成第二滤波器。
第一和第二滤波器是双极低通贝塞尔(Bessel)滤波器，其截止频率为1105赫兹。截止频率的选定与DAS34的取样频率有关，正如业已讨论过的那样，是为了使其具有“抗混杂”的特性。尽管双极贝塞尔滤波器方案是一种最佳方案，但是从以下的介绍中本技术领域普通人都明白，像白托沃斯(Butterworth)或契比雪夫(Tcheby-cheff)等其他低通滤波器也可以采用。
由滤波器子单元74和76构成的第一和第二滤波器也是“可间断”的，当相应的转接开关68或70接通时，它们具有低通滤波器功能，但是，当相应的转接开关68或70断开时，它们却“冻结”了。在其输出端保特着转接开关68或70断开以前所具有的最后的滤波值，并维持其内部数值不随时间而变。
滤波器子单元74和76的输出为SPST分立元件转接开关78和80的“动片”端所接收，它们的“定片”端联在一起，向DAS34提供单一的取样输出信号82，这一信号随周期的变化或者由滤波器子单元74提供，或者由滤波器子单元76提供。转接开关78和80分别相对于转接开关68和70取相反的状态。当转接开关68接通时，转接开关78断开，滤波器子单元74呈滤波状态，此时为周期I；当转接开关68断开时，转接开关78接通，滤波器子单元74呈保持状态，此时为周期II。因此，滤波器子单元74输出端的取样过程可以发生在周期II的任何时刻，并还将给出严格对应于周期I和II之间最后过渡时刻的滤波值。这使得DAS34的取样动作与机架16特定时间或位置的同步比与上面所述的连接波滤波器更容易。这种取样与机架16位置的协调对消除上述利用焦点跳动技术形成图像中的赝象可能是有用的。
另一种方式是当转接开关80允许滤波器子单元76输出端在周期I的任何时间取样时，转接开关70取相反状态。
参看图5和图6，转接开关68，84和80由时钟信号90同步接通，时钟信号90在第一周期I内定时为“通”或高电位，在第二周期II内，定时为“断”。滤波器74在第一周期I内相应为滤波状态，在第二周期II内为保持状态。
相反，转接开关70、86和78由时钟信号92接通，时钟信号92在周期II内定时为通，在周期I内定时为断，滤波器76在周期I内为保持状态，在周期II内为滤波状态。信号90和92是矩形波，其相位设定将保证信号90和92的接通时间不发生重叠，因而在任何时刻，只有信号90或92处于接通状态。
分立元件转接开关的作用是将滤波器子单元74或76与前置放大器66接通，并且通过反馈元件88提供一适当的反馈回路。滤波器子单元74或76具有高直流增益，因此，滤波器19的直流增益将由反馈元件88所确定，而与滤波器子单元74或76哪一个被接通无关，也与通过前置放大器66的正向增益和通过滤波器子单元74或76的正向增益无关。反馈元件88基本上是电阻性的，主要由一个单独电阻所决定，因此，是一个比前述的积分和转储电路的电容器稳定得多的增益控制元件。
滤波器子单元74和76的电路是相同的，为简化起见，将只对滤波器子单元74作详细介绍，通过这一介绍，将会明了滤波器子单元76的工作过程。
参看图7，输入至滤波器子单元74的信号，被接至由输入电阻94形成的积分电路上，输入电阻94又通过转接开关68接至差分放大器96反向输入端。放大器96有一电容器98联接在其输出端和反向输入端，众所周知，电容器98被用来积分流经电阻94的电流。差分放大器96的非反向输入端接地。
因此，当转接开关68接通时，放大器96的反向输入端经过电阻94接收来自前置放大器66的电源，并对该电流积分。但当转接开关68断开时，放大器96的输出端保持着转接开关68断开前所具有的最后的值。
放大器96的输出经电阻100和SPST分立元件转接开关102接至放大器104的反向输入端，其中转接开关102的“动片”端接至电阻100，“定片”端接至放大器104的反向输入端。电容器106接在放大器104的输出端和其反向输入端之间，反馈电阻108接在放大器104的输出端和转接开关102的“动片”端。因此当转接开关102断开时，放大器104的电阻反馈线路是断开的，而放大器104的作用如同放大器96，处于保持状态，保持着转接开关102断开前它的最后输出值。
由于分立元件转接开关102的开关作用而产生的对放大器104的电荷注入将会因为SPST固态转接开关110和112以及114至122诸元件的作用而减至最小。转接开关110的“动片”端和“定片”端皆接至转接开关102的“动片”端上，当转接开关102断开以抵制由其本身而产生的充电时，转接开关110接通。放大器104的非反相端也经过电容器114接地，电容器114被转接开关112旁路，转接开关112与转接开关102同步断开和接通。
放大器104的输出经由一网络接至转接开关84的“动片”端，该网络由一串联电阻116和与之相联的接地电容118组成。放大器104的输出还通过一个同样的网络与转接开关78相联接。这些网络也缓解了由分立元件转解开关84和78产生的充电影响。
反馈元件88由引自转接开关84“定片”端的串联电阻124、与之相联的接地电容126和自电容126与电阻124的结点引向前置放大器66的反相输入端的串联电阻组成。可以知道，反馈电路不仅包括电阻116、124和128，而且也包括电容118和126。这些电容对直流增益无影响，对低频增益的影响也可忽略，因为电容的值和电阻116和124的值都很小，它们大小的选定是使得RC时间常数近似等于0.05微秒。
由前置放大器66、反馈回路中滤波器子单元74或76同反馈元件88组成的贝塞尔滤波器的截止频率的选定原则是防止混杂信号进入将被DAS34取样的信号之中。因此，贝塞尔滤波器的截止频率可以通过改变取样速率加以调节，取样速率一般取决于机架16的速度和探测器单元26的数目和尺寸。
需要注意的是：当信号90很小时，转接开关68和102的作用实际上“冻结”了由前置放大器66和滤波器子单元74构成的滤波器的操作。当信号90很小时，滤波器的中间值，特别是放大器96和104的输出，将不随时间或探测器输入的变化而改变，以电容器98和106的电压表示的滤波器固有的“存贮记忆”也不会改变。这种采用滤波器子单元74或者采用滤波器子单元76的滤波器“保特”状态允许把简单的低通滤波器用于图4(a)所示的不连续探测器信号62的情形。
参看图4(b)，输入到前置放大器66和滤波器子单元74构成的滤波器的有效输入仅仅是那些发生在周期I期间探测器信号62的各分段的联接，这是因为在周期II期间滤波器冻结在保持状态。相反，输入到前置放大器66和滤波器子单元76(未示出)构成的滤波器的有效输入是那些发生在周期II期间探测器信号62的各分段的联接。
参看图4(C)，滤波器19的输出表现为一系列恒定的电压，这些恒定的电压反映出在上一个滤波状态期间每一电流保持状态时滤波器子单元74和76的最后输出。滤波器子单元74和76的相对长的保持状态使得探测器阵列18中所有探测器单元26的取样过程能同步在同样的机架位置，对于积分和转储滤波器，第一和第二滤波器的结构，允许为提高稳定性而引入的阻性反馈。
本最佳实施例的许多改进和改变方案仍然属于本发明的精神和范畴之中，这对那些具有这一领域常识的人们来说是险而易见的。例如，这一滤波器也可以与双能量扫描装置配合使用，此时比如第一周期I是高能X射线周期，则第二周期II是低能X射线周期。另外显而易见的是，该滤波器也可与其他机架和探测器装置配合使用，其中包括所谓的“平动和转动”结构的装置。
Claims (12)
1.一种用于双束CT成象系统的数据采集滤波器，所述的双束CT成象系统具有在第一和第二周期被激活的第一和第二X射线束，该数据采集滤波器接收探测器信号并产生取样输出，包括一个与第一和第二周期同步用以产生时钟信号的时钟；其特征在于包括：
用于接收时钟信号(90、92)的第一和第二可间断滤波器(19)，所述可间断滤波器具有产生取样输出(82)的滤波器输出端，每个滤波器具有滤波状态和保持状态，在滤波状态下每一滤波器的滤波输出只取决于当前和前一个滤波状态下信号(62)的频率分量，而在保持状态下，滤波输出取决于前一个滤波状态的最后输出值；
其中第一可间断滤波器在第一周期处于滤波状态，在第二周期处于保持状态，而第二可间断滤波器在第一周期处于保持状态在第二周期处于滤波状态。
2.如权利要求1的数据采集滤波器，其特征在于包括在第一周期将第一滤波器输出端同取样输出端联接起来，而在第二周期将第二滤波器输出端同取样输出端联接起来的转接开关(68、70)。
3.如权利要求1的数据采集滤波器，其特征在于：第一和第二可间断滤波器包括并共用一个用于放大探测信号的前置放大器装置(66)和一个可转接的反馈元件(72、88)，该可转接的反馈元件用于在第一周期接收来自第一可间断滤波器输出端的反馈信号，在第二周期接收来自第二可间断滤波器输出端的反馈信号，并向前置放大器提供反馈。
4.如权利要求3的数据采集滤波器，其特征在于：反馈元件基本上为电阻性的。
5.如权利要求1至4任何一项的数据采集滤波器，其特征在于：可间断滤波器作为二级低通贝塞尔滤波器，其截止频率基本上等于滤波状态期间取样速率的一半。
6.一种CT成象系统，具有在第一和第二周期激发的第一和第二X射线束，一个探测器(18)接收该第一和第二X射线束，用以产生探测器信号(62)，其特征在于，采用了一个如权利要求1所述的数据采集滤波器。
7.如权利要求6的CT成象系统，其特征在于它包括在第一周期将第一滤波器输出端同取样输出端联接起来，而在第二周期将第二滤波器输出端同取样输出端联接起来的转接开关(68、70)。
8.如权利要求6的CT成象系统，其特征在于：第一和第二可间断滤波器包括并共用一个用于放大探测器信号的前置放大器装置(66)和可转接的反馈元件(72、88)，该可转接的反馈元件用于在第一周期接收来自第一可间断滤波器输出端的反馈信号，在第二周期接收来自第二可间断滤波器输出端的反馈信号，并为前置放大器提供反馈。
9.如权利要求8的CT成象系统，其特征在于：反馈元件基本上为电阻性的。
10.如权利要求6任何一项的CT成象系统，其特征在于：可间断滤波器作为二级低通贝塞尔滤波器，其截止频率基本上等于滤波状态下取样速率的一半。
11.如权利要求6的CT成象系统，其特征在于：第一和第二X射线束在空间相互偏置。
12.如权利要求6至10任何一项的CT成象系统，其特征在于：第一和第二X射线束具有不同的能量。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US562,766 | 1990-08-06 | ||
| US07/562,766 US5218533A (en) | 1990-08-06 | 1990-08-06 | Stable interruptible filter for dual beam computed tomography |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1058710A CN1058710A (zh) | 1992-02-19 |
| CN1056506C true CN1056506C (zh) | 2000-09-20 |
Family
ID=24247683
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN91105322A Expired - Fee Related CN1056506C (zh) | 1990-08-06 | 1991-08-03 | 用于双束计算机ｘ射线断层扫描术的稳定的可间断滤波器 |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US5218533A (zh) |
| EP (1) | EP0473296B1 (zh) |
| JP (1) | JPH0790023B2 (zh) |
| CN (1) | CN1056506C (zh) |
| CA (1) | CA2042226A1 (zh) |
| DE (1) | DE69125668T2 (zh) |
| IL (1) | IL98928A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1509686B (zh) * | 2002-11-27 | 2010-06-23 | Ge医药系统环球科技公司 | 定量组织脂肪含量的方法和装置 |
Families Citing this family (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5301108A (en) * | 1993-02-01 | 1994-04-05 | General Electric Company | Computed tomography system with z-axis correction |
| WO1996042022A1 (en) * | 1993-03-31 | 1996-12-27 | Eg & G Astrophysics Research Corporation | Three-dimensional reconstruction based on a limited number of x-ray projections |
| US5841829A (en) * | 1997-05-13 | 1998-11-24 | Analogic Corporation | Optimal channel filter for CT system with wobbling focal spot |
| JP4550179B2 (ja) * | 1999-02-12 | 2010-09-22 | 株式会社東芝 | Ｘ線コンピュータ断層撮影装置 |
| US6597352B1 (en) * | 1999-05-05 | 2003-07-22 | Nokia Corporation | Method and device for filtered sync detection |
| AU2003252103A1 (en) * | 2002-07-23 | 2004-02-09 | Ge Medical Systems Global Technology Company, Llc | Methods and systems for detecting components of plaque |
| US7187748B2 (en) * | 2003-12-30 | 2007-03-06 | Ge Medical Systems Global Technology Company, Llc | Multidetector CT imaging method and apparatus with reducing radiation scattering |
| US7453987B1 (en) * | 2004-03-04 | 2008-11-18 | Science Applications International Corporation | Method and system for high energy, low radiation power X-ray imaging of the contents of a target |
| WO2006047718A2 (en) * | 2004-10-22 | 2006-05-04 | Scantech Holdings, Llc | Angled-beam detection system for container inspection |
| JP4294610B2 (ja) * | 2005-04-14 | 2009-07-15 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 放射線ｃｔ装置、データ収集装置及びデータ収集方法 |
| JP4299806B2 (ja) * | 2005-05-11 | 2009-07-22 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ｘ線ｃｔ装置 |
| RU2414724C2 (ru) * | 2005-10-28 | 2011-03-20 | Конинклейке Филипс Электроникс, Н.В. | Способ и устройство для спектральной компьютерной томографии |
| US7283609B2 (en) * | 2005-11-10 | 2007-10-16 | General Electric Company | CT detector photodiode having multiple charge storage devices |
| US7789838B2 (en) * | 2006-06-02 | 2010-09-07 | Massachusetts Eye & Ear Infirmary | Audible range oculocometry for assessment of vestibular function |
| RU2465826C2 (ru) * | 2006-12-15 | 2012-11-10 | Конинклейке Филипс Электроникс Н.В. | Устройство формирования рентгеновских изображений со спектральным разрешением |
| FR2947658B1 (fr) * | 2009-07-03 | 2011-07-29 | Norbert Beyrard | Procede et dispositif d'imagerie a rayons x a traitement tridimentionnel |
| US8314394B1 (en) | 2009-11-04 | 2012-11-20 | Science Applications International Corporation | System and method for three-dimensional imaging using scattering from annihilation coincidence photons |
| DE102011075577B4 (de) * | 2011-05-10 | 2013-01-31 | MULTIPOND Wägetechnik GmbH | Signalverarbeitungsverfahren, Vorrichtung zur Signalverarbeitung und Waage mit Vorrichtung zur Signalverarbeitung |
| CN111096761B (zh) * | 2018-10-29 | 2024-03-08 | 上海西门子医疗器械有限公司 | 修正楔形滤波器散射的方法、装置和相关设备 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE2503789C3 (de) * | 1975-01-30 | 1980-10-09 | Philips Patentverwaltung Gmbh, 2000 Hamburg | Einrichtung zum Ermitteln der Absorption einer Strahlung in einer Ebene eines Körpers, mit einer Anordnung einer Strahlenquelle und einer Vielzahl von Strahlendetektoren, die relativ zu dem Körper kontinuierlich gedreht wird |
| US4029963A (en) * | 1976-07-30 | 1977-06-14 | The Board Of Trustees Of Leland Stanford Junior University | X-ray spectral decomposition imaging system |
| US4331877A (en) * | 1980-03-21 | 1982-05-25 | University Patents, Inc. | Transverse tomography with optical filtering |
| US4393352A (en) * | 1980-09-18 | 1983-07-12 | The Perkin-Elmer Corporation | Sample-and-hold hybrid active RC filter |
| US4547893A (en) * | 1981-06-12 | 1985-10-15 | Analogic Corporation | Continuous wave fan beam tomography system having a best-estimating filter |
| US4583240A (en) * | 1983-08-15 | 1986-04-15 | General Electric Company | Data acquisition circuitry for use in computerized tomography system |
| US5115394A (en) * | 1983-11-25 | 1992-05-19 | Technicare Corporation | Dual energy computerized tomography system |
| US4620313A (en) * | 1984-06-20 | 1986-10-28 | Technicare Corporation | Controlled cutoff frequency analog filter for CT scanners |
| US4815118A (en) * | 1987-06-29 | 1989-03-21 | General Electric Company | Data converter for CT data acquisition system |
-
1990
- 1990-08-06 US US07/562,766 patent/US5218533A/en not_active Expired - Fee Related
-
1991
- 1991-05-09 CA CA002042226A patent/CA2042226A1/en not_active Abandoned
- 1991-07-23 IL IL9892891A patent/IL98928A/en not_active IP Right Cessation
- 1991-07-29 JP JP3210504A patent/JPH0790023B2/ja not_active Expired - Lifetime
- 1991-08-01 DE DE69125668T patent/DE69125668T2/de not_active Expired - Fee Related
- 1991-08-01 EP EP91307057A patent/EP0473296B1/en not_active Expired - Lifetime
- 1991-08-03 CN CN91105322A patent/CN1056506C/zh not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1509686B (zh) * | 2002-11-27 | 2010-06-23 | Ge医药系统环球科技公司 | 定量组织脂肪含量的方法和装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| DE69125668T2 (de) | 1997-11-06 |
| IL98928A (en) | 1994-06-24 |
| CA2042226A1 (en) | 1992-02-07 |
| DE69125668D1 (de) | 1997-05-22 |
| US5218533A (en) | 1993-06-08 |
| EP0473296A3 (en) | 1992-10-28 |
| IL98928A0 (en) | 1992-07-15 |
| EP0473296A2 (en) | 1992-03-04 |
| EP0473296B1 (en) | 1997-04-16 |
| JPH04244138A (ja) | 1992-09-01 |
| CN1058710A (zh) | 1992-02-19 |
| JPH0790023B2 (ja) | 1995-10-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN1056506C (zh) | 用于双束计算机ｘ射线断层扫描术的稳定的可间断滤波器 | |
| US5361291A (en) | Deconvolution filter for CT system | |
| US6865247B2 (en) | Three dimensional back projection method and an X-ray CT apparatus | |
| US6166384A (en) | Method and apparatus for minimizing blurring and generating a high resolution image in a radiation imaging system | |
| JPH05324801A (ja) | 円及び線軌道を組み合わせて用いた円錐ビーム再構成 | |
| US20050063611A1 (en) | Super-resolution processor and medical diagnostic imaging apparatus | |
| CN1672637A (zh) | 利用多个偏移x－射线发射点成像的方法和系统 | |
| JPH0775603B2 (ja) | 検知装置 | |
| JPH08336522A (ja) | Ｃｔスキャナおよびｃｔスキャナ画像の生成方法 | |
| JP2004160218A (ja) | Ｘ線コンピュータ断層装置、ｘ線コンピュータ断層装置制御方法、及びｘ線コンピュータ断層撮影プログラム | |
| JPS61209641A (ja) | Ｘ線ｃｔ装置 | |
| CN1261078C (zh) | X射线ct扫描仪和图像处理器 | |
| US4707822A (en) | Tomographic apparatus | |
| IL116257A (en) | Helical interpolative algorithm for image reconstruction in a ct system | |
| JP2000033084A (ja) | コンピュ―タ利用断層撮影装置及びこのコンピュ―タ利用断層撮影装置を動作させる方法 | |
| JP2000139893A (ja) | Ｃｔ装置 | |
| CN110176045A (zh) | 一种由单能ct图像生成双能ct图像的方法 | |
| EP0851393A2 (en) | Methods and apparatus for simplified filtering of scan data in an imaging system | |
| Zeng et al. | Short-scan cone beam algorithm for circular and noncircular detector orbits | |
| CN113017663A (zh) | 一种ct扫描数据采集方法及系统、ct扫描仪 | |
| JP3285944B2 (ja) | コンピュータ断層撮影装置 | |
| JPH10314161A (ja) | 倍密度補間方法およびｘ線ｃｔ装置 | |
| KR100685560B1 (ko) | 데이터 보정 방법 및 ｘ선 ｃｔ 장치 | |
| CA2024475A1 (en) | System and method for collecting and transmitting pictorial information | |
| JP4224425B2 (ja) | Ｃｔ画像生成方法およびｘ線ｃｔ装置 |
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 | ||
| C15 | Extension of patent right duration from 15 to 20 years for appl. with date before 31.12.1992 and still valid on 11.12.2001 (patent law change 1993) | ||
| OR01 | Other related matters | ||
| C19 | Lapse of patent right due to non-payment of the annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |