【youcans 的 OpenCV 例程200篇】203. 伪彩色图像处理

OpenCV 例程200篇 总目录

【youcans 的 OpenCV 例程200篇】203. 伪彩色图像处理

NASA 公布了蟹状星云 （Crab Nebula ）的观测图像。茫茫太空，距离我们几亿光年的宇宙真是这样绚丽迷人吗？

NASA 专家撰文指出，蟹状星云彩色照片实际上人工合成的图像，这是不是暗示照片中的彩色是伪造的？

什么是伪彩色图像？

伪彩色图像是指对单色图像进行处理，结果转换得到颜色分量，构造为彩色效果的图像。伪彩色图像在形式和视觉表现为彩色图像，但其所呈现的颜色并非图像的真实色彩重现，仅仅是各颜色分量的像素值合成的结果。

• （1）一类伪彩色图像是对灰度图像的增强，形成彩色效果。

伪彩色图像增强是指按照特定的准则对灰度图像进行处理，将不同的灰度级按照某种映射关系变换为不同的颜色分量。例如，天气预报中的气象云图，红外测温图像，一般都是伪彩色图像。

人眼的视觉特性，只能分辨 20 级左右的灰度，但可以分辨几千种色调和亮度。通过伪彩色图像增强，可以让灰度图像看起来更清楚，更容易分辨。

将灰度图像进行伪彩色增强，通常是基于一个设定的颜色查找表，将图像像素的灰度值替换为颜色查找表中对应的颜色值。这显然是典型的 LUT 应用场景。

• （2）另一类伪彩色图像，是多光谱图像的编码合成图像。

光谱图像可以是可见光谱，也可以是红外、紫外、X光或其它信号。这些非可见光谱图像实际上既不是彩色图像，也不是灰度图像，“像素值”往往反映某种检测信号的强度。合成为彩色图像进行显示，只是为了便于观察和分析。

这在天文学观测中十分常用。NASA 公布的绚烂迷人的太空照片，并不是真实世界的色彩还原，而是将多种光谱信号图像进行编码组合和特效处理后，构造出的彩色图像。

函数原型

函数 cv.applyColorMap() 根据色彩映射表，将灰度图像变换为伪彩色图像。

cv.applyColorMap(src, colormap[, dst]) → dst
cv.applyColorMap(src, userColor[, dst]) → dst

参数说明：

• src：输入图像，8 位灰度图像或彩色图像，cv_8U
• dst：输出图像，大小和通道数与 src 相同
• colormap：色彩映射表，OpenCV 自带色彩风格类型的颜色查找表
• userColor：用户自定义的色彩映射表， 256个元素

注意事项：

1. 输入图像可以是彩色图像。根据作者的测试，先将输入彩色图像转换为灰度图像，再按色彩映射表进行变换。
2. 用户自定义的色彩映射表，是 256*3 数组。系统提供的 22种色彩映射表一般已经够用，就别折腾了。
3. 例程中没有出现 cv.LUT 查表替换，在 cv.applyColorMap 内核中已经应用了查表替换方法。

色彩映射表：

• OpenCV 提供了 22 种色彩风格类型，类型描述关键字与色彩效果如下图所示。

• OpenCV 提供的色彩风格类型与 Matlab、Matplotlib 热图的颜色映射是类似的，描述关键字也是一致的。
  参见：http://matplotlib.org/examples/color/colormaps_reference.html

例程 14.5：灰度图像转换为伪彩色图像

    # 14.5 灰度图像转换为伪彩色图像    
    gray = cv.imread("../images/Fig0525a.tif", flags=0)  # 读取灰度图像
    h, w = gray.shape[:2]  # 图片的高度, 宽度

    # 伪彩色处理
    pseudo1 = cv.applyColorMap(gray, colormap=cv.COLORMAP_HOT)
    pseudo2 = cv.applyColorMap(gray, colormap=cv.COLORMAP_PINK)
    pseudo3 = cv.applyColorMap(gray, colormap=cv.COLORMAP_RAINBOW)
    pseudo4 = cv.applyColorMap(gray, colormap=cv.COLORMAP_HSV)
    pseudo5 = cv.applyColorMap(gray, colormap=cv.COLORMAP_TURBO)

    plt.figure(figsize=(9, 6))
    plt.subplot(231), plt.axis('off'), plt.title("GRAY"), plt.imshow(gray, cmap='gray')
    plt.subplot(232), plt.axis('off'), plt.title("cv.COLORMAP_HOT")
    plt.imshow(cv.cvtColor(pseudo1, cv.COLOR_BGR2RGB))
    plt.subplot(233), plt.axis('off'), plt.title("cv.COLORMAP_PINK")
    plt.imshow(cv.cvtColor(pseudo2, cv.COLOR_BGR2RGB))
    plt.subplot(234), plt.axis('off'), plt.title("cv.COLORMAP_RAINBOW")
    plt.imshow(cv.cvtColor(pseudo3, cv.COLOR_BGR2RGB))
    plt.subplot(235), plt.axis('off'), plt.title("cv.COLORMAP_HSV")
    plt.imshow(cv.cvtColor(pseudo4, cv.COLOR_BGR2RGB))
    plt.subplot(236), plt.axis('off'), plt.title("cv.COLORMAP_TURBO")
    plt.imshow(cv.cvtColor(pseudo5, cv.COLOR_BGR2RGB))
    plt.tight_layout()
    plt.show()

例程 14.7：多光谱合成的伪彩色图像

蟹状星云 （Crab Nebula ）的观测图像，是由钱德拉 X 射线天文台数据所形成的X射线图、哈勃太空望远镜拍摄的光学图像，和斯必泽空间望远镜拍摄的红外图像合成得到的。

例程结果图中，第一行是 3 种天文观测图像的观测数据图像，图像中的灰度都只是观测的信号强度，而不是实际的物理世界的亮度或灰度。

第二行是分别对 3 种观测图像进行伪彩色变换的结果，如果选择不同的色彩风格可以获得不同的效果。

第三行是将 3 种观测信号合成得到的伪彩色图像，其中左图、中图是简单地将 3 种观测信号分别作为 RGB 的一个色彩通道得到的伪彩色图像，右图是 NASA 公布的多光谱合成伪彩色图像。

虽然 NASA 图像的视觉效果更好，但仍然可以明显地看出，该图也是由第一行 3 种观测数据图像合成的。

    # 14.7 多光谱编码合成的伪彩色图像
    # 蟹状星云 (Crab Nebula), 图片来源：
    # https://www.nasa.gov/mission_pages/chandra/multimedia/photo09-096.html
    composite = cv.imread("../images/CrabNebula.png", flags=1)  # 读取多光谱合成图像
    nebulaOpti = cv.imread("../images/CrabNebula_Optical.jpg", flags=1)  # 读取 Optical
    nebulaXray = cv.imread("../images/CrabNebula_Xray.jpg", flags=1)  # 读取 Xray
    nebulaInfr = cv.imread("../images/CrabNebula_Infrared.jpg", flags=1)  # 读取 Infrared
    h, w = nebulaOpti.shape[:2]  # 图片的高度, 宽度

    grayOpti = cv.cvtColor(nebulaOpti, cv.COLOR_BGR2GRAY)
    grayXray = cv.cvtColor(nebulaXray, cv.COLOR_BGR2GRAY)
    grayInfr = cv.cvtColor(nebulaInfr, cv.COLOR_BGR2GRAY)

    # 伪彩色处理
    pseudoXray = cv.applyColorMap(nebulaXray, colormap=cv.COLORMAP_TURBO)
    pseudoOpti = cv.applyColorMap(nebulaOpti, colormap=cv.COLORMAP_MAGMA)
    pseudoInfr = cv.applyColorMap(nebulaInfr, colormap=cv.COLORMAP_HOT)

    # 多光谱编码合成
    compose1 = np.zeros((h, w, 3), np.uint8)  # 创建黑色图像 BGR=0
    compose1[:, :, 0] = grayOpti  # Optical -> B
    compose1[:, :, 1] = grayXray  # Xray -> G
    compose1[:, :, 2] = grayInfr  # Infrared -> R

    compose2 = np.zeros((h, w, 3), np.uint8)  # 创建黑色图像 BGR=0
    compose2[:, :, 0] = grayXray  # Xray -> B
    compose2[:, :, 1] = grayOpti  # Optical -> G
    compose2[:, :, 2] = grayInfr  # Infrared -> R

    # composite = cv.normalize(nebulaInfr + nebulaOpti, None, 0, 255, cv.NORM_MINMAX)  # 归一化为 [0,255]

    plt.figure(figsize=(9, 9))
    plt.subplot(331), plt.axis('off'), plt.title("CrabNebula-Xray")
    plt.imshow(grayXray, cmap='gray')
    plt.subplot(332), plt.axis('off'), plt.title("CrabNebula-Optical")
    plt.imshow(grayOpti, cmap='gray')
    plt.subplot(333), plt.axis('off'), plt.title("CrabNebula-Infrared")
    plt.imshow(grayInfr, cmap='gray')
    plt.subplot(334), plt.axis('off'), plt.title("Xray Pseudo")
    plt.imshow(cv.cvtColor(pseudoXray, cv.COLOR_BGR2RGB))
    plt.subplot(335), plt.axis('off'), plt.title("Optical Pseudo")
    plt.imshow(cv.cvtColor(pseudoOpti, cv.COLOR_BGR2RGB))
    plt.subplot(336), plt.axis('off'), plt.title("Infrared Pseudo")
    plt.imshow(cv.cvtColor(pseudoInfr, cv.COLOR_BGR2RGB))
    plt.subplot(337), plt.axis('off'), plt.title("Spectral composited")
    plt.imshow(cv.cvtColor(compose1, cv.COLOR_BGR2RGB))
    plt.subplot(338), plt.axis('off'), plt.title("Spectral composited")
    plt.imshow(cv.cvtColor(compose2, cv.COLOR_BGR2RGB))
    plt.subplot(339), plt.axis('off'), plt.title("Composite image")
    plt.imshow(cv.cvtColor(composite, cv.COLOR_BGR2RGB))
    plt.tight_layout()
    plt.show()

---

【本节完】

版权声明：
youcans@xupt 原创作品，转载必须标注原文链接：(https://blog.csdn.net/youcans/article/details/125298385)
Copyright 2022 youcans, XUPT
Crated：2022-6-14
欢迎关注 『youcans 的 OpenCV 例程 200 篇』 系列，持续更新中
欢迎关注 『youcans 的 OpenCV学习课』 系列，持续更新中