直方图统计增强方法

直方图统计增强方法的原理：

直方图统计增强是一种基于像素值分布的图像增强技术，通过调整像素值的分布来增强图像的对比度和细节。其原理是根据图像的直方图信息，将原始像素值映射到一个新的像素值域，从而改变图像的亮度和对比度。
直方图统计增强方法的意义：

这种方法的目的是改善图像的视觉效果，使得图像中的细节更加清晰可见，从而提高图像在视觉上的识别和分析能力。通过将图像的像素值重新分布，使得图像的动态范围更好地覆盖了整个像素值域，凸显了图像中的细节，增强了图像的信息内容。
适用场景：

直方图统计增强方法在以下场景中特别适用：

低对比度图像： 当图像的亮度范围分布不均匀，即图像的对比度较低时，直方图统计增强可以使得图像的像素值分布更加均匀，增强图像的对比度。
局部对比度调整： 有时候图像的一部分区域可能亮度较低，另一部分区域亮度较高。直方图统计增强可以通过针对不同区域使用不同的增强参数，实现局部对比度的调整。
图像细节增强： 当图像中的细节由于亮度过低或过高而不易分辨时，直方图统计增强可以突出这些细节，使得它们更加清晰可见。

公式：

直方图统计增强的基本公式如下：

给定一个原始图像 I ( x , y ) I(x, y) I(x,y)，其像素值为 i i i，进行直方图统计增强后的图像为 I enhanced ( x , y ) I_{\text{enhanced}}(x, y) Ienhanced(x,y)，其像素值为 i ′ i' i′。增强函数可以表示为：
i ′ = T ( i ) = round ( c ⋅ i ) i' = T(i) = \text{round}(c \cdot i) i′=T(i)=round(c⋅i)

其中， c c c是增强常数，可以通过调整来控制增强的强度。一般情况下， c > 1 c > 1 c>1 用于增强亮度， c < 1 c < 1 c<1用于降低亮度。函数 ( \text{round} ) 将结果四舍五入到最接近的整数，确保结果在合法的像素值范围内。

另外，为了避免像素值越界，通常还会使用截断操作将增强后的像素值 ( i' ) 限制在图像像素值域的合法范围内（通常是 $0, 255$ ）：
i ′ = { 0 if c ⋅ i < 0 255 if c ⋅ i > 255 round ( c ⋅ i ) otherwise i' = \begin{cases} 0 & \text{if } c \cdot i < 0 \\ 255 & \text{if } c \cdot i > 255 \\ \text{round}(c \cdot i) & \text{otherwise} \end{cases} i′=⎩ ⎨ ⎧0255round(c⋅i)if c⋅i<0if c⋅i>255otherwise

这样，直方图统计增强就可以通过调整增强常数 c c c 来改变图像的亮度和对比度，从而实现图像的增强，代码如下所示：

bash 复制代码

import cv2
import matplotlib.pyplot as plt
import numpy as np


def enhanceHistStat(image, const):
    enhanced_image = np.clip(image * const, 0, 255).astype(np.uint8)
    return enhanced_image

def enhance_Histstat(image):
    # 读取图像
    image = cv2.imread(image, flags=0)

    if image is None:
        print('Unable to load input_image!')
        return
    else:

        # 计算比例因子
        maxImg = image.max()
        maxROI = image[12:120, 12:120].max()
        const = maxImg / maxROI

        # 调用增强函数
        image_HSE = enhanceHistStat(image, const)
        image_Equ=cv2.equalizeHist(image)

        plt.figure(figsize=(10,6))

        plt.subplot(131),plt.title("Original image"),plt.axis('off')
        plt.imshow(image,cmap='gray',vmin=0,vmax=255)

        plt.subplot(132),plt.title("Global equalize histogram"),plt.axis('off')
        plt.imshow(image_Equ,cmap='gray',vmin=0,vmax=255)

        plt.subplot(133),plt.title("Histogram statistic enhance"),plt.axis('off')
        plt.imshow(image_HSE,cmap='gray',vmin=0,vmax=255)

        plt.show()

imgfile1='./Images/cat.jpg'
enhance_Histstat(imgfile1)