学习OpenCV：滤镜系列（15）——羽化（模糊边缘）

【原文：http://blog.csdn.net/yangtrees/article/details/9210153】

在PHOTOSHOP里，羽化就是使你选定范围的图边缘达到朦胧的效果。 羽化值越大，朦胧范围越宽，羽化值越小，朦胧范围越窄。可根据你想留下图的大小来调节。
算法分析：
1、通过对rgb值增加额外的V值实现朦胧效果
2、通过控制V值的大小实现范围控制。
3、V  = 255 * 当前点Point距中点距离的平方s1 / (顶点距中点的距离平方 *mSize)s2;
4、s1 有根据 ratio 修正 dx dy值。

#include <math.h>
#include <opencv/cv.h>
#include <opencv/highgui.h>
#define MAXSIZE (32768)
using namespace cv;
using namespace std;float mSize = 0.5;int main()
{Mat src = imread("D:/img/arrow04.jpg",1);imshow("src",src);int width=src.cols;int heigh=src.rows;int centerX=width>>1;int centerY=heigh>>1;int maxV=centerX*centerX+centerY*centerY;int minV=(int)(maxV*(1-mSize));int diff= maxV -minV;float ratio = width >heigh ? (float)heigh/(float)width : (float)width/(float)heigh;Mat img;src.copyTo(img);Scalar avg=mean(src);Mat dst(img.size(),CV_8UC3);Mat mask1u[3];float tmp,r;for (int y=0;y<heigh;y++){uchar* imgP=img.ptr<uchar>(y);uchar* dstP=dst.ptr<uchar>(y);for (int x=0;x<width;x++){int b=imgP[3*x];int g=imgP[3*x+1];int r=imgP[3*x+2];float dx=centerX-x;float dy=centerY-y;if(width > heigh)dx= (dx*ratio);elsedy = (dy*ratio);int dstSq = dx*dx + dy*dy;float v = ((float) dstSq / diff)*255;r = (int)(r +v);g = (int)(g +v);b = (int)(b +v);r = (r>255 ? 255 : (r<0? 0 : r));g = (g>255 ? 255 : (g<0? 0 : g));b = (b>255 ? 255 : (b<0? 0 : b));dstP[3*x] = (uchar)b;dstP[3*x+1] = (uchar)g;dstP[3*x+2] = (uchar)r;}}imshow("羽化",dst);waitKey();imwrite("D:/img/羽化.jpg",dst);}

Reference：http://www.cnblogs.com/lipeil/archive/2012/09/21/2696519.html
更加简单的方式：
分析PS的羽化结果可以知道，羽化达成了两个目的：1. 平滑轮廓线 2. 扩宽过渡区域
1.平滑轮廓线：可以采用均值滤波和cvSnakeImage()两种方式，前者维护一个宽度为H的窗口，窗口内均值滤波；而后者是OpenCV的C语言版本函数C++没有包含，其原理是能量最小化，经过测试前者的速度略高于后者，且当H较大时，可以采用窗口加权减一加一的方式来代替每次都求H次加权的方式；
2.扩宽过度区域：采用对mask采用全图均值滤波方法即可，卷积核的半径越大，过渡区域越宽。