【原文:http://blog.csdn.net/yangtrees/article/details/9210153】
#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);
else
dy = (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采用全图均值滤波方法即可,卷积核的半径越大,过渡区域越宽。