[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN105326519A - X射线图像形成、散射成分计算以及重建的方法及装置 - Google Patents

X射线图像形成、散射成分计算以及重建的方法及装置 Download PDF

Info

Publication number
CN105326519A
CN105326519A CN201510836614.8A CN201510836614A CN105326519A CN 105326519 A CN105326519 A CN 105326519A CN 201510836614 A CN201510836614 A CN 201510836614A CN 105326519 A CN105326519 A CN 105326519A
Authority
CN
China
Prior art keywords
image
ray
modulation
collimator
radioscopic
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
CN201510836614.8A
Other languages
English (en)
Other versions
CN105326519B (zh
Inventor
李贵
杨宏成
徐璐
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.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
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 Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201510836614.8A priority Critical patent/CN105326519B/zh
Publication of CN105326519A publication Critical patent/CN105326519A/zh
Priority to US15/361,036 priority patent/US10898154B2/en
Application granted granted Critical
Publication of CN105326519B publication Critical patent/CN105326519B/zh
Priority to US17/156,655 priority patent/US11724128B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048—Monitoring, verifying, controlling systems and methods
    • A61N5/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
    • 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/06—Diaphragms
    • 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/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • 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/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • A61B6/5282—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/70—Denoising; Smoothing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/73—Deblurring; Sharpening
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4035—Arrangements for generating radiation specially adapted for radiation diagnosis the source being combined with a filter or grating
    • 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/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/54—Control of apparatus or devices for radiation diagnosis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048—Monitoring, verifying, controlling systems and methods
    • A61N5/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
    • A61N2005/1054—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam using a portal imaging system
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N2005/1092—Details
    • A61N2005/1095—Elements inserted into the radiation path within the system, e.g. filters or wedges
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1042—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head
    • A61N5/1045—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head using a multi-leaf collimator, e.g. for intensity modulated radiation therapy or IMRT
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10116—X-ray image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20212—Image combination
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20212—Image combination
    • G06T2207/20224—Image subtraction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21K—TECHNIQUES FOR HANDLING PARTICLES OR IONISING RADIATION NOT OTHERWISE PROVIDED FOR; IRRADIATION DEVICES; GAMMA RAY OR X-RAY MICROSCOPES
    • G21K1/00—Arrangements for handling particles or ionising radiation, e.g. focusing or moderating
    • G21K1/02—Arrangements for handling particles or ionising radiation, e.g. focusing or moderating using diaphragms, collimators
    • G21K1/04—Arrangements for handling particles or ionising radiation, e.g. focusing or moderating using diaphragms, collimators using variable diaphragms, shutters, choppers
    • G21K1/043—Arrangements for handling particles or ionising radiation, e.g. focusing or moderating using diaphragms, collimators using variable diaphragms, shutters, choppers changing time structure of beams by mechanical means, e.g. choppers, spinning filter wheels

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Biomedical Technology (AREA)
  • Physics & Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Optics & Photonics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Pulmonology (AREA)
  • Physiology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Quality & Reliability (AREA)

Abstract

本发明公开了一种X射线图像成像、X射线图像散射成分计算、X射线图像重建的方法及装置，X射线图像成像方法包括：采集依次经过准直器和成像对象后的X射线，形成X射线图像组；所述采集是在所述准直器沿第一方向运动时进行的，所述X射线图像组为一组在准直器运动的不同时刻采集得到的X射线图像；提取所述X射线图像组内各X射线图像的子区域图像；沿所述第一方向将所述子区域图像进行组合，形成调制X射线图像。本发明技术方案实现使用准直器对X射线强度分布进行灵活的调制，不增加额外硬件，对X射线图像中的散射成分进行计算，最终减少X射线图像的散射量。

Description

X射线图像形成、散射成分计算以及重建的方法及装置
技术领域
本发明涉及医学图像领域，特别涉及一种X射线图像的形成、散射成分计算以及重建的方法及装置。
背景技术
在X射线图像中，散射线一直都是导致图像质量不理想的一个重要因素。特别是对于有较大探测器面积的X射线系统来说，例如放射治疗设备，用到的锥束CT(conebeamcomputedtomography)采用的是电子射野影像装置(EPID，ElectronicPortalImagingDevice)作为接收X射线的探测器，在其接收到的X射线中散射线占了很大的比例。散射线在图像上引起杯状伪影，降低了图像的对比度，极大降低了重建图像的图像精度。在放射治疗过程中，锥束CT用来实时获得病人图像，对放射治疗进行图像引导，图像精度的降低导致了放射治疗的位置误差，无法达到对肿瘤区域精确放疗的目的。
目前对图像中散射线的矫正方法主要有两种，一种是硬件矫正，是指在X射线影像系统中添加若干硬件工具，减少达到探测器的散射线，从而达到散射矫正目的。常用的硬件矫正工具包括X射线准直器、防散射滤线栅等。第二种是软件矫正方法，是指对已经采集到的X射线投影图像通过对投影图像的分析和被照射物体的估计，得到一个散射分布图，由散射分布图对X射线投影图像进行散射矫正。
发明内容
本发明要解决的问题是提供一种减少散射线、提高X射线图像质量的方法和装置。
为解决上述问题，本发明提供了一种调制X射线图像形成方法，包括：采集依次经过准直器和成像对象后的X射线，形成X射线图像组；所述采集是在所述准直器沿第一方向运动时进行的，所述X射线图像组为一组在准直器运动的不同时刻采集得到的X射线图像；提取所述X射线图像组内各X射线图像的子区域图像；沿所述第一方向将所述子区域图像进行组合，形成调制X射线图像。
优选的，所述采集得到X射线图像组的采集频率和所述准直器的运动速度为正比例关系。
优选的，所述准直器沿第一方向运动包括：所述准直器为多叶准直器，包括一组准直器叶片；所述一组准直器叶片长度相同，共同沿第一方向运动。
为解决上述问题，本发明还提供了一种调制X射线图像的形成方法，其特征在于，包括：采集依次经过准直器和成像对象后的X射线，形成X射线图像；所述采集是在所述准直器沿第一方向运动时进行的；所述准直器包括一对在第一方向上间隔设置的准直器；在所述运动过程中，所述准直器保持间隔距离以周期性速度运动。
优选的，所述以周期性速度运动包括：速度以周期性匀速增加、匀速减少的方式运动。
为解决上述问题，本发明还提供了一种X射线图像散射成分的计算方法，包括：根据上述调制X射线图像的形成方法获得调制X射线图像；对所述调制X射线图像进行低通滤波，得到所述调制X射线图像的低频部分和散射成分；对所述调制X射线图像进行高通滤波，得到所述调制X射线图像的高频部分；对所述调制X射线图像的高频部分解调制后结合权重，计算得到所述调制X射线低频部分的估计值；将经所述低通滤波后得到的调制X射线图像的低频部分和散射成分，减去所述调制X射线低频部分的估计值，计算得到散射成分。
为解决上述问题，本发明还提供了一种X射线图像的重建方法，包括：采集得到成像对象各投影角度下的投影图像；根据上述计算方法计算得到所述投影图像的散射成分；将所述投影图像去掉所述散射成分后进行图像重建，得到重建图像。
为解决上述问题，在本发明调制X射线图像形成方法、X射线图像散射成分计算方法以及X射线图像重建方法的基础上，还提供了对应的装置。
与现有技术相比，本发明技术方案实现使用准直器对X射线强度分布进行灵活的调制，不增加额外硬件，形成明暗相间的条纹状X射线图像，以便对X射线图像中的散射成分进行计算，最终减少X射线图像的散射量。更进一步的，由于准直器安装在治疗头，在成像时机械稳定性好，不会产生因几何位置抖动引起的图像运动模糊。
附图说明
图1是一种放射治疗系统的结构图；
图2是放射治疗系统治疗头的结构图；
图3是放射治疗系统治疗头内次级准直器的结构图；
图4是本发明X射线图像散射成分计算方法的流程图；
图5a-5d是本发明次级准直器运动时的状态示意图；
图6a、6b是本发明由X射线图像组形成调制X射线图像的示意图；
图7是本发明对投影图像进行去散射处理的示意图；
图8是放射治疗系统的多叶准直器的结构图；
图9是本发明X射线图像重建的流程图；
图10是本发明X射线图像重建装置的结构图；
图11是本发明形成调制X射线图像的示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更为明显易懂，下面结合附图对本发明的具体实施方式做详细的说明。在以下描述中阐述了具体细节以便于充分理解本发明。但是本发明能够以多种不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似推广。因此本发明不受下面公开的具体实施方式的限制。
图1是一种放射治疗系统的结构图，如图1所示，放射治疗系统100包括固定部分101和旋转部分102，旋转部分102安装在固定部分101上，旋转部分101可以绕中心轴106进行旋转，从而实现在不同角度对患者进行放射治疗。
旋转部分102的一侧是治疗头103，治疗头103可以产生高能级的X射线(通常为兆伏级)，对在病床105上的患者进行放射治疗。对于同源双束的放射治疗系统，治疗头103还可以产生低能级的X射线(通常为千伏级)，低能级X射线可以用来对患者进行成像，利用得到的患者图像对患者进行图像引导放射治疗(ImageGuidedRadiationTherapy)。
在由低能级X射线进行成像时，治疗头103发出锥束X射线，旋转部分102另一侧的EPID104接收到穿过患者的X射线，形成该角度下的投影图像(projectionimage)。治疗头103在不同角度照射，从而形成多个角度的投影图，再对投影图进行图像重建(imagereconstruction)得到患者的CT(computedtomography)图像。
治疗头103内部结构如图2所示，其通常包括：X射线靶(X-raytarget)、初级准直器(primarycollimator)、均整器(flatteningfilter)、电离室(ionchamber)、次级准直器(secondarycollimator)、楔形块(wedge)等。首先，直线加速器(图中未示出)发出的电子束(electronbeam)打到X射线靶上，产生X射线。X射线经过初级准直器得到方向、范围一定的X射线，之后经过均整器对X射线进行均整，由电离室对X射线的能量进行测量，次级准直器进一步对X射线的强度和/或形状进行调整，最后由楔形块对X射线的能量分布进行调整，最终X射线从治疗头103的出口发出。
图3是沿X射线发射方向看到的次级准直器结构图，其具体由两组次级准直器构成，每组次级准直器又由两块在某一方向上相对设置的次级准直器构成。X射线可以不被遮挡的穿过区域205，或者穿过由准直器覆盖的区域从而射线强度被衰减。其中，一组次级准直器201、次级准直器202相对设置，均可由电机控制沿a方向移动。另一组次级准直器203、次级准直器204也是相对设置，可沿b方向运动。
结合以上描述的结构，下面对本发明X射线图像散射成分的计算方法进行说明，参见图4，其具体包括以下步骤：
S401，获得调制X射线图像。
具体获得调制X射线图像的方法结合图5a-5d和图6a-6b进行说明。图5a是由X射线进行成像的示意图，X射线502由X射线源501(可以认为是X射线靶)发出。在X射线502穿过成像对象(图中未示出，通常为患者，也可以是动物或者水模)前，通过在a方向上运动的次级准直器(201、202)对X射线502的射线强度分布进行调制。X射线502的发射方向和a方向大致垂直，因为X射线502本身有一定的发射角度，而且通常也不用将次级准直器的运动方向做成和X射线发射方向完全垂直。
在本发明技术方案的一种实施方式中，参见图5a、图5b，由于次级准直器结构是由两组次级准直器构成，此时一组次级准直器203、次级准直器204(未在图中示出)全开野。全开野指的是准直器完全打开，不在X射线的发射路径上对X射线产生遮蔽。另一组次级准直器中，次级准直器201全开野，次级准直器202沿a方向远离次级准直器201运动。
在图5a状态下，探测器104上形成的如503所示的X射线图像，其中黑色表示由于次级准直器对X射线有遮挡，图像强度相对较暗。而白色表示由于未被准直器遮挡，图像强度相对较亮。并且由于次级准直器端处206(参见图3)的截面为矩形，因此在X射线图像503内的仅有一条窄亮条。
当次级准直器202运动到图5b中的位置时，探测器104采集到的如504所示的X射线图像。由于次级准直器202沿a方向远离次级准直器201运动了一段距离，因此图像上的亮条的宽度增加。
为了形成最终图像，在次级准直器202沿a方向(也可称为第一方向)运动的过程中，采集得到多张X射线图像形成X射线图像组，每张X射线图像对应次级准直器运动过程的不同时刻采集到的X射线图像。之后在每张X射线图像中对应移动准直器的边缘区域附近，提取出一个子区域图像，子区域图像包含一个亮条纹和一个暗条纹。例如对于图像503，得到子区域图像601，子区域图像601包含一个亮条纹6011和一个暗条纹6012。对于图像504，得到子区域图像602。各子区域图像601、602、603从a方向上进行组合，得到最终的X射线图像600。X射线图像600是明暗相间条纹的图像，其是在普通X射线图像的基础上通过准直器的移动和图像处理方法进行了调制，便于后续对图像的散射成分进行计算，在这里称为调制X射线图像。
采集X射线图像组的采集频率可以根据准直器的移动速度，以及探测器的宽度、最终调制X射线图像内明暗相间条纹的间隔距离来确定。准直器移动速度和采集频率成正比例关系，准直器移动速度越快，则需要的采集频率越高。
在本发明技术方案的另一种实施过程中，如图5c、图5d所示，次级准直器201、次级准直器202在a方向上保持固定距离205，共同沿a方向运动。为了形成调制X射线图像，如图6b所示，从各X射线图像505、X射线图像506中分别提取得到子区域图像604、子区域图像605。将多个子区域图像(604、605、606等)从a方向上进行组合，形成最终的X射线图像600。
在本发明技术方案的另一种实施过程中，次级准直器201、次级准直器202在a方向上保持固定距离205，以周期性速度共同沿a方向运动，沿a方向从左侧运动到右侧，如图11所示。由于单位时间内X射线的的脉冲数目一定，准直器运动速度较慢时，对应的探测器区域所形成的X射线图像上，X射线的脉冲数目多；而准直器运动较快时，对应的探测器区域所形成的X射线图像上，X射线的脉冲数目少。因此在准直器运动过程中进行采集，可以得到一张明暗相间条纹的调制X射线图像600。优选的，速度以周期性的匀速增加、匀速较少的方式。优选的，运动速度还可以为一定比例的匀速进行，例如以速度为3v匀速运动时间t，之后以速度为v匀速运动时间3t，之后再以速度为3v匀速运动时间t，这种周期性的方式运动，形成调制X射线图像。
若采用静态采集方式，与亮条纹投影区域相比，暗条纹投影区域的泊松量子噪声增大，图像信噪比较低。而在准直器运动过程中进行图像采集形成调制X射线图像，图像暗条纹部分信噪比的明显提高。
之后的步骤S402至步骤S405可以参考公开号为US7463712B2的美国专利说明书记载内容。以下结合图7，仅作简单描述说明。
对X射线图像600在计算时，可以对每个明暗条纹的中心线处进行下采样，避免图像中准直器边缘处的边缘效应对图像上X射线强度产生的影响。
S402，对调制X射线图像进行低通滤波，得到经调制X射线图像的低频部分和散射成分。
由于X射线图像中散射成分的特性，步骤S401获取的调制X射线图像中散射成分主要分布在图像数据的低频部分。当对投影图701(这里的投影图是指一定投影角度下的调制X射线图像600)进行低通滤波后，可以得到图像的低频部分和散射成分之和702的图像数据。
S403，对所述调制X射线图像进行高通滤波，得到所述调制X射线图像的高频部分。
而对投影图701进行高通滤波，由于散射成分几乎不包含在调制X射线图像的高频部分，可以得到图像不包含散射成分的高频部分704。
S404，对所述高频部分进行解调制并结合权重，计算得到所述调制X射线低频部分的估计值。
将高频部分数据704解调制后乘以权重系数其中∝为准直器对X射线的透射系数，从而计算得到低频部分的估计值705。
S405，将所述低通滤波后得到的调制X射线图像的低频部分和散射成分，减去所述调制X射线低频部分的估计值，计算得到散射成分。
由步骤S402低通滤波后得到的低频部分和散射成分702之和，减去步骤S404计算得到的低频部分估计，得到散射成分估计706。
在本发明的一种技术方案中，可以采用多叶准直器替代次级准直器。多叶准直器结构如图8所示，多叶准直器800包括多组相对设置的叶片，叶片组801和叶片组802相对设置，当进行对X射线调制时，一侧叶片组的长度可变为相同，从而实现类似次级准直器的作用。
当采用多叶准直器时，由于多叶准直器的透射率几乎为0，因此步骤S404中的权重系数为1。
本发明还提供了一种X射线图像的重建方法，如图9所示，包括：
S901，采集得到各投影角度下的X射线投影图像。
如图1所示的放射治疗设备在进行成像时，旋转机构102旋转，治疗头103发出低能级X射线，探测器104接收多个旋转角度下成像对象的X射线投影图像，此时采集到的投影图像是常用X射线成像方法采集得到的投影图像。
S902，由上述X射线图像散射成分的计算方法计算得到投影图像的散射成分。
由前述的X射线图像散射成分计算方法，首先采集到成像对象的调制X射线图像，之后计算得到调制X射线图像的散射成分，由于成像对象相同，散射情况也相同，所以该方法计算得到的散射成分也是步骤S901得到的X射线投影图像的散射成分。
S903，将所述投影图像去掉所述散射成分后进行图像重建，得到重建图像。
将步骤S901得到的投影图像减去步骤S902计算得到的散射成分，完成对投影图像的散射去除，之后对投影图像进行图像重建，重建得到成像对象的CT图像。具体图像重建过程为本领域技术人员已知，这里不再赘述。
对应于上述的调制X射线图像形成方法本发明还提供了一种调制X射线图像的形成装置，参见图10，调制X射线图像的形成装置1010包括：
第一采集单元1011，用于采集依次经过准直器和成像对象后的X射线，形成X射线图像组；所述采集是在所述准直器沿第一方向运动时进行的，所述X射线图像组为一组在准直器运动的不同时刻采集得到的X射线图像。
提取单元1012，用于提取所述X射线图像组内各X射线图像的子区域图像。
组合单元1013，用于沿所述第一方向将所述子区域图像进行组合，形成调制X射线图像。
在上述X射线图像散射成分计算方法的基础上，本发明还提供了一种X射线图像散射成分计算装置，参见图10，X射线图像散射成分计算装置1020包括：
调制X射线图像形成装置1010，用于获得调制X射线图像；
低通滤波单元1021，用于对所述调制X射线图像进行低通滤波，得到所述调制X射线图像的低频部分和散射成分；
高通滤波单元1022，用于对所述调制X射线图像进行高通滤波，得到所述调制X射线图像的高频部分；
低频部分估计单元1023，用于对所述调制X射线图像的高频部分解调制后结合权重，计算得到所述调制X射线低频部分的估计值；
散射计算单元1024，用于将经所述低通滤波后得到的调制X射线图像的低频部分和散射成分，减去所述调制X射线低频部分的估计值，计算得到散射成分。
在本发明X射线图像重建方法的基础上，本发明还提供了一种X射线图像重建装置，参见图10，X射线图像重建装置1030包括：
第二采集单元1031，用于采集得到成像对象各投影角度下的投影图像；
射线图像散射成分计算单元1020，用于计算所述投影图像的散射成分；
重建单元1032，用于将所述投影图像去掉所述散射成分后进行图像重建，得到重建图像。
本发明调制X射线图像的形成装置、X射线图像散射成分计算装置以及X射线图像重建装置的实施方式可以分别参考对应方法的实施方式，这里不再赘述。
本发明技术方案实现使用准直器对X射线强度分布进行灵活的调制，不增加额外硬件，形成明暗相间的条纹状X射线图像，以便对X射线图像中的散射成分进行计算，最终减少X射线图像的散射量。更进一步的，由于准直器安装在治疗头，在成像时机械稳定性好，不会产生因几何位置抖动引起的图像运动模糊。
本发明虽然已以较佳实施例公开如上，但其并不是用来限定本发明，任何本领域技术人员在不脱离本发明的精神和范围内，都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案的保护范围。

Claims (10)

1.一种调制X射线图像的形成方法，其特征在于，包括：
采集依次经过准直器和成像对象后的X射线，形成X射线图像组；所述采集是在所述准直器沿第一方向运动时进行的，所述X射线图像组为一组在准直器运动的不同时刻采集得到的X射线图像；
提取所述X射线图像组内各X射线图像的子区域图像；
沿所述第一方向将所述子区域图像进行组合，形成调制X射线图像。
2.根据权利要求1所述的形成方法，其特征在于，还包括：所述采集得到X射线图像组的采集频率和所述准直器的运动速度为正比例关系。
3.根据权利要求1所述的形成方法，其特征在于，所述准直器沿第一方向运动包括：所述准直器为多叶准直器，包括一组准直器叶片；所述一组准直器叶片长度相同，共同沿第一方向运动。
4.一种调制X射线图像的形成方法，其特征在于，包括：
采集依次经过准直器和成像对象后的X射线，形成X射线图像；所述采集是在所述准直器沿第一方向运动时进行的；
所述准直器包括一对在第一方向上间隔设置的准直器；
在所述运动过程中，所述准直器保持间隔距离以周期性速度运动。
5.根据权利要求4所述的形成方法，其特征在于，所述以周期性速度运动包括：速度以周期性匀速增加、匀速减少的方式运动。
6.一种X射线图像散射成分的计算方法，其特征在于，包括：
根据权利要求1或4任一项所述的调制X射线图像的形成方法获得调制X射线图像；
对所述调制X射线图像进行低通滤波，得到所述调制X射线图像的低频部分和散射成分；
对所述调制X射线图像进行高通滤波，得到所述调制X射线图像的高频部分；
对所述调制X射线图像的高频部分解调制后结合权重，计算得到所述调制X射线低频部分的估计值；
将经所述低通滤波后得到的调制X射线图像的低频部分和散射成分，减去所述调制X射线低频部分的估计值，计算得到散射成分。
7.一种X射线图像的重建方法，其特征在于，包括：
采集得到成像对象各投影角度下的X射线投影图像；
根据权利要求6所述的计算方法计算得到所述投影图像的散射成分；
将所述投影图像去掉所述散射成分后进行图像重建，得到重建图像。
8.一种调制X射线图像的形成装置，其特征在于，包括：
第一采集单元，用于采集依次经过准直器和成像对象后的X射线，形成X射线图像组；所述采集是在所述准直器沿第一方向运动时进行的，所述X射线图像组为一组在准直器运动的不同时刻采集得到的X射线图像；
提取单元，用于提取所述X射线图像组内各X射线图像的子区域图像；
组合单元，用于沿所述第一方向将所述子区域图像进行组合，形成调制X射线图像。
9.一种X射线图像散射成分计算装置，其特征在于，包括：
权利要求8所述的调制X射线图像形成装置，用于获得调制X射线图像；
低通滤波单元，用于对所述调制X射线图像进行低通滤波，得到所述调制X射线图像的低频部分和散射成分；
高通滤波单元，用于对所述调制X射线图像进行高通滤波，得到所述调制X射线图像的高频部分；
低频部分估计单元，用于对所述调制X射线图像的高频部分解调制后结合权重，计算得到所述调制X射线低频部分的估计值；
散射计算单元，用于将经所述低通滤波后得到的调制X射线图像的低频部分和散射成分，减去所述调制X射线低频部分的估计值，计算得到散射成分。
10.一种X射线图像重建装置，其特征在于，包括：
第二采集单元，用于采集得到成像对象各投影角度下的投影图像；
权利要求9所述的X射线图像散射成分计算单元，用于计算所述投影图像的散射成分；
重建单元，用于将所述投影图像去掉所述散射成分后进行图像重建，得到重建图像。
CN201510836614.8A 2015-11-26 2015-11-26 X射线图像形成、散射成分计算以及重建的方法及装置 Active CN105326519B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201510836614.8A CN105326519B (zh) 2015-11-26 2015-11-26 X射线图像形成、散射成分计算以及重建的方法及装置
US15/361,036 US10898154B2 (en) 2015-11-26 2016-11-24 System and method for x-ray imaging
US17/156,655 US11724128B2 (en) 2015-11-26 2021-01-25 System and method for X-ray imaging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510836614.8A CN105326519B (zh) 2015-11-26 2015-11-26 X射线图像形成、散射成分计算以及重建的方法及装置

Publications (2)

Publication Number Publication Date
CN105326519A true CN105326519A (zh) 2016-02-17
CN105326519B CN105326519B (zh) 2017-03-22

Family

ID=55277444

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510836614.8A Active CN105326519B (zh) 2015-11-26 2015-11-26 X射线图像形成、散射成分计算以及重建的方法及装置

Country Status (2)

Country Link
US (2) US10898154B2 (zh)
CN (1) CN105326519B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106618625A (zh) * 2017-02-08 2017-05-10 李宝生 一种基于衰减板的数字化x线摄影系统
WO2019205924A1 (zh) * 2018-04-28 2019-10-31 西安大医集团有限公司 放射治疗头及放射治疗装置
WO2023179227A1 (zh) * 2022-03-25 2023-09-28 康达洲际医疗器械有限公司 一种基于散射识别的锥束计算机断层扫描成像方法与系统

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105326519B (zh) * 2015-11-26 2017-03-22 上海联影医疗科技有限公司 X射线图像形成、散射成分计算以及重建的方法及装置
US10799717B2 (en) * 2018-09-29 2020-10-13 Uih America, Inc. Systems and methods for reducing a radial distance of a collimator assembly occupying
KR102314902B1 (ko) * 2019-08-20 2021-10-21 한국원자력의학원 진단용 방사선 장치에 장착 가능한 선량 조절 장치 및 이를 포함하는 선량 조절 시스템

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070268997A1 (en) * 2006-05-18 2007-11-22 The Board Of Trustees Of The Leland Stanford Junior University Scatter correction for x-ray imaging using modulation of primary x-ray spatial spectrum
CN102036610A (zh) * 2008-05-21 2011-04-27 皇家飞利浦电子股份有限公司 飞越式扫描期间动态可调的源准直
CN102090900A (zh) * 2009-12-11 2011-06-15 通用电气公司 缓和双能量ct的低信号数据的系统和方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0657057A (ja) 1992-08-11 1994-03-01 Sumitomo Chem Co Ltd ポリプロピレン組成物及びそのフィルム
JPH11511376A (ja) * 1996-03-13 1999-10-05 アナロジック コーポレーション コンピュータ・トモグラフィの自己校正リング抑制フィルタ
US20040120457A1 (en) * 2002-12-20 2004-06-24 University Of Massachusetts Medical Center Scatter reducing device for imaging
JP3942178B2 (ja) * 2003-07-29 2007-07-11 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｘ線ｃｔシステム
US8989469B2 (en) * 2010-12-20 2015-03-24 The Board Of Trustees Of The Leland Stanford Junior University Systems and methods for simultaneous acquisition of scatter and image projection data in computed tomography
DE102011006660A1 (de) * 2011-04-01 2012-10-04 Siemens Aktiengesellschaft Verfahren und Vorrichtung zur Korrektur von Artefakten bei einer Röntgenbilderzeugung, insbesondere Computertomographie, mit bewegtem Modulatorfeld
US8867699B2 (en) * 2011-08-31 2014-10-21 Shimadzu Corporation Radiographic device
CN104783819B (zh) 2014-08-27 2017-10-03 上海联影医疗科技有限公司 散射校正方法及装置
CN105326519B (zh) * 2015-11-26 2017-03-22 上海联影医疗科技有限公司 X射线图像形成、散射成分计算以及重建的方法及装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070268997A1 (en) * 2006-05-18 2007-11-22 The Board Of Trustees Of The Leland Stanford Junior University Scatter correction for x-ray imaging using modulation of primary x-ray spatial spectrum
CN102036610A (zh) * 2008-05-21 2011-04-27 皇家飞利浦电子股份有限公司 飞越式扫描期间动态可调的源准直
CN102090900A (zh) * 2009-12-11 2011-06-15 通用电气公司 缓和双能量ct的低信号数据的系统和方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106618625A (zh) * 2017-02-08 2017-05-10 李宝生 一种基于衰减板的数字化x线摄影系统
WO2019205924A1 (zh) * 2018-04-28 2019-10-31 西安大医集团有限公司 放射治疗头及放射治疗装置
US12036424B2 (en) 2018-04-28 2024-07-16 Our United Corporation Radiation therapy head and radiation therapy apparatus
WO2023179227A1 (zh) * 2022-03-25 2023-09-28 康达洲际医疗器械有限公司 一种基于散射识别的锥束计算机断层扫描成像方法与系统

Also Published As

Publication number Publication date
US10898154B2 (en) 2021-01-26
US11724128B2 (en) 2023-08-15
US20170151444A1 (en) 2017-06-01
CN105326519B (zh) 2017-03-22
US20210137480A1 (en) 2021-05-13

Similar Documents

Publication Publication Date Title
CN105326519A (zh) X射线图像形成、散射成分计算以及重建的方法及装置
CN113164138B (zh) 用于在锥形束计算机断层摄影中进行散射估计的方法和装置
Sun et al. Improved scatter correction using adaptive scatter kernel superposition
EP3400876B1 (en) Automatic estimating and reducing scattering in computed tomography scans
US6950493B2 (en) Dynamic multi-spectral CT imaging
US8270561B2 (en) Motion weighting in computed tomography (CT) with cone angle
US10249063B2 (en) Method and apparatus for proton computerised tomography
Wu et al. Metal artifact correction for x‐ray computed tomography using kV and selective MV imaging
CN102764137B (zh) 一种静态ct扫描仪及其散射x光子校正方法
Dong et al. Low‐dose and scatter‐free cone‐beam CT imaging using a stationary beam blocker in a single scan: phantom studies
Szczykutowicz et al. Experimental realization of fluence field modulated CT using digital beam attenuation
CN106408543A (zh) 一种用于锥束ct图像散射修正的阻挡光栅优化方法及装置
KR101076321B1 (ko) 콘빔 ｃｔ 장치에서의 3차원 영상 획득 방법 및 이를 적용한 콘빔 ｃｔ 장치
Chen et al. Optimization of the geometry and speed of a moving blocker system for cone‐beam computed tomography scatter correction
US9545231B2 (en) Optimized aperture selection imaging computed tomography system and method
KR20190140641A (ko) 이진 이동 변조 차단기를 이용하여 듀얼 에너지 방식으로 동작하는 단일 소스 콘빔 전산화 단층 촬영 시스템
CN115297779A (zh) 用于图像散射校正的噪声和伪影减少
CN105427355A (zh) 罩壳、x射线图像的散射成分计算、重建的方法及装置
Montesdeoca Scattered radiation in projection X-ray mammography and digital breast tomosynthesis
CN104548374A (zh) 基于旋转准直器的cbct图像采集系统及其预处理方法
Ghazi A fluence modulation and scatter shielding apparatus for dedicated breast CT: Theory of operation
Deng et al. Multi‐energy blended CBCT spectral imaging and scatter‐decoupled material decomposition using a spectral modulator with flying focal spot (SMFFS)
CN110840479B (zh) 一种骨成分确定方法、装置、设备及存储介质
Cai Artefact Reduction Methods for Iterative Reconstruction in Full-fan Cone Beam CT Radiotherapy Applications
CN117561545A (zh) 锚定内核散射估计

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
CP01 Change in the name or title of a patent holder

Address after: 201807 No. 2258 Chengbei Road, Jiading Industrial Zone, Jiading District, Shanghai.

Patentee after: Shanghai Lianying Medical Technology Co., Ltd

Address before: 201807 No. 2258 Chengbei Road, Jiading Industrial Zone, Jiading District, Shanghai.

Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.

CP01 Change in the name or title of a patent holder