DR图像亮度自适应调整方法及系统技术方案

本发明专利技术实施例公开了一种DR图像亮度自适应调整方法及系统，所述方法包括：提取组织区域步骤：对待调整图像I进行分割，提取组织区域的二值图像I

【技术实现步骤摘要】
DR图像亮度自适应调整方法及系统
本专利技术涉及图像诊断
，尤其涉及一种DR图像亮度自适应调整方法及系统。
技术介绍
DR(DigitalRadiography，数字X线摄影)系统以X线为信息载体，利用其物理特性，即穿透人体时因组织密度差而发生不同强度的衰减，来获取人体组织状态信息，以供医学诊断。作为医疗影像诊断设备，DR系统输出的图像质量在客观上决定了医学诊断的精度。因此，任何可能影响图像质量的因素，都必须认真考虑。对于医学图像，其质量通常表现在亮度、对比度、锐利度、信噪比等方面。其中，图像亮度是图像信息输出的重要反应，是人眼最直接也是最敏感的方面之一，其适宜与否不仅影响用户观感体验，更对诊断结果起着关键性影响。临床摄片中，曝光射线量的多少以及不同患者体型对X线的衰减强弱，直接决定了DR系统探测器接收到的原始信号强度，从而间接决定了所获图像的明暗程度。因此，要获取一张亮度适宜的DR图像，往往要求拍片医生能够根据患者体型，精准把握曝光剂量。然而在实际中，这种控制很难把握。要么剂量过高，图像偏暗，要么剂量过低，图像偏亮，以致常常需要重复曝光，不仅增加医生工作量，也无形中增加了患者对射线的吸收，造成患者二次伤害。对于DR图像亮度控制，目前常用的有两类方案。一类是借助AEC(AutomaticExposureControl，自动曝光控制)硬件模块，自动控制曝光量，来实现图像亮度稳定。另一类是应用图像处理技术，对图像亮度进行调整，如常用的Gamma亮度校正。从实际应用来看，前种方案基本无需额外人工干预，具有操作简便和图像效果好的优点，然而由于需要额外增加硬件部件，导致经济性方面不够友好。后种方案为纯图像处理技术，不依赖于额外硬件，经济性好，然而就目前通常的实现方式来看，一般采用等力度的亮度调整方式，因此常常需要频繁地手动调整处理参数，增加了医生工作量。申请号为201910683490.2，专利技术名称为《一种医学图像亮度均一化校正方法》的专利技术专利为深度神经网络方案，通过构建GAN生成对抗网络，利用有代表性的已标注的亮度缺陷数据集进行训练，获取缺陷数据到正常数据的特征映射，最终实现对亮度不均图像的亮度校正。该方案对数据集的依赖性大，数据集的制作优劣，包括标注准确性和数据量大小等，直接影响模型的精度性能。另一方面，神经网络模型通常计算量巨大，对硬件资源要求较高，不利于临床应用。
技术实现思路
本专利技术实施例所要解决的技术问题在于，提供一种DR图像亮度自适应调整方法及系统，以实现图像亮度适宜稳定。为了解决上述技术问题，本专利技术实施例提出了一种DR图像亮度自适应调整方法，包括：提取组织区域步骤：利用最大类间方差法，对待调整图像I进行分割，提取组织区域的二值图像Ib；亮度计算步骤：根据组织区域的二值图像Ib计算组织区域归一化灰度均值b，将灰度均值b作为组织区域亮度；调整系数计算步骤：根据组织区域亮度计算亮度调整系数s；映射表计算步骤：根据亮度调整系数s计算亮度调整映射表；亮度调整步骤：根据亮度调整映射表，对待调整图像I进行亮度调整映射变换，调整图像亮度。进一步地，所述亮度计算步骤中，采用下式计算灰度均值b：其中，H和W分别为待调整图像的高度和宽度，a为与图像灰阶相对应的尺度因子。进一步地，a取值如下式：进一步地，调整系数计算步骤中，通过下式计算亮度调整系数s：s＝Cgain·b+Coffset；其中，Cgain和Coffset分别为亮度增益系数和亮度偏移系数。进一步地，所述映射表计算步骤中，通过下式计算亮度调整映射表：其中，k为输入灰度值，T(k)为输出灰度值，Imax和Imin分别为待调整图像I的最大值和最小值。相应地，本专利技术实施例还提供了一种DR图像亮度自适应调整系统，包括：提取组织区域模块：利用最大类间方差法，对待调整图像I进行分割，提取组织区域的二值图像Ib；亮度计算模块：根据组织区域的二值图像Ib计算组织区域归一化灰度均值b，将灰度均值b作为组织区域亮度；调整系数计算模块：根据组织区域亮度计算亮度调整系数s；映射表计算模块：根据亮度调整系数s计算亮度调整映射表；亮度调整模块：根据亮度调整映射表，对待调整图像I进行亮度调整映射变换，调整图像亮度。进一步地，所述亮度计算模块中，采用下式计算灰度均值b：其中，H和W分别为待调整图像的高度和宽度，a为与图像灰阶相对应的尺度因子。进一步地，a取值如下式：进一步地，调整系数计算模块中，通过下式计算亮度调整系数s：s＝Cgain·b+Coffset；其中，Cgain和Coffset分别为亮度增益系数和亮度偏移系数。进一步地，所述映射表计算模块中，通过下式计算亮度调整映射表：其中，k为输入灰度值，T(k)为输出灰度值，Imax和Imin分别为待调整图像I的最大值和最小值。本专利技术的有益效果为：本专利技术首先通过自动提取组织区域，来计算获取组织区域亮度，然后根据组织亮度，自适应计算相应亮度调整系数，最后执行亮度调整，实现图像亮度适宜稳定。作为纯图像处理技术方案，本专利技术无需额外增加硬件部件，经济性好。同时，整个亮度调整过程都是自适应进行，无需人工介入，节省了医生工作量。附图说明图1是本专利技术实施例的DR图像亮度自适应调整方法的流程示意图。图2是本专利技术实施例的DR图像亮度自适应调整系统的结构示意图。图3是本专利技术实施例的DR图像的调整效果对比图。具体实施方式需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互结合，下面结合附图和具体实施例对本专利技术作进一步详细说明。本专利技术实施例中若有方向性指示(诸如上、下、左、右、前、后……)仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。另外，在本专利技术中若涉及“第一”、“第二”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。请参照图1，本专利技术实施例的DR图像亮度自适应调整方法包括提取组织区域步骤、亮度计算步骤、调整系数计算步骤、映射表计算步骤及亮度调整步骤。提取组织区域步骤：利用最大类间方差法，对待调整图像I进行分割，提取组织区域的二值图像Ib：其中Ωt和Ωb分别表示组织区域和背景区域。亮度计算步骤：根据组织区域的二值图像Ib计算组织区域归一化灰度均值b，将灰度均值b作为组织区域亮度。调整系数计算步骤：根据组织区域亮度计算亮度调整系数s。映射表计算步骤：根据亮度调整系数s计算亮度调整映射表。亮度调整步骤：根据亮度调整映射表，对待本文档来自技高网...

【技术保护点】
1.一种DR图像亮度自适应调整方法，其特征在于，包括：/n提取组织区域步骤：利用最大类间方差法，对待调整图像I进行分割，提取组织区域的二值图像I

【技术特征摘要】
1.一种DR图像亮度自适应调整方法，其特征在于，包括：
提取组织区域步骤：利用最大类间方差法，对待调整图像I进行分割，提取组织区域的二值图像Ib；
亮度计算步骤：根据组织区域的二值图像Ib计算组织区域归一化灰度均值b，将灰度均值b作为组织区域亮度；
调整系数计算步骤：根据组织区域亮度计算亮度调整系数s；
映射表计算步骤：根据亮度调整系数s计算亮度调整映射表；
亮度调整步骤：根据亮度调整映射表，对待调整图像I进行亮度调整映射变换，调整图像亮度。


2.如权利要求1所述的DR图像亮度自适应调整方法，其特征在于，所述亮度计算步骤中，采用下式计算灰度均值b：



其中，H和W分别为待调整图像的高度和宽度，a为与图像灰阶相对应的尺度因子。


3.如权利要求2所述的DR图像亮度自适应调整方法，其特征在于，a取值如下式：





4.如权利要求1所述的DR图像亮度自适应调整方法，其特征在于，调整系数计算步骤中，通过下式计算亮度调整系数s：
s＝Cgain·b+Coffset；
其中，Cgain和Coffset分别为亮度增益系数和亮度偏移系数。


5.如权利要求1所述的DR图像亮度自适应调整方法，其特征在于，所述映射表计算步骤中，通过下式计算亮度调整映射表：



其中，k为输入灰度值，T(k)为输出灰度值，Imax和Imin分别为待调整图像I的最大值和最小值。

【专利技术属性】
技术研发人员：王永贞，郭朋，鄢照龙，郑杰，
申请(专利权)人：深圳蓝韵医学影像有限公司，
类型：发明
国别省市：广东;44