数字图像处理(18): 图像灰度变换——线性灰度变换和非线性灰度变换(对数变换与伽马变换)

迷南。 2022-02-05 04:19 1818阅读 0赞

![20190702121722457.png_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly94aW9uZ3lpbWluZy5ibG9nLmNzZG4ubmV0_size_16_color_FFFFFF_t_70][]

**目录**

1 灰度变换简介

2 线性灰度变换—图像反转

3 非线性灰度变换

3.1 对数变换

3.2 伽马变换

参考资料

--------------------

# 1 灰度变换简介 #

灰度变换是图像增强的一种重要手段，用于改善图像显示效果，属于空间域处理方法，它可以使图像动态范围加大，使图像对比度扩展，图像更加清晰，特征更加明显。灰度变换其实质就是按一定的规则修改图像每一个像素的灰度，从而改变图像的灰度范围。常见的灰度变换图像反转，对数变换和伽马变换等。其具体分类如下图所示：

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70][]

--------------------

# 2 线性灰度变换—图像反转 #

灰度线性变换最常见的就是图像反转，在灰度图像灰度级范围![\[0,L-1\]][0_L-1]中，其反转的公式如下所示：

![s=L-1-r][s_L-1-r]

其中，![r][]表示原始图像的灰度级，![s][]表示变换后的灰度级。

下图所示为图像反转的例子，原图像是数字乳房X射线照片，其中显示有一小块病变，通过图像反转就很容易看到病变区域。

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 1][]

代码如下所示：

# -*- coding: utf-8 -*-
    import cv2
    import numpy as np
    import matplotlib.pyplot as plt
    
    #读取原始图像
    img = cv2.imread('zxp.jpg')
    
    #图像灰度转换
    grayImage = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY)
    
    #获取图像高度和宽度
    height = grayImage.shape[0]
    width = grayImage.shape[1]
    
    #创建一幅图像
    result = np.zeros((height, width), np.uint8)
    
    #图像灰度反色变换 s=255-r
    for i in range(height):
        for j in range(width):
            gray = 255 - grayImage[i,j]
            result[i,j] = np.uint8(gray)
    
    #显示图像
    
    cv2.imshow("Gray Image", grayImage)
    cv2.imshow("Result", result)
    
    
    #等待显示
    cv2.waitKey(0)
    cv2.destroyAllWindows()

运行结果如下图所示：

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 2][]

--------------------

# 3 非线性灰度变换 #

## 3.1 对数变换 ##

图像灰度**对数变换**一般表示如下所示：

![s=c\\log (1+r)][s_c_log _1_r]

> 其中，![r][]表示原始图像的灰度级，![s][]表示变换后的灰度级，![c][]为常数。
> 
> 假设![r\\ge 0][r_ge 0]，下图所示的对数曲线的形状表明，改变换将输入中范围较窄的低灰度值映射为输出中较宽范围的灰度值。相反的，对高的输入灰度值也是如此。我们使用这种类型的变换来扩展图像中暗像素的值，同时压缩更高灰度级的值。反对数变换的作用与此相反。
> 
> ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 3][]

代码如下所示：

# -*- coding: utf-8 -*-
    import numpy as np
    import matplotlib.pyplot as plt
    import cv2
    
    #绘制曲线
    def log_plot(c):
        x = np.arange(0, 256, 0.01)
        y = c * np.log(1 + x)
        plt.plot(x, y, 'r', linewidth=1)
        plt.rcParams['font.sans-serif']=['SimHei'] #正常显示中文标签
        plt.title(u'对数变换函数')
        plt.xlim(0, 255), plt.ylim(0, 255)
        plt.show()
    
    #对数变换
    def log(c, img):
        output = c * np.log(1.0 + img)
        output = np.uint8(output + 0.5)
        return output
    
    #读取原始图像
    img = cv2.imread('test8.bmp')
    
    #绘制对数变换曲线
    log_plot(42)
    
    #图像灰度对数变换
    output = log(42, img)
    
    #显示图像
    cv2.imshow('Input', img)
    cv2.imshow('Output', output)
    cv2.waitKey(0)
    cv2.destroyAllWindows()

运行结果如下图所示：

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 4][]

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 5][]

--------------------

## 3.2 伽马变换 ##

**伽玛变换**又称为 **指数变换** 或 **幂次变换**，是另一种常用的灰度非线性变换。图像灰度的伽玛变换一般表示如下所示：

![s=c\{\{r\}^\{\\gamma \}\}][s_c_r_gamma]

> 其中，![r][]表示原始图像的灰度级，![s][]表示变换后的灰度级，![c][]和![\\gamma][gamma]正常数
> 
> 1）当![\\gamma][gamma]>1时，会拉伸图像中灰度级较高的区域，压缩灰度级较低的部分；
> 
> 2）当![\\gamma][gamma]<1时，会拉伸图像中灰度级较低的区域，压缩灰度级较高的部分；
> 
> 3）当![\\gamma][gamma]=1时，该灰度变换是线性的，此时通过线性方式改变原图像。
> 
> 如下图所示，不同![\\gamma][gamma]值的变换曲线：
> 
> ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 6][]

下图所示为图像伽马变换的例子：

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 7][]

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 8][]

代码如下所示：

# -*- coding: utf-8 -*-
    import numpy as np
    import matplotlib.pyplot as plt
    import cv2
    
    #绘制曲线
    def gamma_plot(c, v):
        x = np.arange(0, 256, 0.01)
        y = c*x**v
        plt.plot(x, y, 'r', linewidth=1)
        plt.rcParams['font.sans-serif']=['SimHei'] #正常显示中文标签
        plt.title(u'伽马变换函数')
        plt.xlim([0, 255]), plt.ylim([0, 255])
        plt.show()
    
    #伽玛变换
    def gamma(img, c, v):
        lut = np.zeros(256, dtype=np.float32)
        for i in range(256):
            lut[i] = c * i ** v
        output_img = cv2.LUT(img, lut) #像素灰度值的映射
        output_img = np.uint8(output_img+0.5)
        return output_img
    
    #读取原始图像
    img = cv2.imread('test9.bmp')
    
    #绘制伽玛变换曲线
    gamma_plot(0.00000005, 4.0)
    
    #图像灰度伽玛变换
    output = gamma(img, 0.00000005, 4.0)
    
    #显示图像
    cv2.imshow('Imput', img)
    cv2.imshow('Output', output)
    cv2.waitKey(0)
    cv2.destroyAllWindows()

运行结果如下图所示：

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 9][]

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 10][]

--------------------

# 参考资料 #

\[1\] [https://blog.csdn.net/Eastmount/article/details/88858696][https_blog.csdn.net_Eastmount_article_details_88858696]

\[2\] [https://blog.csdn.net/Eastmount/article/details/88929290][https_blog.csdn.net_Eastmount_article_details_88929290]

\[3\] [冈萨雷斯. 数字图像处理(第三版)][.]

[20190702121722457.png_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly94aW9uZ3lpbWluZy5ibG9nLmNzZG4ubmV0_size_16_color_FFFFFF_t_70]: /images/20220205/155c7d6d6ff648a79de4c512285bf2b4.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70]: /images/20220205/e671e18112ab47108859ca82ae72d58e.png
[0_L-1]: https://private.codecogs.com/gif.latex?%5B0%2CL-1%5D
[s_L-1-r]: https://private.codecogs.com/gif.latex?s%3DL-1-r
[r]: https://private.codecogs.com/gif.latex?r
[s]: https://private.codecogs.com/gif.latex?s
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 1]: /images/20220205/79b3f991f55e428d9a8ef15dbac187b7.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 2]: /images/20220205/987db9f241554742bf4fc0c0d7a18da1.png
[s_c_log _1_r]: https://private.codecogs.com/gif.latex?s%3Dc%5Clog%20%281&plus;r%29
[c]: https://private.codecogs.com/gif.latex?c
[r_ge 0]: https://private.codecogs.com/gif.latex?r%5Cge%200
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 3]: /images/20220205/599691dcc6cd4d348b44c96fbaaef7ac.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 4]: /images/20220205/0c9ad6832888440995e1e389a877de31.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 5]: /images/20220205/9be745914b0a4843b38c2021a2d5a911.png
[s_c_r_gamma]: https://private.codecogs.com/gif.latex?s%3Dc%7B%7Br%7D%5E%7B%5Cgamma%20%7D%7D
[gamma]: https://private.codecogs.com/gif.latex?%5Cgamma
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 6]: /images/20220205/c1f9ed820bc148fab04321d101452c45.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 7]: /images/20220205/5b9e084a78944da786bb503cb42502fc.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 8]: /images/20220205/8768d2fa8c3449e590f878b7eb8cb561.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 9]: /images/20220205/0bd2af7ed9984c7e831d22c1cf69935e.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3phaXNodWl5aWZhbmd4eW0_size_16_color_FFFFFF_t_70 10]: /images/20220205/642a0785579245df99e682fc22aa6285.png
[https_blog.csdn.net_Eastmount_article_details_88858696]: https://blog.csdn.net/Eastmount/article/details/88858696
[https_blog.csdn.net_Eastmount_article_details_88929290]: https://blog.csdn.net/Eastmount/article/details/88929290
[.]: https://download.csdn.net/download/zaishuiyifangxym/10044192