CN1764930A - 成像内部结构 - Google Patents
成像内部结构 Download PDFInfo
- Publication number
- CN1764930A CN1764930A CNA2004800070367A CN200480007036A CN1764930A CN 1764930 A CN1764930 A CN 1764930A CN A2004800070367 A CNA2004800070367 A CN A2004800070367A CN 200480007036 A CN200480007036 A CN 200480007036A CN 1764930 A CN1764930 A CN 1764930A
- Authority
- CN
- China
- Prior art keywords
- image
- equipment
- dimension
- series
- scanner
- 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
- 238000003384 imaging method Methods 0.000 title claims description 6
- 238000000034 method Methods 0.000 claims abstract description 33
- 230000000737 periodic effect Effects 0.000 claims abstract description 11
- 230000029058 respiratory gaseous exchange Effects 0.000 claims abstract description 10
- 238000004458 analytical method Methods 0.000 claims abstract 2
- 230000008859 change Effects 0.000 claims description 9
- 230000005855 radiation Effects 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims description 7
- 230000006835 compression Effects 0.000 claims description 6
- 238000007906 compression Methods 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 4
- 230000004069 differentiation Effects 0.000 claims description 3
- 238000002591 computed tomography Methods 0.000 abstract description 25
- 230000033001 locomotion Effects 0.000 abstract description 7
- 239000002131 composite material Substances 0.000 abstract 1
- 230000002596 correlated effect Effects 0.000 abstract 1
- 230000008569 process Effects 0.000 description 13
- 210000000038 chest Anatomy 0.000 description 9
- 206010028980 Neoplasm Diseases 0.000 description 4
- 230000014509 gene expression Effects 0.000 description 4
- 210000004072 lung Anatomy 0.000 description 4
- 238000002360 preparation method Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 208000037841 lung tumor Diseases 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 230000007306 turnover Effects 0.000 description 1
- 230000004304 visual acuity 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/54—Control of apparatus or devices for radiation diagnosis
- A61B6/541—Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/0816—Measuring devices for examining respiratory frequency
-
- 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
- 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/04—Positioning of patients; Tiltable beds or the like
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Optics & Photonics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physiology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Radiation-Therapy Devices (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Image Processing (AREA)
Abstract
一种设备，方法和软件模块，用于从诸如锥形射束CT扫描仪的扫描仪输出中选取时相相关的图像，其工作方式是，把从二维系列图像中导出的图像压缩成一维，这是对与该一维垂直方向的像素强度求和实现的；从系列图像得到的一维图像组合中产生另一个图像；分析周期性图形的另一个图像；和从系列图像中选取该周期性图形中有类似时相的图像。如果需要，可以导出不同时相的多个重构图像。分析周期性图形的另一个图像可以包括：比较其中的一维图像，用于识别沿该一维的特征运动。这可以(尤其是)准确地确定病人的呼吸循环并利用时相相关的图像提高CT扫描的质量。
Description
技术领域
本发明涉及成像体积内部结构。本发明具体涉及(但不限于)分析CT图像以消除病人呼吸循环造成的人为现象和不准确性。
背景技术
现有的计算机X线断层照相术(CT)扫描仪依靠围绕病人转动的辐射源和检测器，并观察射束在从各个方向透射通过病人之后的衰减。根据这个数据，计算病人内部结构的三维表示。
这些扫描仪分成两种不同的类别。第一，简单的CT扫描仪引导狭窄的“笔型”射束通过病人。每个完整的转动导致通过病人的“切片”或断面形式二维图像。然后，扫描仪(或病人)沿转动轴变动位置，和扫描另一个切片。这些切片可以组合成三维图像。
第二，锥形射束CT扫描仪引导发散的射束指向病人，可以产生(在任何选取的时间)整个有用区(region of interest)的二维投影。当射束转动时，可以从不同的方向获取投影并构造三维图像。
一般地说，锥形射束CT优于简单的CT，因为它在各个方向的分辨率都是相同的。简单的CT在每个切片平面上有高分辨率，但是，垂直于这个切片平面的分辨率受到变动距离的限制。
这两种技术都假设病人是静止的。这对于有生命的病人是一个无效的假设，因为诸如心脏，肺和横膈膜的某些部位总是不可避免地运动。当病人呼吸时，围绕肺和横膈膜的结构就产生运动，这就给获得高质量扫描带来困难。在简单的CT中，图像中出现人为的现象；因为切片是沿病人变动位置，呼吸循环可以移动各种结构进出此时正在被扫描的切片。因此，在重构的体积中出现周期性人为现象。在锥形射束CT中，一些人为现象来自重构过程，但主要的问题是均化过程形成的图像模糊。
例如，在治疗肺部肿瘤时，重要的是知道肿瘤的位置以及在呼吸循环进展中肿瘤的运动。从锥形射束CT中导出的时间平均信息是不合适的，因为它不能区分一个区域中驻留的小肿瘤和短暂通过一个区域的大肿瘤。
至今，在进行胸廓的CT扫描时，要求病人根据外部的刺激控制他们的呼吸，或检测呼吸循环时相(phase)的替代物。已使用过的替代物例子包括围绕鼻孔的局部温度，和胸廓的大小。这些替代物已证明是有帮助的，但是一般不能令人满意。
发明内容
因此，准确地确定病人的呼吸循环是各种用途的重要诊断和治疗计划手段。
所以，本发明提供一种用于成像体积内部结构以展示内部变化的设备，包括：穿透辐射源和用于该辐射的二维检测器，辐射源和检测器安排成产生一系列的体积投影图像；重构装置，用于从选取的系列图像中导出有关体积中三维结构的信息；选择装置，用于从系列中选取有类似时相的图像供给重构装置使用，其中选择装置安排成把从二维系列图像中导出的图像压缩成一维图像，这是沿与该一维垂直的方向对像素强度求和实现的，从系列图像得到的一维图像组合中产生另一个图像，分析图形的另一个图像，和从系列图像中选取该图形中有类似时相的图像。
可以想象到，利用计算机装置实现图像处理，其形式是配置合适程序的通用计算机装置或特别为这个任务设计的专用计算机装置。在这种情况下，所涉及的图像通常是数字形式，其中一个或多个文件中存储的数值对应于该图像中特定位置的强度或其他的图像特征。
因此，当我们讨论图像操作步骤时，我们并不试图推断必须按照这种方式实际操作图像，而是包括利用算术运算处理数值或数组的过程。同样地，计算机装置不必实际显示，描绘或其他方式产生临时图像，例如，以上所涉及的另一个图像，只要实施算术运算或其他过程可以得到相应的数值数据。
若在选择装置工作之前对图像进行预处理，则可以得到提高的准确度。最好是，这种预处理操作包括若干步骤，例如，用于压缩图像强度范围的压缩滤波器和增强图像边缘的导数滤波器。在体积包括病人的情况下，最好是沿病人的颅索轴方向实施取导数操作，因为横膈膜有垂直于该方向的边缘。
预处理操作可以包括应用于图像的掩模，例如，为了选取包含边缘的区域和排除有用区之外的区域。可利用的掩模是这样一种掩模，它排除体积内目标(例如，病人)之外的图像区。可以从应用于图像的阈值导出边缘掩模，该图像经导数滤波器的滤波。
如果需要，可以从多个子集合中导出多个重构图像，每个子集合包含系列图像中时相相关的图像，每个子集合的时相相关不同于其他子集合的时相相关。
最好是，分析周期性图形的另一个图像包括比较其中一维图像的步骤，用于识别沿该一维的特征运动。最好是比较相邻的图像，因为这些图像大体上是相似的，因此，它们之间的差别很可能是最重要的。通过计算不同相对位移图像的图像之间强度差，可以比较这些一维图像。均方根(rms)差通常是最佳的差量度。
在分析周期性图形之前，最好也对另一个图像进行处理。合适的预处理操作包括用于增强另一个图像边缘的导数滤波器，最好是沿垂直于该一维的方向实施导数滤波，和/或选取有用区。特别有利的是，利用导数滤波后图像增强有锐变化的区域，从而选取部分的原始另一个图像。
辐射源和检测器可以相对于该体积转动，因此，系列投射图像可以展示不同取向的体积。
本发明特别适合于周期性变化，但事实上可应用于体积中的任何变化。
一般地说，本发明找到最有利的应用是，周期性变化是活生物体展示的自然变化，例如，呼吸循环。然而，本发明有广泛的应用，并不必限于这种形式。
本发明还提供一种用于从扫描仪输出中选取时相相关图像的方法，包括以下步骤：对于扫描仪产生系列图像的多个图像中每个图像，把二维图像压缩成一维图像，这是沿与该一维垂直的方向对像素强度求和实现的，从系列图像得到的一维图像组合中产生另一个图像，分析周期性图形的另一个图像，和从系列图像中选取该周期性图形中有类似时相的图像。
本发明还提供一种软件模块，用于从扫描仪的输出中选取时相相关的图像，对于扫描仪产生系列图像的多个图像中每个图像，安排成执行以下的步骤：把二维图像压缩成一维图像，这是沿与该一维垂直的方向对像素强度求和实现的，从系列图像得到的一维图像组合中产生另一个图像，分析周期性图形的另一个图像，和从系列图像中选取该周期性图形中有类似时相的图像。
目前，本发明的这些附加特征最适合于锥形射束CT扫描仪，特别是用于呼吸病人图像的CT扫描仪。
应当强调的是，本发明对于分析各种结构的内部排列有潜在的广泛应用，因此，它并不限于呼吸时相的检测。可以检测活体和非活体结构中的其他变化，而且，本发明可适用于锥形射束CT扫描仪之外的其他扫描仪。这个特点源于本发明的性质，它并不试图具体识别与呼吸有关的结构(例如，横膈膜)，而是通过新颖的处理路线仅仅寻求各种变化。
因此，在应用于识别呼吸时相的问题时，本发明的确可以识别横膈膜，但是通过直接寻找该形状的对象进行识别。的确，本发明发现如此做可以起到反作用，在扫描仪围绕病人转动时，横膈膜的投影图像改变形状。
附图说明
现在借助于例子并参照附图描述本发明一个实施例，其中：
图1表示利用简单CT扫描仪得到没有呼吸相关的人胸廓典型图像；
图2表示利用锥形射束CT扫描仪得到没有呼吸相关的人胸廓典型图像；
图3表示从锥形射束CT扫描仪得到系列图像中的人胸廓图像；
图4表示从图3所示图像中导出的病人掩模；
图5和6表示在均衡操作之前和之后的图3所示图像；
图7和8表示在增强边缘的滤波操作之前和之后的图3所示图像；
图9和10表示边缘掩模的制备；
图11和12表示应用梯度规范滤波器；
图13和14表示应用病人掩模和边缘掩模的结果；
图15和16表示把图像压缩成一维图像；
图17表示另一个图像；
图18表示滤波后的另一个图像，滤波操作可以增强与一维垂直方向上的变化；
图19表示选取另一个图像中有用区的方法；
图20表示有用区；
图21表示相邻一维图像的比较结果；
图22表示应用于图像的过程路线概要图；
图23表示得到的重构图像。
具体实施方式
图1表示病人胸廓区的常规CT扫描图。例如，在100，102和104可以看到周期性人为现象，这些人为现象对应于呼吸过程与CT扫描过程相互作用造成的误差。同样地，图2表示病人胸廓的典型锥形射束CT图像，和区域106，108表示病人肺和横膈膜的内部运动造成的严重模糊。
现在描述按照本发明的过程。图22表示大致的过程路线，其中采用的参考数字对应于本申请中以前的相关图号，它说明在该处得到的图像或掩模。因此，在利用锥形射束CT扫描仪或其他类似扫描仪采集图像之后或期间，过程路线A应用于每个图像3。首先，从原始图像3开始，我们得到病人掩模4，它能掩蔽病人之外的部分图像。为此目的，可以设定简单的阈值，而原始图像3中强度超过阈值的像素或像素组可以分配到掩模的已占区。一般地说，试凑方法适合于确定合适的病人掩模，它与所用具体扫描仪的参数有关，但在本发明者使用图像的情况下，我们发现皮肤阈值为20,000是合适的。
再从CT扫描仪得到的原始图像3开始，可以均衡图像以减小对比度的范围。有许多可利用的合适均衡滤波器，但在本发明的例子中，本发明者采用比例对数方法，其中各个强度值被转换成原始强度值对数的1000倍。
从均衡的图像开始，制备导数图像，导数图像是其强度值代表预滤波图像中强度变化率的一种图像。导数滤波器必须有测量变化率的相关方向，而在病人的胸廓区图像情况下，最好是利用颅索(CC)的方向。这可以增强沿该方向的变化。在一般的情况下，这个导数的方向通常应当与以下步骤16中的方向一致。从导数图像开始，可以提供导数掩模10，它是一种用于增强具有高导数值像素的掩模。所以，它增强原始均衡图像中强度变化率高的区域，即，均衡图像6中有边缘的区域。
最后，从均衡的图像开始，制备梯度规范图像12。这是一种均衡导数图像的形式，其中对水平梯度和垂直梯度进行矢量求和，即，水平梯度的平方与垂直梯度的平方相加，再取这个结果的平方根。
然后，导数掩模10和病人掩模4可应用于梯度规范图像12以产生处理后图像14。所以，这个图像包含均衡图像的像素，它们是在病人内和强度变化率高的区域。实质上，这个过程是从病人内的均衡图像中选出边缘，其中略微优选CC对准边缘。这种优选是不完全的，因为虽然CC边缘是部分的掩模，但其他的边缘是以梯度规范形式被部分地引入。应当注意，该算法对于如何提取边缘是故意地不灵敏，因为识别横膈膜的试图往往是失败的。只要相关的特征(例如，横膈膜)是在边缘中，该算法可以检测运动远远超过所有其他边缘的边缘。
处理扫描仪产生系列图像中这个图像的最后步骤是把图像压缩成一维图像，在病人图像的情况下，这个一维方向应当与病人的CC轴对准。可以通过以下方法实现，沿一个水平行像素的总强度求和，并采用该强度或由此导出的强度作为该点沿水平轴的单个像素强度。因此，利用这个过程，把扫描仪得到的每个图像压缩成简单的直线图像，即，一个像素宽和“n”个像素高。
然后，对从扫描仪得到系列图像中的全部图像或基本上全部图像重复这个过程。这可以作为采集整个系列图像之后的后处理步骤，或者，在查询以后图像的同时，可以从早期的边缘开始。一旦获得一维图像的系列，则可以利用并置方法组合这些一维图像以产生组合图像17，其中一个轴对应于病人的CC轴(或一维压缩操作采用的轴)，而另一个轴对应于时间的量度(表示成图像编号)。图17表示得到的图像。虽然它并不包含关于病人整个结构的有用信息，但它确实表示图像随时间变化的有用信息。利用应用于组合图像的系列处理步骤B可以提取这个信息。
首先，应用导数滤波器，导数方向对应于像素编号或时间轴。这产生图18所示的图像，其中可以实际增强随时间快速变化的区域。
把以前图像(图18)的绝对值再映射到CC轴上以形成强度直方图。这个一维信号增强沿病人CC轴的位置，其中有随图像编号或时间的最大变化率。根据强度直方图，可以应用阈值以选取特定的有用区20。各种取阈值方法是已知的，并适合于这种任务，它包括：设置简单的固定阈值，或可变化阈值，它寻找强度直方图的最大值和最小值，并采用二者之间一点的阈值。有用区可以限制于从直方图选取的正确区域，或可以是任何一侧的边缘。
在选取有用区之后，我们选择和分析周期性图形的这个图像区17。从图20中可以看出，提取的有用区中有一个清晰的周期性图形，在扫描结束时重新会聚之前，虽然这个有用区发散成一对周期性图形。这大致对应于病人横膈膜左侧和右侧的图像。
为了克服这个困难和其他的困难，我们提出一个极其简单的比较过程，其中在经受可变位移之后，比较图像中每个单像素“切片”与下一个邻近切片。然后，该系统确定在什么位移下这两个一维切片之间的差最小。如图21所示，一种方法是对两个单像素切片之间强度差平方求和，并识别该和值在什么像素位移下最小。在重复每个相邻对单像素切片的运算之后，可以导出图21所示的曲线，其中位移是随X线图像编号而变化。这是一个简单的周期性曲线，利用已知的方法可以分析该曲线以确定它的频率，从而确定每个图像的时相位置。
然后，该系统可以返回到每个原始图像并给它分配时相信息。收集的图像可以组织成多个时相相关族，而单个相关族内的图像可用于正常的锥形射束CT重构过程。这个结果表示在图23中，并展示较高的精确度和清晰度。可以大大消除呼吸的人为现象。
在本发明的实验应用中，提取三个病人的4个X线系列呼吸循环。我们发现，优化第一系列重构所采用的参数值证明适合于所有的系列。具体地说，三个病人中的一个病人只有一个肺，但由于在X线图像中呼吸循环是可见的，就可以成功地提取呼吸循环。因此，本发明提供一种从顺序X线图像中提取呼吸循环的独立算法，它不需要用户的相互作用，并可在足够快速处理机器上可用的时间内产生高质量重构图像。假设算法所需少量参数的原始初值证明适合于整个测试系列，则该方法是相应健全的。此外，算法的相对简单性和大部分的处理是在逐个图像的基础上，这意味着大多数处理操作可以在X线图像采集时间内完成。在适当快速的处理器上，用于显示呼吸循环频率的组合图像处理不应当超过约1秒钟。
因为呼吸循环是直接从相关X线图像中导出的，而不是根据一些辅助指示，例如，鼻孔温度和胸廓大小，它的不准确性就较少，因为呼吸循环时相与X线图像之间相关是直接的。
应当明白，在不偏离本发明范围的条件下，对以上描述的实施例作各种变化是可能的。
Claims (29)
1.一种用于成像体积内部结构以展示内部变化的设备，包括：
穿透辐射源和用于该辐射的二维检测器，辐射源和检测器安排成产生一系列的体积投影图像；
重构装置，用于从选取的系列图像中导出有关体积中三维结构的信息；
选择装置，用于从系列图像中选取有类似时相的图像供给重构装置使用；
其中选择装置安排成：
把从二维系列图像中导出的图像压缩成一维图像，这是沿与该一维垂直的方向对像素强度求和实现的，
从系列图像得到的一维图像组合中产生另一个图像，
分析图形的另一个图像，和
从系列图像中选取该图形中有类似时相的图像。
2.按照权利要求1的设备，其中辐射源和检测器可以相对于该体积旋转，因此，该系列投影图像展示不同取向的体积。
3.按照权利要求1或权利要求2的设备，其中变化是周期性的。
4.按照权利要求1的设备，其中图像在选择装置运行之前被预处理。
5.按照权利要求4的设备，其中预处理操作包括用于缩小图像强度范围的滤波器。
6.按照权利要求4或权利要求5的设备，其中预处理操作包括用于增强图像边缘的导数滤波器。
7.按照权利要求6的设备，其中该体积包括病人，和取导数操作是沿病人颅索轴的方向实施的。
8.按照权利要求4至7中任何一个的设备，其中预处理操作包括应用于的图像的掩模，用于选取包含边缘的区域。
9.按照权利要求8的设备，其中从应用于图像的阈值导出边缘掩模，该图像经导数滤波器的滤波。
10.按照权利要求4至9中任何一个的设备，其中预处理操作包括加掩模以排除该体积内目标之外的图像区域。
11.按照以上权利要求中任何一个的设备，其中多个重构图像是从多个子集合中导出，每个子集合包含系列图像中时相相关的图像，每个子集合的时相相关不同于其他子集合的时相相关。
12.按照以上权利要求中任何一个的设备，其中分析另一个图像的周期性图形包括比较步骤：比较其中的一维图像，用于识别沿该一维特征的运动。
13.按照权利要求12的设备，其中比较相邻的图像。
14.按照权利要求12或13的设备，其中比较一维图像是通过计算各个图像中不同相对位移的图像之间强度差。
15.按照权利要求14的设备，其中比较方均根差。
16.按照以上权利要求中任何一个的设备，其中在分析周期性图形之前，另一个图像经受处理。
17.按照权利要求16的设备，其中预处理操作包括增强另一个图像边缘的导数滤波器。
18.按照权利要求17的设备，其中取导数操作是在与该一维垂直的方向实施。
19.按照权利要求16至18中任何一个的设备，其中预处理操作包括选取另一个图像的有用区和排除另一种处理操作的其他区。
20.按照权利要求19的设备，其中通过分析包含最高导数的区域，选取有用区。
21.按照以上权利要求中任何一个的设备，其中该变化是活生物体展示的自然变化。
22.按照权利要求21的设备，其中该变化是由呼吸造成的。
23.一种从扫描仪输出中选取时相相关图像的方法，对于扫描仪产生系列图像的多个图像中每个图像，包括以下步骤：
把二维图像压缩成一维图像，这是沿与该一维垂直的方向对像素强度求和实现的，
从系列图像得到的一维图像组合中产生另一个图像，
分析周期性图形的另一个图像，和
从系列图像中选取该周期性图形中有类似时相的图像。
24.按照权利要求23的方法，其中扫描仪是锥形射束CT扫描仪。
25.按照权利要求24的方法，其中图像是呼吸病人的图像。
26.一种从扫描仪输出中选取时相相关图像的软件模块，对于扫描仪产生系列图像的多个图像中每个图像，该软件模块安排成执行以下步骤：
把二维图像压缩成一维图像，这是沿与该一维垂直的方向对像素强度求和实现的，
从系列图像得到的一维图像组合中产生另一个图像，
分析周期性图形的另一个图像，和
从系列图像中选取该周期性图形中有类似时相的图像。
27.按照权利要求26的方法，其中扫描仪是锥形射束CT扫描仪。
28.按照权利要求27的方法，其中图像是呼吸病人的图像。
29.一种用于成像体积内部结构的设备，该设备基本上是此处参照附图所描述的设备。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| GB0301278.8 | 2003-01-21 | ||
| GB0301278A GB2397738B (en) | 2003-01-21 | 2003-01-21 | Computed tomography scanning |
| GB0327675.5 | 2003-11-28 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1764930A true CN1764930A (zh) | 2006-04-26 |
| CN100410970C CN100410970C (zh) | 2008-08-13 |
Family
ID=9951464
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2004800017740A Expired - Lifetime CN100484475C (zh) | 2003-01-21 | 2004-01-20 | 计算层析x射线摄影扫描 |
| CNB2004800070367A Expired - Fee Related CN100410970C (zh) | 2003-01-21 | 2004-01-20 | 成像内部结构 |
Family Applications Before (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2004800017740A Expired - Lifetime CN100484475C (zh) | 2003-01-21 | 2004-01-20 | 计算层析x射线摄影扫描 |
Country Status (6)
| Country | Link |
|---|---|
| EP (1) | EP1589876B2 (zh) |
| JP (1) | JP4499087B2 (zh) |
| CN (2) | CN100484475C (zh) |
| AT (1) | ATE529044T1 (zh) |
| GB (1) | GB2397738B (zh) |
| WO (1) | WO2004064641A1 (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101559260A (zh) * | 2008-04-16 | 2009-10-21 | 西门子公司 | 用于实施照射的装置以及监视该照射的方法 |
| CN101743568B (zh) * | 2007-05-17 | 2012-04-18 | 伊利克塔股份有限公司 | 在线锥形束ct重建 |
| CN101536039B (zh) * | 2006-10-31 | 2013-02-06 | 皇家飞利浦电子股份有限公司 | 用于基于体积图像数据集绘制图像的系统、方法和装置 |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2397738B (en) † | 2003-01-21 | 2007-08-29 | Elekta Ab | Computed tomography scanning |
| GB2416101A (en) * | 2004-07-06 | 2006-01-11 | Elekta Ab | Imaging internal structures and correcting for motion artifacts |
| KR20080044250A (ko) | 2005-07-23 | 2008-05-20 | 토모테라피 인코포레이티드 | 갠트리 및 진료대의 조합된 움직임을 이용하는 방사선치료의 영상화 및 시행 |
| CN100574827C (zh) * | 2005-08-25 | 2009-12-30 | 深圳市海博科技有限公司 | 放射治疗装置 |
| JP2009515586A (ja) * | 2005-11-10 | 2009-04-16 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 画像を生成するための信号処理ユニット |
| CN101161203B (zh) * | 2006-10-11 | 2011-01-26 | 株式会社东芝 | X射线计算机断层及医用图像摄影装置、呼吸指示装置 |
| CN101820948A (zh) | 2007-10-25 | 2010-09-01 | 断层放疗公司 | 用于放疗实施的运动适应性优化的系统和方法 |
| US8467497B2 (en) | 2007-10-25 | 2013-06-18 | Tomotherapy Incorporated | System and method for motion adaptive optimization for radiation therapy delivery |
| US9224202B2 (en) * | 2011-04-14 | 2015-12-29 | Koninklijke Philips N.V. | Device and method for extracting information from characteristic signals |
| US8897527B2 (en) * | 2011-06-07 | 2014-11-25 | Varian Medical Systems, Inc. | Motion-blurred imaging enhancement method and system |
| US20130035588A1 (en) * | 2011-08-03 | 2013-02-07 | Siemens Corporation | Magnetic resonance imaging for therapy planning |
| EP2962309B1 (en) | 2013-02-26 | 2022-02-16 | Accuray, Inc. | Electromagnetically actuated multi-leaf collimator |
| JP6181459B2 (ja) | 2013-08-05 | 2017-08-16 | 株式会社日立製作所 | 放射線治療システム |
Family Cites Families (39)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA1101936A (en) * | 1977-04-01 | 1981-05-26 | Robert H. Wake | Cardiac imaging with ct scanner |
| JPS58206724A (ja) * | 1982-05-27 | 1983-12-02 | 株式会社東芝 | 心拍連動ｘ線ｃｔ装置 |
| JPS62104388A (ja) * | 1985-10-31 | 1987-05-14 | Shimadzu Corp | Ｘ線画像処理装置 |
| US5254948A (en) * | 1988-04-01 | 1993-10-19 | Hitachi, Ltd. | Method and apparatus for body motion corrective imaging |
| US4994965A (en) * | 1988-11-23 | 1991-02-19 | General Electric Company | Method for reducing motion induced image artifacts in projection imaging |
| US5137026A (en) † | 1990-01-04 | 1992-08-11 | Glaxo Australia Pty., Ltd. | Personal spirometer |
| US5251128A (en) * | 1990-11-19 | 1993-10-05 | General Electric Company | Motion artifact reduction in projection imaging |
| JPH05212024A (ja) * | 1992-02-04 | 1993-08-24 | Toshiba Corp | Ｃｔ装置 |
| JP3512875B2 (ja) * | 1993-11-26 | 2004-03-31 | 東芝医用システムエンジニアリング株式会社 | Ｘ線コンピュータ断層撮影装置 |
| EP0734742B1 (en) * | 1995-03-31 | 2005-05-11 | Kabushiki Kaisha Toshiba | Ultrasound therapeutic apparatus |
| US6076005A (en) * | 1998-02-25 | 2000-06-13 | St. Jude Children's Research Hospital | Respiration responsive gating means and apparatus and methods using the same |
| US6937696B1 (en) † | 1998-10-23 | 2005-08-30 | Varian Medical Systems Technologies, Inc. | Method and system for predictive physiological gating |
| JP4266422B2 (ja) * | 1998-12-28 | 2009-05-20 | ジーイー横河メディカルシステム株式会社 | 放射線断層撮影方法および装置 |
| JP3326597B2 (ja) * | 1999-01-14 | 2002-09-24 | 独立行政法人放射線医学総合研究所 | 呼吸同期制御放射線治療機器 |
| JP2000201920A (ja) * | 1999-01-19 | 2000-07-25 | Fuji Photo Film Co Ltd | 撮影画像デ―タ取得方法および撮影画像デ―タ取得装置 |
| US6535570B2 (en) † | 1999-06-17 | 2003-03-18 | Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of National Defence Of Her Majesty's Canadian Government | Method for tracing organ motion and removing artifacts for computed tomography imaging systems |
| EP1061474A1 (de) † | 1999-06-17 | 2000-12-20 | VAMP Verfahren und Apparate der Medizinischen Physik GmbH | Computertomograph mit objektbezogener Bewegungsartefaktreduktion und Extraktion der Objektbewegungsinformation (Kymogramm) |
| AU5561400A (en) * | 1999-07-02 | 2001-01-22 | Ultra-Guide Ltd. | Apparatus and methods for medical interventions |
| US6252924B1 (en) † | 1999-09-30 | 2001-06-26 | General Electric Company | Method and apparatus for motion-free cardiac CT imaging |
| EP2311527B1 (en) † | 2000-02-18 | 2019-08-28 | William Beaumont Hospital | Cone-beam computerized tomography with a flat-panel imager |
| WO2001067960A2 (en) * | 2000-03-16 | 2001-09-20 | Cedara Software Corp. | System and method for processing an image sequence |
| JP4733809B2 (ja) * | 2000-05-23 | 2011-07-27 | 株式会社東芝 | 放射線治療計画装置 |
| AU7046601A (en) * | 2000-05-30 | 2001-12-17 | Patrik Rogalla | Method for the combined representation of morphology and dynamics in split-imageand volume-image methods |
| JP3442346B2 (ja) † | 2000-06-01 | 2003-09-02 | カナガワ アンド カンパニー株式会社 | 画像形成装置及びこれを用いた画像の構成方法 |
| US7031504B1 (en) † | 2000-09-26 | 2006-04-18 | Vital Images, Inc. | Image data based retrospective temporal selection of medical images |
| ATE233517T1 (de) * | 2000-11-22 | 2003-03-15 | Brainlab Ag | Verfahren zur bestimmung der lungenfüllung |
| JP2002210029A (ja) † | 2001-01-19 | 2002-07-30 | Mitsubishi Electric Corp | 放射線治療装置 |
| US6480560B2 (en) † | 2001-03-16 | 2002-11-12 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for motion gating using CT projection data |
| EP1244058A1 (de) † | 2001-03-22 | 2002-09-25 | VAMP Verfahren und Apparate der Medizinischen Physik GmbH | Computertomograph mit merkmalsgewichteter Bild-zu-Volumen Interpolation in Echtzeit |
| DE10117751A1 (de) * | 2001-04-09 | 2002-11-21 | Siemens Ag | Modellierung von Momentanzuständen medizinischer Objekte in Abhängigkeit von zumindest einer Messgröße |
| US7058440B2 (en) * | 2001-06-28 | 2006-06-06 | Koninklijke Philips Electronics N.V. | Dynamic computed tomography imaging using positional state modeling |
| US7006862B2 (en) † | 2001-07-17 | 2006-02-28 | Accuimage Diagnostics Corp. | Graphical user interfaces and methods for retrospectively gating a set of images |
| US7142703B2 (en) † | 2001-07-17 | 2006-11-28 | Cedara Software (Usa) Limited | Methods and software for self-gating a set of images |
| DE10162768A1 (de) † | 2001-12-20 | 2003-07-03 | Philips Intellectual Property | Computertomograph |
| US6879656B2 (en) † | 2002-07-23 | 2005-04-12 | General Electric Company | Method and apparatus for deriving motion information from projection data |
| US6888914B2 (en) † | 2002-11-26 | 2005-05-03 | General Electric Company | Methods and apparatus for computing volumetric perfusion |
| AU2003279484A1 (en) † | 2002-12-13 | 2004-07-09 | Koninklijke Philips Electronics N.V. | System and method for processing a series of image frames representing a cardiac cycle |
| JP4589124B2 (ja) * | 2003-01-09 | 2010-12-01 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | コンピュータ断層撮影のための呼吸モニタ |
| GB2397738B (en) † | 2003-01-21 | 2007-08-29 | Elekta Ab | Computed tomography scanning |
-
2003
- 2003-01-21 GB GB0301278A patent/GB2397738B/en not_active Expired - Lifetime
-
2004
- 2004-01-20 EP EP04703444.2A patent/EP1589876B2/en not_active Expired - Lifetime
- 2004-01-20 CN CNB2004800017740A patent/CN100484475C/zh not_active Expired - Lifetime
- 2004-01-20 AT AT04703444T patent/ATE529044T1/de not_active IP Right Cessation
- 2004-01-20 CN CNB2004800070367A patent/CN100410970C/zh not_active Expired - Fee Related
- 2004-01-20 JP JP2006500227A patent/JP4499087B2/ja not_active Expired - Lifetime
- 2004-01-20 WO PCT/GB2004/000220 patent/WO2004064641A1/en active Application Filing
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101536039B (zh) * | 2006-10-31 | 2013-02-06 | 皇家飞利浦电子股份有限公司 | 用于基于体积图像数据集绘制图像的系统、方法和装置 |
| CN101743568B (zh) * | 2007-05-17 | 2012-04-18 | 伊利克塔股份有限公司 | 在线锥形束ct重建 |
| CN101559260A (zh) * | 2008-04-16 | 2009-10-21 | 西门子公司 | 用于实施照射的装置以及监视该照射的方法 |
| CN101559260B (zh) * | 2008-04-16 | 2014-05-28 | 西门子公司 | 用于实施照射的装置以及监视该照射的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| ATE529044T1 (de) | 2011-11-15 |
| EP1589876B2 (en) | 2021-02-17 |
| CN1722983A (zh) | 2006-01-18 |
| EP1589876B1 (en) | 2011-10-19 |
| CN100484475C (zh) | 2009-05-06 |
| GB2397738B (en) | 2007-08-29 |
| WO2004064641A1 (en) | 2004-08-05 |
| GB0301278D0 (en) | 2003-02-19 |
| WO2004064641A8 (en) | 2004-09-10 |
| JP4499087B2 (ja) | 2010-07-07 |
| CN100410970C (zh) | 2008-08-13 |
| EP1589876A1 (en) | 2005-11-02 |
| GB2397738A (en) | 2004-07-28 |
| JP2006516417A (ja) | 2006-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN100410970C (zh) | 成像内部结构 | |
| US7349564B2 (en) | Imaging internal structures | |
| US8798353B2 (en) | Apparatus and method for two-view tomosynthesis imaging | |
| CN1762305A (zh) | 用具角位移的射线束产生层析图像的方法和对应的ct设备 | |
| CN1516854A (zh) | 数字层析x射线照相组合中的通用滤波反投影重构 | |
| CN1647768A (zh) | 用焦点探测器组合产生周期运动物体断层造影图像的方法 | |
| EP2147410B1 (en) | On-line cone beam CT reconstruction | |
| CN1957847A (zh) | 再现对象的断层造影图像的方法和断层造影设备 | |
| CN1977287A (zh) | 伪像缩减 | |
| CN1565001A (zh) | 从投影数据得出运动信息的方法和设备 | |
| CN1206581A (zh) | 射线层析摄影方法和装置 | |
| CN1515228A (zh) | X射线ct装置 | |
| CN1947151A (zh) | 用于在图像中使用发散梯度场响应进行基于滑降的对象分割的系统和方法 | |
| CN103679706B (zh) | 一种基于图像各向异性边缘检测的ct稀疏角度重建方法 | |
| CN1833613A (zh) | 用于产生对象断层造影图像的方法和计算机断层造影系统 | |
| CN1487478A (zh) | Ct中调节清晰度和噪声的图像回溯滤波及窗口控制滤波 | |
| CN101056584A (zh) | 用于检查周期性运动的对象的ct方法 | |
| CN1742296A (zh) | 锥形光束计算机断层成像的立体像素驱动的螺旋重建方法 | |
| CN101065065A (zh) | X射线ct装置和图像处理装置 | |
| CN1615801A (zh) | 计算机辅助图像获取和诊断系统 | |
| CN1934591A (zh) | 锥束计算机层析成像重构 | |
| JP4982091B2 (ja) | 放射線断層撮影装置 | |
| CN1620988A (zh) | 用于扫描活动对象的时间序列的方法和设备 | |
| CN1644167A (zh) | 用焦点探测器组合产生周期运动物体断层造影图像的方法 | |
| CN1613427A (zh) | Ct图像产生方法和x射线ct设备 |
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 | ||
| C17 | Cessation of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080813 Termination date: 20120120 |