目录
1 金字塔
平时你听到、见到的金字塔是什么样的？
这样？
还是这样？
实际上除了这些，还有 图像金字塔
图像金字塔有什么用？为什么要称作图像金字塔？本文带你研究这些问题。
; 2 什么是图像金字塔？
正如生物视觉系统会处理分层次的尺寸一样，计算机视觉系统实现多分辨率图像处理的基础是 图像金字塔。
考虑这样一个场景：输入系统一幅图像来检测人脸。由于事先并不知道人脸在这张图片中可能的尺寸，所以需要根据输入生成一个不同大小图像组成的金字塔(应用时常用向量存储)，扫描不同层次来检测可能的人脸。
此外，图像金字塔在保持细节的条件下进行图像融合等多尺度编辑操作非常有用。
图像金字塔分为两种：
- 高斯金字塔
高斯金字塔是最基本的图像金字塔，用于获得原图像尺度连续的的降采样序列。高斯金字塔共分O O O组(Octave)，每组分为S S S层(Layer)。组内各层图像分辨率相同但尺度逐层递增；组间图像按隔点降采样取得。
- 拉普拉斯金字塔
拉普拉斯金字塔通过高斯差分得到残差图像序列，用于图像重建、图像融合、特征点检测等。具体做法如图所示：构造高斯金字塔，每组图像内相邻两层相减得到残差图像序列，从而得到拉普拉斯金字塔。所以高斯金字塔每组有s + 3 s+3 s +3层，而拉普拉斯金字塔每组只有s + 2 s+2 s +2层。
3 图像金字塔有什么用？
前面已经提到，图像金字塔的作用主要有：
- 图像重建
- 图像多尺度特征检测
- 图像修复
- 图像融合
- …
下面给出一个基于图像金字塔的图像融合算法。
- 选择待融合图像A与B，以及融合掩码M，并重采样使之具有相同大小；
- 分别取得A、B的拉普拉斯金字塔，以及M的高斯金字塔；
- 使用M的高斯金字塔作为加权，对A、B的每层特征进行融合，得到融合金字塔；
- 融合金字塔仍为拉普拉斯金字塔，对其图像重建得到融合图像。
上面的图可能有点抽象，看看下面的例子。
; 4 OpenCV实战图像金字塔
主函数非常简单易懂
int main()
{
Mat img_1 = imread("1.jpg", 1);
Mat img_2 = imread("2.jpg", 1);
Mat dstImg;
imgFusion(img_1, img_2, dstImg, 0.2);
waitKey(0);
return 0;
}
试试效果
原图
融合效果图
是不是毫无违和感？其中关键的融合函数 imgFusion()具体是如何实现的呢？
void imgFusion(Mat leftImg, Mat rightImg, Mat& dstImg, float threshold)
{
vector<Mat> gaussPyrLeft, gaussPyrRight, laplacePyrLeft, laplacePyrRight;
vector<Mat> maskGaussPyr;
vector<Mat> blendLapPyr;
Mat imgHighest;
Mat mask = Mat::zeros(PYRHEIGHT, PYRWIDTH, CV_32FC1);
mask(Range::all(), Range(mask.cols * threshold, mask.cols)) = 1.0;
cvtColor(mask, mask, COLOR_GRAY2BGR);
buildGaussPyr(mask, maskGaussPyr, 3);
resize(leftImg, leftImg, Size(PYRWIDTH, PYRHEIGHT));
resize(rightImg, rightImg, Size(PYRWIDTH, PYRHEIGHT));
leftImg.convertTo(leftImg, CV_32F);
rightImg.convertTo(rightImg, CV_32F);
buildGaussPyr(leftImg, gaussPyrLeft, 3);
buildGaussPyr(rightImg, gaussPyrRight, 3);
buildLaplacePyr(gaussPyrLeft, laplacePyrLeft, 3);
buildLaplacePyr(gaussPyrRight, laplacePyrRight, 3);
imgHighest = gaussPyrLeft.back().mul(maskGaussPyr.back()) +
((gaussPyrRight.back()).mul(Scalar(1.0, 1.0, 1.0) - maskGaussPyr.back()));
blendLaplacePyr(laplacePyrLeft, laplacePyrRight, maskGaussPyr, blendLapPyr);
dstImg = imgLaplaceBlend(imgHighest, blendLapPyr);
dstImg.convertTo(dstImg, CV_8UC3);
imshow("imgProcess::seamOpt_laplace", dstImg);
}
限于篇幅，其中部分函数的实现无法贴完整，本文的完整工程代码请通过下方名片联系我获取
🚀 计算机视觉基础教程说明
章号 内容
0 色彩空间与数字成像
1 计算机几何基础
2 图像增强、滤波、金字塔
3 图像特征提取
4 图像特征描述
5 图像特征匹配
6 立体视觉
7 项目实战
🔥 更多精彩专栏：
👇源码获取 · 技术交流 · 抱团学习 · 咨询分享 请联系👇
Original: https://blog.csdn.net/FRIGIDWINTER/article/details/123139599
Author: Mr.Winter`
Title: 计算机视觉教程2-5：图像金字塔送她一朵玫瑰(附代码)
原创文章受到原创版权保护。转载请注明出处：https://www.johngo689.com/633145/
转载文章受原作者版权保护。转载请注明原作者出处！