一种图像动态范围调整方法及装置的制造方法
一种图像动态范围调整方法及装置的制造方法
【技术领域】
[0001] 本发明涉及图像处理技术领域，具体涉及一种图像动态范围调整方法及装置。
【背景技术】
[0002] 在图像处理技术领域，通常用动态范围（Dynamic Range)表示图像的最亮部分和 最暗部分之间的相对比值。图像的动态范围越大，其表现的层次越丰富，所包含的色彩空间 也越广。现有技术中，存在一些高动态范围图像，例如数字X光图像，一般可用14bit来区分 图像亮度，包含大量的信息。而计算机在表示图像时，通常是用8bit来区分图像亮度的，动 态范围较小。如何将高动态范围的图像显示在普通的计算机上，并在图像中展示更丰富的 信息成为一个重要的问题。因此，需要对图像进行合理的动态范围调整，以适应显示及观察 的需要。
[0003] 现有技术中，已存在一些图像动态范围调整方法。其中，使用直方图均衡化是调整 图像动态范围的一种常见方法。它根据直方图分布形成灰度变换曲线，这种方法能很容易 使得处理后的图像直方图更平滑，但是这种方法仅考虑直方图本身的密度分布，图像对比 度会被过度增大以至于图像灰度不连续，很难适用于X线图像的处理。
[0004] 现有技术还存在一种方法，使用比较平滑过渡的非线性的S形曲线对图像灰度进 行映射，保证图像灰度连续分布，但是单纯的应用这种曲线会在一定程度上降低图像的对 比度和细节信息。
[0005] 发明人在实现本发明的过程中发现，现有技术存在的方法单纯从直方图分布和曲 线的平滑角度出发，未考虑解剖部位特征，造成某些部位被过度处理，适应性差，不能达到 很好的图像处理效果。
【发明内容】
[0006] 本发明实施例提供了一种图像动态范围调整方法及装置，可以针对不同的解剖部 位进行不同的动态范围调整，可以适应不同图像类型的显示需求，适用性好，具有更佳的显 示效果。
[0007] 为此，本发明实施例提供如下技术方案：
[0008] 第一方面，本发明实施例提供了一种图像动态范围调整方法，所述方法包括：
[0009 ]分析原始图像，获得人体组织前景图像；
[0010] 对原始图像进行多分辨率调节，获得增强图像；
[0011] 根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定非线性LUT曲 线，利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出图像。
[0012] 进一步地，所述分析原始图像，获得人体组织前景图像包括：
[0013] 确定原始图像中的遮挡区域，根据所述遮挡区域确定非遮挡区域；
[0014] 对所述非遮挡区域进行图像分割，得到人体组织前景图像和背景图像。
[0015] 进一步地，所述对原始图像进行多分辨率调节，获得增强图像包括：
[0016] 将原始图像分解成一个低频图像和多个高频图像；
[0017] 对所述低频图像的动态范围进行压缩；
[0018] 对所述高频图像进行强度均衡化处理；
[0019] 将压缩后的低频图像以及处理后的高频图像进行重构处理，得到增强图像。
[0020] 进一步地，所述对所述低频图像的动态范围进行压缩包括：
[0021] 根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应的感兴趣灰 度分布范围；
[0022]根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图像进行压缩。 [0023]进一步地，所述对所述高频图像进行强度均衡化处理包括：
[0024]获取所述尚频图像的最大幅度值；
[0025] 根据所述最大幅度值对所述高频图像进行归一化处理；
[0026] 根据不同的控制参数对所述归一化处理后的高频图像进行频率调制，以扩展小幅 度信号的强度、抑制大幅度信号的强度；
[0027]将处理后的高频图像重新映射到归一化处理前的图像的数据范围。
[0028] 进一步地，所述根据原始图像对应的解剖部位类型以及所述人体组织前景图像确 定非线性LUT曲线包括：
[0029] 对所述人体组织前景图像和所述增强图像进行预处理；
[0030] 根据原始图像对应的解剖部位类型确定参考灰度值；
[0031] 根据所述前景图像、所述参考灰度值、预设的控制参数确定非线性LUT曲线。
[0032] 进一步地，所述利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出 图像包括：
[0033] 对所述增强图像应用非线性LUT曲线，得到处理后的增强图像；
[0034] 获得所述处理后的增强图像的感兴趣区域灰度分布特征；
[0035] 判断所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预设条件；
[0036] 若不满足，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线，并对 所述增强图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的增强 图像的感兴趣区域灰度分布特征满足预设条件；
[0037] 若满足，则输出所述处理后的增强图像。
[0038] 进一步地，所述感兴趣区域灰度分布特征具体为图像动态范围比例值，所述确定 所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预设条件包括：
[0039] 判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比 例值之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足 预设条件。
[0040] 进一步地，所述方法还包括：
[0041] 若判断调整后的控制参数满足第二预设条件，则不再调整控制参数，输出所述处 理后的增强图像。
[0042] 另一方面，本发明实施例还提供了一种图像动态范围调整装置，所述装置包括：
[0043] 前景图像获得单元，用于分析原始图像，获得人体组织前景图像；
[0044] 增强图像获得单元，用于对原始图像进行多分辨率调节，获得增强图像；
[0045] 动态范围调整单元，用于根据原始图像对应的解剖部位类型以及所述人体组织前 景图像确定非线性LUT曲线，利用所述非线性LUT曲线调整所述增强图像的动态范围，获得 输出图像。
[0046] 进一步地，所述前景图像获得单元包括：
[0047] 确定单元，用于确定原始图像中的遮挡区域，根据所述遮挡区域确定非遮挡区域；
[0048] 分割单元，用于对所述非遮挡区域进行图像分割，得到人体组织前景图像和背景 图像。
[0049]进一步地，所述增强图像获得单元包括：
[0050]分解单元，用于将原始图像分解成一个低频图像和多个高频图像；
[0051] 低频压缩单元，用于对所述低频图像的动态范围进行压缩；
[0052] 高频处理单元，用于对所述高频图像进行强度均衡化处理；
[0053]重构单元，用于将压缩后的低频图像以及处理后的高频图像进行重构处理，得到 增强图像。
[0054] 进一步地，所述低频压缩单元具体用于：
[0055] 根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应的感兴趣灰 度分布范围；根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图像进行压 缩。
[0056] 进一步地，所述高频处理单元包括：
[0057]最大幅度值获得单元，用于获取所述高频图像的最大幅度值；
[0058]归一化单元，用于根据所述最大幅度值对所述高频图像进行归一化处理；
[0059] 调制单元，用于根据不同的控制参数对所述归一化处理后的高频图像进行频率调 制，以扩展小幅度信号的强度、抑制大幅度信号的强度；
[0060] 映射单元，用于将处理后的高频图像重新映射到归一化处理前的图像的数据范 围。
[0061 ] 进一步地，所述动态范围调整单元包括曲线确定单元以及调整单元，所述曲线确 定单元包括：
[0062]预处理单元，用于对所述人体组织前景图像和所述增强图像进行预处理；
[0063] 参考灰度值确定单元，用于根据原始图像对应的解剖部位类型确定参考灰度值；
[0064] 曲线生成单元，用于根据所述前景图像、所述参考灰度值、预设的控制参数确定非 线性LUT曲线。
[0065] 进一步地，所述动态范围调整单元包括曲线确定单元以及调整单元，所述调整单 元包括：
[0066] 应用单元，用于对所述增强图像应用非线性LUT曲线，得到处理后的增强图像；
[0067] 特征获得单元，用于获得所述处理后的增强图像的感兴趣区域灰度分布特征；
[0068] 判断单元，用于判断所述处理后的增强图像的感兴趣区域灰度分布特征是否满足 预设条件；
[0069] 迭代调整单元，用于若判断所述处理后的增强图像的感兴趣区域灰度分布特征不 满足预设条件，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线，并对所述 增强图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的增强图像 的感兴趣区域灰度分布特征满足预设条件；
[0070] 第一输出单元，用于若判断所述处理后的增强图像的感兴趣区域灰度分布特征满 足预设条件，输出所述处理后的增强图像。
[0071] 进一步地，所述特征获得单元获得的所述感兴趣区域灰度分布特征具体为图像动 态范围比例值，所述判断单元具体用于：
[0072] 判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比 例值之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足 预设条件。
[0073] 进一步地，所述装置还包括：
[0074] 第二输出单元，用于若判断调整后的控制参数满足第二预设条件，则不再调整控 制参数，输出所述处理后的增强图像。
[0075] 再一方面，本发明还提供了一种图像动态范围调整设备，包括：
[0076] 处理器；
[0077]用于存储所述处理器可执行指令的存储器；
[0078]其中，所述处理器被配置为：
[0079 ]分析原始图像，获得人体组织前景图像；
[0080] 对原始图像进行多分辨率调节，获得增强图像；
[0081] 根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定非线性LUT曲 线，利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出图像。
[0082] 进一步地，所示处理器还被配置用于：
[0083] 确定原始图像中的遮挡区域，根据所述遮挡区域确定非遮挡区域；
[0084] 对所述非遮挡区域进行图像分割，得到人体组织前景图像和背景图像。
[0085] 进一步地，所示处理器还被配置用于：
[0086] 将原始图像分解成一个低频图像和多个高频图像；
[0087] 对所述低频图像的动态范围进行压缩；
[0088]对所述高频图像进行强度均衡化处理；
[0089]将压缩后的低频图像以及处理后的高频图像进行重构处理，得到增强图像。
[0090] 进一步地，所示处理器还被配置用于：
[0091]根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应的感兴趣灰 度分布范围；
[0092]根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图像进行压缩。 [0093] 进一步地，所示处理器还被配置用于：
[0094]获取所述尚频图像的最大幅度值；
[0095] 根据所述最大幅度值对所述高频图像进行归一化处理；
[0096] 根据不同的控制参数对所述归一化处理后的高频图像进行频率调制，以扩展小幅 度信号的强度、抑制大幅度信号的强度；
[0097]将处理后的高频图像重新映射到归一化处理前的图像的数据范围。
[0098 ] 进一步地，所示处理器还被配置用于：
[0099]对所述人体组织前景图像和所述增强图像进行预处理；
[0100] 根据原始图像对应的解剖部位类型确定参考灰度值；
[0101] 根据所述前景图像、所述参考灰度值、预设的控制参数确定非线性LUT曲线。
[0102] 进一步地，所示处理器还被配置用于：
[0103] 对所述增强图像应用非线性LUT曲线，得到处理后的增强图像；
[0104] 获得所述处理后的增强图像的感兴趣区域灰度分布特征；
[0105] 判断所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预设条件；
[0106] 若不满足，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线，并对 所述增强图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的增强 图像的感兴趣区域灰度分布特征满足预设条件；
[0107]若满足，则输出所述处理后的增强图像。
[0108] 进一步地，所示处理器还被配置用于：
[0109] 判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比 例值之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足 预设条件。
[0110] 进一步地，所示处理器还被配置用于：
[0111] 若判断调整后的控制参数满足第二预设条件，则不再调整控制参数，输出所述处 理后的增强图像。
[0112] 本发明实施例提供的图像动态范围调整方法及装置，可以根据不同的解剖部位类 型进行不同的动态范围调整，以突出感兴趣区域的灰度分布范围，降低背景图像的干扰，并 能够提供更佳的图像对比度，显示更丰富的信息，达到更佳的显示效果。
【附图说明】
[0113] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下， 还可以根据这些附图获得其他的附图。
[0114] 图1本发明实施例提供的图像动态范围调整方法流程示意图；
[0115] 图2为本发明实施例多分辨率调节获得增强图像的示意图；
[0116] 图3为本发明实施例动态范围调整流程图；
[0117] 图4为本发明一个实施例中腰椎图像调整过程示意图；
[0118] 图5为本发明实施例提供的增强图像的直方图分布示意图；
[0119] 图6为本发明实施例提供的动态范围调整后的直方图分布示意图；
[0120] 图7为本发明实施例提供的输入前景直方图；
[0121 ]图8为本发明实施例提供的动态范围调整后的前景直方图；
[0122] 图9为本发明实施例提供的胸片增强图像示意图；
[0123] 图10为本发明实施例提供的动态范围调整后的胸片图像示意图；
[0124] 图11为本发明实施例提供的调整得到的胸片动态范围LUT曲线示意图；
[0125] 图12为本发明实施例提供的胸片增强图像直方图；
[0126] 图13为本发明实施例提供的动态范围调整后的胸片图像直方图。
[0127] 图14为本发明实施例提供的图像动态范围调整装置示意图；
[0128] 图15为本发明实施例提供的图像动态范围调整装置的实施例硬件架构图。
【具体实施方式】
[0129] 本发明实施例提供了一种图像动态范围调整方法及装置，可以针对不同的解剖部 位进行不同的动态范围调整，可以适应不同图像类型的显示需求，适用性好，具有更佳的显 示效果。
[0130] 为了使本技术领域的人员更好地理解本发明中的技术方案，下面将结合本发明实 施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施 例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通 技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护 的范围。
[0131] 参见图1，本发明实施例提供的图像动态范围调整方法流程示意图。如图1所示，所 述方法可以包括：
[0132] S101，分析原始图像，获得人体组织前景图像。
[0133] 具体地，对输入的原始图像进行分析，首先确定原始图像中的遮挡区域，根据所述 遮挡区域确定非遮挡区域。一般地，在X线图像中存在一定的遮挡区域，一般为限束器遮挡 区域。X线图像中限束器遮挡区域是指通过铅叶遮挡，避免人体受到多余辐射的非X线直接 照射区域，该区域在图像中的主要特征是灰度偏低，一般情况下遮挡区域边缘与其他部分 存在较大的灰度反差，且边缘为直线。本发明实施例具体实现时，利用灰度、边缘直线检测 相结合的方法确定限束器遮挡区域。举例说明，首先对原始图像进行边缘检测，再对边缘检 测图像检测直线(hough变换），结合限束器遮挡区域灰度与其他区域的关系判定是否属于 限束器遮挡区域。当然，本领域技术人员可以理解的是，还可以采用其他方法识别限束器遮 挡区域，在此不进行限定。
[0134] 在确定了限束器遮挡区域后，根据原始图像和束器遮挡区域，即可以确定非遮挡 区域。然后，对非遮挡区域进行图像分割，分别得到人体组织前景图像和背景图像，其中人 体组织前景图像是指人体有效组织区域，而背景是指人体组织区域外的区域。在进行前景 图像分割时，具体采用的方法可以是灵活多样的，例如基于直方图的波峰波谷检测的方法， 基于区域增长的方法，以及基于纹理特征的方法。发明人经过研究发现，在实际应用过程 中，由于图像部位的多样性，单纯的使用某一特征进行检测往往满足不了具体要求。在本发 明实施例具体实现时，对于前景背景反差较大的部位采用基于直方图波峰波谷检测的方法 进行图像分割;而对于前景与背景灰度反差较小的部位采用基于纹理特征与灰度特征结合 的方法进行图像分割。在判断适用何种方法时，具体可以根据前景图像和背景图像的平均 灰度值进行对比，从而确定采用直方图波峰波谷检测方法，还是采用纹理特征与灰度特征 结合的方法。由于后续的动态范围调整并不是把背景直接进行灰度截止，因此对于分割结 果边缘精度并不高，只要能保证前景背景分割不出错即可。因此，本领域技术人员还可以采 用其他方法实现前景图像与背景图像的分割。需要说明的是，灰度截止是指将图像中像素 值大于或小于一定阈值的点，修改成预先的设定值的方法。
[0135] S102,对原始图像进行多分辨率调节，获得增强图像。
[0136] 参见图2,为本发明实施例多分辨率调节获得增强图像的示意图。具体实现时，将 原始图像分解成一个低频图像和一系列高频图像，压缩低频图像动态范围，对高频图像分 别进行强度均衡化处理，再将这一系列处理后的低频和高频图像进行重构得到增强图像。 具体地，本发明实施例中分解、重构图像的过程，使用拉普拉斯金字塔框架，分解层数N可以 根据需要确定。如果不对低频或高频进行处理，重构后图像与原始图像相同。下面的实施例 中，以分解层数N为7层为例进行说明。下面结合附图2,对具体实现进行详细地阐述。
[0137] S102A，将原始图像分解成一个低频图像和多个高频图像。
[0138] 具体实现时，使用拉普拉斯金字塔框架将原始图像分解成一个低频图像和7个高 频图像。
[0139 ] SI 02B，对所述低频图像的动态范围进行压缩。
[0140] 在本发明实施例中，可以对低频图像预先进行预处理，例如降采样。
[0141] 在对低频图像的动态范围进行压缩时，发明人发现，若对于不同的图像，采用同一 的压缩方式，例如采用相同的压缩强度系数进行压缩，效果并不好。例如，若采用的所述压 缩强度系数太小，则会导致压缩幅度不够;若采用的所述压缩强度系数过大，则会导致对比 度被过度加强导致图像失真。因此本发明是实施例中针对不同的解剖部位图像，进行额外 的压缩补偿。具体地，根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应 的感兴趣灰度分布范围;根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图 像进行压缩。举例说明，参见表1，为各个部位类型对应的感兴趣灰度比例范围表。参照表1 确定的解剖部位类型对应的感兴趣灰度比例范围，可以按照如下三种情况进行低频动态范 围压缩：
[0142] 表1
[0143]
[0144] (1)对于感兴趣灰度分布范围分布在整个动态范围的0-40 %的情况：
[0145]
I )
[0146] 其中，Gs'表示动态范围压缩后的低频图像的像素灰度值，GS表示低频图像的像素 灰度值，BaseGray表示基准灰度，kc为压缩强度系数。??Η为对应高灰度部分的额外补偿系 数，取值范围可以[0.5，1]。
[0147] 其中，基准灰度BaseGray可以通过以下方式得到，将步骤101中得到的前景图像按 照均值降采样成与低频图像大小相同的图像，将降采样的前景部分图像与低频图像进行与 运算可得到低频前景图像，统计该低频前景图像的中间灰度值或平均灰度值即可确定该基 准灰度BaseGray。
[0148] 举例说明，参见表1，对于解剖部位类型为头部、腹部、骨盆，其感兴趣灰度分布比 例为30%，可以参照公式（1)进行压缩。当图像中的像素灰度值小于基准灰度BaseGray，则 采用普通的压缩方式；当图像中的像素灰度值大于等于基准灰度BaseGray，则采用额外的 补偿系数??Η进行压缩补偿。
[0149] (2)而对于感兴趣灰度分布在整个动态范围的60-100%的情况：
[0150]
[0151] 其中，Gs'表示动态范围压缩后的低频图像的像素灰度值，GS表示低频图像的像素 灰度值，BaseGray表示基准灰度，kc为压缩强度系数。fEL为对应低灰度部分的额外补偿系 数，取值范围可以为[0.5，1]。其中，基准灰度BaseGray的获得方式与前面介绍的相同。
[0152] 举例说明，参见表1，对于解剖部位类型为胸部，其感兴趣灰度分布比例为80%，可 以参照公式(2)进行压缩。当图像中的像素灰度值小于基准灰度BaseGray，则采用普通的压 缩方式；当图像中的像素灰度值大于等于基准灰度BaseGray，则采用额外的补偿系数fEL进 行压缩补偿。
[0153] (3)对于感兴趣灰度分布在整个动态范围的40-60%的情况：
[0154] Gs' =kcX (Gs-BaseGray)+BaseGray (3)
[0155] 其中，Gs'表示动态范围压缩后的低频图像的像素灰度值，Gs表示低频图像的像素 灰度值，BaseGray表示基准灰度，kc为压缩强度系数。
[0156] 举例说明，参见表1，对于解剖部位类型为手、足、四肢、股部、膝部，由于均是相对 不太需要关注的区域，因此不采用额外的补偿系数进行补偿，直接采用公式(3)进行压缩即 可。
[0157] 在本发明实施例中，针对不同的解剖部位类型，采用不同的方法或者不同的补偿 系数进行低频压缩处理，从而压缩了背景以及非感兴趣区域的灰度范围，突出了感兴趣区 域的灰度范围，使得感兴趣区域灰度分布更均衡。。
[0158] S102C，对所述高频图像进行强度均衡化处理。
[0159]对所述高频图像进行强度均衡化处理，意即对不同频率上的目标高频信号进行调 制，扩展小幅度信号强度，抑制幅度较大的信号的强度。根据解剖部位分布特点，所需的控 制均衡幅度会有所不同。一般而言，0-3层含有较丰富的信息，4-6层结构和层次信息丰富，7 层以上主要包含大范围亮度信息，通过调整各层次的控制参数可以突出显示感兴趣区域的 频率，从而达到更好的显示效果。
[0160] 各层次的强度均衡化处理过程类似，具体可以包括如下步骤：
[0161] S102C1，获取所述高频图像的最大幅度值。
[016 2]其中，尚频图像的最大幅度值为：
[0163] lapmax=max(abs(Lap))
[0164] 意即获取每个高频图像的灰度最大值。
[0165] S102C2,根据所述最大幅度值对所述高频图像进行归一化处理。
[0166] 根据最大幅度值将高频图像归一化至(_1，1)之间。
[0167] S102C3,根据不同的控制参数对所述归一化处理后的高频图像进行频率调制，以 扩展小幅度信号的强度、抑制大幅度信号的强度。
[0168] 具体地，由于拉普拉斯金字塔分解后噪声主要分布在低层次的高频信息中，例如 1-3层、因此频率处理过程对于1-3层高频图像按照公式(4)进行，对于3层-N层按照公式(5) 进行。N的取值范围在本发明实施例可以取值为10,也可以是大于等于7的任意整数。
[0169]
[0171] 其中，c表示归一化的高频信息，c'表示处理后的归一化高频信息，es为控制参数， 控制频率均衡强度es e [0. 5，1 ]，该值越大均衡强度越大，越小均衡强度越弱。一般的，es取 值为0.8。
[0172] 其中，cn为高频信息中噪声与细节信息之间的阈值。该阈值可根据离线采样平均 灰度与噪声的关系，也可以通过均方差统计出现频率确定。本发明实施例中根据前景平均 灰度确定相对噪声，除以lapmax可以计算得到cn值。
[0173] S102C4,将处理后的高频图像重新映射到归一化处理前的图像的数据范围。
[0174] 具体处理公式如下：
[0175] lapmax' = lapmax X c' (6)
[0176] 由于频率调制前进行了归一化处理，因此，这里需要通过重新映射处理，使得高频 图像重新映射到归一化处理前的数据范围，以使得高频图像的灰度范围与原始高频图像的 灰度范围相同。
[0177] S102D，将压缩后的低频图像以及处理后的高频图像进行重构处理，得到增强图 像。
[0178] 具体的重构方法，可以采用拉普拉斯金字塔框架进行重构。
[0179] S103,根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定非线性 LUT(英文全称为Look Up Table，中文名称为查找表）曲线，利用所述非线性LUT曲线调整所 述增强图像的动态范围，获得输出图像。
[0180] 图像动态范围调整目的主要包括两方面一是于压缩背景灰度，二是调整前景图像 灰度分布。对于背景灰度分布，按照一定比例需要压缩其动态范围，而对于前景区域则依据 解剖部位特征，扩展感兴趣灰度范围，压缩其余的灰度部分，并且避免前景中非感兴趣出现 过明或过暗的情况。基于此，本发明实施例根据原始图像对应的解剖部位类型以及所述人 体组织前景图像确定了非线性LUT曲线，并利用所述非线性LUT曲线调整所述增强图像的动 态范围，获得具有更高对比度的输出图像。
[0181] 具体实现时，根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定 非线性LUT曲线可以包括:对所述人体组织前景图像和所述增强图像进行预处理;根据原始 图像对应的解剖部位类型确定参考灰度值;根据所述前景图像、所述参考灰度值、预设的控 制参数确定非线性LUT曲线。
[0182] 所述利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出图像包括： 对所述增强图像应用非线性LUT曲线，得到处理后的增强图像;获得所述处理后的增强图像 的感兴趣区域灰度分布特征;判断所述处理后的增强图像的感兴趣区域灰度分布特征是否 满足预设条件;若不满足，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线， 并对所述增强图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的 增强图像的感兴趣区域灰度分布特征满足预设条件;若满足，则输出所述处理后的增强图 像。
[0183] 下面结合附图3,对本发明实施例中动态范围调整流程图进行说明。
[0184] S301，对输入信息进行预处理，结合解剖部位类型确定参考灰度。
[0185] 对输入信息的图像预处理过程可以包含增强图像和前景图像的降采样、干扰信息 排除等处理。首先对增强图像S和前景图像Μ降采样。降采样过程为了改善产生LUT效率，而 最终LUT应用在整幅图像上。由于DR图像分辨率较大，为了提升处理效率需要将图像按照一 定采样率抽样，优选地，一个可行的实现方式中，可以设置降采样比率为8:1，得到降采样图 像增强图像S1，降采样前景图Ml。需要注意的是接下来的曲线形成过程都是基于降采样的 S1和Ml进行的，但最终生成处理后的数据是在原始大尺寸图像上应用LUT得到的。
[0186] S302,根据初始参数构造非线性曲线。
[0187] 按照初始设定形成灰度映射曲线，对增强图像S1应用灰度映射曲线，得到处理图 像SP1，LUT曲线可以通过以下方法得到：
[0188] 定义如下几个控制点pl(x0，y0)，p2(xl，yl)，p3(x2，y2)，p4(x3，y3)，通过非线性 曲线拟合可以得到映射曲线，拟合方式可选分段多项式方法，注意x〇，xl，x2，x3和yO，y 1， y 2，y 3都呈升序分布，避免形成异常曲线。
[0189] χ〇 = 〇
[0190] y〇 = 〇
[0191] xl =bcg_gray
[0192] yl=klX(xl-xO)+yO
[0193] x2 = ref_gray
[0194] y2 = k2X(x2-xl)+yl
[0195] x3=max_gray
[0196] y3 = k3X(x3-x2)+y2
[0197] 其中，kl，k2，k3分别表示p2，p3，p4的控制参数，kl，k2，k3取值范围分别是[0_ 0·2]，[0·1-2]，[0·5-2]。
[0198] 其中，bcg_gray表示背景与前景区域分割阈值，增强后的图像前景和背景灰度分 布会存在一定交叉，可以选择该交叉点所对应的灰度值作为该分割阈值，也可以选择直接 在前景图像中，按照一定比例计算得到。其中ref_gray表示参考灰度，参考灰度需要根据具 体图像来确定的，根据大量数据统计，相同的摄影部位和体位在灰度分布存在共同的特点。 因此本发明实施例根据部位类型和前景图像，按照预定比例计算到该灰度值，表2给出了几 种典型部位参考灰度在前景直方图的比例分布。
[0199] 表2不同部位对应的特征灰度分布
[0200]
[0201]
[0202]其中所述按照预定比例计算灰度值，具体的做法是:先计算前景图像直方图，统计 出前景图像中像素灰度总数，乘以设定好的比例值得到累加器门限值，初始化累加器，遍历 直方图所有通道，依次将直方图中各个通道对应的统计值加入累加器，当超过累加器门限 值时，该通道对应的灰度值即是按照一定比例得到的灰度值。
[0203] S303,对增强图像应用非线性曲线得到处理后图像。
[0204]具体地，对于增强图像的每一个像素点p(x，y)，都有p'（x，y)=LUT(p(x，y))，得到 的P'就是应用了 LUT的处理结果。
[0205] S304，统计图像感兴趣区域灰度分布特征。
[0206]具体实现时，所述感兴趣区域灰度分布特征具体为图像动态范围比例值。
[0207]举例说明，步骤304统计处理后图像中感兴趣区域灰度分布特征，主要包括计算处 理后图像直方图中参考点左侧动态范围和右侧动态范围的比例，可以按照如下公式得到：
[0208] (6)
[0209] 其中，dynamic_ratio为图像动态范围比例值，data_pro_ref，data_pro_left， da ta_pr o_r i gh t可以按照步骤3 0 2中给出的按照预定比例计算灰度值的方法确定，为了排 除可能的异常像素的干扰，可以为data_pro_left，data_pro_right设定预定比例值。举例 说明，其中data_pro_left，data_pro_right对应的预定比例分别为1 %，99%，data_pro_ left，data_pro_right对应的预定比例也可以分别是0.05%，99.5%，或者分别是2%、 98 %。需要说明的是，所述预定比例值为根据经验选取，一般地，data_pro_lef t对应的预定 比例一般不大于2%，data_pro_right对应的预定比例一般不小于98%，否则可能会导致图 像失真，例如会损失一些图像细节。本领域技术人员可以根据经验选取所述预定比例值，在 此不进行限定。而data_pr 〇_ref对应的预设比例值则取决于输入的部位类型信息，可以根 据表2得到。
[0210] S305,检查感兴趣区域灰度分布特征是否满足预设条件。如果是，则执行S307,得 到最终处理后的图像。否则执行S306,调整参数构造新的非线性曲线，重复步骤303-305,直 到感兴趣区域灰度分布特征满足预设条件。
[0211] 这一步骤是调整LUT曲线的过程。具体地，判断处理后图像中感兴趣区域灰度分布 特征是否满足需求，如果达到则可以得到最终处理后图像，否则需要调整控制参数后重新 形成LUT曲线。
[0212] 所述确定所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预设条件 包括:判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比例值 之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足预设 条件。
[0213] 举例说明，判断调整后的图像动态范围特征比例值dynamic_ratio和目标动态范 围特征比例值dynamic_target偏差的绝对值是否小于预设门限值threshold，本发明实施 例中threshold取0.1，dynamic_target固定取值为1，如果小于该门限值，表明感兴趣灰度 范围满足需求，否则需要调整参数后重新形成新的LUT曲线。
[0214] S306,调整参数构造新的非线性曲线，重复步骤303-305,直到感兴趣区域灰度分 布特征满足预设条件。
[0215] 当确定不满足预设条件时，则需要调整控制参数。优选地，只需要调整控制参数 k2，kl和k3直接可以使用预设值，不需要自动调整。
[0216] 其中，k2的调整方式如下：
[0217] 如果dynamic_ratio_dynamic_target> threshold，则k2 = k2_k_step，
[0218] 如果dynamic_ratio_dynamic_target<-threshold，则k2 = k2+k_step
[0219] 举例说明，若调整后的图像动态范围特征比例值dynamic_ratio和目标动态范围 特征比例值dynamic_target的差值大于预设门限值threshold，则使用控制参数减去调整 步长k_step，得到新的控制参数。若调整后的图像动态范围特征比例值dynamic_ratio和目 标动态范围特征比例值dynamic_target的差值小于-threshold，则使用控制参数加上调整 步长k_step，得到新的控制参数。
[0220]优选地，若调整后的控制参数满足第二预设条件，则不再调整控制参数，输出所述 处理后的增强图像。具体地，为了避免无意义的迭代，当出现以下情况时，不需要再调整k2 值。第二预设条件具体为：
[0221] (l)k2已超出有效范围
[0222] k2<0.1 或 k2>2
[0223] 需要说明的是，k2的有效取值范围为根据经验选取。经过发明人反复研究发现，k2 在[0.1，2]的取值范围时，图像的调整效果较好。若k2超出有效范围，调整后的图像则会出 现一定程度的失真现象。因此，为了避免图像失真，需要对k2的有效范围进行判断，以达到 更好的图像处理效果。< br>[0224] (2)出现震荡情况。也就是说，迭代前后出现与目标动态范围偏差异常，即同时满 足如下条件：
[0225] dy nam i c_ra t i o (n) -dy nam i c_tar ge t < 0
[0226] dy nam i c_ra t i o (n+1) -dy nam i c_tar ge t < 0
[0227] abs(dynamic_ratio(n)-dynamic_target)>0.1
[0228] abs(dynamic_ratio(n+l)-dynamic_target)>0.1
[0229] 其中，dynamic_ratio(n)表示第n次调整k2值后的计算结果；
[0230] dynamic_ratio(n+l)表示第n+1次调整k2后的计算结果。
[0231] S307,得到最终处理后的输出图像。
[0232] 最后将调整后的LUT应用到增强的全尺寸图像，就可以得到最终动态范围调整后 的图像。
[0233] 下面结合附图4-13对本发明实施例能够达到的有益效果进行说明。
[0234] 参见图4,为本发明一个实施例中腰椎图像调整过程示意图，左侧是多分辨率调整 后的增强图像，右侧是调c整好的输出图像，该调整过程难点就在于如何确定合理的LUT曲 线。
[0235] 当图像输出到显示设备上时，根据亮度对比度设置将灰度映射成亮度信息，因此 能够被很好显示的仅是窗宽范围内的目标，因此感兴趣的灰度能尽量均匀分布在整个动态 范围中间，才会更便于观察病灶。图5是增强图像直方图分布，图中给出了背景、皮肤、肺部、 肌肉、低密度、高密度骨骼等一系列组织灰度分布，背景，肺，皮肤组织占据了很宽的灰度范 围，高密度骨骼组织灰度分布过度则集中，动态范围很窄。而图6是经过动态范围调整后的 直方图分布，背景、肺、皮肤组织灰度都被压缩至很小的范围，高密度和低密度骨骼肌肉灰 度分布范围都不同程度被扩展。因此可以看到图6中，按照前景灰度最小最大范围设置窗宽 窗位进行灰度映射的显示效果中，变换前的图像发白，感兴趣区域对比度不好，而变换后的 图像感兴趣区域对比度合适，更适合观察病灶。图7和图8分别是变换前后前景区域，1024级 直方图分布，二者比较可以看出变换对于灰度分布集中程度的改善。
[0236] 图9-13给出了一种胸片图像的应用实例，图9给出了胸片增强图像，图10是经过动 态范围调整后的胸片图像，图11是经过调整得到的灰度映射曲线，图12和13分别是图像9和 图像10对应的直方图分布，从结果中也可以看到曲线调整后图像对比度改善，从直方图变 化中也能看到肺部灰度动态范围被扩展，更有利于肺部组织观察。
[0237] 在本发明实施例中，能够达到如下有益效果：
[0238] (1)对处理结果进行分析判定，调整曲线生成，更利于达到预期的效果。
[0239] (2)综合考虑解剖部位特征、避免通用基于直方图盲调方法造成的对比过强或灰 度分布不连续等问题，抑制增强带来的尾影。
[0240] (3)能突出显示感兴趣区域组织灰度，又保证其他前景组织被观察到，增加图像可 观测的信息量。
[0241] 本发明实施例经过大量图像样本验证，都可以获得比较满意的图像动态范围。
[0242] 与前述图像动态范围调整方法的实施例相对应，本发明还提供了图像动态范围调 整装置的实施例。
[0243] 参照图14,为本发明实施例提供的图像动态范围调整装置框图。
[0244] 所示图像动态范围调整装置1400可以包括：
[0245] 前景图像获得单元1401，用于分析原始图像，获得人体组织前景图像。
[0246] 增强图像获得单元1402,用于对原始图像进行多分辨率调节，获得增强图像。
[0247] 动态范围调整单元1403,用于根据原始图像对应的解剖部位类型以及所述人体组 织前景图像确定非线性LUT曲线，利用所述非线性LUT曲线调整所述增强图像的动态范围， 获得输出图像。
[0248] 进一步地，所述前景图像获得单元包括：
[0249] 确定单元，用于确定原始图像中的遮挡区域，根据所述遮挡区域确定非遮挡区域；
[0250] 分割单元，用于对所述非遮挡区域进行图像分割，得到人体组织前景图像和背景 图像。
[0251] 进一步地，所述增强图像获得单元包括：
[0252] 分解单元，用于将原始图像分解成一个低频图像和多个高频图像；
[0253] 低频压缩单元，用于对所述低频图像的动态范围进行压缩；
[0254]高频处理单元，用于对所述高频图像进行强度均衡化处理；
[0255] 重构单元，用于将压缩后的低频图像以及处理后的高频图像进行重构处理，得到 增强图像。
[0256] 进一步地，所述低频压缩单元具体用于：
[0257] 根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应的感兴趣灰 度分布范围；根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图像进行压 缩。
[0258] 进一步地，所述高频处理单元包括：
[0259] 最大幅度值获得单元，用于获取所述高频图像的最大幅度值；
[0260] 归一化单元，用于根据所述最大幅度值对所述高频图像进行归一化处理；
[0261] 调制单元，用于根据不同的控制参数对所述归一化处理后的高频图像进行频率调 制，以扩展小幅度信号的强度、抑制大幅度信号的强度；
[0262]映射单元，用于将处理后的高频图像重新映射到归一化处理前图像的数据范围。 [0263] 进一步地，所述动态范围调整单元包括曲线确定单元以及调整单元，所述曲线确 定单元包括：
[0264]预处理单元，用于对所述人体组织前景图像和所述增强图像进行预处理；
[0265] 参考灰度值确定单元，用于根据原始图像对应的解剖部位类型确定参考灰度值；
[0266] 曲线生成单元，用于根据所述前景图像、所述参考灰度值、预设的控制参数确定非 线性LUT曲线。
[0267] 进一步地，所述动态范围调整单元包括曲线确定单元以及调整单元，所述调整单 元包括：
[0268] 应用单元，用于对所述增强图像应用非线性LUT曲线，得到处理后的增强图像；
[0269] 特征获得单元，用于获得所述处理后的增强图像的感兴趣区域灰度分布特征；
[0270] 判断单元，用于判断所述处理后的增强图像的感兴趣区域灰度分布特征是否满足 预设条件；
[0271] 迭代调整单元，用于若判断所述处理后的增强图像的感兴趣区域灰度分布特征不 满足预设条件，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线，并对所述 增强图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的增强图像 的感兴趣区域灰度分布特征满足预设条件；
[0272] 第一输出单元，用于若判断所述处理后的增强图像的感兴趣区域灰度分布特征满 足预设条件，输出所述处理后的增强图像。
[0273] 进一步地，所述特征获得单元获得的所述感兴趣区域灰度分布特征具体为图像动 态范围比例值，所述判断单元具体用于：
[0274] 判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比 例值之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足 预设条件。
[0275] 进一步地，所述装置还包括：
[0276] 第二输出单元，用于若判断调整后的控制参数满足第二预设条件，则不再调整控 制参数，输出所述处理后的增强图像
[0277] 其中，本发明装置各单元或模块的设置可以参照图1所示的方法而实现，在此不赘 述。
[0278] 本发明图像动态范围调整装置的实施例可以应用在DR系统的设备或者任何具有 处理器(特别是图像处理器）的电子设备上，所述电子设备可以是现有的、正在研发的或将 来研发的任何电子设备，包括但不限于:现有的、正在研发的或将来研发的台式计算机、膝 上型计算机、移动终端(包括智能手机、非智能手机、各种平板电脑)等。装置实施例可以通 过软件实现，也可以通过硬件或者软硬件结合的方式实现。以软件实现为例，作为一个逻辑 意义上的装置，是通过其所在DR系统的设备或带有处理器的电子设备的处理器将存储器中 对应的计算机程序指令读取到内存中运行形成的。从硬件层面而言，如图15所示，为本发明 图像动态范围调整装置所在DR系统的设备或带有处理器的电子设备的一种硬件结构图，除 了图15所示的处理器、内存、网络接口、以及存储器之外，实施例中装置所在的DR系统的设 备或带有处理器的电子设备通常根据该设备的实际功能，还可以包括其他硬件，对此不再 赘述。
[0279] 其中，存储器中可以存储有图像动态范围调整方法对应的逻辑指令，该存储器例 如可以是非易失性存储器(non-volatile memory)。处理器可以调用执行存储器中的保存 的逻辑指令，以执行上述的图像动态范围调整方法。
[0280] 图像动态范围调整方法对应的逻辑指令的功能，如果以软件功能单元的形式实现 并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的 理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分 可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指 令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各 个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器 (ROM，Read_Only Memory)、随机存取存储器（RAM，Random Access Memory)、磁碟或者光盘 等各种可以存储程序代码的介质。
[0281] 如图15所示，本发明还提供了 一种图像动态范围调整设备，包括：
[0282] 处理器；
[0283] 用于存储所述处理器可执行指令的存储器；
[0284] 其中，所述处理器被配置为：
[0285] 分析原始图像，获得人体组织前景图像；
[0286] 对原始图像进行多分辨率调节，获得增强图像；
[0287] 根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定非线性LUT曲 线，利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出图像。
[0288] 进一步地，所示处理器还被配置用于：
[0289] 确定原始图像中的遮挡区域，根据所述遮挡区域确定非遮挡区域；
[0290] 对所述非遮挡区域进行图像分割，得到人体组织前景图像和背景图像。
[0291] 进一步地，所示处理器还被配置用于：
[0292] 将原始图像分解成一个低频图像和多个高频图像；
[0293] 对所述低频图像的动态范围进行压缩；
[0294] 对所述高频图像进行强度均衡化处理；
[0295] 将压缩后的低频图像以及处理后的高频图像进行重构处理，得到增强图像。
[0296] 进一步地，所示处理器还被配置用于：
[0297] 根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应的感兴趣灰 度分布范围；
[0298] 根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图像进行压缩。
[0299] 进一步地，所示处理器还被配置用于：
[0300]获取所述尚频图像的最大幅度值；
[0301 ]根据所述最大幅度值对所述高频图像进行归一化处理；
[0302]根据不同的控制参数对所述归一化处理后的高频图像进行频率调制，以扩展小幅 度信号的强度、抑制大幅度信号的强度；
[0303]将处理后的高频图像重新映射到归一化处理前的图像的数据范 围。
[0304] 进一步地，所示处理器还被配置用于：
[0305] 对所述人体组织前景图像和所述增强图像进行预处理；
[0306] 根据原始图像对应的解剖部位类型确定参考灰度值；
[0307]根据所述前景图像、所述参考灰度值、预设的控制参数确定非线性LUT曲线。
[0308] 进一步地，所示处理器还被配置用于：
[0309] 对所述增强图像应用非线性LUT曲线，得到处理后的增强图像；
[0310] 获得所述处理后的增强图像的感兴趣区域灰度分布特征；
[0311] 判断所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预设条件；
[0312] 若不满足，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线，并对 所述增强图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的增强 图像的感兴趣区域灰度分布特征满足预设条件；
[0313] 若满足，则输出所述处理后的增强图像。
[0314] 进一步地，所示处理器还被配置用于：
[0315] 判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比 例值之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足 预设条件。
[0316] 进一步地，所示处理器还被配置用于：
[0317] 若判断调整后的控制参数满足第二预设条件，则不再调整控制参数，输出所述处 理后的增强图像。
[0318] 本领域技术人员可以理解的是，以上对方法和装置实施例进行了示例性说明，以 上不视为对本发明的限制，本领域技术人员在不付出创造性劳动下获得的其他实现方式均 属于本发明的保护范围。
[0319]需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实 体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存 在任何这种实际的关系或者顺序。而且，术语"包括"、"包含"或者其任何其他变体意在涵盖 非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要 素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备 所固有的要素。在没有更多限制的情况下，由语句"包括一个……"限定的要素，并不排除在 包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。本发明可以在由计算 机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执 行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分 布式计算环境中实践本发明，在这些分布式计算环境中，由通过通信网络而被连接的远程 处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地 和远程计算机存储介质中。
[0320]本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部 分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于装置实 施例而言，由于其基本相似于方法实施例，所以描述得比较简单，相关之处参见方法实施例 的部分说明即可。以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明 的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是 物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要 选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出 创造性劳动的情况下，即可以理解并实施。以上所述仅是本发明的【具体实施方式】，应当指 出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干 改进和润饰，这些改进和润饰也应视为本发明的保护范围。
【主权项】
1. 一种图像动态范围调整方法，其特征在于，所述方法包括： 分析原始图像，获得人体组织前景图像； 对原始图像进行多分辨率调节，获得增强图像； 根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定非线性LUT曲线， 利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出图像。2. 根据权利要求1所述的方法，其特征在于，所述分析原始图像，获得人体组织前景图 像包括： 确定原始图像中的遮挡区域，根据所述遮挡区域确定非遮挡区域； 对所述非遮挡区域进行图像分割，得到人体组织前景图像和背景图像。3. 根据权利要求1所述的方法，其特征在于，所述对原始图像进行多分辨率调节，获得 增强图像包括： 将原始图像分解成一个低频图像和多个高频图像； 对所述低频图像的动态范围进行压缩； 对所述高频图像进行强度均衡化处理； 将压缩后的低频图像以及处理后的高频图像进行重构处理，得到增强图像。4. 根据权利要求3所述的方法，其特征在于，所述对所述低频图像的动态范围进行压缩 包括： 根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应的感兴趣灰度分 布范围； 根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图像进行压缩。5. 根据权利要求3所述的方法，其特征在于，所述对所述高频图像进行强度均衡化处理 包括： 获取所述高频图像的最大幅度值； 根据所述最大幅度值对所述高频图像进行归一化处理； 根据不同的控制参数对所述归一化处理后的高频图像进行频率调制，以扩展小幅度信 号的强度、抑制大幅度信号的强度； 将处理后的高频图像重新映射到归一化处理前的图像的数据范围。6. 根据权利要求1所述的方法，其特征在于，所述根据原始图像对应的解剖部位类型以 及所述人体组织前景图像确定非线性LUT曲线包括： 对所述人体组织前景图像和所述增强图像进行预处理； 根据原始图像对应的解剖部位类型确定参考灰度值； 根据所述前景图像、所述参考灰度值、预设的控制参数确定非线性LUT曲线。7. 根据权利要求1或6所述的方法，其特征在于，所述利用所述非线性LUT曲线调整所述 增强图像的动态范围，获得输出图像包括： 对所述增强图像应用非线性LUT曲线，得到处理后的增强图像； 获得所述处理后的增强图像的感兴趣区域灰度分布特征； 判断所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预设条件； 若不满足，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线，并对所述 增强图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的增强图像 的感兴趣区域灰度分布特征满足预设条件； 若满足，则输出所述处理后的增强图像。8. 根据权利要求7所述的方法，其特征在于，所述感兴趣区域灰度分布特征具体为图像 动态范围比例值，所述确定所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预 设条件包括： 判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比例值 之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足预设 条件。9. 根据权利要求7所述的方法，其特征在于，所述方法还包括： 若判断调整后的控制参数满足第二预设条件，则不再调整控制参数，输出所述处理后 的增强图像。10. -种图像动态范围调整装置，其特征在于，所述装置包括： 前景图像获得单元，用于分析原始图像，获得人体组织前景图像； 增强图像获得单元，用于对原始图像进行多分辨率调节，获得增强图像； 动态范围调整单元，用于根据原始图像对应的解剖部位类型以及所述人体组织前景图 像确定非线性LUT曲线，利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出 图像。11. 根据权利要求10所述的装置，其特征在于，所述前景图像获得单元包括： 确定单元，用于确定原始图像中的遮挡区域，根据所述遮挡区域确定非遮挡区域； 分割单元，用于对所述非遮挡区域进行图像分割，得到人体组织前景图像和背景图像。12. 根据权利要求10所述的装置，其特征在于，所述增强图像获得单元包括： 分解单元，用于将原始图像分解成一个低频图像和多个高频图像； 低频压缩单元，用于对所述低频图像的动态范围进行压缩； 高频处理单元，用于对所述高频图像进行强度均衡化处理； 重构单元，用于将压缩后的低频图像以及处理后的高频图像进行重构处理，得到增强 图像。13. 根据权利要求12所述的装置，其特征在于，所述低频压缩单元具体用于： 根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应的感兴趣灰度分 布范围;根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图像进行压缩。14. 根据权利要求12所述的装置，其特征在于，所述高频处理单元包括： 最大幅度值获得单元，用于获取所述高频图像的最大幅度值； 归一化单元，用于根据所述最大幅度值对所述高频图像进行归一化处理； 调制单元，用于根据不同的控制参数对所述归一化处理后的高频图像进行频率调制， 以扩展小幅度信号的强度、抑制大幅度信号的强度； 映射单元，用于将处理后的高频图像重新映射到归一化处理前的图像的数据范围。15. 根据权利要求10所述的装置，其特征在于，所述动态范围调整单元包括曲线确定单 元以及调整单元，所述曲线确定单元包括： 预处理单元，用于对所述人体组织前景图像和所述增强图像进行预处理； 参考灰度值确定单元，用于根据原始图像对应的解剖部位类型确定参考灰度值； 曲线生成单元，用于根据所述前景图像、所述参考灰度值、预设的控制参数确定非线性LUT曲线。16. 根据权利要求10或15所述的装置，其特征在于，所述动态范围调整单元包括曲线确 定单元以及调整单元，所述调整单元包括： 应用单元，用于对所述增强图像应用非线性LUT曲线，得到处理后的增强图像； 特征获得单元，用于获得所述处理后的增强图像的感兴趣区域灰度分布特征； 判断单元，用于判断所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预设 条件； 迭代调整单元，用于若判断所述处理后的增强图像的感兴趣区域灰度分布特征不满足 预设条件，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线，并对所述增强 图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的增强图像的感 兴趣区域灰度分布特征满足预设条件； 第一输出单元，用于若判断所述处理后的增强图像的感兴趣区域灰度分布特征满足预 设条件，输出所述处理后的增强图像。17. 根据权利要求16所述的装置，其特征在于，所述特征获得单元获得的所述感兴趣区 域灰度分布特征具体为图像动态范围比例值，所述判断单元具体用于： 判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比例值 之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足预设 条件。18. 根据权利要求16所述的装置，其特征在于，所述装置还包括： 第二输出单元，用于若判断调整后的控制参数满足第二预设条件，则不再调整控制参 数，输出所述处理后的增强图像。19. 一种图像动态范围调整设备，其特征在于，包括： 处理器； 用于存储所述处理器可执行指令的存储器； 其中，所述处理器被配置为： 分析原始图像，获得人体组织前景图像； 对原始图像进行多分辨率调节，获得增强图像； 根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定非线性LUT曲线， 利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出图像。
【专利摘要】本发明实施例提供一种图像动态范围调整方法和装置，所述方法包括：分析原始图像，获得人体组织前景图像；对原始图像进行多分辨率调节，获得增强图像；根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定非线性LUT曲线，利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出图像。本发明实施例可以根据不同的解剖部位类型进行不同的动态范围调整，以突出感兴趣区域的灰度分布范围，降低背景图像的干扰，并能够提供更佳的图像对比度，显示更丰富的信息，达到更佳的显示效果。
【IPC分类】G06T5/00
【公开号】CN105488765
【申请号】CN201510818564
【发明人】李海春
【申请人】沈阳东软医疗系统有限公司
【公开日】2016年4月13日
【申请日】2015年11月23日
【技术领域】
[0001] 本发明涉及图像处理技术领域，具体涉及一种图像动态范围调整方法及装置。
【背景技术】
[0002] 在图像处理技术领域，通常用动态范围（Dynamic Range)表示图像的最亮部分和 最暗部分之间的相对比值。图像的动态范围越大，其表现的层次越丰富，所包含的色彩空间 也越广。现有技术中，存在一些高动态范围图像，例如数字X光图像，一般可用14bit来区分 图像亮度，包含大量的信息。而计算机在表示图像时，通常是用8bit来区分图像亮度的，动 态范围较小。如何将高动态范围的图像显示在普通的计算机上，并在图像中展示更丰富的 信息成为一个重要的问题。因此，需要对图像进行合理的动态范围调整，以适应显示及观察 的需要。
[0003] 现有技术中，已存在一些图像动态范围调整方法。其中，使用直方图均衡化是调整 图像动态范围的一种常见方法。它根据直方图分布形成灰度变换曲线，这种方法能很容易 使得处理后的图像直方图更平滑，但是这种方法仅考虑直方图本身的密度分布，图像对比 度会被过度增大以至于图像灰度不连续，很难适用于X线图像的处理。
[0004] 现有技术还存在一种方法，使用比较平滑过渡的非线性的S形曲线对图像灰度进 行映射，保证图像灰度连续分布，但是单纯的应用这种曲线会在一定程度上降低图像的对 比度和细节信息。
[0005] 发明人在实现本发明的过程中发现，现有技术存在的方法单纯从直方图分布和曲 线的平滑角度出发，未考虑解剖部位特征，造成某些部位被过度处理，适应性差，不能达到 很好的图像处理效果。
【发明内容】
[0006] 本发明实施例提供了一种图像动态范围调整方法及装置，可以针对不同的解剖部 位进行不同的动态范围调整，可以适应不同图像类型的显示需求，适用性好，具有更佳的显 示效果。
[0007] 为此，本发明实施例提供如下技术方案：
[0008] 第一方面，本发明实施例提供了一种图像动态范围调整方法，所述方法包括：
[0009 ]分析原始图像，获得人体组织前景图像；
[0010] 对原始图像进行多分辨率调节，获得增强图像；
[0011] 根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定非线性LUT曲 线，利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出图像。
[0012] 进一步地，所述分析原始图像，获得人体组织前景图像包括：
[0013] 确定原始图像中的遮挡区域，根据所述遮挡区域确定非遮挡区域；
[0014] 对所述非遮挡区域进行图像分割，得到人体组织前景图像和背景图像。
[0015] 进一步地，所述对原始图像进行多分辨率调节，获得增强图像包括：
[0016] 将原始图像分解成一个低频图像和多个高频图像；
[0017] 对所述低频图像的动态范围进行压缩；
[0018] 对所述高频图像进行强度均衡化处理；
[0019] 将压缩后的低频图像以及处理后的高频图像进行重构处理，得到增强图像。
[0020] 进一步地，所述对所述低频图像的动态范围进行压缩包括：
[0021] 根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应的感兴趣灰 度分布范围；
[0022]根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图像进行压缩。 [0023]进一步地，所述对所述高频图像进行强度均衡化处理包括：
[0024]获取所述尚频图像的最大幅度值；
[0025] 根据所述最大幅度值对所述高频图像进行归一化处理；
[0026] 根据不同的控制参数对所述归一化处理后的高频图像进行频率调制，以扩展小幅 度信号的强度、抑制大幅度信号的强度；
[0027]将处理后的高频图像重新映射到归一化处理前的图像的数据范围。
[0028] 进一步地，所述根据原始图像对应的解剖部位类型以及所述人体组织前景图像确 定非线性LUT曲线包括：
[0029] 对所述人体组织前景图像和所述增强图像进行预处理；
[0030] 根据原始图像对应的解剖部位类型确定参考灰度值；
[0031] 根据所述前景图像、所述参考灰度值、预设的控制参数确定非线性LUT曲线。
[0032] 进一步地，所述利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出 图像包括：
[0033] 对所述增强图像应用非线性LUT曲线，得到处理后的增强图像；
[0034] 获得所述处理后的增强图像的感兴趣区域灰度分布特征；
[0035] 判断所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预设条件；
[0036] 若不满足，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线，并对 所述增强图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的增强 图像的感兴趣区域灰度分布特征满足预设条件；
[0037] 若满足，则输出所述处理后的增强图像。
[0038] 进一步地，所述感兴趣区域灰度分布特征具体为图像动态范围比例值，所述确定 所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预设条件包括：
[0039] 判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比 例值之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足 预设条件。
[0040] 进一步地，所述方法还包括：
[0041] 若判断调整后的控制参数满足第二预设条件，则不再调整控制参数，输出所述处 理后的增强图像。
[0042] 另一方面，本发明实施例还提供了一种图像动态范围调整装置，所述装置包括：
[0043] 前景图像获得单元，用于分析原始图像，获得人体组织前景图像；
[0044] 增强图像获得单元，用于对原始图像进行多分辨率调节，获得增强图像；
[0045] 动态范围调整单元，用于根据原始图像对应的解剖部位类型以及所述人体组织前 景图像确定非线性LUT曲线，利用所述非线性LUT曲线调整所述增强图像的动态范围，获得 输出图像。
[0046] 进一步地，所述前景图像获得单元包括：
[0047] 确定单元，用于确定原始图像中的遮挡区域，根据所述遮挡区域确定非遮挡区域；
[0048] 分割单元，用于对所述非遮挡区域进行图像分割，得到人体组织前景图像和背景 图像。
[0049]进一步地，所述增强图像获得单元包括：
[0050]分解单元，用于将原始图像分解成一个低频图像和多个高频图像；
[0051] 低频压缩单元，用于对所述低频图像的动态范围进行压缩；
[0052] 高频处理单元，用于对所述高频图像进行强度均衡化处理；
[0053]重构单元，用于将压缩后的低频图像以及处理后的高频图像进行重构处理，得到 增强图像。
[0054] 进一步地，所述低频压缩单元具体用于：
[0055] 根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应的感兴趣灰 度分布范围；根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图像进行压 缩。
[0056] 进一步地，所述高频处理单元包括：
[0057]最大幅度值获得单元，用于获取所述高频图像的最大幅度值；
[0058]归一化单元，用于根据所述最大幅度值对所述高频图像进行归一化处理；
[0059] 调制单元，用于根据不同的控制参数对所述归一化处理后的高频图像进行频率调 制，以扩展小幅度信号的强度、抑制大幅度信号的强度；
[0060] 映射单元，用于将处理后的高频图像重新映射到归一化处理前的图像的数据范 围。
[0061 ] 进一步地，所述动态范围调整单元包括曲线确定单元以及调整单元，所述曲线确 定单元包括：
[0062]预处理单元，用于对所述人体组织前景图像和所述增强图像进行预处理；
[0063] 参考灰度值确定单元，用于根据原始图像对应的解剖部位类型确定参考灰度值；
[0064] 曲线生成单元，用于根据所述前景图像、所述参考灰度值、预设的控制参数确定非 线性LUT曲线。
[0065] 进一步地，所述动态范围调整单元包括曲线确定单元以及调整单元，所述调整单 元包括：
[0066] 应用单元，用于对所述增强图像应用非线性LUT曲线，得到处理后的增强图像；
[0067] 特征获得单元，用于获得所述处理后的增强图像的感兴趣区域灰度分布特征；
[0068] 判断单元，用于判断所述处理后的增强图像的感兴趣区域灰度分布特征是否满足 预设条件；
[0069] 迭代调整单元，用于若判断所述处理后的增强图像的感兴趣区域灰度分布特征不 满足预设条件，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线，并对所述 增强图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的增强图像 的感兴趣区域灰度分布特征满足预设条件；
[0070] 第一输出单元，用于若判断所述处理后的增强图像的感兴趣区域灰度分布特征满 足预设条件，输出所述处理后的增强图像。
[0071] 进一步地，所述特征获得单元获得的所述感兴趣区域灰度分布特征具体为图像动 态范围比例值，所述判断单元具体用于：
[0072] 判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比 例值之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足 预设条件。
[0073] 进一步地，所述装置还包括：
[0074] 第二输出单元，用于若判断调整后的控制参数满足第二预设条件，则不再调整控 制参数，输出所述处理后的增强图像。
[0075] 再一方面，本发明还提供了一种图像动态范围调整设备，包括：
[0076] 处理器；
[0077]用于存储所述处理器可执行指令的存储器；
[0078]其中，所述处理器被配置为：
[0079 ]分析原始图像，获得人体组织前景图像；
[0080] 对原始图像进行多分辨率调节，获得增强图像；
[0081] 根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定非线性LUT曲 线，利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出图像。
[0082] 进一步地，所示处理器还被配置用于：
[0083] 确定原始图像中的遮挡区域，根据所述遮挡区域确定非遮挡区域；
[0084] 对所述非遮挡区域进行图像分割，得到人体组织前景图像和背景图像。
[0085] 进一步地，所示处理器还被配置用于：
[0086] 将原始图像分解成一个低频图像和多个高频图像；
[0087] 对所述低频图像的动态范围进行压缩；
[0088]对所述高频图像进行强度均衡化处理；
[0089]将压缩后的低频图像以及处理后的高频图像进行重构处理，得到增强图像。
[0090] 进一步地，所示处理器还被配置用于：
[0091]根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应的感兴趣灰 度分布范围；
[0092]根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图像进行压缩。 [0093] 进一步地，所示处理器还被配置用于：
[0094]获取所述尚频图像的最大幅度值；
[0095] 根据所述最大幅度值对所述高频图像进行归一化处理；
[0096] 根据不同的控制参数对所述归一化处理后的高频图像进行频率调制，以扩展小幅 度信号的强度、抑制大幅度信号的强度；
[0097]将处理后的高频图像重新映射到归一化处理前的图像的数据范围。
[0098 ] 进一步地，所示处理器还被配置用于：
[0099]对所述人体组织前景图像和所述增强图像进行预处理；
[0100] 根据原始图像对应的解剖部位类型确定参考灰度值；
[0101] 根据所述前景图像、所述参考灰度值、预设的控制参数确定非线性LUT曲线。
[0102] 进一步地，所示处理器还被配置用于：
[0103] 对所述增强图像应用非线性LUT曲线，得到处理后的增强图像；
[0104] 获得所述处理后的增强图像的感兴趣区域灰度分布特征；
[0105] 判断所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预设条件；
[0106] 若不满足，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线，并对 所述增强图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的增强 图像的感兴趣区域灰度分布特征满足预设条件；
[0107]若满足，则输出所述处理后的增强图像。
[0108] 进一步地，所示处理器还被配置用于：
[0109] 判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比 例值之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足 预设条件。
[0110] 进一步地，所示处理器还被配置用于：
[0111] 若判断调整后的控制参数满足第二预设条件，则不再调整控制参数，输出所述处 理后的增强图像。
[0112] 本发明实施例提供的图像动态范围调整方法及装置，可以根据不同的解剖部位类 型进行不同的动态范围调整，以突出感兴趣区域的灰度分布范围，降低背景图像的干扰，并 能够提供更佳的图像对比度，显示更丰富的信息，达到更佳的显示效果。
【附图说明】
[0113] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下， 还可以根据这些附图获得其他的附图。
[0114] 图1本发明实施例提供的图像动态范围调整方法流程示意图；
[0115] 图2为本发明实施例多分辨率调节获得增强图像的示意图；
[0116] 图3为本发明实施例动态范围调整流程图；
[0117] 图4为本发明一个实施例中腰椎图像调整过程示意图；
[0118] 图5为本发明实施例提供的增强图像的直方图分布示意图；
[0119] 图6为本发明实施例提供的动态范围调整后的直方图分布示意图；
[0120] 图7为本发明实施例提供的输入前景直方图；
[0121 ]图8为本发明实施例提供的动态范围调整后的前景直方图；
[0122] 图9为本发明实施例提供的胸片增强图像示意图；
[0123] 图10为本发明实施例提供的动态范围调整后的胸片图像示意图；
[0124] 图11为本发明实施例提供的调整得到的胸片动态范围LUT曲线示意图；
[0125] 图12为本发明实施例提供的胸片增强图像直方图；
[0126] 图13为本发明实施例提供的动态范围调整后的胸片图像直方图。
[0127] 图14为本发明实施例提供的图像动态范围调整装置示意图；
[0128] 图15为本发明实施例提供的图像动态范围调整装置的实施例硬件架构图。
【具体实施方式】
[0129] 本发明实施例提供了一种图像动态范围调整方法及装置，可以针对不同的解剖部 位进行不同的动态范围调整，可以适应不同图像类型的显示需求，适用性好，具有更佳的显 示效果。
[0130] 为了使本技术领域的人员更好地理解本发明中的技术方案，下面将结合本发明实 施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施 例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通 技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护 的范围。
[0131] 参见图1，本发明实施例提供的图像动态范围调整方法流程示意图。如图1所示，所 述方法可以包括：
[0132] S101，分析原始图像，获得人体组织前景图像。
[0133] 具体地，对输入的原始图像进行分析，首先确定原始图像中的遮挡区域，根据所述 遮挡区域确定非遮挡区域。一般地，在X线图像中存在一定的遮挡区域，一般为限束器遮挡 区域。X线图像中限束器遮挡区域是指通过铅叶遮挡，避免人体受到多余辐射的非X线直接 照射区域，该区域在图像中的主要特征是灰度偏低，一般情况下遮挡区域边缘与其他部分 存在较大的灰度反差，且边缘为直线。本发明实施例具体实现时，利用灰度、边缘直线检测 相结合的方法确定限束器遮挡区域。举例说明，首先对原始图像进行边缘检测，再对边缘检 测图像检测直线(hough变换），结合限束器遮挡区域灰度与其他区域的关系判定是否属于 限束器遮挡区域。当然，本领域技术人员可以理解的是，还可以采用其他方法识别限束器遮 挡区域，在此不进行限定。
[0134] 在确定了限束器遮挡区域后，根据原始图像和束器遮挡区域，即可以确定非遮挡 区域。然后，对非遮挡区域进行图像分割，分别得到人体组织前景图像和背景图像，其中人 体组织前景图像是指人体有效组织区域，而背景是指人体组织区域外的区域。在进行前景 图像分割时，具体采用的方法可以是灵活多样的，例如基于直方图的波峰波谷检测的方法， 基于区域增长的方法，以及基于纹理特征的方法。发明人经过研究发现，在实际应用过程 中，由于图像部位的多样性，单纯的使用某一特征进行检测往往满足不了具体要求。在本发 明实施例具体实现时，对于前景背景反差较大的部位采用基于直方图波峰波谷检测的方法 进行图像分割;而对于前景与背景灰度反差较小的部位采用基于纹理特征与灰度特征结合 的方法进行图像分割。在判断适用何种方法时，具体可以根据前景图像和背景图像的平均 灰度值进行对比，从而确定采用直方图波峰波谷检测方法，还是采用纹理特征与灰度特征 结合的方法。由于后续的动态范围调整并不是把背景直接进行灰度截止，因此对于分割结 果边缘精度并不高，只要能保证前景背景分割不出错即可。因此，本领域技术人员还可以采 用其他方法实现前景图像与背景图像的分割。需要说明的是，灰度截止是指将图像中像素 值大于或小于一定阈值的点，修改成预先的设定值的方法。
[0135] S102,对原始图像进行多分辨率调节，获得增强图像。
[0136] 参见图2,为本发明实施例多分辨率调节获得增强图像的示意图。具体实现时，将 原始图像分解成一个低频图像和一系列高频图像，压缩低频图像动态范围，对高频图像分 别进行强度均衡化处理，再将这一系列处理后的低频和高频图像进行重构得到增强图像。 具体地，本发明实施例中分解、重构图像的过程，使用拉普拉斯金字塔框架，分解层数N可以 根据需要确定。如果不对低频或高频进行处理，重构后图像与原始图像相同。下面的实施例 中，以分解层数N为7层为例进行说明。下面结合附图2,对具体实现进行详细地阐述。
[0137] S102A，将原始图像分解成一个低频图像和多个高频图像。
[0138] 具体实现时，使用拉普拉斯金字塔框架将原始图像分解成一个低频图像和7个高 频图像。
[0139 ] SI 02B，对所述低频图像的动态范围进行压缩。
[0140] 在本发明实施例中，可以对低频图像预先进行预处理，例如降采样。
[0141] 在对低频图像的动态范围进行压缩时，发明人发现，若对于不同的图像，采用同一 的压缩方式，例如采用相同的压缩强度系数进行压缩，效果并不好。例如，若采用的所述压 缩强度系数太小，则会导致压缩幅度不够;若采用的所述压缩强度系数过大，则会导致对比 度被过度加强导致图像失真。因此本发明是实施例中针对不同的解剖部位图像，进行额外 的压缩补偿。具体地，根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应 的感兴趣灰度分布范围;根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图 像进行压缩。举例说明，参见表1，为各个部位类型对应的感兴趣灰度比例范围表。参照表1 确定的解剖部位类型对应的感兴趣灰度比例范围，可以按照如下三种情况进行低频动态范 围压缩：
[0142] 表1
[0143]
[0144] (1)对于感兴趣灰度分布范围分布在整个动态范围的0-40 %的情况：
[0145]
I )
[0146] 其中，Gs'表示动态范围压缩后的低频图像的像素灰度值，GS表示低频图像的像素 灰度值，BaseGray表示基准灰度，kc为压缩强度系数。??Η为对应高灰度部分的额外补偿系 数，取值范围可以[0.5，1]。
[0147] 其中，基准灰度BaseGray可以通过以下方式得到，将步骤101中得到的前景图像按 照均值降采样成与低频图像大小相同的图像，将降采样的前景部分图像与低频图像进行与 运算可得到低频前景图像，统计该低频前景图像的中间灰度值或平均灰度值即可确定该基 准灰度BaseGray。
[0148] 举例说明，参见表1，对于解剖部位类型为头部、腹部、骨盆，其感兴趣灰度分布比 例为30%，可以参照公式（1)进行压缩。当图像中的像素灰度值小于基准灰度BaseGray，则 采用普通的压缩方式；当图像中的像素灰度值大于等于基准灰度BaseGray，则采用额外的 补偿系数??Η进行压缩补偿。
[0149] (2)而对于感兴趣灰度分布在整个动态范围的60-100%的情况：
[0150]
[0151] 其中，Gs'表示动态范围压缩后的低频图像的像素灰度值，GS表示低频图像的像素 灰度值，BaseGray表示基准灰度，kc为压缩强度系数。fEL为对应低灰度部分的额外补偿系 数，取值范围可以为[0.5，1]。其中，基准灰度BaseGray的获得方式与前面介绍的相同。
[0152] 举例说明，参见表1，对于解剖部位类型为胸部，其感兴趣灰度分布比例为80%，可 以参照公式(2)进行压缩。当图像中的像素灰度值小于基准灰度BaseGray，则采用普通的压 缩方式；当图像中的像素灰度值大于等于基准灰度BaseGray，则采用额外的补偿系数fEL进 行压缩补偿。
[0153] (3)对于感兴趣灰度分布在整个动态范围的40-60%的情况：
[0154] Gs' =kcX (Gs-BaseGray)+BaseGray (3)
[0155] 其中，Gs'表示动态范围压缩后的低频图像的像素灰度值，Gs表示低频图像的像素 灰度值，BaseGray表示基准灰度，kc为压缩强度系数。
[0156] 举例说明，参见表1，对于解剖部位类型为手、足、四肢、股部、膝部，由于均是相对 不太需要关注的区域，因此不采用额外的补偿系数进行补偿，直接采用公式(3)进行压缩即 可。
[0157] 在本发明实施例中，针对不同的解剖部位类型，采用不同的方法或者不同的补偿 系数进行低频压缩处理，从而压缩了背景以及非感兴趣区域的灰度范围，突出了感兴趣区 域的灰度范围，使得感兴趣区域灰度分布更均衡。。
[0158] S102C，对所述高频图像进行强度均衡化处理。
[0159]对所述高频图像进行强度均衡化处理，意即对不同频率上的目标高频信号进行调 制，扩展小幅度信号强度，抑制幅度较大的信号的强度。根据解剖部位分布特点，所需的控 制均衡幅度会有所不同。一般而言，0-3层含有较丰富的信息，4-6层结构和层次信息丰富，7 层以上主要包含大范围亮度信息，通过调整各层次的控制参数可以突出显示感兴趣区域的 频率，从而达到更好的显示效果。
[0160] 各层次的强度均衡化处理过程类似，具体可以包括如下步骤：
[0161] S102C1，获取所述高频图像的最大幅度值。
[016 2]其中，尚频图像的最大幅度值为：
[0163] lapmax=max(abs(Lap))
[0164] 意即获取每个高频图像的灰度最大值。
[0165] S102C2,根据所述最大幅度值对所述高频图像进行归一化处理。
[0166] 根据最大幅度值将高频图像归一化至(_1，1)之间。
[0167] S102C3,根据不同的控制参数对所述归一化处理后的高频图像进行频率调制，以 扩展小幅度信号的强度、抑制大幅度信号的强度。
[0168] 具体地，由于拉普拉斯金字塔分解后噪声主要分布在低层次的高频信息中，例如 1-3层、因此频率处理过程对于1-3层高频图像按照公式(4)进行，对于3层-N层按照公式(5) 进行。N的取值范围在本发明实施例可以取值为10,也可以是大于等于7的任意整数。
[0169]
[0171] 其中，c表示归一化的高频信息，c'表示处理后的归一化高频信息，es为控制参数， 控制频率均衡强度es e [0. 5，1 ]，该值越大均衡强度越大，越小均衡强度越弱。一般的，es取 值为0.8。
[0172] 其中，cn为高频信息中噪声与细节信息之间的阈值。该阈值可根据离线采样平均 灰度与噪声的关系，也可以通过均方差统计出现频率确定。本发明实施例中根据前景平均 灰度确定相对噪声，除以lapmax可以计算得到cn值。
[0173] S102C4,将处理后的高频图像重新映射到归一化处理前的图像的数据范围。
[0174] 具体处理公式如下：
[0175] lapmax' = lapmax X c' (6)
[0176] 由于频率调制前进行了归一化处理，因此，这里需要通过重新映射处理，使得高频 图像重新映射到归一化处理前的数据范围，以使得高频图像的灰度范围与原始高频图像的 灰度范围相同。
[0177] S102D，将压缩后的低频图像以及处理后的高频图像进行重构处理，得到增强图 像。
[0178] 具体的重构方法，可以采用拉普拉斯金字塔框架进行重构。
[0179] S103,根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定非线性 LUT(英文全称为Look Up Table，中文名称为查找表）曲线，利用所述非线性LUT曲线调整所 述增强图像的动态范围，获得输出图像。
[0180] 图像动态范围调整目的主要包括两方面一是于压缩背景灰度，二是调整前景图像 灰度分布。对于背景灰度分布，按照一定比例需要压缩其动态范围，而对于前景区域则依据 解剖部位特征，扩展感兴趣灰度范围，压缩其余的灰度部分，并且避免前景中非感兴趣出现 过明或过暗的情况。基于此，本发明实施例根据原始图像对应的解剖部位类型以及所述人 体组织前景图像确定了非线性LUT曲线，并利用所述非线性LUT曲线调整所述增强图像的动 态范围，获得具有更高对比度的输出图像。
[0181] 具体实现时，根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定 非线性LUT曲线可以包括:对所述人体组织前景图像和所述增强图像进行预处理;根据原始 图像对应的解剖部位类型确定参考灰度值;根据所述前景图像、所述参考灰度值、预设的控 制参数确定非线性LUT曲线。
[0182] 所述利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出图像包括： 对所述增强图像应用非线性LUT曲线，得到处理后的增强图像;获得所述处理后的增强图像 的感兴趣区域灰度分布特征;判断所述处理后的增强图像的感兴趣区域灰度分布特征是否 满足预设条件;若不满足，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线， 并对所述增强图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的 增强图像的感兴趣区域灰度分布特征满足预设条件;若满足，则输出所述处理后的增强图 像。
[0183] 下面结合附图3,对本发明实施例中动态范围调整流程图进行说明。
[0184] S301，对输入信息进行预处理，结合解剖部位类型确定参考灰度。
[0185] 对输入信息的图像预处理过程可以包含增强图像和前景图像的降采样、干扰信息 排除等处理。首先对增强图像S和前景图像Μ降采样。降采样过程为了改善产生LUT效率，而 最终LUT应用在整幅图像上。由于DR图像分辨率较大，为了提升处理效率需要将图像按照一 定采样率抽样，优选地，一个可行的实现方式中，可以设置降采样比率为8:1，得到降采样图 像增强图像S1，降采样前景图Ml。需要注意的是接下来的曲线形成过程都是基于降采样的 S1和Ml进行的，但最终生成处理后的数据是在原始大尺寸图像上应用LUT得到的。
[0186] S302,根据初始参数构造非线性曲线。
[0187] 按照初始设定形成灰度映射曲线，对增强图像S1应用灰度映射曲线，得到处理图 像SP1，LUT曲线可以通过以下方法得到：
[0188] 定义如下几个控制点pl(x0，y0)，p2(xl，yl)，p3(x2，y2)，p4(x3，y3)，通过非线性 曲线拟合可以得到映射曲线，拟合方式可选分段多项式方法，注意x〇，xl，x2，x3和yO，y 1， y 2，y 3都呈升序分布，避免形成异常曲线。
[0189] χ〇 = 〇
[0190] y〇 = 〇
[0191] xl =bcg_gray
[0192] yl=klX(xl-xO)+yO
[0193] x2 = ref_gray
[0194] y2 = k2X(x2-xl)+yl
[0195] x3=max_gray
[0196] y3 = k3X(x3-x2)+y2
[0197] 其中，kl，k2，k3分别表示p2，p3，p4的控制参数，kl，k2，k3取值范围分别是[0_ 0·2]，[0·1-2]，[0·5-2]。
[0198] 其中，bcg_gray表示背景与前景区域分割阈值，增强后的图像前景和背景灰度分 布会存在一定交叉，可以选择该交叉点所对应的灰度值作为该分割阈值，也可以选择直接 在前景图像中，按照一定比例计算得到。其中ref_gray表示参考灰度，参考灰度需要根据具 体图像来确定的，根据大量数据统计，相同的摄影部位和体位在灰度分布存在共同的特点。 因此本发明实施例根据部位类型和前景图像，按照预定比例计算到该灰度值，表2给出了几 种典型部位参考灰度在前景直方图的比例分布。
[0199] 表2不同部位对应的特征灰度分布
[0200]
[0201]
[0202]其中所述按照预定比例计算灰度值，具体的做法是:先计算前景图像直方图，统计 出前景图像中像素灰度总数，乘以设定好的比例值得到累加器门限值，初始化累加器，遍历 直方图所有通道，依次将直方图中各个通道对应的统计值加入累加器，当超过累加器门限 值时，该通道对应的灰度值即是按照一定比例得到的灰度值。
[0203] S303,对增强图像应用非线性曲线得到处理后图像。
[0204]具体地，对于增强图像的每一个像素点p(x，y)，都有p'（x，y)=LUT(p(x，y))，得到 的P'就是应用了 LUT的处理结果。
[0205] S304，统计图像感兴趣区域灰度分布特征。
[0206]具体实现时，所述感兴趣区域灰度分布特征具体为图像动态范围比例值。
[0207]举例说明，步骤304统计处理后图像中感兴趣区域灰度分布特征，主要包括计算处 理后图像直方图中参考点左侧动态范围和右侧动态范围的比例，可以按照如下公式得到：
[0208] (6)
[0209] 其中，dynamic_ratio为图像动态范围比例值，data_pro_ref，data_pro_left， da ta_pr o_r i gh t可以按照步骤3 0 2中给出的按照预定比例计算灰度值的方法确定，为了排 除可能的异常像素的干扰，可以为data_pro_left，data_pro_right设定预定比例值。举例 说明，其中data_pro_left，data_pro_right对应的预定比例分别为1 %，99%，data_pro_ left，data_pro_right对应的预定比例也可以分别是0.05%，99.5%，或者分别是2%、 98 %。需要说明的是，所述预定比例值为根据经验选取，一般地，data_pro_lef t对应的预定 比例一般不大于2%，data_pro_right对应的预定比例一般不小于98%，否则可能会导致图 像失真，例如会损失一些图像细节。本领域技术人员可以根据经验选取所述预定比例值，在 此不进行限定。而data_pr 〇_ref对应的预设比例值则取决于输入的部位类型信息，可以根 据表2得到。
[0210] S305,检查感兴趣区域灰度分布特征是否满足预设条件。如果是，则执行S307,得 到最终处理后的图像。否则执行S306,调整参数构造新的非线性曲线，重复步骤303-305,直 到感兴趣区域灰度分布特征满足预设条件。
[0211] 这一步骤是调整LUT曲线的过程。具体地，判断处理后图像中感兴趣区域灰度分布 特征是否满足需求，如果达到则可以得到最终处理后图像，否则需要调整控制参数后重新 形成LUT曲线。
[0212] 所述确定所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预设条件 包括:判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比例值 之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足预设 条件。
[0213] 举例说明，判断调整后的图像动态范围特征比例值dynamic_ratio和目标动态范 围特征比例值dynamic_target偏差的绝对值是否小于预设门限值threshold，本发明实施 例中threshold取0.1，dynamic_target固定取值为1，如果小于该门限值，表明感兴趣灰度 范围满足需求，否则需要调整参数后重新形成新的LUT曲线。
[0214] S306,调整参数构造新的非线性曲线，重复步骤303-305,直到感兴趣区域灰度分 布特征满足预设条件。
[0215] 当确定不满足预设条件时，则需要调整控制参数。优选地，只需要调整控制参数 k2，kl和k3直接可以使用预设值，不需要自动调整。
[0216] 其中，k2的调整方式如下：
[0217] 如果dynamic_ratio_dynamic_target> threshold，则k2 = k2_k_step，
[0218] 如果dynamic_ratio_dynamic_target<-threshold，则k2 = k2+k_step
[0219] 举例说明，若调整后的图像动态范围特征比例值dynamic_ratio和目标动态范围 特征比例值dynamic_target的差值大于预设门限值threshold，则使用控制参数减去调整 步长k_step，得到新的控制参数。若调整后的图像动态范围特征比例值dynamic_ratio和目 标动态范围特征比例值dynamic_target的差值小于-threshold，则使用控制参数加上调整 步长k_step，得到新的控制参数。
[0220]优选地，若调整后的控制参数满足第二预设条件，则不再调整控制参数，输出所述 处理后的增强图像。具体地，为了避免无意义的迭代，当出现以下情况时，不需要再调整k2 值。第二预设条件具体为：
[0221] (l)k2已超出有效范围
[0222] k2<0.1 或 k2>2
[0223] 需要说明的是，k2的有效取值范围为根据经验选取。经过发明人反复研究发现，k2 在[0.1，2]的取值范围时，图像的调整效果较好。若k2超出有效范围，调整后的图像则会出 现一定程度的失真现象。因此，为了避免图像失真，需要对k2的有效范围进行判断，以达到 更好的图像处理效果。< br>[0224] (2)出现震荡情况。也就是说，迭代前后出现与目标动态范围偏差异常，即同时满 足如下条件：
[0225] dy nam i c_ra t i o (n) -dy nam i c_tar ge t < 0
[0226] dy nam i c_ra t i o (n+1) -dy nam i c_tar ge t < 0
[0227] abs(dynamic_ratio(n)-dynamic_target)>0.1
[0228] abs(dynamic_ratio(n+l)-dynamic_target)>0.1
[0229] 其中，dynamic_ratio(n)表示第n次调整k2值后的计算结果；
[0230] dynamic_ratio(n+l)表示第n+1次调整k2后的计算结果。
[0231] S307,得到最终处理后的输出图像。
[0232] 最后将调整后的LUT应用到增强的全尺寸图像，就可以得到最终动态范围调整后 的图像。
[0233] 下面结合附图4-13对本发明实施例能够达到的有益效果进行说明。
[0234] 参见图4,为本发明一个实施例中腰椎图像调整过程示意图，左侧是多分辨率调整 后的增强图像，右侧是调c整好的输出图像，该调整过程难点就在于如何确定合理的LUT曲 线。
[0235] 当图像输出到显示设备上时，根据亮度对比度设置将灰度映射成亮度信息，因此 能够被很好显示的仅是窗宽范围内的目标，因此感兴趣的灰度能尽量均匀分布在整个动态 范围中间，才会更便于观察病灶。图5是增强图像直方图分布，图中给出了背景、皮肤、肺部、 肌肉、低密度、高密度骨骼等一系列组织灰度分布，背景，肺，皮肤组织占据了很宽的灰度范 围，高密度骨骼组织灰度分布过度则集中，动态范围很窄。而图6是经过动态范围调整后的 直方图分布，背景、肺、皮肤组织灰度都被压缩至很小的范围，高密度和低密度骨骼肌肉灰 度分布范围都不同程度被扩展。因此可以看到图6中，按照前景灰度最小最大范围设置窗宽 窗位进行灰度映射的显示效果中，变换前的图像发白，感兴趣区域对比度不好，而变换后的 图像感兴趣区域对比度合适，更适合观察病灶。图7和图8分别是变换前后前景区域，1024级 直方图分布，二者比较可以看出变换对于灰度分布集中程度的改善。
[0236] 图9-13给出了一种胸片图像的应用实例，图9给出了胸片增强图像，图10是经过动 态范围调整后的胸片图像，图11是经过调整得到的灰度映射曲线，图12和13分别是图像9和 图像10对应的直方图分布，从结果中也可以看到曲线调整后图像对比度改善，从直方图变 化中也能看到肺部灰度动态范围被扩展，更有利于肺部组织观察。
[0237] 在本发明实施例中，能够达到如下有益效果：
[0238] (1)对处理结果进行分析判定，调整曲线生成，更利于达到预期的效果。
[0239] (2)综合考虑解剖部位特征、避免通用基于直方图盲调方法造成的对比过强或灰 度分布不连续等问题，抑制增强带来的尾影。
[0240] (3)能突出显示感兴趣区域组织灰度，又保证其他前景组织被观察到，增加图像可 观测的信息量。
[0241] 本发明实施例经过大量图像样本验证，都可以获得比较满意的图像动态范围。
[0242] 与前述图像动态范围调整方法的实施例相对应，本发明还提供了图像动态范围调 整装置的实施例。
[0243] 参照图14,为本发明实施例提供的图像动态范围调整装置框图。
[0244] 所示图像动态范围调整装置1400可以包括：
[0245] 前景图像获得单元1401，用于分析原始图像，获得人体组织前景图像。
[0246] 增强图像获得单元1402,用于对原始图像进行多分辨率调节，获得增强图像。
[0247] 动态范围调整单元1403,用于根据原始图像对应的解剖部位类型以及所述人体组 织前景图像确定非线性LUT曲线，利用所述非线性LUT曲线调整所述增强图像的动态范围， 获得输出图像。
[0248] 进一步地，所述前景图像获得单元包括：
[0249] 确定单元，用于确定原始图像中的遮挡区域，根据所述遮挡区域确定非遮挡区域；
[0250] 分割单元，用于对所述非遮挡区域进行图像分割，得到人体组织前景图像和背景 图像。
[0251] 进一步地，所述增强图像获得单元包括：
[0252] 分解单元，用于将原始图像分解成一个低频图像和多个高频图像；
[0253] 低频压缩单元，用于对所述低频图像的动态范围进行压缩；
[0254]高频处理单元，用于对所述高频图像进行强度均衡化处理；
[0255] 重构单元，用于将压缩后的低频图像以及处理后的高频图像进行重构处理，得到 增强图像。
[0256] 进一步地，所述低频压缩单元具体用于：
[0257] 根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应的感兴趣灰 度分布范围；根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图像进行压 缩。
[0258] 进一步地，所述高频处理单元包括：
[0259] 最大幅度值获得单元，用于获取所述高频图像的最大幅度值；
[0260] 归一化单元，用于根据所述最大幅度值对所述高频图像进行归一化处理；
[0261] 调制单元，用于根据不同的控制参数对所述归一化处理后的高频图像进行频率调 制，以扩展小幅度信号的强度、抑制大幅度信号的强度；
[0262]映射单元，用于将处理后的高频图像重新映射到归一化处理前图像的数据范围。 [0263] 进一步地，所述动态范围调整单元包括曲线确定单元以及调整单元，所述曲线确 定单元包括：
[0264]预处理单元，用于对所述人体组织前景图像和所述增强图像进行预处理；
[0265] 参考灰度值确定单元，用于根据原始图像对应的解剖部位类型确定参考灰度值；
[0266] 曲线生成单元，用于根据所述前景图像、所述参考灰度值、预设的控制参数确定非 线性LUT曲线。
[0267] 进一步地，所述动态范围调整单元包括曲线确定单元以及调整单元，所述调整单 元包括：
[0268] 应用单元，用于对所述增强图像应用非线性LUT曲线，得到处理后的增强图像；
[0269] 特征获得单元，用于获得所述处理后的增强图像的感兴趣区域灰度分布特征；
[0270] 判断单元，用于判断所述处理后的增强图像的感兴趣区域灰度分布特征是否满足 预设条件；
[0271] 迭代调整单元，用于若判断所述处理后的增强图像的感兴趣区域灰度分布特征不 满足预设条件，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线，并对所述 增强图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的增强图像 的感兴趣区域灰度分布特征满足预设条件；
[0272] 第一输出单元，用于若判断所述处理后的增强图像的感兴趣区域灰度分布特征满 足预设条件，输出所述处理后的增强图像。
[0273] 进一步地，所述特征获得单元获得的所述感兴趣区域灰度分布特征具体为图像动 态范围比例值，所述判断单元具体用于：
[0274] 判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比 例值之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足 预设条件。
[0275] 进一步地，所述装置还包括：
[0276] 第二输出单元，用于若判断调整后的控制参数满足第二预设条件，则不再调整控 制参数，输出所述处理后的增强图像
[0277] 其中，本发明装置各单元或模块的设置可以参照图1所示的方法而实现，在此不赘 述。
[0278] 本发明图像动态范围调整装置的实施例可以应用在DR系统的设备或者任何具有 处理器(特别是图像处理器）的电子设备上，所述电子设备可以是现有的、正在研发的或将 来研发的任何电子设备，包括但不限于:现有的、正在研发的或将来研发的台式计算机、膝 上型计算机、移动终端(包括智能手机、非智能手机、各种平板电脑)等。装置实施例可以通 过软件实现，也可以通过硬件或者软硬件结合的方式实现。以软件实现为例，作为一个逻辑 意义上的装置，是通过其所在DR系统的设备或带有处理器的电子设备的处理器将存储器中 对应的计算机程序指令读取到内存中运行形成的。从硬件层面而言，如图15所示，为本发明 图像动态范围调整装置所在DR系统的设备或带有处理器的电子设备的一种硬件结构图，除 了图15所示的处理器、内存、网络接口、以及存储器之外，实施例中装置所在的DR系统的设 备或带有处理器的电子设备通常根据该设备的实际功能，还可以包括其他硬件，对此不再 赘述。
[0279] 其中，存储器中可以存储有图像动态范围调整方法对应的逻辑指令，该存储器例 如可以是非易失性存储器(non-volatile memory)。处理器可以调用执行存储器中的保存 的逻辑指令，以执行上述的图像动态范围调整方法。
[0280] 图像动态范围调整方法对应的逻辑指令的功能，如果以软件功能单元的形式实现 并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的 理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分 可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指 令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各 个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器 (ROM，Read_Only Memory)、随机存取存储器（RAM，Random Access Memory)、磁碟或者光盘 等各种可以存储程序代码的介质。
[0281] 如图15所示，本发明还提供了 一种图像动态范围调整设备，包括：
[0282] 处理器；
[0283] 用于存储所述处理器可执行指令的存储器；
[0284] 其中，所述处理器被配置为：
[0285] 分析原始图像，获得人体组织前景图像；
[0286] 对原始图像进行多分辨率调节，获得增强图像；
[0287] 根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定非线性LUT曲 线，利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出图像。
[0288] 进一步地，所示处理器还被配置用于：
[0289] 确定原始图像中的遮挡区域，根据所述遮挡区域确定非遮挡区域；
[0290] 对所述非遮挡区域进行图像分割，得到人体组织前景图像和背景图像。
[0291] 进一步地，所示处理器还被配置用于：
[0292] 将原始图像分解成一个低频图像和多个高频图像；
[0293] 对所述低频图像的动态范围进行压缩；
[0294] 对所述高频图像进行强度均衡化处理；
[0295] 将压缩后的低频图像以及处理后的高频图像进行重构处理，得到增强图像。
[0296] 进一步地，所示处理器还被配置用于：
[0297] 根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应的感兴趣灰 度分布范围；
[0298] 根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图像进行压缩。
[0299] 进一步地，所示处理器还被配置用于：
[0300]获取所述尚频图像的最大幅度值；
[0301 ]根据所述最大幅度值对所述高频图像进行归一化处理；
[0302]根据不同的控制参数对所述归一化处理后的高频图像进行频率调制，以扩展小幅 度信号的强度、抑制大幅度信号的强度；
[0303]将处理后的高频图像重新映射到归一化处理前的图像的数据范 围。
[0304] 进一步地，所示处理器还被配置用于：
[0305] 对所述人体组织前景图像和所述增强图像进行预处理；
[0306] 根据原始图像对应的解剖部位类型确定参考灰度值；
[0307]根据所述前景图像、所述参考灰度值、预设的控制参数确定非线性LUT曲线。
[0308] 进一步地，所示处理器还被配置用于：
[0309] 对所述增强图像应用非线性LUT曲线，得到处理后的增强图像；
[0310] 获得所述处理后的增强图像的感兴趣区域灰度分布特征；
[0311] 判断所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预设条件；
[0312] 若不满足，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线，并对 所述增强图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的增强 图像的感兴趣区域灰度分布特征满足预设条件；
[0313] 若满足，则输出所述处理后的增强图像。
[0314] 进一步地，所示处理器还被配置用于：
[0315] 判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比 例值之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足 预设条件。
[0316] 进一步地，所示处理器还被配置用于：
[0317] 若判断调整后的控制参数满足第二预设条件，则不再调整控制参数，输出所述处 理后的增强图像。
[0318] 本领域技术人员可以理解的是，以上对方法和装置实施例进行了示例性说明，以 上不视为对本发明的限制，本领域技术人员在不付出创造性劳动下获得的其他实现方式均 属于本发明的保护范围。
[0319]需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实 体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存 在任何这种实际的关系或者顺序。而且，术语"包括"、"包含"或者其任何其他变体意在涵盖 非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要 素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备 所固有的要素。在没有更多限制的情况下，由语句"包括一个……"限定的要素，并不排除在 包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。本发明可以在由计算 机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执 行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分 布式计算环境中实践本发明，在这些分布式计算环境中，由通过通信网络而被连接的远程 处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地 和远程计算机存储介质中。
[0320]本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部 分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于装置实 施例而言，由于其基本相似于方法实施例，所以描述得比较简单，相关之处参见方法实施例 的部分说明即可。以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明 的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是 物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要 选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出 创造性劳动的情况下，即可以理解并实施。以上所述仅是本发明的【具体实施方式】，应当指 出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干 改进和润饰，这些改进和润饰也应视为本发明的保护范围。
【主权项】
1. 一种图像动态范围调整方法，其特征在于，所述方法包括： 分析原始图像，获得人体组织前景图像； 对原始图像进行多分辨率调节，获得增强图像； 根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定非线性LUT曲线， 利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出图像。2. 根据权利要求1所述的方法，其特征在于，所述分析原始图像，获得人体组织前景图 像包括： 确定原始图像中的遮挡区域，根据所述遮挡区域确定非遮挡区域； 对所述非遮挡区域进行图像分割，得到人体组织前景图像和背景图像。3. 根据权利要求1所述的方法，其特征在于，所述对原始图像进行多分辨率调节，获得 增强图像包括： 将原始图像分解成一个低频图像和多个高频图像； 对所述低频图像的动态范围进行压缩； 对所述高频图像进行强度均衡化处理； 将压缩后的低频图像以及处理后的高频图像进行重构处理，得到增强图像。4. 根据权利要求3所述的方法，其特征在于，所述对所述低频图像的动态范围进行压缩 包括： 根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应的感兴趣灰度分 布范围； 根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图像进行压缩。5. 根据权利要求3所述的方法，其特征在于，所述对所述高频图像进行强度均衡化处理 包括： 获取所述高频图像的最大幅度值； 根据所述最大幅度值对所述高频图像进行归一化处理； 根据不同的控制参数对所述归一化处理后的高频图像进行频率调制，以扩展小幅度信 号的强度、抑制大幅度信号的强度； 将处理后的高频图像重新映射到归一化处理前的图像的数据范围。6. 根据权利要求1所述的方法，其特征在于，所述根据原始图像对应的解剖部位类型以 及所述人体组织前景图像确定非线性LUT曲线包括： 对所述人体组织前景图像和所述增强图像进行预处理； 根据原始图像对应的解剖部位类型确定参考灰度值； 根据所述前景图像、所述参考灰度值、预设的控制参数确定非线性LUT曲线。7. 根据权利要求1或6所述的方法，其特征在于，所述利用所述非线性LUT曲线调整所述 增强图像的动态范围，获得输出图像包括： 对所述增强图像应用非线性LUT曲线，得到处理后的增强图像； 获得所述处理后的增强图像的感兴趣区域灰度分布特征； 判断所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预设条件； 若不满足，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线，并对所述 增强图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的增强图像 的感兴趣区域灰度分布特征满足预设条件； 若满足，则输出所述处理后的增强图像。8. 根据权利要求7所述的方法，其特征在于，所述感兴趣区域灰度分布特征具体为图像 动态范围比例值，所述确定所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预 设条件包括： 判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比例值 之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足预设 条件。9. 根据权利要求7所述的方法，其特征在于，所述方法还包括： 若判断调整后的控制参数满足第二预设条件，则不再调整控制参数，输出所述处理后 的增强图像。10. -种图像动态范围调整装置，其特征在于，所述装置包括： 前景图像获得单元，用于分析原始图像，获得人体组织前景图像； 增强图像获得单元，用于对原始图像进行多分辨率调节，获得增强图像； 动态范围调整单元，用于根据原始图像对应的解剖部位类型以及所述人体组织前景图 像确定非线性LUT曲线，利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出 图像。11. 根据权利要求10所述的装置，其特征在于，所述前景图像获得单元包括： 确定单元，用于确定原始图像中的遮挡区域，根据所述遮挡区域确定非遮挡区域； 分割单元，用于对所述非遮挡区域进行图像分割，得到人体组织前景图像和背景图像。12. 根据权利要求10所述的装置，其特征在于，所述增强图像获得单元包括： 分解单元，用于将原始图像分解成一个低频图像和多个高频图像； 低频压缩单元，用于对所述低频图像的动态范围进行压缩； 高频处理单元，用于对所述高频图像进行强度均衡化处理； 重构单元，用于将压缩后的低频图像以及处理后的高频图像进行重构处理，得到增强 图像。13. 根据权利要求12所述的装置，其特征在于，所述低频压缩单元具体用于： 根据所述图像对应的解剖部位的类型，确定与所述解剖部位类型对应的感兴趣灰度分 布范围;根据所述感兴趣灰度分布范围采用对应的补偿系数对所述低频图像进行压缩。14. 根据权利要求12所述的装置，其特征在于，所述高频处理单元包括： 最大幅度值获得单元，用于获取所述高频图像的最大幅度值； 归一化单元，用于根据所述最大幅度值对所述高频图像进行归一化处理； 调制单元，用于根据不同的控制参数对所述归一化处理后的高频图像进行频率调制， 以扩展小幅度信号的强度、抑制大幅度信号的强度； 映射单元，用于将处理后的高频图像重新映射到归一化处理前的图像的数据范围。15. 根据权利要求10所述的装置，其特征在于，所述动态范围调整单元包括曲线确定单 元以及调整单元，所述曲线确定单元包括： 预处理单元，用于对所述人体组织前景图像和所述增强图像进行预处理； 参考灰度值确定单元，用于根据原始图像对应的解剖部位类型确定参考灰度值； 曲线生成单元，用于根据所述前景图像、所述参考灰度值、预设的控制参数确定非线性LUT曲线。16. 根据权利要求10或15所述的装置，其特征在于，所述动态范围调整单元包括曲线确 定单元以及调整单元，所述调整单元包括： 应用单元，用于对所述增强图像应用非线性LUT曲线，得到处理后的增强图像； 特征获得单元，用于获得所述处理后的增强图像的感兴趣区域灰度分布特征； 判断单元，用于判断所述处理后的增强图像的感兴趣区域灰度分布特征是否满足预设 条件； 迭代调整单元，用于若判断所述处理后的增强图像的感兴趣区域灰度分布特征不满足 预设条件，调整控制参数，根据调整后的控制参数重新生成非线性LUT曲线，并对所述增强 图像应用重新生成的非线性LUT曲线获得处理后的增强图像，直到处理后的增强图像的感 兴趣区域灰度分布特征满足预设条件； 第一输出单元，用于若判断所述处理后的增强图像的感兴趣区域灰度分布特征满足预 设条件，输出所述处理后的增强图像。17. 根据权利要求16所述的装置，其特征在于，所述特征获得单元获得的所述感兴趣区 域灰度分布特征具体为图像动态范围比例值，所述判断单元具体用于： 判断所述处理后的增强图像的图像动态范围比例值与预设目标图像动态范围比例值 之间的差值是否小于预设门限值;若小于，确定满足预设条件;若不小于，确定不满足预设 条件。18. 根据权利要求16所述的装置，其特征在于，所述装置还包括： 第二输出单元，用于若判断调整后的控制参数满足第二预设条件，则不再调整控制参 数，输出所述处理后的增强图像。19. 一种图像动态范围调整设备，其特征在于，包括： 处理器； 用于存储所述处理器可执行指令的存储器； 其中，所述处理器被配置为： 分析原始图像，获得人体组织前景图像； 对原始图像进行多分辨率调节，获得增强图像； 根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定非线性LUT曲线， 利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出图像。
【专利摘要】本发明实施例提供一种图像动态范围调整方法和装置，所述方法包括：分析原始图像，获得人体组织前景图像；对原始图像进行多分辨率调节，获得增强图像；根据原始图像对应的解剖部位类型以及所述人体组织前景图像确定非线性LUT曲线，利用所述非线性LUT曲线调整所述增强图像的动态范围，获得输出图像。本发明实施例可以根据不同的解剖部位类型进行不同的动态范围调整，以突出感兴趣区域的灰度分布范围，降低背景图像的干扰，并能够提供更佳的图像对比度，显示更丰富的信息，达到更佳的显示效果。
【IPC分类】G06T5/00
【公开号】CN105488765
【申请号】CN201510818564
【发明人】李海春
【申请人】沈阳东软医疗系统有限公司
【公开日】2016年4月13日
【申请日】2015年11月23日
最新回复(0)