CN101478919B - 改善图像分辨率的系统和方法 - Google Patents
改善图像分辨率的系统和方法 Download PDFInfo
- Publication number
- CN101478919B CN101478919B CN2007800234587A CN200780023458A CN101478919B CN 101478919 B CN101478919 B CN 101478919B CN 2007800234587 A CN2007800234587 A CN 2007800234587A CN 200780023458 A CN200780023458 A CN 200780023458A CN 101478919 B CN101478919 B CN 101478919B
- Authority
- CN
- China
- Prior art keywords
- image
- voxel
- image voxel
- resolution
- detector
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 88
- 238000003384 imaging method Methods 0.000 claims abstract description 18
- 238000001514 detection method Methods 0.000 claims description 32
- 230000000007 visual effect Effects 0.000 claims description 25
- 230000008859 change Effects 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 20
- 238000002591 computed tomography Methods 0.000 claims description 19
- 230000006872 improvement Effects 0.000 claims description 6
- 238000013459 approach Methods 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 239000000470 constituent Substances 0.000 claims 4
- 230000004304 visual acuity Effects 0.000 claims 1
- 238000005303 weighing Methods 0.000 claims 1
- 230000006870 function Effects 0.000 description 36
- 238000012360 testing method Methods 0.000 description 35
- 238000013170 computed tomography imaging Methods 0.000 description 9
- 210000003484 anatomy Anatomy 0.000 description 8
- 238000005259 measurement Methods 0.000 description 7
- 230000006866 deterioration Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 238000002600 positron emission tomography Methods 0.000 description 6
- 238000001914 filtration Methods 0.000 description 5
- 238000003491 array Methods 0.000 description 4
- 230000014509 gene expression Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 1
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007796 conventional method 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
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 238000012633 nuclear imaging Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 238000001959 radiotherapy Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 238000000411 transmission spectrum Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
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
发明背景
本发明总体上涉及成像系统,更具体地,涉及改善图像分辨率的系统和方法。
计算机断层摄影(CT)成像系统通常包括x射线源,x射线源投射扇形x射线束透过患者到达辐射探测器阵列。x射线束经准直落在通常称为“成像平面”的xy平面内。探测器阵列接收到的x射线束的辐射强度取决于患者使x射线束的衰减。单独采集探测器阵列多个探测器单元的衰减测量结果而生成透射谱。
在台架内围绕受检患者转动x射线源和探测器,以使x射线束与患者相交所成的投射角不断变化。将一组x射线衰减测量结果称为“视角(view)”,所述测量结果是来自成一种台架角度或一种投射角度的探测器阵列的模拟投影数据。对患者的“扫描”包括在x射线源和探测器旋转一周的过程中以不同的投射角形成的一组视角。
为缩短获取多个切片所用的总扫描时间,可进行螺旋扫描。螺旋扫描技术允许使用单光子源以较快的速度进行大体积扫描。为进行螺旋扫描,使患者静躺于其上的平台沿台架旋转所围绕的z轴移动,同时采集规定数量切片的模拟投影数据。螺旋扫描生成单螺线。x射线束绘出的螺线形成模拟投影数据,可根据所述模拟投影数据,重建各规定切片中的图像。除缩短扫描时间以外,螺旋扫描还具有其它优势,例如更好地使用注入对比度(injected contrast)、任意位置处改善的图像重建和较好的三维图像。螺旋扫描的实例包括多切片螺旋扫描。在多切片螺旋扫描中,探测器阵列沿z轴延伸。通常,在多切片螺旋扫描中,探测器阵列包括多行,每一行相应于沿z轴不同的位置和不同的测量切片。在轴向扫描中,处理模拟投影数据以构建图像,该图像相应于透过患者所取的二维切片。对于离散切片,可进行全视场迭代重建,以提高图像质量。
对于连续扫描,形成这样的扫描图案,其中患者沿z轴的部位随着台架的旋转而线性变化。在采集数据过程中,对连续扫描图案进行量化,生成患者周围有限个x射线源位置的一组离散的投影视角。常规的直接图像重建技术,例如二维或三维滤波反投影,根据投影数据如下重建图像体素:在投影数据中内插要素以将各投影角的贡献累加到多个图像体素中,从而通过对投影数据的单遍扫描形成图像或图像体(image volume)。应用滤波反投射生成的图像的典型分辨率基于探测器阵列的尺寸、焦斑尺寸、数据采集系统(DAS)采集模拟投影数据的采样速度和滤波反投射过程中过滤投影数据的滤波器的kernel。通常情况下,典型分辨率不高于CT成像系统等角点处各探测单元的投影尺寸。根据尼奎斯特定理,无需以大于两倍的极限典型分辨率采样。然而,通过常规的直接图像重建技术生成的图像体通常不具备高的空间分辨率。
发明内容
一方面,描述断层摄影成像系统中重建图像的方法。该方法包括通过迭代重建图像改善图像的空间分辨率。
另一方面,描述权衡(trade-off)断层摄影成像系统性能的迭代重建方法。该方法包括调整一部分迭代重建图像的噪声和分辨率。
再一方面,描述调整断层摄影成像系统分辨率的迭代重建方法。该方法包括调整迭代重建图像的面内分辨率(in-plane resolution),以及调整迭代重建图像的断面分辨率(cross-plane resolution)。
再一方面,描述改善图像空间分辨率的方法。该方法包括迭代重建图像并改善图像的空间分辨率。该方法通过建立具有平滑曲线的正投影函数和建立具有平滑曲线的图像中的至少一种来改善图像的空间分辨率。正投影的建立包括通过确定一系列投影值来建立正投影。图像的建立包括通过确定所述一系列投影值的逆(inverse)来重建图像。
再一方面,描述改善图像空间分辨率的处理器。配置该处理器以迭代重建图像并改善图像的空间分辨率。配置该处理器,以通过建立具有平滑曲线的正投影函数和建立具有平滑曲线的图像中的至少一种来改善图像的空间分辨率。正投影的建立包括通过确定一系列投影值来建立正投影。图像的建立包括通过确定所述一系列投影值的逆来重建图像。
附图说明
图1是实现图像分辨率改善方法的多切片CT成像系统的实施方式的等视图。
图2是图1的CT成像系统的方框图。
图3是图像分辨率改善方法的实施方式的示意图。
图4是图3的方法的示意图。
图5是图像分辨率改善方法的实施方式的流程图。
图6是图5的流程图的续图。
图7示出多个图像的实施方式,显示应用图像分辨率改善方法的实施方式的效果。
具体实施方式
迭代重建包括如下形成图像的方法:根据实际投影数据和预期投影数据之间匹配的质量基于图像的先验估计值(current estimate)反复调整现有的估计值。对图像多个特性的考虑,例如图像的平滑度和/或预建模型的满意度,也可能影响匹配质量。进行多次迭代,从而形成基于预定标准与实际投影数据匹配最佳的图像。由于使一组重建图像形成患者的三维重现,并且各图像像素或像点对应于三维重建中的单个体素或体元(volumeelement),因而将整组重建图像称为三维重建。
参考图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代表图像体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以加权平方形式表示为:
其中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)中 的值减小。图像重建器44根据计算正投影函数F(r)和正则化函数U(r)之前的前一次迭代中计算得到的 值,在本次迭代中计算正投影函数F(r)和正则化函数U(r)。在本次迭代中,图像重建器44以 代替等式(1)中的r。在前一次迭代中,图像重建器44由根据实际投影数据p滤波反投射重建的图像体318生成 在替换性实施方式中,在前一次迭代中,图像重建器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)中的 并相对于滤波反投影生成的典型分辨率改善空间分辨率
在另一替换性实施方式中,图像重建器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和与其邻近的至少一个图像体素之间尺寸大小之差的函数。例如,U(r)=
其中α为度量常数如正实数,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 (26)
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.一种权衡断层摄影成像系统的性能的迭代重建方法,所述方法包括调整一部分迭代重建图像的噪声和分辨率,其中所述迭代重建包括应用正则化函数,其中所述图像包括图像体素,所述方法还包括模型化所述图像体素,所述图像体素具有平行于第一坐标轴的第一维度和平行于第二坐标轴的第二维度,其中所述第一维度不同于所述第二维度,并且其中所述第一坐标轴不同于所述第二坐标轴。
15.一种调整断层摄影成像系统的分辨率的迭代重建方法,所述方法包括:
调整迭代重建图像的面内分辨率;和
调整迭代重建图像的断面分辨率;
其中所述迭代重建包括应用正则化函数,其中所述图像包括图像体素,所述方法还包括模型化所述图像体素,所述图像体素具有平行于第一坐标轴的第一维度和平行于第二坐标轴的第二维度,其中所述第一维度不同于所述第二维度,并且其中所述第一坐标轴不同于所述第二坐标轴。
16.权利要求15的方法,其中所述面内分辨率包括径向分辨率和方位角分辨率。
17.权利要求15的方法,其中所述面内分辨率等于所述断面分辨率。
18.一种改善图像的空间分辨率的方法,所述方法包括:
迭代重建所述图像,其中所述迭代重建图像包括应用正则化函数;和
通过下述中的至少一种改善所述图像的空间分辨率:
建立具有平滑曲线的正投影函数,其中所述建立正投影包括通过确定一组投影值重建正投影;和
建立具有光滑曲线的图像,其中所述建立图像包括通过确定所述一组投影值的逆重建图像,其中所述图像包括图像体素,所述方法还包括模型化所述图像体素,所述图像体素具有平行于第一坐标轴的第一维度和平行于第二坐标轴的第二维度,其中所述第一维度不同于所述第二维度,并且其中所述第一坐标轴不同于所述第二坐标轴。
19.权利要求18的方法,其中所述改善空间分辨率包括沿探测器信道轴和探测器行轴中的至少一个调整所述图像的图像体素的尺寸。
20.权利要求18的方法,其中所述改善空间分辨率包括沿探测器信道轴和探测器行轴中的至少一个调整由图像体素正投影的探测元的数目,所述图像体素是图像的一部分并以一定的投影角形成。
21.权利要求18的方法,其中所述改善空间分辨率包括沿探测器信道轴和探测器行轴中的至少一个调整探测元的数目,其中所述探测元的中心为图像体素的正投影,所述图像体素是图像的一部分并以一定的投影角形成。
22.权利要求18的方法,其中所述改善空间分辨率包括生成图像体素的正投影,其中所述生成正投影包括将图像体素正投影在至少一个探测元上,其中所述正投影沿探测器信道轴和探测器行轴中至少一个的尺寸与所述图像体素的尺寸成正比。
23.权利要求18的方法,其中所述改善空间分辨率包括沿探测器信道轴和探测器行轴中的至少一个调整正投影区的尺寸,所述正投影区通过正投影图像的图像体素而生成。
24.权利要求18的方法,其中所述改善空间分辨率包括沿探测器信道轴和探测器行轴中的至少一个调整正投影区的形状,所述正投影区通过正投影图像的图像体素而生成。
25.权利要求18的方法,其中所述改善空间分辨率包括在图像所呈现的视场内实现均匀的空间分辨率,其中所述均匀的空间分辨率如下实现:改变多个正投影沿探测器信道轴和探测器行轴中至少一个的多个尺寸,所述多个正投影通过正投影图像的多个图像体素而形成。
26.权利要求18的方法,其中所述改善空间分辨率包括在图像所呈现的视场内实现均匀的空间分辨率,其中所述均匀的空间分辨率如下实现:改变多个正投影沿探测器信道轴和探测器行轴中至少一个的多个形状,所述多个正投影通过正投影图像的多个图像体素而形成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310233382.8A CN103400401B (zh) | 2006-06-22 | 2007-06-22 | 改善图像分辨率的系统和方法 |
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 | ||
PCT/US2007/071917 WO2007150037A2 (en) | 2006-06-22 | 2007-06-22 | Systems and methods for improving a resolution of an image |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310233382.8A Division CN103400401B (zh) | 2006-06-22 | 2007-06-22 | 改善图像分辨率的系统和方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101478919A CN101478919A (zh) | 2009-07-08 |
CN101478919B true CN101478919B (zh) | 2013-07-17 |
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 After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310233382.8A Active CN103400401B (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) |
Families Citing this family (36)
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 | 株式会社東芝 | X線コンピュータ断層撮影装置、医用画像処理装置、及び医用画像処理プログラム |
CN102576467B (zh) * | 2009-07-14 | 2015-11-25 | 皇家飞利浦电子股份有限公司 | 包括移变模糊补偿的图像重建 |
US8655033B2 (en) * | 2009-10-28 | 2014-02-18 | General Electric Company | Iterative reconstruction |
AU2010246365B2 (en) * | 2009-11-27 | 2014-04-24 | Ge Sensing & Inspection Technologies Gmbh | Computed tomography method, computer software, computing device and computed tomography system for determining a volumetric representation of a sample |
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 |
EP2591431B1 (en) * | 2010-07-09 | 2020-01-01 | GE Sensing & Inspection Technologies GmbH | Computed tomography method, computer program, computing device and computed tomography system |
DE102010041191B4 (de) * | 2010-09-22 | 2016-02-18 | Siemens Aktiengesellschaft | Erstellung von MR-Bilddaten mit paralleler Schichtanregung und teilweiser Überlappung der Schichten im Frequenzbereich |
JP5978429B2 (ja) * | 2011-04-28 | 2016-08-24 | 株式会社日立製作所 | 医用画像処理装置、医用画像処理方法 |
US8483471B2 (en) * | 2011-06-30 | 2013-07-09 | General Electric Company | Method and system for scatter correction in X-ray imaging |
US9824467B2 (en) * | 2011-06-30 | 2017-11-21 | Analogic Corporation | Iterative image reconstruction |
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 | 삼성전자주식회사 | 영상 처리를 위한 투사 및 역투사 방법 및 그 영상 처리 장치 |
CN105283132B (zh) * | 2013-05-27 | 2018-10-12 | 东芝医疗系统株式会社 | X射线ct装置以及图像诊断装置 |
DE102014212943B4 (de) | 2014-07-03 | 2016-11-24 | Siemens Healthcare Gmbh | Magnetresonanz-Bildgebung unter Berücksichtigung von unterschiedlichen Frequenzkodiermustern |
US10685452B2 (en) * | 2015-06-30 | 2020-06-16 | 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 |
EP3245634B1 (en) * | 2016-03-31 | 2020-10-14 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image reconstruction |
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 |
CN107495976B (zh) * | 2016-06-14 | 2021-01-01 | 上海联影医疗科技股份有限公司 | 图像重建中的最大值和灰度值图像的获取方法及装置 |
CN106934856B (zh) * | 2017-03-20 | 2020-04-21 | 广东电网有限责任公司电力科学研究院 | 基于x射线检测技术的三维断层重建及切片显示方法 |
CN108064396B (zh) * | 2017-05-27 | 2021-04-30 | 上海联影医疗科技股份有限公司 | 一种在图像引导放射治疗中补偿诊察台下沉的系统和方法 |
CN107492131B (zh) * | 2017-07-01 | 2020-06-12 | 武汉斗鱼网络科技有限公司 | 用于安卓电视的倒影生成方法、存储介质、设备及系统 |
US10482634B2 (en) * | 2017-10-24 | 2019-11-19 | General Electric Company | Systems and methods for imaging with anisotropic voxels |
CN113450287B (zh) * | 2021-07-27 | 2024-07-26 | 北京与光科技有限公司 | 高分辨率光谱恢复方法、装置和电子设备 |
CN115631232B (zh) * | 2022-11-02 | 2023-07-25 | 佛山读图科技有限公司 | 一种确定双探头探测器径向位置方法 |
Citations (3)
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 |
CN1268034A (zh) * | 1997-07-01 | 2000-09-27 | 模拟技术公司 | 迭代锥束ct再现 |
US6266388B1 (en) * | 1999-09-07 | 2001-07-24 | General Electric Company | Methods and apparatus for two-pass cone beam image reconstruction |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5909476A (en) * | 1997-09-22 | 1999-06-01 | University Of Iowa Research Foundation | Iterative process for reconstructing cone-beam tomographic images |
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 |
US20060052821A1 (en) * | 2001-09-06 | 2006-03-09 | Ovalis, Inc. | Systems and methods for treating septal defects |
US6702835B2 (en) * | 2001-09-07 | 2004-03-09 | Core Medical, Inc. | Needle apparatus for closing septal defects and methods for using such apparatus |
US6776784B2 (en) * | 2001-09-06 | 2004-08-17 | Core Medical, Inc. | Clip apparatus for closing septal defects and methods of use |
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 |
-
2006
- 2006-06-22 US US11/472,572 patent/US7583780B2/en active Active
-
2007
- 2007-06-22 CN CN2007800234587A patent/CN101478919B/zh active Active
- 2007-06-22 WO PCT/US2007/071917 patent/WO2007150037A2/en active Application Filing
- 2007-06-22 JP JP2009516750A patent/JP5198443B2/ja active Active
- 2007-06-22 DE DE112007001451.2T patent/DE112007001451B4/de active Active
- 2007-06-22 CN CN201310233382.8A patent/CN103400401B/zh active Active
Patent Citations (3)
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 |
CN1268034A (zh) * | 1997-07-01 | 2000-09-27 | 模拟技术公司 | 迭代锥束ct再现 |
US6266388B1 (en) * | 1999-09-07 | 2001-07-24 | General Electric Company | Methods and apparatus for two-pass cone beam image reconstruction |
Also Published As
Publication number | Publication date |
---|---|
US20070297660A1 (en) | 2007-12-27 |
US7583780B2 (en) | 2009-09-01 |
JP2009540966A (ja) | 2009-11-26 |
DE112007001451T5 (de) | 2009-05-20 |
DE112007001451B4 (de) | 2016-05-25 |
CN101478919A (zh) | 2009-07-08 |
WO2007150037A3 (en) | 2008-10-23 |
CN103400401A (zh) | 2013-11-20 |
CN103400401B (zh) | 2017-04-12 |
WO2007150037A2 (en) | 2007-12-27 |
JP5198443B2 (ja) | 2013-05-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101478919B (zh) | 改善图像分辨率的系统和方法 | |
US4149248A (en) | Apparatus and method for reconstructing data | |
JP5340600B2 (ja) | ラドンデータから(n+1)次元イメージ関数を再構成する方法および装置 | |
CN101849247B (zh) | 用于生成衰减图的设备和方法 | |
CN102655810B (zh) | 计算机断层摄影装置 | |
EP0520778B1 (en) | Tomographic image reconstruction using cross-plane rays | |
CN102160085B (zh) | 成像装置 | |
GB1577014A (en) | Tomographic apparatus and method for reconstructing planar slices from non-absorbed and non-scattered radiation | |
US20140219417A1 (en) | Filter for Tomographic Reconstructions | |
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 | |
US6490335B1 (en) | Helical segment image reconstruction | |
WO1993016435A1 (en) | Method of recovering signals in tomographic measurement systems | |
US11947053B2 (en) | Weighting detector configurations in SPECT imaging | |
Zeng et al. | Image reconstruction algorithm for a rotating slat collimator | |
Huang et al. | Reconstruction from uniformly attenuated SPECT projection data using the DBH method | |
US9508164B2 (en) | Fast iterative image reconstruction method for 3D computed tomography | |
Guillemaud et al. | A multifocal collimator with circularly distributed focal points for SPECT imaging | |
Hsieh et al. | Optimization of system design for 64-slice cone-beam computed tomography | |
Schmidt et al. | MTF analysis of a prototype table-top inverse-geometry volumetric CT system | |
Guillemaud et al. | A circular multifocal collimator for 3D SPECT imaging | |
Siemens | Tomographic apparatus for reconstructing planar slices from non-absorbed and non-scattered radiation |
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 |