c++编写的自相关超分辨率重建算法

#include <opencv2\highgui\highgui.hpp> #include <opencv2\imgproc\imgproc.hpp> #include <opencv2\core\core.hpp> #include <iostream> #include <stdlib.h> #include <vector> using namespace std; using namespace cv; int main(int argc, char* argv[]) { Mat test=Mat::eye(16,12,CV_32F); Mat t,d; Mat testout=Mat::ones(12,16,CV_32F); int i,j; Scalar sum1; int dirScale; double xIn,w,yOut; //cout<<test<<endl; //copyMakeBorder(test,test,2,2,2,2,BORDER_REPLICATE); //sum1=sum(test-20); //testout(Rect(0,0,3,3)).copyTo(test(Rect(0,0,3,3))); //cout<<test<<endl; //test(Rect(0,0,3,3))=test(Rect(0,0,3,3))+1; Mat filt=Mat::ones(1,5,CV_32F); //cout<<testout<<endl; Mat tm(32,36,CV_32F,Scalar(0)); copyMakeBorder(testout,testout,10,10,10,10,BORDER_CONSTANT); //cout<<testout<<endl; //cout<<tm<<endl; //cout<<filt<<endl; for(i=0;i<testout.rows;i++) { filter2D(testout.row(i),tm.row(i),CV_32F,filt); } cout<<tm<<endl; Mat out=tm(Rect(10,10,tm.cols-20,tm.rows-20)); cout<<out<<endl; //vector<double> vec[9]; //double syn5[4][11]= {-0.0280753,-0.0528051, +0.0609643, +0.9250912, +0.3045216, +0.0075195, +0.0146154, +0.0000000, +0.0000000, +0.0000000, +0.0000000, // +0.0000000, +0.0004889, +0.0382805, -0.0861971, +0.8620317, +0.5205558, -0.1288526, +0.0618612, +0.0000000, +0.0000000, +0.0000000, // +0.0000000, +0.0000000, +0.0000000, +0.0618612, -0.1288526, +0.5205558, +0.8620317, -0.0861971, +0.0382805, +0.0004889, +0.0000000, // +0.0000000, +0.0000000, +0.0000000, +0.0000000, +0.0146154, +0.0075195, +0.3045216, +0.9250912, +0.0609643, -0.0528051, -0.0280753}; //double syn4[3][10]={-0.0420048, -0.0414416, +0.0696145, +0.9278507, +0.3652888, -0.0085849, +0.0300705, +0.0000000 +0.0000000, +0.0000000, // +0.0000000, +0.0000000, +0.0293302, -0.0568660, +0.7267425, +0.7267425, -0.0568660, +0.0293302 +0.0000000, +0.0000000, // +0.0000000, +0.0000000, +0.0000000, +0.0300705, -0.0085849, +0.3652888, +0.9278507, +0.0696145 -0.0414416, -0.0420048}; //double syn3[2][7]={1,2,3,4,5,6,7, // -0.1005021, +0.1197507, +0.9271653, +0.6005021, -0.0469160,+0.0000000, +0.0000000}; //for(i=0;i<4;i++) //{ // for(j=0;j<11;j++) // { // vec[i+5].push_back(syn5[i][j]); // } //} //for(i=0;i<3;i++) //{ // for(j=0;j<11;j++) // { // vec[i+2].push_back(syn4[i][j]); // } //} //for(i=0;i<2;i++) //{ // for(j=0;j<7;j++) // { // vec[i].push_back(syn3[i][j]); // } //} //xIn=0.7; //w=abs(xIn); // if(w<1&&w>=0) //{ // yOut=1-2*w*w+w*w*w; //} // else if( w>=1&&w<2) //{ // yOut=4-8*w+5*w*w-w*w*w; //} // else yOut=0; //cout<<yOut<<endl; //xIn=1.3; //w=abs(xIn); // if(w<1&&w>=0) //{ // yOut=1-2*w*w+w*w*w; //} // else if( w>=1&&w<2) //{ // yOut=4-8*w+5*w*w-w*w*w; //} // else yOut=0; //cout<<yOut<<endl; //xIn=0; //w=abs(xIn); // if(w<1&&w>=0) //{ // yOut=1-2*w*w+w*w*w; //} // else if( w>=1&&w<2) //{ // yOut=4-8*w+5*w*w-w*w*w; //} // else yOut=0; //cout<<yOut<<endl; //test.col(5)=test.col(5)+2; //test.col(7)=test.col(7)+3; //test.row(2)=test.row(2)+4; //test.row(1)=test.row(1)+1; //cout<<test<<endl; ////test.row(0).copyTo(testout.col(1)); //transpose(test,test); //transpose(test,test); //cout<<test<<endl; //copyMakeBorder(test,t,2,2,2,2,BORDER_CONSTANT,0); //cout<<t<<endl; //d=t(Rect(2,2,test.cols,test.rows)); //cout<<d<<endl; /*test.row(0).copyTo(t); cout<<t<<endl; cout<<vec[1][0]<<vec[1][1]<<vec[1][2]<<vec[1][3]<<vec[1][4]<<vec[1][5]<<vec[1][6]<<endl; filter2D(t,d,CV_32F,vec[1]); cout<<d<<endl;*/ //t=Mat::eye(test.Size,CV_32F); //cout<<t<<endl; /*int a=3; if(a==5) { double syn[4][11]= {-0.0280753,-0.0528051, +0.0609643, +0.9250912, +0.3045216, +0.0075195, +0.0146154, +0.0000000, +0.0000000, +0.0000000, +0.0000000, +0.0000000, +0.0004889, +0.0382805, -0.0861971, +0.8620317, +0.5205558, -0.1288526, +0.0618612, +0.0000000, +0.0000000, +0.0000000, +0.0000000, +0.0000000, +0.0000000, +0.0618612, -0.1288526, +0.5205558, +0.8620317, -0.0861971, +0.0382805, +0.0004889, +0.0000000, +0.0000000, +0.0000000, +0.0000000, +0.0000000, +0.0146154, +0.0075195, +0.3045216, +0.9250912, +0.0609643, -0.0528051, -0.0280753}; } else if(a==4) { double syn[3][10]={-0.0420048, -0.0414416, +0.0696145, +0.9278507, +0.3652888, -0.0085849, +0.0300705, +0.0000000 +0.0000000, +0.0000000, +0.0000000, +0.0000000, +0.0293302, -0.0568660, +0.7267425, +0.7267425, -0.0568660, +0.0293302 +0.0000000, +0.0000000, +0.0000000, +0.0000000, +0.0000000, +0.0300705, -0.0085849, +0.3652888, +0.9278507, +0.0696145 -0.0414416, -0.0420048}; } else if(a==3) { double syn[2][7]={-0.1005021, +0.1197507, +0.9271653, +0.6005021, -0.0469160, +0.0000000, +0.0000000, +0.0000000, +0.0000000, -0.0469160, +0.6005021, +0.9271653, +0.1197507, -0.1005021}; } cout<<syn[1][6]<<endl; */ //value[0]=1; //value[1]=2; //value[2]=3; //bgr[0].create(10,10,CV_8UC1); //bgr[1].create(10,10,CV_8UC1); //bgr[2].create(10,10,CV_8UC1); ////i=bgr[0].rows; ////j=bgr[0].cols; ////cout<<i<<j<<endl; //for (i=0;i<5;i++) // for (j=0;j<5;j++) // { // bgr[0].at<uchar>(i,j)=i; // bgr[1].at<uchar>(i,j)=i; // bgr[2].at<uchar>(i,j)=i; // } // bgr[0]=bgr[0]*5+15; //merge(bgr,3,src); // cout<<"out"<<src<<endl; // // //namedWindow("test",WINDOW_AUTOSIZE); //while(true) //{ // c=waitKey(); // if(c=='a') // { // copyMakeBorder(src,dst,2,2,2,2, BORDER_REFLECT101); // cout<<dst<<endl; // } // else if(c=='b') // { // copyMakeBorder(src,dst,2,2,2,2,BORDER_CONSTANT,value); // cout<<dst<<endl; // } // else if(c=='d') // { // break; // } // //} ///*src.convertTo(src,CV_32F); //src=src/2; //cout<<"out"<<src<<endl; //split(src,bgr); //cout<<"out"<<bgr[0]<<endl;*/ // ////Mat inp=Mat::eye (10,10,CV_8UC1); ////Mat out(19,15,CV_8UC1,0); ////out=inp(Rect(1,1,4,8)); ////for(i=0;i<8;i++) // //for(j=0;j<4;j++) // //{ // //out.at<uchar>(i,j)=i+j; // //} ////cout<<"out"<<src(Rect(0,0,2,2)); getchar(); return 0; }