CN103400401A - 改善图像分辨率的系统和方法 - Google Patents

改善图像分辨率的系统和方法 Download PDF

Info

Publication number
CN103400401A
CN103400401A CN2013102333828A CN201310233382A CN103400401A CN 103400401 A CN103400401 A CN 103400401A CN 2013102333828 A CN2013102333828 A CN 2013102333828A CN 201310233382 A CN201310233382 A CN 201310233382A CN 103400401 A CN103400401 A CN 103400401A
Authority
CN
China
Prior art keywords
image
voxel
estimated value
function
projection
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2013102333828A
Other languages
English (en)
Other versions
CN103400401B (zh
Inventor
吉昂.西
琼-巴普蒂斯特.蒂鲍尔特
肯.D.索尔
查尔斯.A.鲍曼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
General Electric Co
University of Notre Dame
Purdue Research Foundation
Original Assignee
General Electric Co
University of Notre Dame
Purdue Research Foundation
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by General Electric Co, University of Notre Dame , Purdue Research Foundation filed Critical General Electric Co
Publication of CN103400401A publication Critical patent/CN103400401A/zh
Application granted granted Critical
Publication of CN103400401B publication Critical patent/CN103400401B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • 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/58—Testing, adjusting or calibrating thereof
    • A61B6/582—Calibration
    • A61B6/583—Calibration using calibration phantoms
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/424—Iterative
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S378/00—X-ray or gamma ray systems or devices
    • Y10S378/901—Computer tomography program or processor

Landscapes

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

Abstract

本发明描述了一种断层摄影成像系统中重建图像的方法。该方法包括通过迭代重建图像改善图像的空间分辨率。

Description

改善图像分辨率的系统和方法
本申请是基于申请日为2007年6月22日，优先权日为2006年6月22日，申请号为200780023458.7（PCT/US2007/071917），发明名称为：“改善图像分辨率的系统和方法”的专利申请的分案申请。
发明背景
本发明总体上涉及成像系统，更具体地，涉及改善图像分辨率的系统和方法。
计算机断层摄影(CT)成像系统通常包括x射线源，x射线源投射扇形x射线束透过患者到达辐射探测器阵列。x射线束经准直落在通常称为“成像平面”的xy平面内。探测器阵列接收到的x射线束的辐射强度取决于患者使x射线束的衰减。单独采集探测器阵列多个探测器单元的衰减测量结果而生成透射谱。
在台架内围绕受检患者转动x射线源和探测器，以使x射线束与患者相交所成的投射角不断变化。将一组x射线衰减测量结果称为“视角(view)”，所述测量结果是来自成一种台架角度或一种投射角度的探测器阵列的模拟投影数据。对患者的“扫描”包括在x射线源和探测器旋转一周的过程中以不同的投射角形成的一组视角。
为缩短获取多个切片所用的总扫描时间，可进行螺旋扫描。螺旋扫描技术允许使用单光子源以较快的速度进行大体积扫描。为进行螺旋扫描，使患者静躺于其上的平台沿台架旋转所围绕的z轴移动，同时采集规定数量切片的模拟投影数据。螺旋扫描生成单螺线。x射线束绘出的螺线形成模拟投影数据，可根据所述模拟投影数据，重建各规定切片中的图像。除缩短扫描时间以外，螺旋扫描还具有其它优势，例如更好地使用注入对比度(injectedcontrast)、任意位置处改善的图像重建和较好的三维图像。螺旋扫描的实例包括多切片螺旋扫描。在多切片螺旋扫描中，探测器阵列沿z轴延伸。通常，在多切片螺旋扫描中，探测器阵列包括多行，每一行相应于沿z轴不同的位置和不同的测量切片。在轴向扫描中，处理模拟投影数据以构建图像，该图像相应于透过患者所取的二维切片。对于离散切片，可进行全视场迭代重建，以提高图像质量。
对于连续扫描，形成这样的扫描图案，其中患者沿z轴的部位随着台架的旋转而线性变化。在采集数据过程中，对连续扫描图案进行量化，生成患者周围有限个x射线源位置的一组离散的投影视角。常规的直接图像重建技术，例如二维或三维滤波反投影，根据投影数据如下重建图像体素：在投影数据中内插要素以将各投影角的贡献累加到多个图像体素中，从而通过对投影数据的单遍扫描形成图像或图像体(image volume)。应用滤波反投射生成的图像的典型分辨率基于探测器阵列的尺寸、焦斑尺寸、数据采集系统(DAS)采集模拟投影数据的采样速度和滤波反投射过程中过滤投影数据的滤波器的kernel。通常情况下，典型分辨率不高于CT成像系统等角点处各探测单元的投影尺寸。根据尼奎斯特定理，无需以大于两倍的极限典型分辨率采样。然而，通过常规的直接图像重建技术生成的图像体通常不具备高的空间分辨率。
发明内容
一方面，描述断层摄影成像系统中重建图像的方法。该方法包括通过迭代重建图像改善图像的空间分辨率。
另一方面，描述权衡(trade-off)断层摄影成像系统性能的迭代重建方法。该方法包括调整一部分迭代重建图像的噪声和分辨率。
再一方面，描述调整断层摄影成像系统分辨率的迭代重建方法。该方法包括调整迭代重建图像的面内分辨率(in-plane resolution)，以及调整迭代重建图像的断面分辨率(cross-plane resolution)。
再一方面，描述改善图像空间分辨率的方法。该方法包括迭代重建图像并改善图像的空间分辨率。该方法通过建立具有平滑曲线的正投影函数和建立具有平滑曲线的图像中的至少一种来改善图像的空间分辨率。正投影的建立包括通过确定一系列投影值来建立正投影。图像的建立包括通过确定所述一系列投影值的逆(inverse)来重建图像。
再一方面，描述改善图像空间分辨率的处理器。配置该处理器以迭代重建图像并改善图像的空间分辨率。配置该处理器，以通过建立具有平滑曲线的正投影函数和建立具有平滑曲线的图像中的至少一种来改善图像的空间分辨率。正投影的建立包括通过确定一系列投影值来建立正投影。图像的建立包括通过确定所述一系列投影值的逆来重建图像。
本发明包括以下内容：
实施方式1.一种在断层摄影成像系统中重建图像的方法，所述方法包括通过迭代重建图像改善图像的空间分辨率。
实施方式2.实施方式1的方法，其中所述断层摄影成像系统包括透射计算机断层摄影系统和发射计算机断层摄影系统之一。
实施方式3.实施方式1的方法，其中所述迭代重建图像包括生成分辨率高于典型分辨率的图像。
实施方式4.实施方式1的方法，其中所述迭代重建图像包括根据x射线多切片螺旋数据迭代重建图像，其中所述图像的空间分辨率高于典型分辨率。
实施方式5.实施方式1的方法，其中所述迭代重建图像包括应用正则化函数。
实施方式6.实施方式1的方法，其中所述空间分辨率包括各向同性空间分辨率。
实施方式7.实施方式1的方法，还包括生成正模型，所述正模型表现出小于典型分辨率的图像体素尺寸，其中所述图像体素尺寸包括图像的图像体素的尺寸。
实施方式8.实施方式1的方法，其中所述图像包括图像体素，所述方法还包括模型化所述图像体素，所述图像体素具有平行于第一坐标轴的第一维度和平行于第二坐标轴的第二维度，其中所述第一维度不同于所述第二维度，并且其中所述第一坐标轴不同于所述第二坐标轴。
实施方式9.实施方式1的方法，其中所述图像包括图像体素，并且其中所述图像体素为立方形或除立方形以外的形状。
实施方式10.实施方式1的方法，其中所述图像包括多个图像体素，其中所述多个图像体素中的第一个和第二个之间的间距不同于所述多个图像体素中第一个的尺寸。
实施方式11.实施方式1的方法，还包括由所述图像的图像体素生成正模型，其中所述正模型结合下述中的至少一种：
所述图像体素的几何形状、位置和取向；
所述正模型内包含的多个探测元的几何形状、位置和取向；
用于生成所述正模型的焦斑的几何形状、位置和取向；
时间采样，其取决于扫描的触发频率；和
成像系统的平台的运动。
实施方式12.实施方式1的方法，其中所述图像包括图像体素，所述方法还包括由所述图像体素生成正模型，其中所述正模型包括平滑函数。
实施方式13.实施方式1的方法，其中所述图像包括图像体素，所述方法还包括：
由所述图像体素生成正模型；和
通过调整所述正模型使断层摄影成像系统的组元表现为大于该组元的物理尺寸。
实施方式14.实施方式1的方法，其中所述图像包括图像体素，所述方法还包括：
由所述图像体素生成正模型；和
通过调整所述正模型使断层摄影成像系统的组元表现为小于该组元的物理尺寸。
实施方式15.实施方式1的方法，其中所述图像包括图像体素，所述方法还包括由所述图像体素生成正模型，其中所述正模型的置信度取决于所述图像体素的位置。
实施方式16.一种权衡断层摄影成像系统的性能的迭代重建方法，所述方法包括调整一部分迭代重建图像的噪声和分辨率。
实施方式17.一种调整断层摄影成像系统的分辨率的迭代重建方法，所述方法包括：
调整迭代重建图像的面内分辨率；和
调整迭代重建图像的断面分辨率。
实施方式18.实施方式17的方法，其中所述面内分辨率包括径向分辨率和方位角分辨率。
实施方式19.实施方式17的方法，其中所述面内分辨率等于所述断面分辨率。
实施方式20.一种改善图像的空间分辨率的方法，所述方法包括：
迭代重建所述图像；和
通过下述中的至少一种改善所述图像的空间分辨率：
建立具有平滑曲线的正投影函数，其中所述建立正投影包括通过确定一组投影值重建正投影；和
建立具有光滑曲线的图像，其中所述建立图像包括通过确定所述一组投影值的逆重建图像。
实施方式21.实施方式20的方法，其中所述改善空间分辨率包括沿探测器信道轴和探测器行轴中的至少一个调整所述图像的图像体素的尺寸。
实施方式22.实施方式20的方法，其中所述改善空间分辨率包括沿探测器信道轴和探测器行轴中的至少一个调整由图像体素正投影的探测元的数目，所述图像体素是图像的一部分并以一定的投影角形成。
实施方式23.实施方式20的方法，其中所述改善空间分辨率包括沿探测器信道轴和探测器行轴中的至少一个调整探测元的数目，其中所述探测元的中心为图像体素的正投影，所述图像体素是图像的一部分并以一定的投影角形成。
实施方式24.实施方式20的方法，其中所述改善空间分辨率包括生成图像体素的正投影，其中所述生成正投影包括将图像体素正投影在至少一个探测元上，其中所述正投影沿探测器信道轴和探测器行轴中至少一个的尺寸与所述图像体素的尺寸成正比。
实施方式25.实施方式20的方法，其中所述改善空间分辨率包括沿探测器信道轴和探测器行轴中的至少一个调整正投影区的尺寸，所述正投影区通过正投影图像的图像体素而生成。
实施方式26.实施方式20的方法，其中所述改善空间分辨率包括沿探测器信道轴和探测器行轴中的至少一个调整正投影区的形状，所述正投影区通过正投影图像的图像体素而生成。
实施方式27.实施方式20的方法，其中所述改善空间分辨率包括在图像所呈现的视场内实现均匀的空间分辨率，其中所述均匀的空间分辨率如下实现：改变多个正投影沿探测器信道轴和探测器行轴中至少一个的多个尺寸，所述多个正投影通过正投影图像的多个图像体素而形成。
实施方式28.实施方式20的方法，其中所述改善空间分辨率包括在图像所呈现的视场内实现均匀的空间分辨率，其中所述均匀的空间分辨率如下实现：改变多个正投影沿探测器信道轴和探测器行轴中至少一个的多个形状，所述多个正投影通过正投影图像的多个图像体素而形成。
实施方式29.一种改善图像空间分辨率的处理器，所述处理器经配置：
迭代重建图像；和
通过下述中的至少一种改善所述图像的空间分辨率：
建立具有平滑曲线的正投影函数，其中所述建立正投影包括通过确定一组投影值建立正投影；和
建立具有光滑曲线的图像，其中所述建立图像包括通过确定所述一组投影值的逆重建图像。
实施方式30.实施方式29的处理器，其中配置所述处理器，以通过沿探测器信道轴和探测器行轴中的至少一个调整所述图像的图像体素的尺寸，来改善空间分辨率。
本申请还包括以下内容：
实施方案1.使用断层摄影成像系统重建图像的方法，所述方法包括：
生成图像第一估计值；
在生成图像第一估计值之后，生成图像第二估计值；
比较实际投影数据和图像第二估计值的预期投影数据；
使用所述比较调整所述第一估计值；
生成对应于断层摄影成像系统的真实X射线源的虚拟X射线源；
生成包括多个对应于所述断层摄影成像系统的真实X射线探测器阵列的虚拟探测元的虚拟X射线探测面，所述真实X射线探测器阵列包括多个虚拟探测元所对应于的多个探测器单元；
使用虚拟X射线源和虚拟X射线探测器确定包括表示图像体素的X射线密度的值的图像体数据；
赋予实际投影数据的第一部分和第一正投影函数第一权重，其中实际投影数据的第一部分和第一正投影函数是当虚拟X射线源处于第一位置时所生成的并用于重建图像体素；
赋予实际投影数据的第二部分和第二正投影函数第二权重，其中实际投影数据的第二部分和第二正投影函数是当虚拟X射线源处于第二位置时所生成的并用于重建图像体素，所述第二权重小于所述第一权重，且所述第一位置比所述第二位置更接近所述图像体素；
使用图像体数据调整所述第一估计值；和
基于图像第一估计值输出图像，所述输出图像的分辨率高于所述第一估计值的分辨率。
实施方案2.实施方案1的方法，还包括使用一系列图像的后来的估计值迭代调整所述第一估计值。
实施方案3.实施方案2的方法，其中反复调整所述第一估计值还包括用现有的估计值替代先前的估计值以确定之后的估计值。
实施方案4.实施方案1的方法，其中使用图像体数据调整第一估计值还包括使用图像体数据的函数调整第一估计值，其中所述函数使用实际投影数据，正投影函数，正则化函数，和失真量度中的至少一个来调整所述第一估计值。
实施方案5.实施方案4的方法，其中使用图像体数据的函数还包括使用
Figure BDA00003339665000071
作为函数，其中
Figure BDA00003339665000072
是调整的第一估计值，r是图像体数据，D是失真量度，p是实际投影数据，F(r)是正投影函数，以及U(r)是正则化函数。
实施方案6.实施方案5的方法，还包括，当一系列估计值用于调整所述第一估计值时，用来自前一次迭代的
Figure BDA00003339665000073
替代本次迭代中的r以确定之后的
Figure BDA00003339665000074
实施方案7.实施方案4的方法，还包括使用失真量度，该失真量度为实际投影数据和图像第一估计值的正投影函数之间失配的量度。
实施方案8.实施方案4的方法，还包括使用失真量度，该失真量度是非负凸函数、概率密度函数的负对数和加权平方函数之一。
实施方案9.实施方案4的方法，还包括使用正则化函数，该正则化函数根据相邻图像体素中至少一个的值补偿图像体素的值。
实施方案10.实施方案4的方法，还包括使用正则化函数，该正则化函数根据与该图像体素不相邻并与该图像体素相隔预定间距的至少一个图像体素的值补偿该图像体素的值。
实施方案11.实施方案4的方法，还包括调整邻近图像体素并补偿图像体素的值的图像体素的数目以达到以下的至少一种：相对于图像体素的第一空间分辨率提高图像体素的空间分辨率，在图像的视场内获得均匀的空间分辨率，在图像的视场内获得各向同性的空间分辨率，和在图像的视场内获得预定的空间分辨率。
实施方案12.实施方案1的方法，其中使用图像体数据调整所述第一估计值还包括对所述图像体素应用正投影函数以生成包括至少一个虚拟探测面的虚拟探测元的正投影区。
实施方案13.实施方案12的方法，其中对图像体素应用正投影函数还包括沿虚拟探测面的虚拟探测器行和虚拟探测面的虚拟探测器信道中至少一个调整图像体素相对于正投影区的尺寸。
实施方案14.实施方案13的方法，其中调整图像体素的尺寸还包括调整图像体素的尺寸以提高图像体素的空间分辨率，这是与所述调整第一估计值的前一次迭代过程中确定的图像体素的之前空间分辨率相比而言的。
实施方案15.实施方案1的方法，其中使用图像体数据调整所述第一估计值还包括通过从虚拟X射线源以一定的投影角将图像体素正投影在虚拟探测面上，确定贡献于图像体素的分辨率的一组虚拟探测元。
实施方案16.实施方案15的方法，其中确定一组虚拟探测元还包括：
将图像体素正投影在虚拟探测面上以在虚拟探测面上生成图像体素的阴影；和
在生成的阴影内选择一组虚拟探测元作为确定组的虚拟探测元。
实施方案17.实施方案15的方法，其中确定一组虚拟探测元还包括：
通过使多条射线从虚拟X射线源延伸到虚拟探测面上以在虚拟探测面上产生正投影表面，将图像体素正投影在虚拟探测面上；和
在正投影表面内选择虚拟探测元作为确定组的虚拟探测元。
实施方案18.实施方案1的方法，其中使用图像体数据调整所述第一估计值还包括将虚拟探测元之一的尺寸调整为以下之一：大于和小于真实X射线探测器阵列的探测器单元的相应尺寸。
实施方案19.实施方案1的方法，其中使用图像体数据调整所述第一估计值还包括确定在虚拟探测面内并且用于重建图像体数据的一组虚拟探测元的多个值。
实施方案20.实施方案19的方法，其中确定一组虚拟探测元的多个值还包括给虚拟探测面的预定部分内的每个虚拟探测元赋值。
实施方案21.实施方案19的方法，其中确定一组虚拟探测元的多个值还包括根据沿虚拟探测器行轴和虚拟探测器信道轴之一的每个虚拟探测元离虚拟探测面上预定点的距离给每个虚拟探测元赋值。
实施方案22.实施方案19的方法，其中确定一组虚拟探测元的多个值还包括在虚拟探测面的预定部分内根据预定曲线函数给每个虚拟探测元赋值。
实施方案23.实施方案19的方法，其中确定一组虚拟探测元的多个值还包括在虚拟探测面内根据多个预定条件给每个虚拟探测元赋值。
实施方案24.实施方案1的方法，其中使用图像体数据调整所述第一估计值还包括在虚拟探测面内确定多个虚拟探测元以在图像的视场内获得均匀的空间分辨率。
实施方案25.实施方案1的方法，其中赋予第一权重还包括使用第一失真量度以赋予第一权重，并且权重赋予第二权重还包括使用第二失真量度以赋予第二权重，所述第二失真量度不同于所述第一失真量度。
实施方案26.断层摄影成像系统，其包括：
X射线源，其配置为生成射线束；
探测器阵列，其配置为接收射线束；和
与所述X射线源和所述探测器阵列相通的图像重建器用于接收图像数据，所述图像重建器配置为：
生成图像第一估计值；
在生成图像第一估计值之后，生成图像第二估计值；
比较实际投影数据和图像第二估计值的预期投影数据；
使用所述比较调整所述第一估计值；
生成对应于断层摄影成像系统的真实X射线源的虚拟X射线源；
生成包括多个对应于所述断层摄影成像系统的真实X射线探测器阵列的虚拟探测元的虚拟X射线探测面，所述真实X射线探测器阵列包括多个虚拟探测元所对应于的多个探测器单元；
使用虚拟X射线源和虚拟X射线探测器确定包括表示图像体素的X射线密度的值的图像体数据；和
如下调整第一估计值：
对正投影函数应用图像体素以生成包括虚拟探测面的至少一个虚拟探测元的正投影区；
通过从虚拟X射线源以一定的投影角将图像体素正投影在虚拟探测面上，确定贡献于图像体素的分辨率的一组虚拟探测元；
将虚拟探测元之一的尺寸调整为以下之一：大于和小于真实X射线探测器阵列的探测器单元的相应尺寸；
确定在虚拟探测面内且用于重建图像体数据的一组虚拟探测元的多个值；
赋予实际投影数据的第一部分和第一正投影函数第一权重，其中实际投影数据的第一部分和第一正投影函数是当虚拟X射线源处于第一位置时所生成的并用于重建图像体素；和
赋予实际投影数据的第二部分和第二正投影函数第二权重，其中实际投影数据的第二部分和第二正投影函数是当虚拟X射线源处于第二位置时所生成的并用于重建图像体素，所述第二权重小于所述第一权重，且所述第一位置比所述第二位置更接近图像体素；和
基于调节的图像第一估计值输出图像，所述输出图像的分辨率高于所述第一估计值的分辨率。
实施方案27.使用断层摄影成像系统重建图像的方法，所述方法包括：
在视场内生成图像第一估计值；
在生成图像第一估计值之后，在视场内生成图像第二估计值；
比较实际投影数据和图像第二估计值的预期投影数据；
使用第二估计值的比较调整所述第一估计值以在视场内生成图像第三估计值；
在调整第一估计值以生成第三估计值之后，在视场内生成图像第四估计值；
比较实际投影数据和图像第四估计值的预期投影数据；
使用第四估计值的比较调整所述第三估计值以在视场内生成图像第五估计值；和
基于图像第五估计值输出图像，所述输出图像的分辨率高于所述第一估计值的分辨率。
附图说明
图1是实现图像分辨率改善方法的多切片CT成像系统的实施方式的等视图。
图2是图1的CT成像系统的方框图。
图3是图像分辨率改善方法的实施方式的示意图。
图4是图3的方法的示意图。
图5是图像分辨率改善方法的实施方式的流程图。
图6是图5的流程图的续图。
图7示出多个图像的实施方式，显示应用图像分辨率改善方法的实施方式的效果。
具体实施方式
迭代重建包括如下形成图像的方法：根据实际投影数据和预期投影数据之间匹配的质量基于图像的先验估计值(current estimate)反复调整现有的估计值。对图像多个特性的考虑，例如图像的平滑度和/或预建模型的满意度，也可能影响匹配质量。进行多次迭代，从而形成基于预定标准与实际投影数据匹配最佳的图像。由于使一组重建图像形成患者的三维重现，并且各图像像素或像点对应于三维重建中的单个体素或体元(volume element)，因而将整组重建图像称为三维重建。
参考图1和图2，示出了多切片透射计算机断层摄影(CT)成像系统10的实施方式，该系统应用重建受检体12(如医学上的患者或影像)的解剖部位图像的方法。CT成像系统10包括台架14，台架14具有包括x射线源18和探测器阵列20的旋转内部16。x射线源18和探测器阵列20随台架14的旋转而转动。x射线源18朝向探测器阵列20投射x射线束32。x射线源18和探测器阵列20围绕置于可操作平移的平台22上的受检体12旋转。平台22沿平行于z方向的z轴在源18和探测器阵列20之间平移来进行解剖部位的螺旋扫描，或者在解剖部位的整个轴向扫描过程中停留在沿z轴的同一位置。在探测器阵列20处探测透过患者孔24内受检体12之后的x射线束32，从而建立用于形成解剖部位CT图像的模拟投影数据。
x射线源18和探测器阵列20围绕平行于z轴的中心轴30旋转。束32由探测器阵列20的多个探测器行中的多个探测器单元34接收。探测器阵列20包括多行探测器单元34和多道探测器单元34。探测器信道(detectorchannel)平行于信道轴，信道轴平行于台架24的平面。探测器行平行于行轴(row axis)，行轴平行于z轴。各探测器行沿着台架24旋转所围绕的z轴在z方向上偏离所有其它探测器行。探测器单元34生成模拟投影数据，所述模拟投影数据表现为相应于束32强度的电信号。束32由于穿过受检体12而衰减。台架14的旋转和源18的操作由控制机构36控制。控制机构36包括向x射线源18提供动力和定时信号的x射线控制器38以及控制台架14的转速和位置的台架电动机控制器40。数据采集系统(DAS)42从探测器单元34采集模拟投影数据并将模拟投影数据从模拟形式转换为数字信号，从而生成采样数字化投影数据，即实际投影数据。图像重建器44从DAS 42接收实际投影数据并进行图像重建，例如改善图像分辨率的方法，从而生成CT图像。主控制器46将CT图像存储在大容量存储装置48中。大容量存储装置48的实例包括非易失性存储器如只读存储器(ROM)和易失性存储器如随机存取存储器(RAM)。大容量存储装置48的其它实例包括软盘、紧致盘——ROM(CD-ROM)、磁光盘(MOD)和数字多功能盘(DVD)。
主控制器46还经由操作员控制台50接收来自操作员的指令和扫描参数。监视器52使操作员能够观察到CT图像和来自主控制器46的其它数据。监视器52可以是阴极射线管(CRT)或液晶显示器(LCD)。主控制器46将操作员发出的指令和参数用于DAS 42、x射线控制器38和台架电动机控制器40的操控。另外，主控制器46操控平台电动机控制器54，平台电动机控制器54平移平台22，从而将解剖部位置于台架14中。
x射线控制器38、台架电动机控制器40、图像重建器44、主控制器46和平台电动机控制器54均不仅仅限于本领域称为控制器的集成电路，而是广泛地表示计算机、处理器、微控制器、微型计算机、可编程逻辑控制器、专用集成电路和/或任意其它可编程电路。x射线控制器38、台架电动机控制器40、图像重建器44、主控制器46和平台电动机控制器54可以是中央控制器的一部分或者各自为如图所示的独立部件。
尽管上述具体实施方式是指第三代CT成像系统10，但迭代重建图像的方法等同地应用于具有固定探测器和旋转x射线源的第四代CT系统、具有固定探测器和电子束偏转x射线源的第五代CT系统、包括多个x射线源和/或探测器的未来各代CT系统、以及发射CT系统如单光子发射CT系统(SPECT)或正电子发射断层摄影系统(PET)。
另外，尽管就医疗设施对迭代重建图像的方法进行描述，但可以预期所述方法的多个技术效果可由非医疗成像系统实现，如通常应用于工业设施或运输设施的系统，例如但不限于，非破坏性检测系统，用于机场、其它运输中心、政府大楼和办公大楼的行李扫描系统。所述技术效果还可由微PET和CT系统实现，裁定所述系统的尺寸以研究实验用动物体而非人体。
图3和4为示意图，图5和6为流程图，示出了迭代重建图像的方法的实施方式。源302示例x射线源18的虚拟像，探测面304示例探测器阵列20面的虚拟像。探测器阵列20面朝向x射线源18。探测面304包括多个探测元306、308、310、312和314，探测面304的各探测元是探测器单元34的虚拟像。图像重建器44根据实验数据例如线影的变化获取x射线源18的位置，并生成x射线源18的虚拟像。根据实验数据获取源302的位置变化。例如，通过在x射线源18和探测器阵列20之间放置线来确定源302的位置。线相对x射线源18是固定的，相对探测器阵列20是不固定的。形成在探测器阵列20上的线影的位置变化提供了源302的位置变化。
图像重建器44还从位置编码器获取探测器阵列20的位置，位置编码器相对于x射线源18探测位置。可根据x射线源18和探测器阵列20的位置(相对图1所示的xyz坐标系)，通过图像重建器44生成x射线源18和探测器阵列20的虚拟像。源302的位置改变以配合x射线源18的位置变化，探测面304的位置改变以配合探测器阵列20面的位置变化。源302和探测面304之间的距离与x射线源18和探测器阵列20面之间的距离成比例。
对于特定的视角，迭代重建图像方法的优化问题表达为：
r ^ ( r ) = arg min r { D ( p - F ( r ) ) + U ( r ) } - - - ( 1 )
其中r代表图像体318，该图像体318通过应用改善图像分辨率的方法而重建并代表受检体12的解剖部位。图像体318是图像体r的实例。图像体318包括表示x射线密度的值，例如以Hounsfield为单位测量的图像体素320的CT值。尽管图像体素320为立方形，但在替换性实施方式中，图像体318内的图像体素具有其它形状，例如球形、椭球形、柱形、平行六面体形、梯形或多面体形。在替换性实施方式中，图像体318包括多个表示代表解剖部位的图像体素320和多个图像体素322、324、326、330、332的x射线密度的值。在替换性实施方式中，沿x、y和z轴之一平行于图像体素322设置图像体素320。在另一替换性实施方式中，可将图像体素322设置在图像体318周围。在等式(1)中，p表示实际投影数据。在等式(1)中，F(r)表示探测面304上r的正投影函数或正投影模型。正投影F(r)以模拟CT成像系统10的方式转换r。
在等式(1)中，U(r)为补偿局部图像体素差的正则化函数。U(r)的实例包括凸函数，例如图像体素320的CT值和与其相邻的图像体素328的CT值之差的平方。
另外，在等式(1)中，D是实际投影数据p和图像估计值的正投影函数F(r)之间失配的失真量度(distortion measure)。失真量度D的实例包括非负凸函数、概率密度函数的负对数和其它补偿函数。再例如，在等式(2)中失真量度D以加权平方形式表示为：
r ^ ( r ) = arg min r { ( p - F ( r ) ) T W ( p - F ( r ) ) + U ( r ) } - - - ( 2 )
其中W为加权函数，T为转置运算。加权函数的实例包括具有多个对角元的对角矩阵，所述多个对角元为实际投影数据的逆。在替换性实施方式中，等式(1)和(2)中任意一个均不包括U(r)。在另一替换性实施方式中，等式(2)不包括加权函数W并包括(p-F(r))而不是W(p-F(r))。
r和p之间关系的实例包括线性关系，例如p=Ar+n，其中A为探测面304的探测元贡献值的矩阵，n表示一组至少一个噪声值，该噪声值代表p围绕实际投射数据平均值的随机涨落，例如方差。在等式(1)的每次迭代中，图像重建器44计算r的微扰，随着每次迭代r的微扰使等式(1)中
Figure BDA00003339665000142
的值减小。图像重建器44根据计算正投影函数F(r)和正则化函数U(r)之前的前一次迭代中计算得到的
Figure BDA00003339665000143
值，在本次迭代中计算正投影函数F(r)和正则化函数U(r)。在本次迭代中，图像重建器44以代替等式(1)中的r。在前一次迭代中，图像重建器44由根据实际投影数据p滤波反投射重建的图像体318生成
Figure BDA00003339665000145
在替换性实施方式中，在前一次迭代中，图像重建器44可由具有操作员经由操作员控制台50选择的值的图像体318生成的质量强烈依赖于正投影函数F(r)反映CT成像系统10的物理现实的程度。
图像重建器44将正投影函数F(r)应用于图像体素320，以生成包括探测面304的探测元308的正投影区。例如，由尺寸例如面积小于典型分辨率的正投影图像体素320形成正投影区。在替换性实施方式中，正投影区包括多个探测元，例如探测元306、308、310和312中的至少两个。在一种实施方式中，图像重建器44沿探测器行轴329和探测器信道轴331中的至少一个调整502例如增大或者减小图像体素320的尺寸。探测器行轴329平行于z轴，探测器信道轴331平行于台架14平面。探测元如探测元308和310平行于探测器信道轴331。探测元如探测元306和308平行于探测器行轴329。图像重现器44通过沿x轴和y轴中的至少一个调整图像体素320的尺寸，沿探测器信道轴331调整图像体素320的尺寸。图像重建器44沿探测器行轴329和探测器信道轴331中的至少一个相对于正投影区域的尺寸调整例如减小图像体素320的尺寸。例如，图像重建器44将沿探测器行轴329测量的图像体素320的深度减小为沿探测器行轴329测量的正投影区的长度的一半。再例如，图像重建器44将沿探测器信道轴331测量的图像体素320的宽度减小为沿探测器信道轴331测量的正投影区的宽度的三分之一。
图像重建器44在本次迭代中调整图像体素320的尺寸，从而相对于等式(1)的前一次迭代中图像体素320的分辨率以及相对于常规重建技术如滤波反投射，改善图像体素320的空间分辨率，例如断面分辨率和面内分辨率中的至少一个。例如，图像重建器44沿探测器行轴329调整图像体素320的尺寸，以改善图像体素320的断面分辨率。再例如，图像重建器44沿探测器信道轴331调整图像体素320的尺寸，以改善图像体素320的面内分辨率。再例如，图像重建器44沿探测器行轴329和探测器信道轴331中的至少一个调整图像体素320的尺寸，使得图像体素320的面内分辨率等于图像体素320的断面分辨率。面内分辨率是在平面内沿径向方向以及在平面内成一定方位角的平面内分辨率。图像重建器44调整图像体素320的尺寸，使得图像体素320的所有尺寸不同于图像体素320和332之间的距离。
在替换性实施方式中，图像重建器44通过从源302以一定的投影角将图像体素320正投影在探测面304上，来确定或绘出贡献于图像体素320的分辨率的大量探测元。例如，图像重建器44将图像体素320的中心321正投影在探测面304上以生成正投影点334，并沿探测器行轴329和探测器信道轴331中的至少一个在探测面304内选择预定数目例如两个或三个探测元。在一个实例中，经由操作员控制器50获取操作员发出的预定数目。在一个实例中，所选择的探测元与包括正投影点334的探测元308相邻。在一个实例中，替换性地，所选择的探测元不与包括正投影点334的探测元308相邻，而是通过探测面304的大量探测元与探测元308相隔。再例如，图像重建器44将图像体素320的中心321正投影在探测面304上以生成正投影点334，并选择位于正投影区内的大量探测元，所述探测元沿探测器信道轴331和探测器行轴329中的至少一个设置并且尺寸与图像体素320的尺寸成正比。在一个实例中，图像重建器44所选择的并且与图像体素320的尺寸成正比的探测元的尺寸包括沿探测器行轴329的探测元深度。在一个实例中，该深度为图像体素320沿探测器行轴329的深度的两倍或三倍。正投影区内探测元数目的变化改变正投影区的尺寸和正投影区的形状中的至少一个。再例如，图像重建器44将图像体素320正投影在探测面304上以在探测面304上形成图像体素320的阴影，并沿探测器行轴329和探测器信道轴331中的至少一个选择落在图像体素320的阴影内的大量探测元。图像体素320的阴影沿探测器行轴329和探测器信道轴331具有与图像体素320底面336的尺寸相同的尺寸。底面336面向探测面304。或者，图像体素320沿探测器行轴329的阴影尺寸至少等于底面336沿探测器行轴329的尺寸。在另一替换性实施方式中，图像体素320的阴影尺寸至少等于底面336沿探测器信道轴331的尺寸。在一个实例中，图像重建器44选择位于探测面304内并落在图像体素320的阴影内的探测元。再例如，通过使多条射线从源302经由图像体素320的多条边339、341、343和345延伸到探测面304上，图像重建器44将图像体素320正投影在投影面304上，从而在探测面304上生成正投影表面。在一个实例中，图像重建器沿探测器行轴329和探测器信道轴331中的至少一个选择落在探测面304上正投影表面内的大量探测元。再例如，通过使多条射线从源302经由图像体素320的边339、341、343和345上的多个点延伸到探测面304上，图像重建器44将图像体素320正投影在探测面304上，从而在探测面304上生成正投影带。在一个实例中，图像重建器44沿探测器行轴329和探测器信道轴331中的至少一个选择落在探测面304上正投影带内的大量探测元。正投影点334是正投影区的中心。在替换性实施方式中，正投影点334不是正投影区的中心。
图像重建器44通过选择探测面304的探测元，选择沿探测器行轴329的多列探测元和沿探测器信道轴331的多道探测元中至少之一。图像重建器44将探测元308的尺寸调整为大于探测器单元34的尺寸。例如，图像重建器44将探测元308的体积增大为大于探测器单元34的体积。或者，图像重建器44将探测元308的尺寸调整为小于探测器单元34的尺寸。例如，图像重建器44将探测元308的体积调整为小于探测器单元34的体积。图像重建器44应用所选择的并且贡献于图像体素320的分辨率的大量探测元，重建等式(1)中的并相对于滤波反投影生成的典型分辨率改善空间分辨率
Figure BDA00003339665000171
在另一替换性实施方式中，图像重建器44确定落在探测面304内并用于重建图像体318的探测元的多个值。例如，图像重建器44为探测面304的正投影区、正投影表面和正投影带内各探测元随机赋值例如正值或负值。再例如，图像重建器44根据图像体素320的中心321和源302之间的距离333以及源302和图像体素322的中心337之间的距离335，为正投影区内各探测元和正投影面内各探测元赋值。当距离333大于距离335时，图像重建器44为图像体素322的正投影面内大量探测元赋值，所述探测元的数目大于图像体素320的正投影区内探测元的数目。为图像体素322的正投影面内数目较大的探测元赋值使得图像体素322的分辨率大于图像体素320的分辨率。图像体素322的分辨率提高使得图像体素322内的噪声增加，图像体素320的分辨率下降使得图像体素320内的噪声降低。
或者，当距离335大于距离333时，图像重建器44为图像体素320的正投影区内大量探测元赋值，所述探测元的数目大于图像体素322的正投影面内探测元的数目。为图像体素320的正投影区内数目较大的探测元赋值使得图像体素320的分辨率大于图像体素322的分辨率。图像体素320的分辨率提高使得图像体素320内的噪声增加，图像体素322的分辨率下降使得图像体素322内的噪声降低。
图像重建器44通过将射线从源302经由图像体素322正投影在探测面304上而在探测面304上形成正投影面。图像重建器44通过应用改善图像分辨率的系统和方法从正投影面迭代重建图像体素322。例如，在改变正投影面内探测元的数目时图像重建器44从正投影面迭代重建图像体素322，以改善图像体素322的分辨率。正投影面与正投影区相分离并且不包括正投影区。源302和中心321之间的距离决定图像体素320在x射线源18的视场内的位置。视场包括图像体318。另外，源302和图像体素322的中心337之间的距离335决定图像体素322在视场内的位置。
图像重建器44继续确定探测面304内用于重建图像体318的多个探测元的值。例如，图像重建器44如下确定正投影区内探测元的值：根据探测元与正投影点334沿探测器行轴329和探测器信道轴331中至少一个的距离，为探测元赋值。例如，当图像重建器44确定探测元306和正投影点334之间沿探测器行轴329的距离大于探测元308和正投影点334之间沿探测器行轴329的距离时，图像重建器44赋予探测元306的值小于赋予探测元308的值。再例如，当图像重建器44确定探测元310和正投影点334之间沿探测器信道轴331的距离大于探测元308和正投影点334之间沿探测器信道轴331的距离时，图像重建器44赋予探测元310的值小于赋予探测元308的值。
图像重建器44进一步继续确定探测面304内用于重建图像体318的多个探测元的值。例如，根据图像体素320的中心321和正投影区的正投影点334之间的距离338以及图像体素322的中心337和正投影面的正投影点342之间的距离340，图像重建器44为正投影区内各探测元和正投影面内各探测元赋值。当距离338大于距离340时，图像重建器44赋予正投影区内探测元的值大于赋予正投影面内探测元的值。或者，当距离340大于距离338时，图像重建器44赋予正投影面内探测元的值大于赋予正投影区内探测元的值。正投影点342为正投影面的中心。在替换性实施方式中，正投影点342不是正投影面的中心。中心321和正投影点334之间的距离338决定中心321的位置，中心337和正投影面的正投影点342之间的距离340决定中心337的位置。再例如，图像重建器44根据预定曲线为正投影区内多个探测元赋值。图像重建器44经由操作员控制台50接收来自操作员的预定曲线。预定曲线可以是下述中至少一个的函数：正投影区内探测元的数目，正投影区沿探测器信道轴331和探测器行轴329中至少一个的尺寸，图像体素320沿x、y和z维度中至少一个的尺寸。x维度平行于x轴，y维度平行于y轴，z维度平行于z轴。y轴垂直于z轴，x轴同时垂直于y轴和z轴。图像重建器44赋予探测面304内与预定曲线重合的探测元预定曲线上的值。
在替换性实施方式中，图像重建器44为探测面304内的探测元赋值以满足至少两个条件。第一个条件是正投影函数F(r)为不包括高频分量如脉冲的平滑曲线。高频分量的另一实例包括不符合正投影函数F(r)通常预期的多个缓慢变化的多个快速局部变化。平滑曲线是预定曲线的一种实例。正投影函数F(r)的实例为Ar+n。当图像重建器44设定r=1且n=0时，图像重建器44确定A的值以生成平滑曲线形式的正投影函数F(r)。第二个条件是图像体318为平滑曲线。当图像重建器44设定n=0时，图像重建器44确定A-1的值以生成不包括高频分量的平滑曲线形式的图像体318。
在另一替换性实施方式中，图像重建器44确定506探测面304内探测元的数量，以在视场内沿探测器信道轴331和探测器行轴329中的至少一个获得均匀的空间分辨率。例如，图像重建器44确定正投影区内探测元的数目等于正投影面内探测元的数目。再例如，图像重建器44确定正投影区内沿探测器行轴329的探测元数目等于正投影面内沿探测器行轴329的探测元数目。再例如，图像重建器44确定正投影区内沿探测器信道轴331的探测元数目等于正投影面内沿探测器行轴329的探测元数目。图像重建器44经由操作员控制台50接收操作员发出的均匀空间分辨率的请求。
在再一替换性实施方式中，图像重建器44根据CT成像系统10进行扫描的触发频率改变正投影区内探测元的数目。触发频率决定x射线源18围绕受检体12旋转一周在受检体12周围采集的投影数据视角的离散值。例如，当采集用于重建图像体素320的实际投影数据的触发频率增大时，正投影区内所选择的探测元的数目减少。再例如，当采集实际投影数据的触发频率减小时，正投影区内所选择的探测元的数目增加。
正则化函数U(r)根据邻近图像体素320的图像体素的值约束图像体素320。例如，正则化函数U(r)根据相邻图像体素324、326、328和332中至少一个的值补偿图像体素320的值。再例如，U(r)根据与图像体素320不相邻并与图像体素320相隔一定数目如一个或两个图像体素的图像体素的值补偿图像体素320的值。根据邻近图像体素320的图像体素的值补偿图像体素320的值的实例包括图像体素320的值不大于接近图像体素320的图像体素的平均值。根据邻近图像体素320的图像体素的值补偿图像体素320的值的另一实例包括图像体素320的值不大于邻近图像体素320的图像体素的中间值。
在再一替换性实施方式中，图像重建器44调整508沿探测器行轴329和探测器信道轴331中至少一个邻近图像体素320并补偿图像体素320的值的图像体素的数目，以改善图像体素320的空间分辨率。例如，图像重建器44调整沿探测器行轴329的图像体素的数目以补偿图像体素320的值。再例如，图像重建器44调整沿探测器信道轴331的图像体素的数目以补偿图像体素320的值。再例如，图像重建器44根据下述函数中的至少一个改变U(r)：视场内图像体素320的位置的函数以及沿x轴、y轴和z轴中至少一个测量的图像体素320的尺寸的函数。再例如，图像重建器44根据下述函数改变U(r)：图像体素320和与其邻近的至少一个图像体素之间尺寸大小之差的函数。例如，
Figure BDA00003339665000202
U ( r ) = a Σ c b ij | r i - r j | p - - - ( 3 )
其中α为度量常数如正实数，p为指数参数如正实数，bij为定向权重系数如正实数，ri表示图像体318内位置i处的图像体素320的值，rj表示图像体318内邻近图像体素320的图像体素j的值，C是选择用于补偿图像体素320的邻近图像体素320的图像体素的邻域。例如，图像重建器44根据视场内图像体素320的位置改变α和bij。
在一种实施方式中，图像重建器44调整邻近图像体素320并补偿图像体素320的值的图像体素的数目，并调整邻近图像体素322并补偿图像体素322的值的图像体素的数目，以沿视场内径向方向、视场内方位角方向和视场内z方向中的至少一个获得均匀的空间分辨率。例如，图像重建器44确定邻近图像体素320并补偿图像体素320的值的图像体素的数目等于邻近图像体素322并补偿图像体素322的值的图像体素的数目。在替换性实施方式中，图像重建器44调整邻近图像体素320并补偿图像体素320的值的图像体素的数目，并调整邻近图像体素322并补偿图像体素322的值的图像体素的数目，以沿视场内径向方向、方位角方向和z方向中至少两个获得各向同性的空间分辨率。图像重建器44经由操作员控制台50接收操作员发出的沿视场内径向方向、视场内方位角方向、视场内z方向中至少两个的多个所需的空间分辨率值。在另一替换性实施方式中，图像重建器44调整邻近图像体素320并补偿图像体素320的值的图像体素的数目，并调整邻近图像体素322并补偿图像体素322的值的图像体素的数目，以沿视场内径向方向、方位角方向和z方向中的至少一个获得所需的空间分辨率。图像重建器44经由操作员控制台50接收操作员发出的所需的空间分辨率。
图像重建器44还可改变正投影区内探测元的数目和接近图像体素320并补偿图像体素320的值的图像体素的数目，以沿视场内径向方向、视场内方位角方向和视场内z方向中至少一个方向获得空间分辨率，例如均匀的空间分辨率、所需的分辨率和各向同性的分辨率。
在一种实施方式中，在对解剖部位进行半扫描或全扫描过程中，图像重建器44赋予510实际投影数据的一部分和正投影函数F(r)大的置信度(confidence)或权重，其中所述实际投影数据的一部分和正投影函数F(r)是源302处于第一位置时所生成的并用于迭代重建图像体素320。所述大的置信度大于图像重建器44赋予实际投影数据的一部分和正投影函数F(r)的置信度，其中所述实际投影数据的一部分和正投影函数F(r)是在源302处于第二位置时生成的并用于迭代重建图像体素320。例如，当确定在x射线源18处于0度投影角时由第一组实际投影数据迭代重建图像体素320以及在源302处于90度投影角时由第二组实际投影数据迭代重建图像体素320时，并且当确定源302处于0度投影角时比其处于90度投影角时更接近图像体素320时，图像重建器44赋予第一组实际投影数据的置信度高于其赋予第二组实际投影数据的置信度。第一位置比第二位置更接近图像体素320。图像体素320相对源302的位置随源302的投影角而改变。图像重建器44改变应用于实际投影数据p和正投影函数F(r)的失真量度D，以改变置信度。例如，图像重建器44将失真量度D从(p-F(r))的二次方提高到(p-F(r))的三次方，以增大置信度。
在替换性实施方式中，图像重建器44调整例如增加正投影区内探测元的数目，以改善由于采集用于迭代重建图像体素320的实际投影数据过程中焦斑尺寸的改变而劣化的图像体素320的分辨率。焦斑的尺寸包括x射线源18的面积，束32源于该面积。此外，在另一替换性实施方式中，图像重建器44沿探测器行轴329和探测器信道轴331中的至少一个移动正投影点334的位置，以改善由于采集用于迭代重建图像体素320的实际投影数据过程中焦斑位置的改变而劣化的图像体素320的分辨率。
在再一替换性实施方式中，图像重建器44沿z轴移动正投影点334的位置，以改善由于采集用于迭代重建图像体素320的实际投影数据过程中平台22沿z轴的位置测量误差而劣化的图像体素320的分辨率。在一种实施方式中，图像重建器44移动正投影点334的位置和/或调整例如增加正投影区内沿z轴的探测元数目，以改善由于平台22下沉而劣化的图像体素320的分辨率。
在替换性实施方式中，图像重建器44随着时间的变化和/或源302的投影角的变化沿探测器信道轴331移动正投影区的正投影点334的位置，以改善由于采集用于重建图像体素320的一部分实际投影数据过程中出现的台架14振动而劣化的图像体素320的分辨率。源302的投影角随x射线源18的位置的改变而改变。在另一替换性实施方式中，图像重建器44调整正投影区内探测元的数目和正投影区的正投影点334的位置中至少之一，以改善由于采集用于重建图像体素320的实际投影数据过程中出现的焦斑尺寸的改变、焦斑位置的改变、平台22沿z轴的位置改变、平台22的下沉、台架14的多个振动中的至少一种而劣化的图像体素320的分辨率。
图7示出了应用常规方法和改善图像分辨率的方法形成的多个图像602、604、606和608的实施方式。通过低分辨率探测装置进行的同一扫描重建图像602和608。通过高分辨率探测装置进行的两个独立扫描重建图像604和606。通过应用Feldkamp类重建重建图像602、604和606，通过应用改善图像分辨率的方法重建图像608。显然图像608的分辨率优于图像602的分辨率并且与图像604和606的分辨率相近。
尽管利用多切片CT成像系统10对改善图像分辨率的系统和方法进行了描述，但所述系统和方法能够用于多种成像系统，包括但不限于，磁共振成像(MRI)系统、光学扫描系统、CT系统、放射治疗系统、x射线成像系统、超声系统、核成像系统、磁共振波谱系统和正电子发射断层摄影(PET)成像系统。还应当指出的是，图像重建器44确定图像体318内任意图像体素相对于源302和探测面304的位置。还应当指出的是，在一种实施方式中，操作员经由操作员控制台50控制图像重建器44以改善图像的分辨率并确定分辨率是否改善。
迭代重建图像的系统和方法的技术效果包括改善图像的分辨率。其它技术效果包括通过实现均匀的空间分辨率降低图像体素320的空间分辨率的差异。通过应用迭代重建图像的系统和方法而迭代重建的图像的分辨率的提高使得图像的噪声降低。图像体318的空间分辨率的差异可以是台架14的等角点和图像体素320、322中任意一个之间距离的函数。例如，图像体素320的空间分辨率随图像体素320和322中任意一个与等角点之间距离的增大而降低。等角点和图像体素320、322中任意一个之间的距离由于束32的投影角、DAS 42的有限积分时间、heel效应和其它因素中的至少一个因素而改变。在有限积分时间内，DAS 42对以一定的投影角采集的实际投影数据的一部分进行积分。然而，实际上，实际投影数据中进行积分的部分不是以精确的投影角采集到的，因而，在图像体318内产生方位角模糊。图像重建器44通过根据DAS触发频率选择正投影区内多个探测元来减轻方位角模糊。
尽管针对不同的具体实施方式对本发明进行了描述，但本领域技术人员应当理解的是，可在不脱离权利要求的构思和范围进行改进的情况下实现本发明。

Claims (10)

1.使用断层摄影成像系统重建图像的方法，所述方法包括：
生成图像第一估计值；
在生成图像第一估计值之后，生成图像第二估计值；
比较实际投影数据和图像第二估计值的预期投影数据；
使用所述比较调整所述第一估计值；
生成对应于断层摄影成像系统的真实X射线源的虚拟X射线源；
生成包括多个对应于所述断层摄影成像系统的真实X射线探测器阵列的虚拟探测元的虚拟X射线探测面，所述真实X射线探测器阵列包括多个虚拟探测元所对应于的多个探测器单元；
使用虚拟X射线源和虚拟X射线探测器确定包括表示图像体素的X射线密度的值的图像体数据；
赋予实际投影数据的第一部分和第一正投影函数第一权重，其中实际投影数据的第一部分和第一正投影函数是当虚拟X射线源处于第一位置时所生成的并用于重建图像体素；
赋予实际投影数据的第二部分和第二正投影函数第二权重，其中实际投影数据的第二部分和第二正投影函数是当虚拟X射线源处于第二位置时所生成的并用于重建图像体素，所述第二权重小于所述第一权重，且所述第一位置比所述第二位置更接近所述图像体素；
使用图像体数据调整所述第一估计值；和
基于图像第一估计值输出图像，所述输出图像的分辨率高于所述第一估计值的分辨率。
2.权利要求1的方法，还包括使用一系列图像的后来的估计值迭代调整所述第一估计值。
3.权利要求2的方法，其中反复调整所述第一估计值还包括用现有的估计值替代先前的估计值以确定之后的估计值。
4.权利要求1的方法，其中使用图像体数据调整第一估计值还包括使用图像体数据的函数调整第一估计值，其中所述函数使用实际投影数据，正投影函数，正则化函数，和失真量度中的至少一个来调整所述第一估计值。
5.权利要求4的方法，其中使用图像体数据的函数还包括使用作为函数，其中
Figure FDA00003339664900022
是调整的第一估计值，r是图像体数据，D是失真量度，p是实际投影数据，F(r)是正投影函数，以及U(r)是正则化函数。
6.权利要求5的方法，还包括，当一系列估计值用于调整所述第一估计值时，用来自前一次迭代的
Figure FDA00003339664900023
替代本次迭代中的r以确定之后的
Figure FDA00003339664900024
7.权利要求4的方法，还包括使用失真量度，该失真量度为实际投影数据和图像第一估计值的正投影函数之间失配的量度。
8.权利要求4的方法，还包括使用失真量度，该失真量度是非负凸函数、概率密度函数的负对数和加权平方函数之一。
9.权利要求4的方法，还包括使用正则化函数，该正则化函数根据相邻图像体素中至少一个的值补偿图像体素的值。
10.权利要求4的方法，还包括使用正则化函数，该正则化函数根据与该图像体素不相邻并与该图像体素相隔预定间距的至少一个图像体素的值补偿该图像体素的值。
CN201310233382.8A 2006-06-22 2007-06-22 改善图像分辨率的系统和方法 Active CN103400401B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US11/472,572 US7583780B2 (en) 2006-06-22 2006-06-22 Systems and methods for improving a resolution of an image
US11/472,572 2006-06-22
CN2007800234587A CN101478919B (zh) 2006-06-22 2007-06-22 改善图像分辨率的系统和方法

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CN2007800234587A Division CN101478919B (zh) 2006-06-22 2007-06-22 改善图像分辨率的系统和方法

Publications (2)

Publication Number Publication Date
CN103400401A true CN103400401A (zh) 2013-11-20
CN103400401B CN103400401B (zh) 2017-04-12

Family

ID=38834424

Family Applications (2)

Application Number Title Priority Date Filing Date
CN2007800234587A Active CN101478919B (zh) 2006-06-22 2007-06-22 改善图像分辨率的系统和方法
CN201310233382.8A Active CN103400401B (zh) 2006-06-22 2007-06-22 改善图像分辨率的系统和方法

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN2007800234587A Active CN101478919B (zh) 2006-06-22 2007-06-22 改善图像分辨率的系统和方法

Country Status (5)

Country Link
US (1) US7583780B2 (zh)
JP (1) JP5198443B2 (zh)
CN (2) CN101478919B (zh)
DE (1) DE112007001451B4 (zh)
WO (1) WO2007150037A2 (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107492131A (zh) * 2017-07-01 2017-12-19 武汉斗鱼网络科技有限公司 用于安卓电视的倒影生成方法、存储介质、设备及系统
CN107495976A (zh) * 2016-06-14 2017-12-22 上海联影医疗科技有限公司 图像重建中的最大值和灰度值图像的获取方法及装置
CN108352077A (zh) * 2016-03-31 2018-07-31 上海联影医疗科技有限公司 图像重建的系统和方法
CN113450287A (zh) * 2021-07-27 2021-09-28 北京与光科技有限公司 高分辨率光谱恢复方法、装置和电子设备

Families Citing this family (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7949172B2 (en) * 2007-04-27 2011-05-24 Siemens Medical Solutions Usa, Inc. Iterative image processing
US8233682B2 (en) * 2007-06-05 2012-07-31 General Electric Company Methods and systems for improving spatial and temporal resolution of computed images of moving objects
KR20100133950A (ko) * 2008-01-30 2010-12-22 더 리전츠 오브 더 유니버시티 오브 캘리포니아 동적인 제약들에 따른 오브젝트 주변의 사용을 통한 단층 촬영에 있어서의 양 감소 및 이미지 강화
US9420974B2 (en) * 2008-06-06 2016-08-23 Koninklijke Philips N.V. Method and apparatus for attenuation correction
JP5537132B2 (ja) * 2008-12-11 2014-07-02 株式会社東芝 Ｘ線コンピュータ断層撮影装置、医用画像処理装置、及び医用画像処理プログラム
RU2543946C2 (ru) * 2009-07-14 2015-03-10 Конинклейке Филипс Электроникс Н.В. Реконструкция изображения, включающая в себя коррекцию изменяющегося со смещением размытия
US8655033B2 (en) * 2009-10-28 2014-02-18 General Electric Company Iterative reconstruction
JP5681924B2 (ja) 2009-11-27 2015-03-11 ジーイー センシング アンド インスペクション テクノロジーズ ゲ−エムベーハー サンプルの立体表現を決定するためのコンピュータ断層撮影法、コンピュータソフトウェア、計算装置およびコンピュータ断層撮影システム
US8731266B2 (en) * 2009-12-17 2014-05-20 General Electric Company Method and system for correcting artifacts in image reconstruction
US8447089B2 (en) * 2010-04-27 2013-05-21 Chunlei Liu Systems and methods for susceptibility tensor imaging
JP5675978B2 (ja) * 2010-07-09 2015-02-25 ジーイー センシング アンド インスペクション テクノロジーズ ゲ−エムベーハー コンピュータ断層撮影方法、コンピュータプログラム、コンピュータデバイスおよびコンピュータ断層撮影システム
DE102010041191B4 (de) 2010-09-22 2016-02-18 Siemens Aktiengesellschaft Erstellung von MR-Bilddaten mit paralleler Schichtanregung und teilweiser Überlappung der Schichten im Frequenzbereich
US9123098B2 (en) * 2011-04-28 2015-09-01 Hitachi Medical Corporation Medical image processing device and medical image processing method, applying weighted penalty term in iterative approximation
WO2013002805A1 (en) * 2011-06-30 2013-01-03 Analogic Corporation Iterative image reconstruction
US8483471B2 (en) * 2011-06-30 2013-07-09 General Electric Company Method and system for scatter correction in X-ray imaging
WO2013008702A1 (ja) * 2011-07-08 2013-01-17 株式会社 日立メディコ 画像再構成装置及び画像再構成方法
US8837797B2 (en) * 2012-01-10 2014-09-16 Kabushiki Kaisha Toshiba Spatial resolution improvement in computer tomography (CT) using iterative reconstruction
US8861827B2 (en) 2012-02-28 2014-10-14 General Electric Company System and method for determining confidence measurements of single volume elements in computer tomography
DE102012217940A1 (de) * 2012-10-01 2014-04-03 Siemens Aktiengesellschaft Rekonstruktion von Bilddaten
US10679104B2 (en) * 2012-10-26 2020-06-09 Varian Medical Systems, Inc. Forward projection cylinder volume method
DE102012222149B4 (de) 2012-12-04 2014-09-18 Siemens Aktiengesellschaft Verfahren zur Erzeugung von Magnetresonanzmessdaten
KR102060659B1 (ko) * 2013-03-20 2019-12-30 삼성전자주식회사 영상 처리를 위한 투사 및 역투사 방법 및 그 영상 처리 장치
JP6466085B2 (ja) * 2013-05-27 2019-02-06 キヤノンメディカルシステムズ株式会社 Ｘ線ｃｔ装置及び画像診断装置
DE102014212943B4 (de) 2014-07-03 2016-11-24 Siemens Healthcare Gmbh Magnetresonanz-Bildgebung unter Berücksichtigung von unterschiedlichen Frequenzkodiermustern
WO2017001210A1 (en) * 2015-06-30 2017-01-05 Koninklijke Philips N.V. Registration apparatus for registering images
DE102015224176A1 (de) * 2015-12-03 2017-06-08 Siemens Healthcare Gmbh Tomografieanlage und Verfahren zum Erzeugen einer Abfolge von Volumenbildern eines Gefäßsystems
US10115211B2 (en) * 2016-03-25 2018-10-30 L3 Security & Detection Systems, Inc. Systems and methods for reconstructing projection images from computed tomography volumes
US10055860B2 (en) * 2016-05-27 2018-08-21 Toshiba Medical Systems Corporation Computed tomography apparatus and empirical pre-weighting method for decreasing image noise nonuniformity
CN106934856B (zh) * 2017-03-20 2020-04-21 广东电网有限责任公司电力科学研究院 基于x射线检测技术的三维断层重建及切片显示方法
CN108064396B (zh) * 2017-05-27 2021-04-30 上海联影医疗科技股份有限公司 一种在图像引导放射治疗中补偿诊察台下沉的系统和方法
US10482634B2 (en) * 2017-10-24 2019-11-19 General Electric Company Systems and methods for imaging with anisotropic voxels
US20230021332A1 (en) * 2021-07-26 2023-01-26 GE Precision Healthcare LLC Methods and system for dynamically annotating medical images
CN115631232B (zh) * 2022-11-02 2023-07-25 佛山读图科技有限公司 一种确定双探头探测器径向位置方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5404293A (en) * 1991-06-11 1995-04-04 The University Of Utah Cone beam reconstruction using helical data collection paths
US5907594A (en) * 1997-07-01 1999-05-25 Analogic Corporation Reconstruction of volumetric images by successive approximation in cone-beam computed tomography systems
US5909476A (en) * 1997-09-22 1999-06-01 University Of Iowa Research Foundation Iterative process for reconstructing cone-beam tomographic images
US6266388B1 (en) * 1999-09-07 2001-07-24 General Electric Company Methods and apparatus for two-pass cone beam image reconstruction
US6507633B1 (en) * 2001-02-15 2003-01-14 The Regents Of The University Of Michigan Method for statistically reconstructing a polyenergetic X-ray computed tomography image and image reconstructor apparatus utilizing the method
US6776784B2 (en) * 2001-09-06 2004-08-17 Core Medical, Inc. Clip apparatus for closing septal defects and methods of use
US6702835B2 (en) * 2001-09-07 2004-03-09 Core Medical, Inc. Needle apparatus for closing septal defects and methods for using such apparatus
US20060052821A1 (en) * 2001-09-06 2006-03-09 Ovalis, Inc. Systems and methods for treating septal defects
US6907102B1 (en) * 2002-12-16 2005-06-14 Ken Sauer Iterative reconstruction methods for multi-slice computed tomography
US6768782B1 (en) * 2002-12-16 2004-07-27 University Of Notre Dame Du Lac Iterative method for region-of-interest reconstruction

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108352077A (zh) * 2016-03-31 2018-07-31 上海联影医疗科技有限公司 图像重建的系统和方法
CN108352077B (zh) * 2016-03-31 2022-03-01 上海联影医疗科技股份有限公司 图像重建的系统和方法
CN107495976A (zh) * 2016-06-14 2017-12-22 上海联影医疗科技有限公司 图像重建中的最大值和灰度值图像的获取方法及装置
CN107492131A (zh) * 2017-07-01 2017-12-19 武汉斗鱼网络科技有限公司 用于安卓电视的倒影生成方法、存储介质、设备及系统
CN107492131B (zh) * 2017-07-01 2020-06-12 武汉斗鱼网络科技有限公司 用于安卓电视的倒影生成方法、存储介质、设备及系统
CN113450287A (zh) * 2021-07-27 2021-09-28 北京与光科技有限公司 高分辨率光谱恢复方法、装置和电子设备

Also Published As

Publication number Publication date
DE112007001451B4 (de) 2016-05-25
JP2009540966A (ja) 2009-11-26
JP5198443B2 (ja) 2013-05-15
US20070297660A1 (en) 2007-12-27
WO2007150037A3 (en) 2008-10-23
CN101478919B (zh) 2013-07-17
CN103400401B (zh) 2017-04-12
DE112007001451T5 (de) 2009-05-20
WO2007150037A2 (en) 2007-12-27
CN101478919A (zh) 2009-07-08
US7583780B2 (en) 2009-09-01

Similar Documents

Publication Publication Date Title
CN101478919B (zh) 改善图像分辨率的系统和方法
CN101849247B (zh) 用于生成衰减图的设备和方法
CN105452852B (zh) 通过使用相衬ct的单色衰减对比图像生成
US6339652B1 (en) Source-assisted attenuation correction for emission computed tomography
RU2606561C2 (ru) Компенсация усечения для итерационной реконструкции в компьютерной томографии (кт) с коническим пучком в комбинированных системах офэкт/кт
US4149248A (en) Apparatus and method for reconstructing data
JP5340600B2 (ja) ラドンデータから（ｎ＋１）次元イメージ関数を再構成する方法および装置
CN102655810B (zh) 计算机断层摄影装置
CN102160085B (zh) 成像装置
EP0520778B1 (en) Tomographic image reconstruction using cross-plane rays
GB1577014A (en) Tomographic apparatus and method for reconstructing planar slices from non-absorbed and non-scattered radiation
CN101473348A (zh) 用于误差补偿的方法和系统
JP4468352B2 (ja) コンピュータトモグラフィにおける局所的患者線量の再構成
US8389943B2 (en) Modeling of the point-spread-function in single-pinhole and multi-pinhole spect reconstruction
US6539103B1 (en) Method and apparatus for image reconstruction using a knowledge set
US11947053B2 (en) Weighting detector configurations in SPECT imaging
Zeng et al. Image reconstruction algorithm for a rotating slat collimator
Seeram et al. Computed tomography: Physical principles, instrumentation, and quality control
Huang et al. Reconstruction from uniformly attenuated SPECT projection data using the DBH method
Ghani et al. An empirical method for geometric calibration of a photon counting detector-based cone beam CT system
US20230309935A1 (en) Methods and systems for contrast-to-noise evaluation of computed tomography systems
Welch et al. Accurate attenuation correction in SPECT without transmission measurements
Kärnä Quality Assurance and Patient Dose Monitoring Methods in Computed Tomography
Natterer et al. Accurate attenuation correction in SPECT without transmission measurements
Hamilton et al. Single Photon Tomographic Imaging

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant