OpenCV——图像金字塔与拉普拉斯金字塔

前言

这篇文章主要讲的是图像金字塔，拉普拉斯金字塔。在说明什么是图像金字塔的过程之前，我们需要需要了解一个概念：尺度
尺度，顾名思义就是说图像的尺寸和分辨率。在我们进行图像处理的时候，会经常对源图像的尺寸进行放大或者缩小的变换，进而转换为我们指定尺寸的目标图像。在对图像进行放大和缩小的变换的这个过程，我们称为尺度调整。
而图像金字塔则是图像多尺度调整表达的一种重要的方式，图像金字塔方法的原理是：将参加融合的的每幅图像分解为多尺度的金字塔图像序列，将低分辨率的图像在上层，高分辨率的图像在下层，上层图像的大小为前一层图像大小的1/4。层数为0,1,2……N。将所有图像的金字塔在相应层上以一定的规则融合，就可得到合成金字塔，再将该合成金字塔按照金字塔生成的逆过程进行重构，得到融合金字塔。这个总的思路就是一下所有基于金字塔融合的算法过程，不同点就在于分解构造的金字塔不同，每层的融合规则不一样，重构的方法不同而已。

原理

高斯金字塔

高斯金字塔是最基本的图像塔。首先将原图像作为最底层图像G0（高斯金字塔的第0层），利用高斯核（5*5）对其进行卷积(卷积的解释)，然后对卷积后的图像进行下采样（去除偶数行和列）得到上一层图像G1，将此图像作为输入，重复卷积和下采样操作得到更上一层图像，反复迭代多次，形成一个金字塔形的图像数据结构，即高斯金字塔。

Opencv中使用pyrdown函数就可以获得高斯金字塔。

拉普拉斯金字塔

在高斯金字塔的运算过程中，图像经过卷积和下采样操作会丢失部分高频细节信息。为描述这些高频信息，人们定义了拉普拉斯金字塔(Laplacian Pyramid， LP)。用高斯金字塔的每一层图像减去其上一层图像上采样并高斯卷积之后的预测图像，得到一系列的差值图像即为 LP 分解图像。
在使用拉普拉斯金字塔操作中，实际上也就实现了上采样的过程
①：首先，将图像在每个方向扩大为原来的两倍，新增的行和列以0填充(0)
②：使用先前同样的内核(乘以4)与放大后的图像卷积，获得 “新增像素” 的近似值。

正文

PyrDown降采样

效果图

函数：pyrDown(src[, dst[, dstsize[, borderType]]])

src参数表示输入图像。
dst参数表示输出图像，它与src类型、大小相同,如果不进行制定就有返回值。
dstsize参数表示降采样之后的目标图像的大小。它是有默认值的，如果我们调用函数的时候不指定第三个参数，那么这个值是按照 Size((src.cols+1)/2, (src.rows+1)/2) 计算的。而且不管你自己如何指定这个参数，一定必须保证满足以下关系式：|dstsize.width * 2 - src.cols| ≤ 2; |dstsize.height * 2 - src.rows| ≤ 2。也就是说降采样的意思其实是把图像的尺寸缩减一半，行和列同时缩减一半。
borderType参数表示表示图像边界的处理方式。

方法

1. 先复制原始图像为temp
2. 使用cv.pyrDown进行降采样。
3. 不断循环。

code

def pyramid_demo(image):
    level = 4
    temp = image.copy()
    pyramid_images = []
    for i in range(level):
        dst = cv.pyrDown(temp)
        pyramid_images.append(dst)
        cv.imshow("pyramid_down_"+str(i+1),dst)
        temp = dst.copy()
    return pyramid_images

pyrUp函数升采样

效果图：

函数：pyrUp(src[, dst[, dstsize[, borderType]]])

src参数表示输入图像。
dst参数表示输出图像，它与src类型、大小相同。
dstsize参数表示降采样之后的目标图像的大小。在默认的情况下，这个尺寸大小是按照 Size(src.cols2, (src.rows2) 来计算的。如果你自己要指定大小，那么一定要满足下面的条件：
|dstsize.width - src.cols * 2| ≤ (dstsize.width mod 2); //如果width是偶数，那么必须dstsize.width是src.cols的2倍
|dstsize.height - src.rows * 2| ≤ (dstsize.height mod 2);
borderType参数表示表示图像边界的处理方式。

方法：

1. 先用pyramid_demo进行降采样，得到各个降采样的图片。
2. 接下来，对这些降采样的图像进行pyrUp扩展。然后，用原始图像进去该扩展后的图像。
3. 就得到了那部分在降采样时被搞掉的那部分像素。

code

def lapalas_demo(image):
    pyramid_images = pyramid_demo(image)
    level = len(pyramid_images)

    for i in range(level-1,-1,-1):
        if i-1<0:
            expand = cv.pyrUp(pyramid_images[i],dstsize=image.shape[:2])
            lpls = cv.subtract(image,expand)
            cv.imshow("laplace_demo"+str(i),lpls)
        else:
            expand = cv.pyrUp(pyramid_images[i], dstsize=pyramid_images[i - 1].shape[:2])
            lpls = cv.subtract(pyramid_images[i - 1], expand)
            cv.imshow("laplace_demo" + str(i), lpls)

参考文献

1. CV学习笔记（十二）：图像金字塔