學習OpenCV：濾鏡系列（11）——高反差保留 （6.30修改版）

查志強發表於2014-11-25
OpenCV

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

==============================================

版權所有：小熊不去實驗室CSDN部落格

==============================================


高反差保留就是高通濾波

r=(pix[x,y]-avg(R))/128

pix[x,y]*r+128*(1-r)

[cpp] view plaincopyprint?
  1. #include <math.h>  
  2. #include <opencv/cv.h>  
  3. #include <opencv/highgui.h>  
  4.   
  5. using namespace cv;  
  6. using namespace std;  
  7.   
  8. int R=5;  
  9.   
  10. int main()  
  11. {  
  12.     Mat src = imread("D:/10.jpg",1);  
  13.     int width=src.cols;  
  14.     int heigh=src.rows;  
  15.     Mat img;  
  16.     src.copyTo(img);  
  17.     Mat avg;  
  18.     //GaussianBlur(img,avg,Size(R,R),0.0);  
  19.     blur(img,avg,Size(R,R));  
  20.     Mat dst(img.size(),CV_8UC3);  
  21.     float tmp;  
  22.     for (int y=0;y<heigh;y++)  
  23.     {  
  24.         uchar* imgP=img.ptr<uchar>(y);  
  25.         uchar* avgP=avg.ptr<uchar>(y);  
  26.         uchar* dstP=dst.ptr<uchar>(y);  
  27.         for (int x=0;x<width;x++)  
  28.         {  
  29.             float r0 = abs((float)imgP[3*x]-(float)avgP[3*x])/128;    
  30.             tmp = abs( ((float)imgP[3*x]  )*r0 + 128*(1-r0) );  
  31.             tmp=tmp>255?255:tmp;  
  32.             tmp=tmp<0?0:tmp;  
  33.             dstP[3*x]=(uchar)(tmp);  
  34.   
  35.             float r1 = abs((float)imgP[3*x+1]-(float)avgP[3*x+1])/128;  
  36.             tmp = (uchar)abs( ((float)imgP[3*x+1])*r1 + 128*(1-r1) );  
  37.             tmp=tmp>255?255:tmp;  
  38.             tmp=tmp<0?0:tmp;  
  39.             dstP[3*x+1]=(uchar)(tmp);  
  40.               
  41.             float r2 = abs((float)imgP[3*x+2]-(float)avgP[3*x+2])/128;  
  42.             tmp = (uchar)abs( ((float)imgP[3*x+2])*r2 + 128*(1-r2) );  
  43.             tmp=tmp>255?255:tmp;  
  44.             tmp=tmp<0?0:tmp;  
  45.             dstP[3*x+2]=(uchar)(tmp);  
  46.         }  
  47.     }  
  48.     imshow("high",dst);  
  49.       
  50.     //高通濾波測試  
  51.   
  52.     Mat kern = (Mat_<char>(3,3) <<  -1, -1,  -1,  
  53.                                     -1,  5, -1,  
  54.                                     -1, -1,  -1);  
  55.     Mat dstF;  
  56.     filter2D(img,dstF,img.depth(),kern);  
  57.     imshow("kernel",dstF);  
  58.   
  59.     waitKey();  
  60.     imwrite("D:/高反差保留.jpg",dst);  
  61.     imwrite("D:/高通濾波.jpg",dstF);  
  62. }  


原圖：



高反差保留：



高通濾波器：




修改版 高反差保留 對比PS的效果，更加接近。

    r    = (pix[x,y]-avg(R))

dst[x,y] = 128+|r|*r/(2*R)


[cpp] view plaincopyprint?
  1. #include <math.h>  
  2. #include <opencv/cv.h>  
  3. #include <opencv/highgui.h>  
  4.   
  5. using namespace cv;  
  6. using namespace std;  
  7.   
  8. int R=11;  
  9.   
  10. int main()  
  11. {  
  12.     Mat src = imread("D:/img/liushishi02.jpg",1);  
  13.     imshow("src",src);  
  14.     int width=src.cols;  
  15.     int heigh=src.rows;  
  16.     Mat img;  
  17.     src.copyTo(img);  
  18.     Mat avg;  
  19.     //GaussianBlur(img,avg,Size(R,R),0.0);  
  20.     blur(img,avg,Size(R,R));  
  21.     Mat dst(img.size(),CV_8UC3);  
  22.     float tmp;  
  23.     for (int y=0;y<heigh;y++)  
  24.     {  
  25.         uchar* imgP=img.ptr<uchar>(y);  
  26.         uchar* avgP=avg.ptr<uchar>(y);  
  27.         uchar* dstP=dst.ptr<uchar>(y);  
  28.         for (int x=0;x<width;x++)  
  29.         {  
  30.   
  31.             float r0 = ((float)imgP[3*x]-(float)avgP[3*x]);    
  32.             tmp = 128+abs(r0)*r0/(2*R);  
  33.             tmp=tmp>255?255:tmp;  
  34.             tmp=tmp<0?0:tmp;  
  35.             dstP[3*x]=(uchar)(tmp);  
  36.   
  37.             float r1 = ((float)imgP[3*x+1]-(float)avgP[3*x+1]);  
  38.             tmp = 128+abs(r1)*r1/(2*R);  
  39.             tmp=tmp>255?255:tmp;  
  40.             tmp=tmp<0?0:tmp;  
  41.             dstP[3*x+1]=(uchar)(tmp);  
  42.               
  43.             float r2 = ((float)imgP[3*x+2]-(float)avgP[3*x+2]);  
  44.             tmp = 128+abs(r2)*r2/(2*R);  
  45.             tmp=tmp>255?255:tmp;  
  46.             tmp=tmp<0?0:tmp;  
  47.             dstP[3*x+2]=(uchar)(tmp);  
  48.         }  
  49.     }  
  50.     imshow("high",dst);  
  51.       
  52.     //高通濾波測試  
  53.   
  54.     Mat kern = (Mat_<char>(3,3) <<  0, -1,  0,  
  55.                                     -1,  5, -1,  
  56.                                     0, -1,  0);  
  57.     Mat dstF;  
  58.     filter2D(img,dstF,img.depth(),kern);  
  59.     imshow("kernel",dstF);  
  60.   
  61.     waitKey();  
  62.     //imwrite("D:/高反差保留.jpg",dst);  
  63. //  imwrite("D:/高通濾波.jpg",dstF);  
  64. }  

相關文章