CN102067177B - 优化剂量控制的图像生成装置 - Google Patents

优化剂量控制的图像生成装置 Download PDF

Info

Publication number
CN102067177B
CN102067177B CN200980123925.2A CN200980123925A CN102067177B CN 102067177 B CN102067177 B CN 102067177B CN 200980123925 A CN200980123925 A CN 200980123925A CN 102067177 B CN102067177 B CN 102067177B
Authority
CN
China
Prior art keywords
measurement data
interest
noise
generation device
image
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
Application number
CN200980123925.2A
Other languages
English (en)
Other versions
CN102067177A (zh
Inventor
J·维格特
M·贝尔特拉姆
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
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 Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN102067177A publication Critical patent/CN102067177A/zh
Application granted granted Critical
Publication of CN102067177B publication Critical patent/CN102067177B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computerised tomographs
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/54Control of apparatus or devices for radiation diagnosis
    • A61B6/542Control of apparatus or devices for radiation diagnosis involving control of exposure

Abstract

本发明涉及用于根据测量数据生成图像的图像生成装置,其中,针对感兴趣区域优化图像质量,并且涉及包括这种图像生成装置的成像系统。所述图像生成装置包括用于确定感兴趣区域的投影域中的噪声分布的噪声确定单元,以及用于基于所述确定的噪声分布使用噪声传播算法确定针对所述图像生成装置的辐射源(2)的剂量分布图的剂量控制单元(32)。因此,可以改善重建体积的信噪比而不敏感地依赖于选定的感兴趣区域。

Description

优化剂量控制的图像生成装置
技术领域
本发明涉及图像生成装置、图像生成方法和用于根据测量数据生成图像的计算机程序,其中,针对感兴趣区域优化图像质量。本发明还涉及包括图像生成装置以及相应的成像方法和计算机程序的成像系统。
背景技术
测定体积的成像功能,诸如计算机断层摄影(CT)、介入式CT、CT荧光透视术、3D血管成像等,在过去的几十年间已经成为越来越有价值的工具。X射线CT是一种通过将X射线穿过对象并测量穿过对象的X射线衰减来确定对象内部构造的技术。在这种技术中,对象被细分为许多体素,体素是用于成像目的的基本测定体积的单元。使用基于X射线的测定体积的成像技术的放射程序的数目已经相应增加。然而,这些技术占据放射学中总应用剂量的大部分,并且因此,应用的患者剂量已经成为关键性问题。为了应对日益增加的患者剂量,需要有效的剂量节约技术。
EP 1 172 069 A1公开了一种剂量优化的CT,所述剂量优化通过实时设置经优化的管电流(自动曝光控制)、管电流调制(剂量最小化)以及据此通过自适应3D滤波器的后处理(降噪)实现。基于在探测器中心测量的衰减来计算用于采集投影数据的剂量分布图。对于重建的数据而言,这将导致以下方式的体素-噪声的分布:即,使重建对象中心的噪声最优。尽管这种技术已经在CT中得到广泛认可,但由于若干原因,该程序并不是最优的。第一原因是,仅使用测量的衰减作为独立投影对重建体积中总噪声贡献的估计的基础仅仅是近似的。源于例如射束硬化的频谱效应和散射辐射的影响被忽略。第二原因是,如果在重建之前使用散射偏差校正技术,在计算最优剂量分布图时需要充分考虑表示测量信号重要部分的散射辐射以及散射减小自身的影响。第三原因是,如果重建体积的中心区域之外的区域是感兴趣区域(例如,对于心脏CT而言),针对重建图像中心和感兴趣区域的最优剂量分布图可能显著不同。针对重建体积中心的剂量分布图甚至可能降低感兴趣区域的对比噪声比。因此,关键是要具体针对感兴趣区域进行噪声/剂量优化。
US2004/0062341A1公开了一种用于CT扫描中的自动曝光控制的方法和设备,其中,根据绕受检者的X射线聚焦的回转的第一半期间受检者的切片(slice)计算实际衰减分布图,并调节诸如管电流的X射线源的操作参数,从而根据推断计算的衰减分布图,在回转的第二半期间修改由X射线源发射的辐射剂量,以实现针对当前正在被扫描的切片在图像中预确定目标像素噪声。
发明内容
本发明的目的是提供图像生成装置、图像生成方法和用于根据测量数据生成图像的计算机程序,其提供了增强的剂量控制。
根据本发明的第一方面,提出了用于根据测量数据生成图像的图像生成装置,其中,针对感兴趣区域优化图像质量,其中,该图像生成装置包括:
-噪声确定单元,其用于确定感兴趣区域的投影域中的噪声分布;以及
-剂量控制单元,其用于基于所述确定的噪声分布使用噪声传播算法确定针对所述图像生成装置的辐射源的剂量分布图。
本发明是基于这样的思想:即基于投影中的噪声分布在重建体积的独立体素中确定噪声。使用该投影噪声分布,可以基于噪声的传播计算最优剂量分布图。特别针对重建体积的选定感兴趣区域改善信噪比。这种解决方案提供了多个优点,特别是对于不对称的患者几何结构,包括通过使用相同的剂量得到更高的图像质量、通过保持图像质量降低剂量以及噪声诱发地伪影图案的均匀度。
在噪声传播中,可以考虑对象中X射线的散射以及诸如射束硬化的频谱效应。
在优选实施例中,可以由用户指定感兴趣区域。
在另一优选实施例中,可以利用自动或半自动的器官分割选取感兴趣区域。
在另一实施例中,可以基于之前采集的投影数据确定剂量分布图。在对相同或相似的几何结构重复采集的情况下,可以从先前的扫描中获取所需的信息。
在另一优选实施例中,可以基于包括所述感兴趣区域的对象模型来确定剂量分布图。因此,可以使用预先计算的数据以获取剂量分布图。
在另一优选实施例中,可以基于追踪扫描确定剂量分布图,追踪扫描通常被用作低剂量和/或低分辨率的扫描,以便确定待成像的FOV。
在又一实施例中,噪声传播算法可以特别考虑频谱效应(例如射束硬化)和散射辐射中的至少一个。
优选地,基于低于所述图像的最终分辨率的经降低的空间分辨率和经简化的噪声传播算法中的至少一个确定剂量分布图,从而加速对最优剂量调制的确定。
根据本发明的另一方面,提出了成像系统,其中,该成像系统包括:
-测量数据生成单元,其用于生成图像的测量数据,以及
-图像生成装置,其用于根据权利要求1中限定的测量数据生成所述图像。
根据本发明的另一方面,提出了用于根据测量数据生成图像的图像生成方法,其中,针对感兴趣区域优化图像质量,并且其中,该图像生成方法包括以下步骤:
-确定感兴趣区域的投影域中的噪声分布;以及
-基于所述确定的噪声分布通过使用噪声传播算法确定针对辐射的剂量分布图。
根据本发明的另一方面,成像方法包括以下步骤:
-通过测量数据生成单元生成测量数据,以及
-根据权利要求1中限定的步骤通过图像生成装置根据测量数据生成图像。
根据本发明的另一方面,提出了用于根据测量数据生成图像的计算机程序,其中,该计算机程序包括:当该计算机程序在控制图像生成装置的计算机上运行时,用于令如权利要求1中限定的图像生成装置执行如权利要求11中限定的检测方法的步骤的程序代码模块。
根据本发明的另一方面,提出了用于成像的计算机程序,其中,该计算机程序包括:当该计算机程序在控制成像系统的计算机上运行时,用于令如权利要求10中限定的成像系统执行如权利要求12中限定的成像方法的步骤的程序代码模块。
应当理解,图像生成装置、成像系统、成像生成方法、成像方法和计算机程序具有如从属权利要求中限定的相似和/或相同的优选实施例。
应当理解,本发明的优选实施例还可以是从属权利要求与各个独立权利要求的任意组合。
附图说明
本发明的这些方面和其他方面将从下文描述的实施例中变得显而易见,并将参考这些实施例加以阐述。在以下附图中:
图1示意性和示例性地示出了用于对感兴趣区域成像、包括根据本发明的增强剂量控制的成像系统的表示;
图2示出了根据本发明的实施例的剂量控制过程;
图3示出了根据本发明用于确定噪声贡献的过程;以及
图4图示说明了本发明中使用的几何结构和术语。
具体实施方式
图1示意性和示例性地示出了是计算机断层摄影系统的成像系统。计算机断层摄影系统包括机架1,机架1能够绕平行于z方向延伸的旋转轴R旋转。辐射源2被安装在机架1上,在本实施例中,辐射源2是X射线管。辐射源2配备有准直器3,在本实施例中,准直器3根据辐射源2生成的辐射形成圆锥形辐射束4。辐射穿透检查区域中诸如患者的对象(未示出)和优选位于对象内的感兴趣区域(ROI),在本实施例中,检查区域是圆柱形的。辐射束4穿透检查区域5之后入射到检测装置6上,检测装置6包括二维检测表面。检测装置6安装在机架1上。该计算机断层摄影系统可能适于利用针对整个患者而不仅仅是针对感兴趣区域测量的数据来生成“正常”的计算机断层摄影扫描,同时通过以下方式执行剂量控制:即针对包括整个患者(并且不仅仅是ROI)的重建图像内的(用户指定的)ROI优化图像质量。
计算机断层摄影系统包括两个电动机7、8。由电动机7以优选恒定但可调节的角速度驱动机架1。电动机8被用于移动对象,例如平行于旋转轴R或z轴方向在检查区域5中的患者台上布置的患者。这些电动机7、8受控制单元9控制,例如,使得辐射源和检查区域5以及因此检查区域内的感兴趣区域相对彼此沿着螺旋形轨线移动。然而,还可能不移动对象或检查区域5,而仅旋转辐射源2,即辐射源沿着相对于对象或检查区域5的圆形轨线移动。此外,在另一实施例中,准直器3可能适于形成另一射束形状,具体为扇形射束,并且检测装置6可以包括检测表面,该检测表面对应于另一射束形状成形,具体为扇形射束。
在辐射源和检查区域5的相对运动期间,检测装置6根据检测装置6的检测表面上的辐射入射生成测量值。因此,辐射源2、用于使辐射源2相对于检查区域移动的元件,具体为电动机7、8和机架1,以及检测装置6形成用于生成测量数据的测量数据生成单元。
在本实施例中,测量数据为投影数据,该测量数据被提供给图像生成和显示装置10,用于根据测量数据、即根据投影数据生成和显示包括感兴趣区域的图像。感兴趣区域位于检查区域内,并且优选包含对象或对象的一部分。图像生成和显示装置10包括用于根据测量数据生成和处理图像的图像生成单元12,以及具有用于显示生成的图像的显示屏幕14的图像显示单元13。
图像生成和显示装置10也可以优选由控制单元9控制。在另一实施例中,控制单元9还可以执行对图像生成单元12或其各部分的控制。
此外,提供了剂量调制单元11(例如剂量调制处理器等),剂量调制单元11控制辐射源2的供电电源(未示出),或与辐射源2的供电电源通信,从而控制断层摄影成像期间由辐射源2生成的辐射强度。在图示说明的实施例中,辐射源2是X射线管,并且剂量调制单元11可能适于控制X射线管灯丝或阴极电流,从而调制由辐射源2生成的X射线强度。在其他涵盖的实施例中,可以通过以下方式产生剂量调制单元11的剂量调制:通过关闭或过滤辐射束、通过调制管电压、通过调制维纳尔圆柱电极(Wehneltcylinder)上的偏置电压等。
当在感兴趣区域中采集断层摄影成像数据时,剂量调制单元11的剂量调制通常降低由辐射源2生成的辐射强度。
在所有的本发明的实施例中,本发明的基本程序是采用噪声传播模型,该噪声传播模型允许基于投影数据中得到的噪声分布空间分辨地计算重建图像中的体素噪声。具体而言,噪声传播模型使得能够精确地计算源自每个独立投影方向的噪声贡献。使用这些来自每个投影方向的噪声贡献,所提出的程序使得能够计算所需的剂量/管电流分布图,用以实现针对用户选定的感兴趣区域(ROI),即在可重建视野(FOV)中的任何选定位置,在对比噪声比方面改善图像质量。
倘若是散射模型或是散射估计程序,该程序或算法可以明确地考虑散射辐射的影响,还提高了基于电流衰减方法的剂量效率。
图2示出了所提出的剂量优化处理的实例的流程图,该剂量优化处理可以在控制单元9或附加的分立剂量控制单元或处理器(未示出)中执行,以控制剂量调制单元11。在第一控制阶段31中,确定了投影域中预期的噪声分布,其被称为“投影噪声分布”。这可以通过不同的备选方法实现,例如,可以基于之前采集的投影数据ACQ、基于患者模型MOD或基于追踪扫描SS获取预期的噪声分布。这些备选方法可以单独或组合应用。使用获取的投影噪声分布,在第二控制阶段32中可以使用噪声传播算法计算最优剂量分布图。对投影中的噪声分布的确定或计算和对噪声传播算法两者而言,可以基于散射模型SM或散射估计算法将投影中包含的散射辐射的量考虑进来。针对待重建体积中用户可选的位置或ROI具体计算最优噪声分布。最后,在第三控制阶段33中,将最优剂量分布(管电流分布图)传送到剂量调制单元11。
图3示出了表示例如常规CT系统或基于C型臂的平板检测器锥形射束CT系统的断层摄影系统的经简化成像链的流程图。
通过二次误差传播一般可以示出:不相关非平稳噪声分布σ1 2(x)通过线性移不变系统gG(x)传播为
σ2 2(x)=gG 2(x)*σ1 2(x)       (1)
滤波反向投影类型的重建算法可以被视为线性系统,并且因此,可以基于等式1开发通过这种算法的噪声传播。
通常,这种系统的成像链包括与检测器相关的图像形成构成41,诸如检测器点扩散函数(PSF)以及基于投影的处理(组合、低通滤波(未示出))。此外,该成像链包括与重建算法相关的图像形成构成42,如重建滤波器,以及涉及插值和反向投影的构成43。
对于图3中给出的实例,一个投影i对重建体素j的总噪声的噪声贡献可以计算如下:
σ 3 , i 2 ( x j ) = a ij 2 · ( g G 2 * σ 1 2 ) ( x 0 , ij ) + b ij 2 · ( g G 2 * σ 1 2 ) ( x 0 , ij + 1 ) + 2 · a ij 2 · b ij 2 · ( ( g G ( x ) g G ( x + 1 ) ) * σ 1 2 ( x ) ) ( x 0 , ij ) - - - ( 2 )
在这一等式中,gG(x)表示联合传递函数gG(x),其包括检测器点扩散函数gS(x)和重建滤波器的传递函数gF(x)。为了简化,对于等式2中的插值步骤而言,已经假定了扇形射束几何结构(两个插值权值),而直接推出锥形射束的扩展情况(一般四个插值权值)。
图4图示说明了此处假定为扇形射束几何结构所使用的几何结构和术语。对于投影i而言,来自焦斑F、通过位置xj处的体素j中心的射线确定检测器上的插值位置以及相应的插值权值aij和bij
使用等式2,计算每个独立投影i对体素j的噪声方差的贡献。根据下式通过加和来自所有投影的独立贡献计算体素j的总噪声方差:
σ 3 2 ( x j ) = Σ i σ 3 , i 2 ( x j ) - - - ( 3 )
可以以线性方式通过加和所有的西格玛(sigma)值,或者也可以使用均方(squared average),即建立均方根(RMS)值,来计算整个感兴趣区域的体素噪声的代表值。例如,后者可以如下实现。作为对ROI内体素噪声的测量值,RMS值可以根据下式计算:
σ RMS ( ROI ) = 1 N j Σ j σ 3 2 ( x j ) = 1 N j Σ j Σ i σ 3 , i 2 ( x j ) - - - ( 4 )
这样,大值(“噪声”离群值)对总噪声有较大影响,并且因此将对下文将讨论的电流调制具有更强的影响。
可以通过以下方式通过在辐射源2处调制管电流来减小感兴趣区域中的平均体素噪声σRMS(ROI),即对于具有大噪声贡献的投影,增大管电流,而对于具有小噪声贡献的投影,使用减小的管电流。作为任选的条件,总共施加的辐射可以保持恒定。
例如通过剂量调制单元11,可以利用相对管电流调制因子Ci将管电流调制整合到上述实施例中。投影的管电流被计算为Ii=Ci·I平均,其中,I平均是未经调制使用的平均管电流。由于管电流调制因子,独立方差贡献与调制因子的倒数成比例,使得此时将一个体素的体素噪声计算为
σ 3 2 ( x j ) = Σ i 1 C i · σ 3 , i 2 ( x j ) . - - - ( 5 )
使用条件
1 N i Σ i C i = ! 1 - - - ( 6 )
然后,最小化任务可以通过下式限定
σ RMS 2 ( ROI ) = 1 N j Σ i 1 C i Σ j σ 3 , i 2 ( x j ) = Σ i 1 C i σ pro , i 2 = ! min . - - - ( 7 )
在以上等式中,源于每个投影对总的贡献已经集中到了单个量中。
使用偏导数和等式6限定的条件,等式7给出的最优任务的解为
C i = 1 1 N i Σ i ′ σ pro , i ′ · σ pro , i - - - ( 8 )
即,等式8限定了最优相对管电流调制因子。
可以基于对输入噪声方差的描述实现噪声传播的确定和最优相对管电流调制因子的计算。因为在传播通过患者期间,多能射束频谱持续变化,测量的能量汇集的检测器信号可能与到达检测器的量子数目不完全成比例。因为量子的数目确定量子噪声,所以考虑频谱变化是很重要的,频谱变化可根据下面的等式9收集。
此外,为了考虑根据测量的X射线分布图的重建过程中涉及的对数变换,等式9包含通过归一化的初级辐射的倒数值计算的噪声增强因子。该因子表示对数函数的局部斜率,并且可以由1阶泰勒(taylor)展开式导出。
在该等式中,nP(E)和n0(E)是撞击到检测装置6上的频谱光子计数率密度,N0是在患者体内无衰减地到达检测器单元的光子的数目,<EP>和<E0>是在患者体内衰减前和衰减后的频谱的平均能量,以及是在患者体内衰减后的频谱的均方能量。对于等式9的实际实施,σP(IP,频谱)的值可以利用查找表(LUT)或利用参数计算方案进行存储。
如果在重建之前使用散射偏差校正技术或测量,可以充分考虑散射辐射。在存在散射辐射的情况下,输入噪声可以类似于等式9根据下式计算:
同样对于这种情况,可以利用LUT预先计算和存储函数σPS(IP,IS,频谱)的值或通过参数模型计算函数σPS(IP,IS,频谱)的值。
在存在散射辐射的情况下,IP和IS可以根据测量的检测器信号I测量=IP+IS确定。这可以通过简单的散射模型或基于散射估计算法(例如,与校正用于重建的投影数据所使用的相同的算法)实现。
应当认识到,正确考虑在测量的检测器信号中得到的散射辐射的量有利于获取真正最优的剂量调制曲线。如在例如EP 1 172 069 A1中描述的已知技术提出的、简单地使用基于衰减的量以导出输入噪声的分布可能导致次优的结果。
可通过简化以上噪声传播公式加速剂量的计算或确定。在函数gG具有短横向范围(类似于delta峰值)的情况下,这在典型CT和基于平面检测器的锥形射束CT系统中是常见的情况,噪声传播公式可以近似为:
&sigma; 3 , i 2 ( x j ) = a ij 2 &CenterDot; G G 0 2 &CenterDot; &sigma; 1 2 ( x 0 , ij ) + b ij 2 &CenterDot; G G 0 2 &sigma; 1 2 ( x 0 , ij + 1 ) + a ij 2 &CenterDot; b ij 2 &CenterDot; G G 1 2 &CenterDot; ( &sigma; 1 2 ( x 0 , ij ) + &sigma; 1 2 ( x 0 , ij + 1 ) ) - - - ( 11 )
使用
G G 0 2 = &Sigma; x = - &infin; &infin; g G 2 ( x ) G G 1 2 = &Sigma; x = - &infin; &infin; g G ( x ) &CenterDot; g G ( x + 1 ) 使得不再需要卷积滤波。可以预先计算
然后仅有对应于ROI投影(在极端情况下,每个投影仅有两个检测器单元)的那些投影值需要包括到计算中。
利用适当的比例因子,可以使用比图像数据的最终拟重建更为粗糙的空间分辨率完成噪声传播,大大减少了输入数据的量,并加速了最优剂量调制的计算。
在下文中,基于图2指示的备选噪声分布源描述了本发明的各种实施例。
在使用追踪扫描SS的情况下,对于CT中的扫描计划而言,通常对患者执行低剂量低分辨率的追踪扫描,以便确定例如待成像的FOV。使用该追踪扫描投影数据和在重建区域内由用户指定的ROI,可以使用等式2、8和9、使用散射校正算法和使用用于加速的低分辨率数据来计算最优剂量分布图。
使用具有诸如身体尺寸、体重、成像区域的参数的患者模型PM,预计算的数据可以用于等式8中的分母然后,可以使用等式2、8和10计算最优剂量分布图,在此,对于独立投影i的管值计算而言,可以使用来自目前采集的所有投影i-n(n=1...10)的投影的输入数据。在这种情况下,可以通过使用经减小的空间分辨率和经简化的噪声传播公式两者来加速计算。
对相同的几何结构进行重复采集的情况下使用先前的扫描数据ACQ,可以从先前的扫描中获取所需的信息。然后可以使用与追踪扫描SS相同的程序。
总之,已经描述了根据测量数据生成图像的图像生成装置和方法,其中,针对感兴趣区域优化图像质量,并且其中,在感兴趣区域的投影域中确定噪声分布,以及基于确定的噪声分布通过使用噪声分布算法确定针对辐射源的剂量分布图。因此,可以特别针对选定的感兴趣区域改善重建体积的信噪比。
所提出的剂量控制程序可以应用于螺旋和轴向CT两者以及基于平板检测器的锥形射束CT系统,诸如三维螺旋血管造影术(3DRA)或三维螺旋X射线(3DRX)产品。
在以上实施例中,成像系统为计算机断层摄影系统,但在其他实施例中,成像系统还可以是任何其他成像模式,例如核成像系统。
本发明可以用于生成活体对象的图像,如患者,但本发明还可以用于生成技术对象的图像。
本领域的技术人员根据对附图、公开和所附权利要求书的研究,在实践权利要求的发明的过程中,可以理解并实现所公开实施例的其他变型。
在权利要求中,“包括”一词并不排除其他元件或步骤,并且不定冠词“一”或“一个”并不排除多个。
单个单元或装置可以实现权利要求中叙述的若干项功能。基本的事实是,在相互不同的从属权利要求中叙述的特定测量并不表明不可以有利地使用这些测量的组合。
根据上述程序对剂量调制的计算、确定和/或控制可以作为计算机程序的程序代码模块和/或作为专用硬件来实施。
计算机程序可以存储/发布在合适的介质上,诸如光存储介质或固态介质,与其他硬件一起或作为其他硬件的一部分提供,但计算机程序还可以以其他形式发布,诸如经由因特网或其他有线或无线电信系统。
权利要求书中的任何附图标记不应解释为对范围的限制。

Claims (14)

1.一种用于根据测量数据生成图像的图像生成装置,其中,针对要被重建的体积中的感兴趣区域优化所述图像的图像质量,所述图像生成装置包括:
-噪声确定单元,其用于在投影域中基于所述感兴趣区域的附加数据确定针对所述感兴趣区域的独立投影的预期的噪声分布;以及
-剂量控制单元,其用于基于所述确定的所述独立投影的噪声分布使用噪声传播算法确定针对所述感兴趣区域的剂量分布图,
其中,所述剂量控制单元基于所述剂量分布图控制辐射源的供电电源,以便控制在断层摄影成像期间由所述辐射源生成的辐射强度,所述断层摄影成像利用成像系统生成所述测量数据;并且
其中,所述图像生成装置处理所述测量数据并根据所述测量数据生成所述图像。
2.如权利要求1所述的图像生成装置,其中,由用户指定所述感兴趣区域。
3.如权利要求1所述的图像生成装置,其中,利用自动或半自动的器官分割选取所述感兴趣区域。
4.如权利要求1所述的图像生成装置,其中,所述附加数据为之前采集的投影数据。
5.如权利要求1所述的图像生成装置,其中,所述附加数据为包括所述感兴趣区域的对象的模型。
6.如权利要求1所述的图像生成装置,其中,所述附加数据为追踪扫描。
7.如权利要求6所述的图像生成装置,其中,通过利用低剂量和低分辨率中的至少一个的扫描操作获取所述追踪扫描。
8.如权利要求1所述的图像生成装置,其中,所述噪声传播算法考虑了频谱效应和散射辐射中的至少一个。
9.如权利要求1所述的图像生成装置,其中,所述噪声确定单元适于基于低于所述图像的最终分辨率的经降低的空间分辨率和经简化的噪声传播算法中的至少一个确定所述噪声分布。
10.一种成像系统,包括:
-测量数据生成单元,其用于生成所述图像的所述测量数据,以及
-图像生成装置,如权利要求1中限定的、其用于根据所述测量数据生成所述图像。
11.一种用于根据测量数据生成图像的图像生成方法,其中,针对要被重建的体积中的感兴趣区域优化所述图像的图像质量,所述图像生成方法包括以下步骤:
-在投影域中基于所述感兴趣区域的附加数据确定针对所述感兴趣区域的独立投影的预期的噪声分布;
-基于所述确定的所述独立投影的噪声分布使用噪声传播算法确定针对所述感兴趣区域的剂量分布图;
-基于所述剂量分布图控制在生成所述测量数据的断层摄影成像期间生成的辐射强度;并且
-处理所述测量数据以根据所述测量数据生成所述图像。
12.如权利要求11所述的用于根据测量数据生成图像的图像生成方法,其中,所述测量数据是通过测量数据生成单元生成的。
13.一种用于根据测量数据生成图像的装置,其中,针对要被重建的体积中的感兴趣区域优化所述图像的图像质量,所述装置包括:
-用于在投影域中基于所述感兴趣区域的附加数据确定针对所述感兴趣区域的独立投影的预期的噪声分布的模块;
-用于基于所述确定的所述独立投影的噪声分布使用噪声传播算法确定针对所述感兴趣区域的剂量分布图的模块;
-用于基于所述剂量分布图控制在生成所述测量数据的断层摄影成像期间生成的辐射强度的模块;以及
-用于处理所述测量数据以根据所述测量数据生成所述图像的模块。
14.如权利要求13所述的用于根据测量数据生成图像的装置,其中,所述测量数据是通过测量数据生成单元生成的。
CN200980123925.2A 2008-06-25 2009-06-23 优化剂量控制的图像生成装置 Active CN102067177B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP08159000 2008-06-25
EP08159000.2 2008-06-25
PCT/IB2009/052686 WO2009156943A1 (en) 2008-06-25 2009-06-23 Image generation device with optimized dose control

Publications (2)

Publication Number Publication Date
CN102067177A CN102067177A (zh) 2011-05-18
CN102067177B true CN102067177B (zh) 2015-05-20

Family

ID=41119938

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200980123925.2A Active CN102067177B (zh) 2008-06-25 2009-06-23 优化剂量控制的图像生成装置

Country Status (5)

Country Link
US (1) US8712121B2 (zh)
EP (1) EP2308028B1 (zh)
CN (1) CN102067177B (zh)
BR (1) BRPI0910207A2 (zh)
WO (1) WO2009156943A1 (zh)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8708561B2 (en) 2009-03-20 2014-04-29 Orthoscan, Inc. Mobile imaging apparatus
EP2583241A1 (en) * 2010-06-21 2013-04-24 Koninklijke Philips Electronics N.V. Method and system for noise reduction in low dose computed tomography
CN102456228B (zh) * 2010-10-29 2015-11-25 Ge医疗系统环球技术有限公司 图像重建方法和装置及ct机
WO2012082799A1 (en) 2010-12-13 2012-06-21 Orthoscan, Inc. Mobile fluoroscopic imaging system
EP2661736B1 (en) * 2011-01-06 2014-04-30 Koninklijke Philips N.V. Imaging system for imaging an object
KR102060659B1 (ko) * 2013-03-20 2019-12-30 삼성전자주식회사 영상 처리를 위한 투사 및 역투사 방법 및 그 영상 처리 장치
US9406128B2 (en) * 2013-04-24 2016-08-02 Koninklijke Philips N.V. X-ray dose distribution calculation for a computed tomography examination
EP3077993A1 (en) * 2013-12-04 2016-10-12 Koninklijke Philips N.V. Image data processing
US10278666B2 (en) 2014-06-13 2019-05-07 Palodex Group Oy Systems and methods of automated dose control in x-ray imaging
CN104116518B (zh) * 2014-06-23 2016-06-01 沈阳东软医疗系统有限公司 一种剂量优化扫描方法及装置
US10932730B2 (en) * 2015-12-17 2021-03-02 Koninklijke Philips N.V. Method for estimating the radiation dose received by an organ during a computed tomography scan
US11276208B2 (en) 2017-06-12 2022-03-15 The Research Foundation For The State University Of New York System, method and computer-accessible medium for ultralow dose computed tomography image reconstruction
EP3420908A1 (en) * 2017-06-27 2019-01-02 Koninklijke Philips N.V. Scattered x-ray detection to detect misuse and prevent harm
WO2019090533A1 (en) * 2017-11-08 2019-05-16 Shenzhen United Imaging Healthcare Co., Ltd. System and method for diagnositic and treatment
DK3773213T3 (da) * 2018-04-10 2023-03-27 Bayer Healthcare Llc Fleksibel dosisestimering med brugerdefinerede volumener
CN109077746B (zh) * 2018-06-29 2022-08-02 上海联影医疗科技股份有限公司 一种确定辐射剂量调制线的方法、系统和装置
US10888296B2 (en) 2018-06-29 2021-01-12 Shanghai United Imaging Healthcare Co., Ltd. Methods and systems for modulating radiation dose
CN110796693B (zh) * 2019-09-11 2023-03-21 重庆大学 一种工业ct切片图像直接生成二维有限元模型的方法
CN115778412B (zh) * 2023-02-09 2023-04-28 之江实验室 X光声成像中造影剂剂量的优化方法、装置及存储介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1496713A (zh) * 2002-09-26 2004-05-19 西门子公司 用于计算机断层扫描中自动曝射控制的方法和设备

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5396531A (en) * 1992-11-05 1995-03-07 General Electric Company Method of achieving reduced dose X-ray fluoroscopy by employing statistical estimation of poisson noise
US5379333A (en) 1993-11-19 1995-01-03 General Electric Company Variable dose application by modulation of x-ray tube current during CT scanning
US5400378A (en) * 1993-11-19 1995-03-21 General Electric Company Dynamic dose control in multi-slice CT scan
EP1172069A1 (de) 2000-07-14 2002-01-16 VAMP Verfahren und Apparate der Medizinischen Physik GmbH Computertomograph mit Dosisoptimierung durch Festlegung der optimalen Wahl des Röhrenstroms in Echtzeit (Belichtungsautomatik), der Röhrenstrommodulation (Dosisminimierung) und darauf aufbauender Nachverarbeitung durch 3D adaptive Filter (Rauschreduzierung
ATE373287T1 (de) * 2001-05-16 2007-09-15 Koninkl Philips Electronics Nv Computertomographiesystem mit filterung
US6493416B1 (en) 2001-11-21 2002-12-10 Ge Medical Systems Global Technology Company, Llc Method and apparatus for noise reduction in computed tomographic systems
SE527138C2 (sv) * 2003-07-08 2005-12-27 Xcounter Ab Skanningsbaserad detektering av joniserande strålning för tomosyntes
JP2007529738A (ja) * 2004-03-17 2007-10-25 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 干渉性散乱ctのビーム硬化補正および減衰補正
SE528234C2 (sv) * 2004-03-30 2006-09-26 Xcounter Ab Anordning och metod för att erhålla tomosyntesdata
EP1782094A2 (en) * 2004-07-23 2007-05-09 ANGELSEN, Bjorn A. J. Ultrasound imaging using non-linear manipulation of forward propagation properties of a pulse
US7570732B2 (en) * 2005-11-09 2009-08-04 Dexela Limited Methods and apparatus for obtaining low-dose imaging
US20070147579A1 (en) * 2005-12-23 2007-06-28 De Man Bruno K B Method and system for radiographic imaging with organ-based radiation profile prescription

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1496713A (zh) * 2002-09-26 2004-05-19 西门子公司 用于计算机断层扫描中自动曝射控制的方法和设备

Also Published As

Publication number Publication date
CN102067177A (zh) 2011-05-18
WO2009156943A1 (en) 2009-12-30
EP2308028B1 (en) 2017-04-26
US8712121B2 (en) 2014-04-29
EP2308028A1 (en) 2011-04-13
US20110110573A1 (en) 2011-05-12
BRPI0910207A2 (pt) 2015-09-29

Similar Documents

Publication Publication Date Title
CN102067177B (zh) 优化剂量控制的图像生成装置
CN101158653B (zh) 一种锥束ct系统的散射测定和校正方法
US9265475B2 (en) Methods and apparatus for scatter correction for CBCT system and cone-beam image reconstruction
Speidel et al. Scanning‐beam digital x‐ray (SBDX) technology for interventional and diagnostic cardiac angiography
US7298882B2 (en) Generalized measure of image quality in medical X-ray imaging
US7734076B2 (en) Material decomposition image noise reduction
CN101346102B (zh) X射线ct装置
US7826585B2 (en) Stereo tube computed tomography
Maltz et al. Algorithm for X-ray scatter, beam-hardening, and beam profile correction in diagnostic (kilovoltage) and treatment (megavoltage) cone beam CT
US8054940B2 (en) Image capturing system for medical use, image processing apparatus, and image processing method
US7668361B2 (en) System and method for the correction of temporal artifacts in tomographic images
US9261467B2 (en) System and method of iterative image reconstruction for computed tomography
EP3671648A1 (en) System and method for statistical iterative reconstruction and material decomposition
CN102697516B (zh) 用于建立断层造影图像显示的方法和计算机断层造影系统
US9589373B2 (en) Monte carlo modeling of field angle-dependent spectra for radiographic imaging systems
CN108471996B (zh) 一种用于估计由器官接收的辐射剂量的方法和成像设备
US11288775B2 (en) Methods and systems for parametric noise modulation in x-ray imaging
CN114387359A (zh) 一种三维x射线低剂量成像方法及装置
US20100054562A1 (en) Systems and methods for adaptive filtering
US10552992B2 (en) Poly-energetic reconstruction method for metal artifacts reduction
Wiegert Scattered radiation in cone beam computed tomography: analysis, quantification and compensation
Bernhardt et al. Spatial frequency-dependent signal-to-noise ratio as a generalized measure of image quality
Bornefalk et al. Dual-energy imaging using a photon counting detector with electronic spectrum-splitting
Saha et al. Single photon emission computed tomography
Forthmann et al. Adaptive two‐pass cone‐beam artifact correction using a FOV‐preserving two‐source geometry: A simulation study

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