OpenCV——PS滤镜算法之 Ellipsoid (凹陷)
2016-07-02 09:05
267 查看
// define head function #ifndef PS_ALGORITHM_H_INCLUDED #define PS_ALGORITHM_H_INCLUDED #include <iostream> #include <string> #include "cv.h" #include "highgui.h" #include "cxmat.hpp" #include "cxcore.hpp" #include "math.h" using namespace std; using namespace cv; void Show_Image(Mat&, const string &); #endif // PS_ALGORITHM_H_INCLUDED #include "PS_Algorithm.h" #include <time.h> using namespace std; using namespace cv; #define pi 3.1415926 int main() { string Img_name("4.jpg"); Mat Img; Img=imread(Img_name); Mat Img_out(Img.size(), CV_8UC3); Img.copyTo(Img_out); int width=Img.cols; int height=Img.rows; float e; float a,b, a0, b0, a1,b1; float alpha=1.0; float K=pi/2; a=width/2; b=height/2; e=(float)width/(float)height; Point2i Center(width/2, height/2); float new_x, new_y; float p,q,x1,y1,x0,y0; float theta; for (int y=0; y<height; y++) { for (int x=0; x<width; x++) { y0=Center.y-y; x0=x-Center.x; theta=atan(y0*e/(x0+0.0001)); if(x0<0) theta=theta+pi; a0=x0/cos(theta); b0=y0/sin(theta+0.0001); if(a0>a || b0>b) continue; a1=a*sin(a0/a*K); b1=b*sin(b0/b*K); a1=(a1-a0)*(alpha)+a0; b1=(b1-b0)*(alpha)+b0; new_x=a1*cos(theta); new_y=b1*sin(theta); new_x=Center.x+new_x; new_y=Center.y-new_y; if(new_x<0) new_x=0; if(new_x>=width-1) new_x=width-2; if(new_y<0) new_y=0; if(new_y>=height-1) new_y=height-2; x1=(int)new_x; y1=(int)new_y; p=new_x-x1; q=new_y-y1; for (int k=0; k<3; k++) { Img_out.at<Vec3b>(y, x)[k]=(1-p)*(1-q)*Img.at<Vec3b>(y1, x1)[k]+ (p)*(1-q)*Img.at<Vec3b>(y1,x1+1)[k]+ (1-p)*(q)*Img.at<Vec3b>(y1+1,x1)[k]+ (p)*(q)*Img.at<Vec3b>(y1+1,x1+1)[k]; } } } Show_Image(Img_out, "out"); imwrite("Out.jpg", Img_out); waitKey(); } // define the show image #include "PS_Algorithm.h" #include <iostream> #include <string> using namespace std; using namespace cv; void Show_Image(Mat& Image, const string& str) { namedWindow(str.c_str(),CV_WINDOW_AUTOSIZE); imshow(str.c_str(), Image); }
原图
效果图
相关文章推荐
- OpenCV——PS滤镜算法之 Ellipsoid (凸出)
- 欢迎使用CSDN-markdown编辑器
- 转:反向代理服务器的工作原理
- Linux折腾记
- 通过strace 监控 fdatasync
- 用OpenLayers解析GML文件
- Expression Add Operators
- 一步步带你,如何网站架构
- 绝杀!Debug 大法,让运维不再尴尬
- 详解linux运维工程师入门级必备技能
- uva 861 Little Bishops 组合数学 (布棋)
- Hadoop伪分布式集群的搭建
- Centos6.7 升级Glib2.14后问题
- 浅析Linux开发工具之Makefile
- linux 新增硬盘分区格式化
- 《三分热度》:mac中意外的Operation not permitted
- Nginx反向代理服务器(博客架构优化)
- Linux下通过SSH远程连接BAE
- Linux文件的三种时间属性
- Linux文件操作