OpenCV图像金字塔：高斯金字塔、拉普拉斯金字塔与图片尺寸缩放
这篇已经写得很好，真心给作者点个赞。题目都是直接转过来的，直接去看吧。
Reference Link ： http://blog.csdn.net/poem_qianmo/article/details/26157633
In case：
【OpenCV入门教程之十三】OpenCV图像金字塔：高斯金字塔、拉普拉斯金字塔与图片尺寸缩放
这篇文章里，我们将一起探讨图像金字塔的一些基本概念，如何使用OpenCV函数 pyrUp 和 pyrDown 对图像进行向上和向下采样，以及了解了专门用于缩放图像尺寸的resize函数的用法。此博文一共有四个配套的简短的示例程序，其详细注释过的代码都在文中贴出，且文章最后提供了综合示例程序的下载。
先尝鲜一下其中一个示例程序的运行截图：
我们经常会将某种尺寸的图像转换为其他尺寸的图像，如果放大或者缩小图片的尺寸，笼统来说的话，可以使用OpenCV为我们提供的如下两种方式：
<1>resize函数。这是最直接的方式，
<2>pyrUp( )、pyrDown( )函数。即图像金字塔相关的两个函数，对图像进行向上采样，向下采样的操作。
pyrUp、pyrDown其实和专门用作放大缩小图像尺寸的resize在功能上差不多，披着图像金字塔的皮，说白了还是在对图像进行放大和缩小操作。另外需要指出的是，pyrUp、pyrDown在OpenCV的imgproc模块中的Image Filtering子模块里。而resize在imgproc 模块的Geometric Image Transformations子模块里。
这篇文章中，我们将先介绍图像金字塔的原理，接着介绍resize函数，然后是pyrUp和pyrDown函数，最后是一个综合示例程序。
二、关于图像金字塔
图像金字塔是图像中多尺度表达的一种，最主要用于图像的分割，是一种以多分辨率来解释图像的有效但概念简单的结构。
图像金字塔最初用于机器视觉和图像压缩，一幅图像的金字塔是一系列以金字塔形状排列的分辨率逐步降低，且来源于同一张原始图的图像集合。其通过梯次向下采样获得，直到达到某个终止条件才停止采样。
金字塔的底部是待处理图像的高分辨率表示，而顶部是低分辨率的近似。
我们将一层一层的图像比喻成金字塔，层级越高，则图像越小，分辨率越低。
素材图是帅气的美剧《绿箭侠》里面的绿箭侠Oliver Queen。
一般情况下有两种类型的图像金字塔常常出现在文献和以及实际运用中。他们分别是:
- 高斯金字塔Gaussianpyramid): 用来向下采样，主要的图像金字塔
- 拉普拉斯金字塔(Laplacianpyramid): 用来从金字塔低层图像重建上层未采样图像，在数字图像处理中也即是预测残差，可以对图像进行最大程度的还原，配合高斯金字塔一起使用。
两者的简要区别：高斯金字塔用来向下降采样图像，而拉普拉斯金字塔则用来从金字塔底层图像中向上采样重建一个图像。
要从金字塔第i层生成第i+1层（我们表示第i+1层为G_i+1），我们先要用高斯核对G_1进行卷积，然后删除所有偶数行和偶数列。当然的是，新得到图像面积会变为源图像的四分之一。按上述过程对输入图像G_0执行操作就可产生出整个金字塔。
当图像向金字塔的上层移动时，尺寸和分辨率就降低。OpenCV中，从金字塔中上一级图像生成下一级图像的可以用PryDown。而通过PryUp将现有的图像在每个维度都放大两遍。
图像金字塔中的向上和向下采样分别通过OpenCV函数 pyrUp 和 pyrDown 实现。
概括起来就是：
- 对图像向上采样：pyrUp函数
- 对图像向下采样：pyrDown函数
这里的向下与向上采样，是对图像的尺寸而言的（和金字塔的方向相反），向上就是图像尺寸加倍，向下就是图像尺寸减半。而如果我们按上图中演示的金字塔方向来理解，金字塔向上图像其实在缩小，这样刚好是反过来了。
但需要注意的是，PryUp和PryDown不是互逆的，即PryUp不是降采样的逆操作。这种情况下，图像首先在每个维度上扩大为原来的两倍，新增的行（偶数行）以0填充。然后给指定的滤波器进行卷积（实际上是一个在每个维度都扩大为原来两倍的过滤器）去估计“丢失”像素的近似值。
PryDown( )是一个会丢失信息的函数。为了恢复原来更高的分辨率的图像，我们要获得由降采样操作丢失的信息，这些数据就和拉普拉斯金字塔有关系了。
高斯金字塔
高斯金字塔是通过高斯平滑和亚采样获得一些列下采样图像，也就是说第K层高斯金字塔通过平滑、亚采样就可以获得K+1层高斯图像，高斯金字塔包含了一系列低通滤波器，其截至频率从上一层到下一层是以因子2逐渐增加，所以高斯金字塔可以跨越很大的频率范围。金字塔的图像如下：
这里增加一个注解： 在数字图像处理教材当中是从塔尖到塔底定义为0,1,2...J 层，而opencv是相反的。
另外，每一层都按从下到上的次序编号， 层级 G_i+1 (表示为 G_i+1尺寸小于第i层G_i)。
2.1.1 对图像的向下取样
为了获取层级为 G_i+1 的金字塔图像，我们采用如下方法:
<1>对图像G_i进行高斯内核卷积
<2>将所有偶数行和列去除
得到的图像即为G_i+1的图像，显而易见，结果图像只有原图的四分之一。通过对输入图像G_i(原始图像)不停迭代以上步骤就会得到整个金字塔。同时我们也可以看到，向下取样会逐渐丢失图像的信息。
以上就是对图像的向下取样操作，即缩小图像。
2.1.2 对图像的向上取样
如果想放大图像，则需要通过向上取样操作得到，具体做法如下：
<1>将图像在每个方向扩大为原来的两倍，新增的行和列以0填充
<2>使用先前同样的内核(乘以4)与放大后的图像卷积，获得 “新增像素”的近似值 此处乘以4是因为在Mask（5*5）中只有6个位置是有值的，*4是一种近似，接近25个值。
得到的图像即为放大后的图像，但是与原来的图像相比会发觉比较模糊，因为在缩放的过程中已经丢失了一些信息，如果想在缩小和放大整个过程中减少信息的丢失，这些数据形成了拉普拉斯金字塔。
那么，我们接下来一起看一看拉普拉斯金字塔的概念吧。
拉普拉斯金字塔
下式是拉普拉斯金字塔第i层的数学定义：
式中的
表示第i层的图像。而UP（）操作是将源图像中位置为(x,y)的像素映射到目标图像的(2x+1,2y+1)位置，即在进行向上取样。符号
表示卷积，
为5x5的高斯内核。
我们下文将要介绍的pryUp，就是在进行上面这个式子的运算。
因此，我们可以直接用OpenCV进行拉普拉斯运算：
也就是说，拉普拉斯金字塔是通过源图像减去先缩小后再放大的图像的一系列图像构成的。
整个拉普拉斯金字塔运算过程可以通过下图来概括：
所以，我们可以将拉普拉斯金字塔理解为高斯金字塔的逆形式。
另外再提一点，关于图像金字塔非常重要的一个应用就是实现图像分割。图像分割的话，先要建立一个图像金字塔，然后在G_i和G_i+1的像素直接依照对应的关系，建立起”父与子“关系。而快速初始分割可以先在金字塔高层的低分辨率图像上完成，然后逐层对分割加以优化。
三、resize( )函数剖析
resize( )为OpenCV中专职调整图像大小的函数。
此函数将源图像精确地转换为指定尺寸的目标图像。如果源图像中设置了ROI（Region Of Interest ，感兴趣区域），那么resize( )函数会对源图像的ROI区域进行调整图像尺寸的操作，来输出到目标图像中。若目标图像中已经设置ROI区域，不难理解resize( )将会对源图像进行尺寸调整并填充到目标图像的ROI中。
很多时候，我们并不用考虑第二个参数dst的初始图像尺寸和类型（即直接定义一个Mat类型，不用对其初始化），因为其尺寸和类型可以由src,dsize,fx和fy这其他的几个参数来确定。
看一下它的函数原型：
- resize(InputArray src,OutputArray dst, Size dsize, double fx=0, double fy=0, interpolation=INTER_LINEAR )
C++: void resize(InputArray src,OutputArray dst, Size dsize, double fx=0, double fy=0, int interpolation=INTER_LINEAR )
- 第一个参数，InputArray类型的src，输入图像，即源图像，填Mat类的对象即可。
- 第二个参数，OutputArray类型的dst，输出图像，当其非零时，有着dsize（第三个参数）的尺寸，或者由src.size()计算出来。
- 第三个参数，Size类型的dsize，输出图像的大小;如果它等于零，由下式进行计算：
其中，dsize，fx，fy都不能为0。
- 第四个参数，double类型的fx，沿水平轴的缩放系数，有默认值0，且当其等于0时，由下式进行计算：
- 第五个参数，double类型的fy，沿垂直轴的缩放系数，有默认值0，且当其等于0时，由下式进行计算：
- 第六个参数，int类型的interpolation，用于指定插值方式，默认为INTER_LINEAR（线性插值）。
可选的插值方式如下：
- INTER_NEAREST - 最近邻插值
- INTER_LINEAR - 线性插值（默认值）
- INTER_AREA - 区域插值（利用像素区域关系的重采样插值）
- INTER_CUBIC –三次样条插值（超过4×4像素邻域内的双三次插值）
- INTER_LANCZOS4 -Lanczos插值（超过8×8像素邻域的Lanczos插值）
若要缩小图像，一般情况下最好用CV_INTER_AREA来插值，
而若要放大图像，一般情况下最好用CV_INTER_CUBIC（效率不高，慢，不推荐使用）或CV_INTER_LINEAR（效率较高，速度较快，推荐使用）。
关于插值，我们看几张图就能更好地理解。先看原图：
当进行6次图像缩小接着6次图像放大操作后，两种不同的插值方式得到的效果图：
效果很明显，第一张全是一个个的像素，非常影响美观。另外一张却有雾化的朦胧美感，所以插值方式的选择，对经过多次放大缩小的图片最终得到的效果是有很大影响的。
接着我们来看两种resize的调用范例。
方式一，调用范例： Mat dst=Mat::zeros(512 ,512, CV_8UC3 );//新建一张512x512尺寸的图片 Mat src=imread(“1.jpg”); //显式指定dsize=dst.size（）,那么fx和fy会其计算出来，不用额外指定。 resize(src, dst, dst.size()); Mat dst=Mat::zeros(512 ,512, CV_8UC3 );//新建一张512x512尺寸的图片 Mat src=imread(“1.jpg”); //显式指定dsize=dst.size（）,那么fx和fy会其计算出来，不用额外指定。 resize(src, dst, dst.size()); 方式二、调用范例： Mat dst; Mat src=imread(“1.jpg”) //指定fx和fy，让函数计算出目标图像的大小。 resize(src, dst, Size(), 0.5, 0.5); Mat dst; Mat src=imread(“1.jpg”) //指定fx和fy，让函数计算出目标图像的大小。 resize(src, dst, Size(), 0.5, 0.5); 接着我们看看完整的示例程序： //-----------------------------------【头文件包含部分】--------------------------------------- // 描述：包含程序所依赖的头文件 //---------------------------------------------------------------------------------------------- #include <opencv2/opencv.hpp> #include <opencv2/imgproc/imgproc.hpp> //-----------------------------------【命名空间声明部分】--------------------------------------- // 描述：包含程序所使用的命名空间 //----------------------------------------------------------------------------------------------- usingnamespace //-----------------------------------【main( )函数】-------------------------------------------- // 描述：控制台应用程序的入口函数，我们的程序从这里开始 //----------------------------------------------------------------------------------------------- main( ) //载入原始图 Mat srcImage = imread("1.jpg"//工程目录下应该有一张名为1.jpg的素材图 Mat tmpImage,dstImage1,dstImage2;//临时变量和目标图的定义 tmpImage=srcImage;//将原始图赋给临时变量 //显示原始图 imshow("【原始图】", srcImage); //进行尺寸调整操作 resize(tmpImage,dstImage1,Size( tmpImage.cols/2, tmpImage.rows/2 ),(0,0),(0,0),3); resize(tmpImage,dstImage2,Size( tmpImage.cols*2, tmpImage.rows*2 ),(0,0),(0,0),3); //显示效果图 imshow("【效果图】之一", dstImage1); imshow("【效果图】之二", dstImage2); waitKey(0); return //-----------------------------------【头文件包含部分】--------------------------------------- // 描述：包含程序所依赖的头文件 //---------------------------------------------------------------------------------------------- #include <opencv2/opencv.hpp> #include <opencv2/imgproc/imgproc.hpp> //-----------------------------------【命名空间声明部分】--------------------------------------- // 描述：包含程序所使用的命名空间 //----------------------------------------------------------------------------------------------- using namespace cv; //-----------------------------------【main( )函数】-------------------------------------------- // 描述：控制台应用程序的入口函数，我们的程序从这里开始 //----------------------------------------------------------------------------------------------- int main( ) { //载入原始图 Mat srcImage = imread("1.jpg"); //工程目录下应该有一张名为1.jpg的素材图 Mat tmpImage,dstImage1,dstImage2;//临时变量和目标图的定义 tmpImage=srcImage;//将原始图赋给临时变量 //显示原始图 imshow("【原始图】", srcImage); //进行尺寸调整操作 resize(tmpImage,dstImage1,Size( tmpImage.cols/2, tmpImage.rows/2 ),(0,0),(0,0),3); resize(tmpImage,dstImage2,Size( tmpImage.cols*2, tmpImage.rows*2 ),(0,0),(0,0),3); //显示效果图 imshow("【效果图】之一", dstImage1); imshow("【效果图】之二", dstImage2); waitKey(0); return 0; }
程序运行截图：
四、pyrUp()函数剖析
pyrUp( )函数的作用是向上采样并模糊一张图像，说白了就是放大一张图片。
- pyrUp(InputArray src, OutputArraydst, const Size& dstsize=Size(), borderType=BORDER_DEFAULT )
C++: void pyrUp(InputArray src, OutputArraydst, const Size& dstsize=Size(), int borderType=BORDER_DEFAULT )
- 第一个参数，InputArray类型的src，输入图像，即源图像，填Mat类的对象即可。
- 第二个参数，OutputArray类型的dst，输出图像，和源图片有一样的尺寸和类型。
- 第三个参数，const Size&类型的dstsize，输出图像的大小;有默认值Size()，即默认情况下，由Size（src.cols*2，src.rows*2）来进行计算，且一直需要满足下列条件：
- 第四个参数，int类型的borderType，又来了，边界模式，一般我们不用去管它。
pyrUp函数执行高斯金字塔的采样操作，其实它也可以用于拉普拉斯金字塔的。
首先，它通过插入可为零的行与列，对源图像进行向上取样操作，然后将结果与pyrDown（）乘以4的内核做卷积，就是这样。
直接看完整的示例程序：
//-----------------------------------【头文件包含部分】--------------------------------------- // 描述：包含程序所依赖的头文件 //---------------------------------------------------------------------------------------------- #include <opencv2/opencv.hpp> #include <opencv2/imgproc/imgproc.hpp> //-----------------------------------【命名空间声明部分】--------------------------------------- // 描述：包含程序所使用的命名空间 //----------------------------------------------------------------------------------------------- usingnamespace //-----------------------------------【main( )函数】-------------------------------------------- // 描述：控制台应用程序的入口函数，我们的程序从这里开始 //----------------------------------------------------------------------------------------------- main( ) //载入原始图 Mat srcImage = imread("1.jpg"//工程目录下应该有一张名为1.jpg的素材图 Mat tmpImage,dstImage;//临时变量和目标图的定义 tmpImage=srcImage;//将原始图赋给临时变量 //显示原始图 imshow("【原始图】", srcImage); //进行向上取样操作 pyrUp( tmpImage, dstImage, Size( tmpImage.cols*2, tmpImage.rows*2 ) ); //显示效果图 imshow("【效果图】", dstImage); waitKey(0); return //-----------------------------------【头文件包含部分】--------------------------------------- // 描述：包含程序所依赖的头文件 //---------------------------------------------------------------------------------------------- #include <opencv2/opencv.hpp> #include <opencv2/imgproc/imgproc.hpp> //-----------------------------------【命名空间声明部分】--------------------------------------- // 描述：包含程序所使用的命名空间 //----------------------------------------------------------------------------------------------- using namespace cv; //-----------------------------------【main( )函数】-------------------------------------------- // 描述：控制台应用程序的入口函数，我们的程序从这里开始 //----------------------------------------------------------------------------------------------- int main( ) { //载入原始图 Mat srcImage = imread("1.jpg"); //工程目录下应该有一张名为1.jpg的素材图 Mat tmpImage,dstImage;//临时变量和目标图的定义 tmpImage=srcImage;//将原始图赋给临时变量 //显示原始图 imshow("【原始图】", srcImage); //进行向上取样操作 pyrUp( tmpImage, dstImage, Size( tmpImage.cols*2, tmpImage.rows*2 ) ); //显示效果图 imshow("【效果图】", dstImage); waitKey(0); return 0; }
程序运行截图：
五、pyrDown()函数剖析
pyrDown( )函数的作用是向下采样并模糊一张图片，说白了就是缩小一张图片。
- pyrDown(InputArray src,OutputArray dst, const Size& dstsize=Size(), borderType=BORDER_DEFAULT)
C++: void pyrDown(InputArray src,OutputArray dst, const Size& dstsize=Size(), int borderType=BORDER_DEFAULT)
- 第一个参数，InputArray类型的src，输入图像，即源图像，填Mat类的对象即可。
- 第二个参数，OutputArray类型的dst，输出图像，和源图片有一样的尺寸和类型。
- 第三个参数，const Size&类型的dstsize，输出图像的大小;有默认值Size()，即默认情况下，由Size Size((src.cols+1)/2, (src.rows+1)/2)来进行计算，且一直需要满足下列条件：
该pyrDown函数执行了高斯金字塔建造的向下采样的步骤。首先，它将源图像与如下内核做卷积运算：
接着，它便通过对图像的偶数行和列做插值来进行向下采样操作。
依然是看看完整的示例程序：
程序运行截图：
六、综合示例篇——在实战中熟稔
依然是每篇文章都会配给大家的一个详细注释的博文配套示例程序，把这篇文章中介绍的知识点以代码为载体，展现给大家。
这个示例程序中，分别演示了用resize，pryUp，pryDown来让源图像进行放大缩小的操作，分别用键盘按键1、2、3、4、A、D、W、S来控制图片的放大与缩小：
OK，上详细注释的代码吧：
放一些程序运行截图。
经过多次按键后的效果图：
另外，还可以放大图像到很大的尺寸，上图的话会很凶残而且不美观。所以就不放出超过原图尺寸的截图了。
好的，就放出这些效果图，具体更多的运行效果大家就自己下载示例程序回去玩~
本篇文章的配套源代码请点击这里下载：
------------------------------- 下面附上python 的图像金字塔代码
# -*- coding: utf-8 -*- import cv2 import numpy as np A = cv2.imread('../burn image/case3-burn-to-head-showing-healing-0.jpg',0) # generate Gaussian pyramid for A G = A.copy() gpA = [G] lpA = [] for i in xrange(3): cv2.imshow('1',G) image_o = G G = cv2.pyrDown(G) GE = cv2.pyrUp(G) L = np.subtract(image_o,GE) gpA.append(G) lpA.append(L) cv2.imshow('lpA',L ) cv2.imshow('上采样+拉普拉斯',L+GE ) cv2.imshow('upsample',GE ) MSE(image_o,L+GE) cv2.waitKey(0)
注：拉普拉斯金子塔是每层的原图减去该图的下采样后的上采样图(一般用高斯下采样)， 所以也叫预测残差金字塔。 通过预测残差金字塔可以100%还原图像，对么？
