CN107845070B - 一种数字x射线透视剂量全自动跟踪的方法 - Google Patents

一种数字x射线透视剂量全自动跟踪的方法 Download PDF

Info

Publication number
CN107845070B
CN107845070B CN201710347095.8A CN201710347095A CN107845070B CN 107845070 B CN107845070 B CN 107845070B CN 201710347095 A CN201710347095 A CN 201710347095A CN 107845070 B CN107845070 B CN 107845070B
Authority
CN
China
Prior art keywords
perspective
mean
value
frame image
calculating
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
CN201710347095.8A
Other languages
English (en)
Other versions
CN107845070A (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.)
Shenzhen Lanying Medical Technology Co.,Ltd.
Original Assignee
Shenzhen Lanyun Medical Image 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 Shenzhen Lanyun Medical Image Co ltd filed Critical Shenzhen Lanyun Medical Image Co ltd
Priority to CN201710347095.8A priority Critical patent/CN107845070B/zh
Publication of CN107845070A publication Critical patent/CN107845070A/zh
Application granted granted Critical
Publication of CN107845070B publication Critical patent/CN107845070B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • 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/008Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction

Landscapes

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

Abstract

本发明公开了一种数字X射线透视剂量全自动跟踪的方法,包括以下步骤:获取首帧图像,并采集与计算首帧图像的灰度值和透视参数;采用首帧图像的灰度值和透视参数计算当前第二帧的透视参数,并采集与计算第二帧的灰度均值;利用当前第二帧图像的灰度值和透视参数计算下一帧的透视参数,并采集与计算下一帧的灰度值,依次循环,直到透视结束。通过上述方式,本发明可降低设备的成本,并且可以避免控制电路的伪影出现在图像上,提高透视图像的质量;此外本发明的计算速度快,结果精度高,可减少患者的透视时间,降低患者的辐射剂量。

Description

一种数字X射线透视剂量全自动跟踪的方法
技术领域
本发明涉及医学图像处理技术领域,特别是涉及一种数字X射线透视剂量全自动跟踪的方法。
背景技术
X射线是一种电磁波,具有很强的穿透能力,甚至能够穿透一定厚度的钢板,因此常被用于对物体内部的透视成像。自1895年被德国物理学家伦琴发现后,X射线首先被运用到医学成像与诊断上,在随后的一百多年当中,在医学、安全检查、无损检测、工业探伤等领域中均发挥了巨大作用。X射线透视成像技术是违禁品检查领域最基本的、也是最早得到广泛应用的方法。
在医学X射线透视中,由于拍摄人体部位可能会动态的发生变化,为了使透视图像有个稳定、适当的亮度(亮度值用BT表示)显示,一般会根据当前拍摄到的影像,对拍摄剂量自动做相应的调整,即图像亮度变暗,增加剂量,反之降低剂量,但是,对于数字X线透视检查,患者需要长时间暴露在X射线中,因此透视剂量是非常重要的,应在获得高质量图像的前提下尽可能的降低透视剂量,因此需要一种数字X射线透视剂量全自动跟踪方法来控制透视剂量,缩短透视时间,降低患者的辐射剂量。
现有技术中通常采用动态调剂量的方法,其一般是逐渐提升剂量,每提升一次剂量,检查一次亮度值BT是否达到预定目标亮度值(目标亮度值用DBT表示)范围,如果没有达到,则逐渐增加剂量,直到满足条件为止。这个逐渐调整剂量的周期长,在调整好剂量前往往部分图像是无效的,不能为临床诊断提供有用的信息,并且增加了透视时间,使得被透视的人吸收了不必要的剂量。
目前调节剂量的方式有两种,一种是调节KV值,另一种是调节MAS值,其中,KV表示X线球管的管电压;mAs:毫安秒,表示X线球管的输出功率。在透视过程中,因为KV可以调节图像亮度BT的同时,还可以调节射线的穿透力,所以一般固定MAS调KV值比较适合。
为了实现动态调剂量及控制透视剂量,公开号为CN105405121A的中国发明专利申请说明书中公开了一种调整X射线透视图像亮度的方法,包括:根据人体不同部位的厚度,采集X射线图像的管电压值及对应的亮度值,并绘制对应的KV-BT拟合曲线;修正所述KV-BT拟合曲线,获得KV-BT修正曲线;根据当前图像的管电压值及平均亮度值,获取对应厚度的KV-BT修正曲线;根据所述对应厚度的KV-BT修正曲线,获取与目标亮度值对应的目标管电压值。然而,此方法仅能获取管电压值,而无法获取透视的另外一个参数----管电流值;此外,需要获取受照体的厚度,增加机械的设计难度。
公开号为CN 2569512Y的中国实用新型专利说明书中揭露了一种X光机亮度自动控制装置,该装置有亮度采集电路和亮度自动控制电路两部分组成,由采样板亮度信号LIGHT,通过亮度自动控制处理板模块的隔离、放大、叠加,得到PLIGHT,同时,由限束器输出到亮度自动控制装置的横向、纵向限束器位置反馈信号分别通过亮度自动控制处理板模块的隔离、调零、放大,得到PHF和PVF,PLIGHT、PHF、PVF3个信号输入到亮度自动控制处理板的模块IV,通过模块IV的运算得到ABS,再经过模块V的隔离得到亮度自动控制装置最终的输出ABS-OUT。此方法需要光电二极管和控制电路,会增加成本,同时也存在图像上出现电路伪影的情况,影响图像质量;另外,需要与限束器进行通信,也就是说需要特制的限束器。
发明内容
本发明的目的是克服现有技术中的不足之处,提供一种成本低、速度快、数据准确性高的一种数字X射线透视剂量全自动跟踪方法,不仅可以获得高质量的透视图像,而且可大幅度减少患者的透视时间,降低患者的辐射剂量。
本发明的目的是通过以下技术方案来实现的:
一种数字X射线透视剂量全自动跟踪的方法,包括以下步骤:
S1、获取首帧图像,并采集与计算首帧图像的灰度值和透视参数;
S2、采用首帧图像的灰度值和透视参数计算当前第二帧的透视参数,并采集与计算第二帧的灰度均值;
S3、利用当前第二帧图像的灰度值和透视参数计算下一帧的透视参数,并采集与计算下一帧的灰度值,依次循环,直到透视结束。具体的,所述步骤S1的具体步骤为:
固定X线球管的输出功率mA0值,将X线球管的管电压KV值从低到高依次曝光,计算每个KV值对应图像的灰度值,建立KV值与灰度值的拟合曲线L;
固定首帧图像的透视参数进行透视,获取首帧图像I1,并计算首帧图像I1的灰度值Mean1
优选的,所述mA0=4.0,所述首帧图像的透视参数中,首帧图像的透视参数的X线球管的管电压KV1=60,首帧图像的透视参数的X线球管的输出功率mA1=mA0=4.0。
具体的,所述步骤S2的具体步骤为:
S201、若首帧图像I1的灰度值Mean1<ThrMin,则认为是致密组织,按照曲线L调高X线球管的管电压KV值,调整后KV值后的获取的图像其等效灰度值记为Mean1D;若Mean1>ThrMax,则认为是疏松组织,按照曲线L适当调小X线球管的管电压KV值,调整后KV值后的获取的图像等效灰度值记为Mean1D;若ThrMin<Mean1<ThrMax,则认为是正常组织,不用调整KV,记为此时的等效灰度值记为Mean1D=Mean1;其中ThrMin和ThrMax分别为正常组织密度的上下阈值;
S202、采用公式mA2=mAO*(MeanT/Mean1D)计算第二帧图像的透视的X线球管的输出功率mA2,其中MeanT为透视图像的目标灰度值,Mean1D为当前帧图像的灰度值;
S203、采用调整后的KV值和mA2值继续透视,直至当前帧图像的灰度值Mean1D满足ThrMin<Mean1D<ThrMax,并获取此时满足条件后的第二帧图像I2
S204、计算图像I2的灰度均值Mean2
具体的,所述利用当前第二帧图像的灰度值和透视参数计算下一帧的透视参数具体为:
在与当前第二帧图像相同KV值的条件下,根据KV值与灰度值的拟合曲线L的线性关系,计算mA0时对应的图像灰度值Mean2D
重复步骤S201-S203的透视循环方式,计算下一帧图像的透视参数。
本发明相比现有技术具有以下优点及有益效果:区别于现有技术,本发明通过首帧图像的灰度值和透视参数计算第二帧的透视参数,然后利用第二帧图像的灰度值和透视参数计算后一帧的透视参数,依次循环,直到透视结束,在整个全自动跟踪过程中不需要二极管和控制电路等额外电子元器件,可降低设备的成本,并且可以避免控制电路的伪影出现在图像上,提高透视图像的质量;此外本发明的计算速度快,结果精度高,可减少患者的透视时间,降低患者的辐射剂量。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明方法实施例的步骤示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然所描述的实施例仅是本发明一个实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明最关键的构思在于:根据首帧图像的灰度值和透视参数计算第二帧的透视参数,然后利用第二帧图像的灰度值和透视参数计算后一帧的透视参数,依次循环,直到透视结束。
如图1,本发明实施例提供一种数字X射线透视剂量全自动跟踪的方法,包括如下步骤:
Y1:采用固定X线球管的输出功率mA0=4.0,KV从低到高依次曝光,依次计算每个KV对应图像的灰度值,建立KV与灰度值的拟合曲线L;
Y2:采用固定首帧图像的透视参数KV1=60,mA1=4.0进行透视,从而获取首帧图像I1
Y3:根据图像I1计算灰度均值Mean1
Y4:假如Mean1<ThrMin,则认为是致密组织,按照曲线L适当调高KV,调整后KV记为KV2,等效灰度值记为Mean1D;假如Mean1>ThrMax,则认为是疏松组织,按照曲线L适当调小KV,调整后KV记为KV2,等效灰度值记为Mean1D;假如ThrMin<Mean1<ThrMax,则认为是正常组织,不用调整KV,KV2=KV,Mean1D=Mean1;其中ThrMin和ThrMax分别为正常组织密度的上下阈值,由临床经验值确定;
Y5:采用公式mA2=mAO*(MeanT/Mean1D)计算第二帧图像的透视的mA2,其中MeanT为透视图像的目标灰度值,由临床经验值确定;
Y6:采用KV2和mA2继续透视,获取第二帧图像I2
Y7:根据图像I2计算灰度均值Mean2
Y8:根据相同KV下,KV值与灰度值的拟合曲线L的线性关系,计算mA0时对应的图像灰度值Mean2D
Y9:重复步骤:Y4),计算第三帧图像的透视参数KV3
Y10:重复步骤Y5),计算第三帧图像的透视参数mA3
Y11:采用KV3和mA3继续透视,获取第三帧图像I3
Y12:循环步骤Y6)——步骤Y11),直到透视结束。
本发明利用前一帧图像的灰度值和透视参数计算当前帧的透视参数,然后利用当前帧图像的灰度值和透视参数计算后一帧的透视参数,依次循环,直到透视结束。
相对现有技术,由于本发明不需要光电二极管和控制电路等额外电子元器件,可适当降低设备的成本,并且可以避免控制电路的伪影出现在图像上,提高透视图像的质量;此外本发明的计算速度快,结果精度高,可减少患者的透视时间,降低患者的辐射剂量。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。

Claims (2)

1.一种数字X射线透视剂量全自动跟踪的方法,其特征在于,包括以下步骤:
S1、获取首帧图像,并采集与计算首帧图像的灰度值和透视参数;
S2、采用首帧图像的灰度值和透视参数计算当前第二帧的透视参数,并采集与计算第二帧的灰度均值;
S3、利用当前第二帧图像的灰度值和透视参数计算下一帧的透视参数,并采集与计算下一帧的灰度值,依次循环,直到透视结束;
所述步骤S1的具体步骤为:
固定X线球管的输出功率mA0值,将X线球管的管电压KV值从低到高依次曝光,计算每个KV值对应图像的灰度值,建立KV值与灰度值的拟合曲线L;
固定首帧图像的透视参数进行透视,获取首帧图像I1,并计算首帧图像I1的灰度值Mean1
所述mA0=4.0,所述首帧图像的透视参数中,首帧图像的透视参数的X线球管的管电压KV1=60,首帧图像的透视参数的X线球管的输出功率mA1=mA0=4.0;
所述步骤S2的具体步骤为:
S201、若首帧图像I1的灰度值Mean1<ThrMin,则认为是致密组织,按照曲线L调高X线球管的管电压KV值,调整后KV值后的获取的图像其等效灰度值记为Mean1D;若Mean1>ThrMax,则认为是疏松组织,按照曲线L适当调小X线球管的管电压KV值,调整后KV值后的获取的图像等效灰度值记为Mean1D;若ThrMin<Mean1<ThrMax,则认为是正常组织,不用调整KV,记为此时的等效灰度值记为Mean1D=Mean1;其中ThrMin和ThrMax分别为正常组织密度的上下阈值;
S202、采用公式mA2=mAO*(MeanT/Mean1D)计算第二帧图像的透视的X线球管的输出功率mA2,其中MeanT为透视图像的目标灰度值,Mean1D为当前帧图像的灰度值;
S203、采用调整后的KV值和mA2值继续透视,直至当前帧图像的灰度值Mean1D满足ThrMin<Mean1D<ThrMax,并获取此时满足条件后的第二帧图像I2
S204、计算图像I2的灰度均值Mean2
2.根据权利要求1所述一种数字X射线透视剂量全自动跟踪的方法,其特征在于,所述利用当前第二帧图像的灰度值和透视参数计算下一帧的透视参数具体为:
在与当前第二帧图像相同KV值的条件下,根据KV值与灰度值的拟合曲线L的线性关系,计算mA0时对应的图像灰度值Mean2D
重复步骤S201-S203的透视循环方式,计算下一帧图像的透视参数。
CN201710347095.8A 2017-05-17 2017-05-17 一种数字x射线透视剂量全自动跟踪的方法 Active CN107845070B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710347095.8A CN107845070B (zh) 2017-05-17 2017-05-17 一种数字x射线透视剂量全自动跟踪的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710347095.8A CN107845070B (zh) 2017-05-17 2017-05-17 一种数字x射线透视剂量全自动跟踪的方法

Publications (2)

Publication Number Publication Date
CN107845070A CN107845070A (zh) 2018-03-27
CN107845070B true CN107845070B (zh) 2021-08-03

Family

ID=61682858

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710347095.8A Active CN107845070B (zh) 2017-05-17 2017-05-17 一种数字x射线透视剂量全自动跟踪的方法

Country Status (1)

Country Link
CN (1) CN107845070B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108924432B (zh) * 2018-06-12 2021-01-05 东软医疗系统股份有限公司 一种信号偏差的校正方法、装置及设备
CN111781633B (zh) * 2019-04-03 2023-12-22 苏州博思得电气有限公司 利用校正自动调节x射线机亮度方法、系统、存储介质、设备
CN113096081B (zh) * 2021-03-30 2024-06-07 海辉医学(北京)科技有限公司 一种x射线曝光亮度控制方法、装置、设备及存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105405121A (zh) * 2015-10-21 2016-03-16 深圳市安健科技有限公司 调整x射线透视图像亮度的方法及系统
CN105528764A (zh) * 2015-12-01 2016-04-27 沈阳东软医疗系统有限公司 一种图像亮度调节方法、装置及设备

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101732056B (zh) * 2009-12-01 2011-07-20 西安集智医疗器械科技有限公司 一种医用x射线影像增强处理用灰度等级线性处理方法
CN104367331B (zh) * 2013-08-15 2017-02-15 深圳市蓝韵实业有限公司 用于数字乳腺x射线机的全数字自动曝光方法
CN104156917A (zh) * 2014-07-30 2014-11-19 天津大学 基于双能谱的x射线ct图像增强方法
CN106308827A (zh) * 2015-07-01 2017-01-11 南京普爱医疗设备股份有限公司 一种三维c形臂自动剂量控制方法
CN105636327B (zh) * 2015-12-28 2017-08-25 上海奕瑞光电子科技有限公司 一种基于透射可见光的曝光控制系统及方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105405121A (zh) * 2015-10-21 2016-03-16 深圳市安健科技有限公司 调整x射线透视图像亮度的方法及系统
CN105528764A (zh) * 2015-12-01 2016-04-27 沈阳东软医疗系统有限公司 一种图像亮度调节方法、装置及设备

Also Published As

Publication number Publication date
CN107845070A (zh) 2018-03-27

Similar Documents

Publication Publication Date Title
US6459765B1 (en) Automatic exposure control and optimization in digital x-ray radiography
CN107845070B (zh) 一种数字x射线透视剂量全自动跟踪的方法
US20130202086A1 (en) Radiation imaging apparatus and control method thereof, and radiation imaging system
JP5042887B2 (ja) 放射線画像撮影装置
US20150139382A1 (en) X-ray imaging apparatus and method of controlling the same
CN104105444B (zh) 放射线摄影装置
US20160302752A1 (en) Radiation image capturing system and body thickness estimation method
KR20150107335A (ko) 디지털 엑스레이 영상 시스템, 엑스레이 조사 조절 장치 및 그 방법
JP6083990B2 (ja) 放射線撮影装置、その制御方法及びプログラム
US9031190B2 (en) X-ray imaging apparatus
US7457452B2 (en) Radiographic image diagnosis device
CN101721222B (zh) 一种修正床板和摆位辅助装置对图像质量影响的方法
CN105310704B (zh) 放射成像的控制系统及控制方法
CN111781633B (zh) 利用校正自动调节x射线机亮度方法、系统、存储介质、设备
JP2018149166A (ja) 放射線画像処理装置
US7558412B2 (en) System and method for compensation of scintillator hysteresis in x-ray detectors
JP2009219529A (ja) 放射線画像撮影装置
EP3708082A1 (en) Image processing device, image processing method, radiation imaging device, and method and program for controlling radiation imaging device
US11690587B2 (en) Apparatus comprising data obtaining unit and image processing unit and method for processing X-ray image
US11950948B2 (en) Systems and methods for mitigating imaging artifacts
EP0974064B1 (de) Verfahren zur bilderzeugung bei der digitalen dentalen radiographie
KR102556513B1 (ko) 엑스선 영상 처리 장치 및 방법
CN108924432B (zh) 一种信号偏差的校正方法、装置及设备
WO2020012520A1 (ja) 医用x線画像処理装置およびx線画像撮影装置
US8300913B2 (en) Apparatus and method for processing radiation image

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP03 Change of name, title or address

Address after: No.103, baguang District Service Center, No.2 BaiShaWan Road, baguang community, Kuiyong street, Dapeng New District, Shenzhen, Guangdong 518000

Patentee after: Shenzhen Lanying Medical Technology Co.,Ltd.

Address before: 518108 floor 1, building 2, Honglong high tech Industrial Park (Lanyun Industrial Park), Liyuan Industrial Zone, Zhuanchang village, Langxin neighborhood committee, Shiyan street, Bao'an District, Shenzhen, Guangdong Province

Patentee before: SHENZHEN LANYUN MEDICAL IMAGE CO.,LTD.

CP03 Change of name, title or address