CN105193436A - 数字乳腺x射线成像系统曝光剂量的评价方法 - Google Patents

数字乳腺x射线成像系统曝光剂量的评价方法 Download PDF

Info

Publication number
CN105193436A
CN105193436A CN201410280940.0A CN201410280940A CN105193436A CN 105193436 A CN105193436 A CN 105193436A CN 201410280940 A CN201410280940 A CN 201410280940A CN 105193436 A CN105193436 A CN 105193436A
Authority
CN
China
Prior art keywords
exposure
index
imaging system
ray imaging
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.)
Granted
Application number
CN201410280940.0A
Other languages
English (en)
Other versions
CN105193436B (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 Landwind Industry 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 Landwind Industry Co Ltd filed Critical Shenzhen Landwind Industry Co Ltd
Priority to CN201410280940.0A priority Critical patent/CN105193436B/zh
Publication of CN105193436A publication Critical patent/CN105193436A/zh
Application granted granted Critical
Publication of CN105193436B publication Critical patent/CN105193436B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Measurement Of Radiation (AREA)

Abstract

本发明公开一种数字乳腺X射线成像系统曝光剂量的评价方法,包括如下步骤:制定目标曝光指数;计算偏差指数DI。使用本发明的技术方案,根据曝光指数与目标曝光指数的偏差来评价乳腺图像曝光剂量高低的方法,从而方便医生后续调整曝光参数,获得曝光剂量适当的图像,算法简单,使用方便。

Description

数字乳腺X射线成像系统曝光剂量的评价方法
技术领域
本发明涉及一种医学影像处理技术,特别涉及一种医学影像曝光剂量的评价方法。
背景技术
乳腺作为女性健康的重要课题,已经越来越被人们所重视。近年来,乳腺癌的发病率在我国逐年上升,在许多大城市中乳腺癌已经成为女性最常见的恶性肿瘤,严重威胁着女性的身体健康,研究表明乳腺癌如果能早期发现,那么将大大提高其治愈率,所以乳腺癌的早期诊断具有非常重要的临床价值。
实际临床发现乳腺癌的早期表现为乳腺内存在微小的钙化点,这些微小的钙化点在使用超声成像时很难发现,而使用数字乳腺X射线成像系统对乳腺进行检查是早期发现钙化点的有效方式。对于X射线检查,辐射剂量是非常重要的,应在获得高质量图像的前提下尽可能的降低曝光参数。因此需要一种评价乳腺图像曝光剂量的方法和相关标准。目前能够搜索到的近似相关专利有:
《一种乳腺机X射线剂量控制方法及系统》,申请人:北京国药恒瑞美联信息技术有限公司。此方法通过针对不同性质的目标乳腺进行预曝光获得预曝光图像,并对预曝光图像进行分析计算从而获得目标乳腺在不同投照角度的X射线剂量值,保证获取最佳的拍摄图像。此方法只能控制当次曝光的曝光参数,而无法评价曝光剂量的高低。
《乳腺X光机及其实现全自动曝光的方法》,申请人:深圳圣诺医疗有限公司。此方法首先根据压迫厚度确定预曝光参数,然后利用预曝光图像的剂量信息确定第二次曝光的剂量参数。此方法可以自动的计算合适的曝光参数,但是无法评价曝光剂量的高低。
发明内容
为了解决上述的技术问题,本发明提出一种数字乳腺X射线成像系统曝光剂量的评价方法,本发明采用如下技术方案来解决上述技术问题:
一种数字乳腺X射线成像系统曝光剂量的评价方法,包括如下步骤:
a.制定目标曝光指数;
b.计算偏差指数DI。
优选的,在上述的一种数字乳腺X射线成像系统曝光剂量的评价方法,中所述步骤a包括如下步骤:
a1.收集临床病例中对应体位的100幅原始图像;
a2.计算每幅图像曝光指数EI;
a3.将所有图像的EI值按照从小到大排序得到EI值的序列LEI
a4.从序列LEI中取中间的10个EI值,求平均值,这个平均值即为目标曝光指数EIT
a5.将EIT写入配置文件或数据库,供计算偏差指数DI时调用。
优选的,在上述一种数字乳腺X射线成像系统曝光剂量的评价方法中,所述步骤a2包括如下步骤:
a21.校正曝光指数,获取图像斜率Slop和偏移量Offset;
a22.读入原始图像;
a23.分割出乳腺组织区域ROI,计算ROI的像素均值Mean和像素点个数Count;
a24.对乳腺组织区域进行二次分割,分割出组织区域中的致密区域secondROI,计算secondROI的像素均值secondMean和像素点个数secondCount;
a25.利用公式计算乳腺组织区域的感兴趣值Value,所述公式为:
Value=(secondMean*secondCount+Mean*(Count-secondCount))/Count;
a26.利用公式EI=10*(Slop*Value+Offset)计算曝光指数EI,并将EI写入Dicom文件。
优选的,在上述的一种数字乳腺X射线成像系统曝光剂量的评价方法中,步骤a21具体包括如下步骤:
A.将剂量仪放置于平板探测器的胸壁侧,确定RQM8辐射质量;
B.分别采集6不同mAs的图像,并记录每幅图像对应的剂量Dose;
C.计算每幅图像的均值Mean;
D.计算Mean与Dose间的线性关系,通过线性拟合公式Dose=Slop*Mean+Offset确定斜率Slop和偏移量Offset,将Slop和Offset写入配置文件或数据库,供计算曝光指数EI时调用。
优选的,在上述的一种数字乳腺X射线成像系统曝光剂量的评价方法中,
所述步骤b为:使用公式计算偏差指数DI,所述公式为将生成的DI写入Dicom文件。
不同的乳腺中致密腺体区域所占整个腺体区域的比例不同,本发明计算乳腺组织区域的感兴趣值时考虑了致密腺体区域占整个腺体区域的比例,提高了相关数据的准确性。使用本发明的技术方案,根据曝光指数与目标曝光指数的偏差来评价乳腺图像曝光剂量高低的方法,从而方便医生后续调整曝光参数,获得曝光剂量适当的图像,算法简单,使用方便。
附图说明
图1为本发明流程图;
图2为本发明曝光指数校正流程图;
图3为计算曝光指数流程图;
图4为制定目标曝光指数流程图。
具体实施方式
下面结合附图对本发明作进一步详细说明,但是本发明不仅仅局限于以下实施例。本发明所述的一种数字乳腺X射线成像系统曝光剂量的评价方法,具体包括以下步骤:
1.曝光指数校正:
因不同平板探测器对剂量的响应不同,故使用不同类型的平板探测器前均需要首先进行曝光指数的校正。
将剂量仪放置于平板探测器的胸壁侧即患者拍摄时靠近身体的一侧,确定RQM8辐射质量。RQM8辐射质量是本领域研究人员所熟知的相关标准,在此不做描述。
分别采集6不同mAs的图像,并记录每幅图像对应的剂量Dose。
计算每幅图像的均值Mean。采用OTSU分割算法分割掉剂量仪区域;为了克服边界的影响,可适当裁剪掉图像的部分边界,本例中胸壁侧不裁剪,其余三个边界均裁剪掉20个像素的宽度。OTSU分割算法是本领域研究人员所熟知的相关技术,故在此不做描述。
计算Mean与Dose间的线性关系。通过线性拟合公式Dose=Slop*Mean+Offset确定斜率Slop和偏移量Offset。可将Slop和Offset写入配置文件或数据库,供计算曝光指数EI时调用。
2.计算曝光指数EI:
在得到Slop和Offset的前提下计算曝光指数EI。该数据表征拍摄图像时对患者的曝光程度。
读入原始图像O。图像允许相关的校正处理,但是不允许进行增强处理。
利用OTSU分割算法分割出乳腺组织区域ROI。计算ROI的像素均值Mean和像素点个数Count。OTSU分割算法是本领域研究人员所熟知的相关技术,故在此不做描述。
利用OTSU分割算法分割对乳腺组织区域进行二次分割,分割出组织区域中的致密区域secondROI,计算secondROI的像素均值secondMean和像素点个数secondCount。
利用公式Value=(secondMean*secondCount+Mean*(Count-secondCount))/Count计算乳腺组织区域的感兴趣值Value。
利用公式EI=10*(Slop*Value+Offset)计算曝光指数EI,并将EI写入Dicom文件。
3.目标曝光指数的制定:
不同体位的目标曝光指数不同,本实施例以RCC位即右侧正位为例,介绍相关的流程。
收集100副RCC位的原始图像。
根据步骤2计算每幅图像的曝光指数EI,利用快速排序法对所有图像的EI值由小到大排序,得到EI值的序列LEI。快速排序法是本领域研究人员所熟知的相关技术,故在此不做描述。
从序列LEI中取中间的10个EI值,这10个EI值的平均值就作为RCC位的目标曝光指数EIT
将EIT写入配置文件或数据库,供计算偏差指数DI时调用。
4.计算偏差指数DI:
DI表示曝光剂量的高低,DI为正时表示曝光剂量偏高,DI为负时表示曝光剂量偏低。由下式计算:
DI = 10 * log 10 ( EI EI T )
将DI写入Dicom文件,方便以后查看图像曝光剂量的高低。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (5)

1.一种数字乳腺X射线成像系统曝光剂量的评价方法,其特征在于,包括如下步骤:
a.制定目标曝光指数;
b.计算偏差指数DI。
2.根据权利要求1所述的一种数字乳腺X射线成像系统曝光剂量的评价方法,其特征在于,所述步骤a包括如下步骤:
a1.收集临床病例对应体位的100幅原始图像;计算每幅图像曝光指数EI;
a2.将所有图像的EI值按照从小到大排序得到EI值的序列LEI
a3.从序列LEI中取中间的10个EI值,求平均值,这个平均值即为目标曝光指数EIT
a4.将EIT写入配置文件或数据库,供计算偏差指数DI时调用。
3.根据权利要求2所述的一种数字乳腺X射线成像系统曝光剂量的评价方法,其特征在于,所述步骤a2包括如下步骤:
a21.校正曝光指数,获取图像斜率Slop和偏移量Offset;
a22.读入原始图像;
a23.分割出乳腺组织区域ROI,计算ROI的像素均值Mean和像素点个数Count;
a24.对乳腺组织区域进行二次分割,分割出组织区域中的致密区域secondROI,计算secondROI的像素均值secondMean和像素点个数secondCount;
a25.利用公式计算乳腺组织区域的感兴趣值Value,所述公式为:
Value=(secondMean*secondCount+Mean*(Count-secondCount))/Count;
a26.利用公式EI=10*(Slop*Value+Offset)计算曝光指数EI,并将EI写入Dicom文件。
4.根据权利要求3所述的一种数字乳腺X射线成像系统曝光剂量的评价方法,其特征在于,步骤a21具体包括如下步骤:
A.将剂量仪放置于平板探测器的胸壁侧,确定RQM8辐射质量;
B.分别采集6不同mAs的图像,并记录每幅图像对应的剂量Dose;
C.计算每幅图像的均值Mean;
D.计算Mean与Dose间的线性关系,通过线性拟合公式Dose=Slop*Mean+Offset确定斜率Slop和偏移量Offset,将Slop和Offset写入配置文件或数据库,供计算曝光指数EI时调用。
5.根据权利要求1所述的一种数字乳腺X射线成像系统曝光剂量的评价方法,其特征在于,所述步骤b为:使用公式计算偏差指数DI,所述公式为
将生成的DI与入Dicom文件。
CN201410280940.0A 2014-06-20 2014-06-20 数字乳腺x射线成像系统曝光剂量的评价方法 Active CN105193436B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410280940.0A CN105193436B (zh) 2014-06-20 2014-06-20 数字乳腺x射线成像系统曝光剂量的评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410280940.0A CN105193436B (zh) 2014-06-20 2014-06-20 数字乳腺x射线成像系统曝光剂量的评价方法

Publications (2)

Publication Number Publication Date
CN105193436A true CN105193436A (zh) 2015-12-30
CN105193436B CN105193436B (zh) 2018-04-20

Family

ID=54941674

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410280940.0A Active CN105193436B (zh) 2014-06-20 2014-06-20 数字乳腺x射线成像系统曝光剂量的评价方法

Country Status (1)

Country Link
CN (1) CN105193436B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107049345A (zh) * 2017-05-23 2017-08-18 深圳艾砾拓科技有限公司 数字x射线乳腺摄影系统中曝光指数及偏差指数计算方法
CN108742663A (zh) * 2018-04-03 2018-11-06 深圳蓝韵医学影像有限公司 曝光剂量评价方法、装置及计算机可读存储介质
CN108846847A (zh) * 2018-05-08 2018-11-20 嘉兴全景医疗器械有限公司 基于平板状结构形状滤波器的乳腺图像分割方法
CN109239755A (zh) * 2018-09-26 2019-01-18 深圳蓝韵医学影像有限公司 计算乳腺x线系统平均腺体剂量方法、系统、设备及介质
CN113114952A (zh) * 2021-04-07 2021-07-13 华南理工大学 一种评价乳腺随访图像质量的方法、系统及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102196772A (zh) * 2008-09-04 2011-09-21 霍罗吉克公司 集成多模式乳房x射线照相术/体层合成x射线系统和方法
CN103826538A (zh) * 2011-09-27 2014-05-28 富士胶片株式会社 放射线成像系统及其操作方法,以及放射线图像检测装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102196772A (zh) * 2008-09-04 2011-09-21 霍罗吉克公司 集成多模式乳房x射线照相术/体层合成x射线系统和方法
CN103826538A (zh) * 2011-09-27 2014-05-28 富士胶片株式会社 放射线成像系统及其操作方法,以及放射线图像检测装置

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107049345A (zh) * 2017-05-23 2017-08-18 深圳艾砾拓科技有限公司 数字x射线乳腺摄影系统中曝光指数及偏差指数计算方法
CN108742663A (zh) * 2018-04-03 2018-11-06 深圳蓝韵医学影像有限公司 曝光剂量评价方法、装置及计算机可读存储介质
CN108846847A (zh) * 2018-05-08 2018-11-20 嘉兴全景医疗器械有限公司 基于平板状结构形状滤波器的乳腺图像分割方法
CN108846847B (zh) * 2018-05-08 2021-12-03 王磊 基于平板状结构形状滤波器的乳腺图像分割方法
CN109239755A (zh) * 2018-09-26 2019-01-18 深圳蓝韵医学影像有限公司 计算乳腺x线系统平均腺体剂量方法、系统、设备及介质
CN113114952A (zh) * 2021-04-07 2021-07-13 华南理工大学 一种评价乳腺随访图像质量的方法、系统及装置
CN113114952B (zh) * 2021-04-07 2022-02-18 华南理工大学 一种评价乳腺随访图像质量的方法、系统及装置

Also Published As

Publication number Publication date
CN105193436B (zh) 2018-04-20

Similar Documents

Publication Publication Date Title
CN104460181A (zh) 数字放射成像曝光剂量的评价方法
CN105193436A (zh) 数字乳腺x射线成像系统曝光剂量的评价方法
CN103347447B (zh) 乳房密度评估
CN103886576B (zh) 一种腺体组织特征灰度检测方法及装置
CN107049345A (zh) 数字x射线乳腺摄影系统中曝光指数及偏差指数计算方法
US20160140721A1 (en) Radiographic image analysis device and method, and recording medium having program recorded therein
US8655039B2 (en) Method of segmenting lesions in images
CN103460254B (zh) 基于图像内容的自动亮度检测方法与装置
US9734574B2 (en) Image processor, treatment system, and image processing method
JP2012035068A (ja) 放射線画像処理装置および方法並びにプログラム
US10430930B2 (en) Image processing apparatus, image processing method, and image processing program for performing dynamic range compression process
CN113034389A (zh) 图像处理方法、装置、计算机设备和存储介质
CN107007294A (zh) X射线图像诊断装置以及骨密度测量方法
US9763636B2 (en) Method and system for spine position detection
JP2016190012A (ja) ノイズ抑制処理装置および方法並びにプログラム
WO2015040547A1 (en) Method and system for spine position detection
US11000239B2 (en) Methods, systems, and apparatus for determining radiation doses
US10089728B2 (en) Radiation-image processing device and method
CN106651768A (zh) 图像校正方法及装置、x射线摄影设备
CN106204548A (zh) 图像分辨方法及装置
CN108742663A (zh) 曝光剂量评价方法、装置及计算机可读存储介质
US20200015767A1 (en) Sensitivity optimized patient positioning system for dark-field x-ray imaging
CN105286960B (zh) 乳腺病灶穿刺位置的定位方法及系统
CN109813259B (zh) 高动态x射线成像方法、存储介质和装置
KR101750173B1 (ko) 방사선 유효 피폭량 자동 계산 시스템 및 방법

Legal Events

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

Effective date of registration: 20180731

Address after: 518000 Shenzhen, Baoan District, Guangdong Shiyan street, Lang Xin neighborhood committee, brick factory village, Liyuan Industrial Zone, hung lung hi tech Industrial Park (blue Yun Industrial Park), 1 1 buildings.

Patentee after: Shenzhen blue ribbon medical imaging Co., Ltd.

Address before: 518000 room 408-413, E garden, Jingjing Road, Jingtian Road, Futian District, Shenzhen, Guangdong.

Patentee before: Landwind Co., Ltd.

CP03 Change of name, title or address
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: 518000 Shenzhen, Baoan District, Guangdong Shiyan street, Lang Xin neighborhood committee, brick factory village, Liyuan Industrial Zone, hung lung hi tech Industrial Park (blue Yun Industrial Park), 1 1 buildings.

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