Add PS模糊滤镜之水波效果算法.cpp

Author: BBuf

PhotoShop Algorithm/PS模糊滤镜之水波效果算法.cpp

@@ -0,0 +1,40 @@
+const double pi = acos(-1.0);
+Mat WaveFilter(Mat src) {
+	int row = src.rows;
+	int col = src.cols;
+	float A = 7;
+	float B = 2.5;
+	Point Center(col / 2, row / 2);
+	Mat dst(row, col, CV_8UC3);
+	for (int y = 0; y < row; y++) {
+		for (int x = 0; x < col; x++) {
+			float y0 = Center.y - y;
+			float x0 = x - Center.x;
+			float theta = atan(y0 / (x0 + 0.00001));
+			if (x0 < 0)
+				theta = theta + pi;
+			float r0 = sqrt(x0*x0 + y0*y0);
+			float r1 = r0 + A * col * 0.01 * sin(B * 0.1 * r0);
+			
+			float new_x = r1 * cos(theta);
+			float new_y = r1 * sin(theta);
+
+			new_x = Center.x + new_x;
+			new_y = Center.y + new_y;
+
+			if (new_x < 0) new_x = 0;
+			if (new_x >= col - 1) new_x = col - 2;
+			if (new_y < 0) new_y = 0;
+			if (new_y >= row - 1) new_y = row - 2;
+			int x1 = (int)new_x;
+			int y1 = (int)new_y;
+			float p = new_x - x1;
+			float q = new_y - y1;
+			for (int k = 0; k < 3; k++) {
+				dst.at<Vec3b>(y, x)[k] = (1 - p)*(1 - q)*src.at<Vec3b>(y1, x1)[k] + (p)*(1 - q)*src.at<Vec3b>(y1, x1 + 1)[k] + 
+					(1 - p)*(q)*src.at<Vec3b>(y1 + 1, x1)[k] + (p)*(q)*src.at<Vec3b>(y1 + 1, x1 + 1)[k];
+			}
+		}
+	}
+	return dst;
+}

PhotoShop Algorithm/README.md

@@ -33,4 +33,5 @@
 - PS色彩平衡(偏色校正).cpp 根据用户指定的R、G、B三个色彩的调整分量，分别附加到对应的色彩分量上，从而改变原始图像的色彩。
 - PS图像素描算法.cpp 原理就是图像去色后转置，然后接高斯滤波，最后对去色的图和模糊的图进行图层混合。
 - PS模糊滤镜之表面模糊算法.cpp 其作用是在保留图像边缘的情况下，对图像的表面进行模糊处理，效果比双边滤波好。
+- PS模糊滤镜之水波效果算法.cpp 这个算法原理不清楚，最先出处是这个博客：https://blog.csdn.net/matrix_space/article/details/42396829 。但他也没说明原理，还需研究下。