fork download
copy 
  1. #include<vector>
  2. #include<opencv\cv.h>
  3. #include<opencv\highgui.h>
  4. #include<opencv2/imgproc/imgproc.hpp>
  5. #include<cstdlib>
  6. #include<iostream>
  7. #include<cmath>
  8. #include<omp.h>
  9.  
  10. using namespace cv;
  11. using namespace std;
  12.  
  13. Scalar meanpixel(Mat matrx)
  14. {
  15.     Scalar meanpix;
  16.  
  17. 	#pragma omp parallel for
  18. 	for(int i=0;i<matrx.rows;i++)
  19. 	{
  20. 		for(int j=0;j<matrx.cols;j++)
  21. 		{
  22. 			meanpix[0]=meanpix[0]+double(matrx.at<Vec3b>(i,j)[0]);
  23. 			meanpix[1]=meanpix[1]+double(matrx.at<Vec3b>(i,j)[1]);
  24. 			meanpix[2]=meanpix[2]+double(matrx.at<Vec3b>(i,j)[2]);
  25. 		}
  26. 	}
  27. 	meanpix[0]=floor((meanpix[0]/(matrx.rows*matrx.cols))+0.5);
  28. 	meanpix[1]=floor((meanpix[1]/(matrx.rows*matrx.cols))+0.5);
  29. 	meanpix[2]=floor((meanpix[2]/(matrx.rows*matrx.cols))+0.5);
  30. 	return meanpix;
  31. }
  32.  
  33. Mat covarmatrixcal(Mat matrx1, Scalar meanscal)
  34. {
  35.  
  36. 	double rr=0,rg=0,rb=0,gr=0,gg=0,gb=0,br=0,bg=0,bb=0;
  37. 	int k=0;
  38. 	k=(matrx1.rows*matrx1.cols);
  39. 	Vec3b temp1;
  40. 	#pragma omp parallel for
  41. 	for(int i = 0; i < matrx1.rows;i++)
  42. 	{
  43. 		for(int j=0;j<matrx1.cols;j++)
  44. 		{
  45. 			temp1=matrx1.at<Vec3b>(i,j);
  46. 			rr=rr+((double(temp1[0])-meanscal[0])*(double(temp1[0])-meanscal[0]));
  47. 			rg=rg+((double(temp1[0])-meanscal[0])*(double(temp1[1])-meanscal[1]));
  48. 			rb=rb+((double(temp1[0])-meanscal[0])*(double(temp1[2])-meanscal[2]));
  49. 			gg=gg+((double(temp1[1])-meanscal[1])*(double(temp1[1])-meanscal[1]));
  50. 			gb=gb+((double(temp1[1])-meanscal[1])*(double(temp1[2])-meanscal[2]));
  51. 			bb=bb+((double(temp1[2])-meanscal[2])*(double(temp1[2])-meanscal[2]));
  52. 		}
  53. 	}
  54. 	rr=rr/(k-1);
  55. 	rg=rg/(k-1);
  56. 	rb=rb/(k-1);
  57. 	gr=rg;
  58. 	gg=gg/(k-1);
  59. 	gb=gb/(k-1);
  60. 	br=rb;
  61. 	bg=gb;
  62. 	bb=bb/(k-1);
  63.  
  64. 	double m[3][3]={{rr,rg,rb},{gr,gg,gb},{br,bg,bb}};
  65.  
  66. 	Mat covmatrix(3,3,CV_64F,m);
  67. 	return covmatrix.clone();
  68.  
  69. }
  70.  
  71. double mahadistance(Mat invcov, Scalar meanvec, Vec3b patternvec)
  72. {
  73. 	double distmeasure=0;
  74. 	Scalar x1,x2;
  75.  
  76. 	//mahalanobis distance is equal to (x-mean)^T*inv(cov)*(x-mean)
  77. 	x1[0]=(patternvec[0]-meanvec[0]);
  78. 	x1[1]=(patternvec[1]-meanvec[1]);
  79. 	x1[2]=(patternvec[2]-meanvec[2]);
  80.  
  81. 	#pragma omp parallel for
  82. 	for(int k = 0; k < 3;k++)
  83. 	{
  84. 		x2[0]=x2[0]+(x1[k]*invcov.at<double>(k,0));
  85. 		x2[1]=x2[1]+(x1[k]*invcov.at<double>(k,1));  
  86. 		x2[2]=x2[2]+(x1[k]*invcov.at<double>(k,2));
  87. 	}
  88.  
  89. 	distmeasure=((x2[0]*x1[0])+(x2[1]*x1[1])+(x2[2]*x1[2]));
  90. 	return distmeasure;
  91. }
  92.  
  93. Mat inversemat(Mat matrx)
  94. {
  95. 	double determ;
  96. 	determ = (matrx.at<double>(0,0)*((matrx.at<double>(1,1)*matrx.at<double>(2,2))-(matrx.at<double>(1,2)*matrx.at<double>(2,1))))-(matrx.at<double>(0,1)*((matrx.at<double>(1,0)*matrx.at<double>(2,2))-(matrx.at<double>(2,0)*matrx.at<double>(1,2))))+(matrx.at<double>(0,2)*((matrx.at<double>(1,0)*matrx.at<double>(2,1))-(matrx.at<double>(1,1)*matrx.at<double>(2,0))));
  97. 	Mat mt=matrx.t();
  98. 	double m11,m12,m13,m21,m22,m23,m31,m32,m33;
  99. 	m11=(mt.at<double>(1,1)*mt.at<double>(2,2))-(mt.at<double>(2,1)*mt.at<double>(1,2));
  100. 	m12=(mt.at<double>(1,0)*mt.at<double>(2,2))-(mt.at<double>(2,0)*mt.at<double>(1,2));
  101. 	m13=(mt.at<double>(1,0)*mt.at<double>(2,1))-(mt.at<double>(2,0)*mt.at<double>(1,1));
  102. 	m21=(mt.at<double>(0,1)*mt.at<double>(2,2))-(mt.at<double>(2,1)*mt.at<double>(0,2));
  103. 	m22=(mt.at<double>(0,0)*mt.at<double>(2,2))-(mt.at<double>(2,0)*mt.at<double>(0,2));
  104. 	m23=(mt.at<double>(0,0)*mt.at<double>(2,1))-(mt.at<double>(2,0)*mt.at<double>(0,1));
  105. 	m31=(mt.at<double>(0,1)*mt.at<double>(1,2))-(mt.at<double>(1,1)*mt.at<double>(0,2));
  106. 	m32=(mt.at<double>(0,0)*mt.at<double>(1,2))-(mt.at<double>(1,0)*mt.at<double>(0,2));
  107. 	m33=(mt.at<double>(0,0)*mt.at<double>(1,1))-(mt.at<double>(1,0)*mt.at<double>(0,1));
  108.  
  109. 	double imat[3][3]={{(m11/determ),(-m12/determ),(m13/determ)},{(-m21/determ),(m22/determ),(-m23/determ)},{(m31/determ),(-m32/determ),(m33/determ)}};
  110. 	Mat invmat(3,3,CV_64F,imat);
  111. 	return invmat.clone();	//as only the address to where the data is stored is passed to the function the object is destroyed so an empty matrix is returned with initial values -9.23842e+61
  112. }
  113.  
  114. Mat imagesegmentation(Mat image, Mat icovar, Scalar meanmat,double distance, Mat mask, Mat element)
  115. {	
  116.    #pragma omp parallel for
  117.     for(int i = 0; i < image.rows;i++)
  118.     {
  119.         for(int j=0; j<image.cols;j++)	//in case it fails changed it from i=1 to i=0
  120. 		{
  121. 			Vec3b pixel= image.at<Vec3b>(i,j);	//prints wierd characters
  122. 			double mdist=mahadistance(icovar,meanmat,pixel);
  123. 			if(mdist<distance)
  124. 				mask.at<uchar>(i,j)=255;
  125. 			else
  126. 				mask.at<uchar>(i,j)=0;
  127. 		}
  128.     }
  129. 	Mat mask1=mask;
  130. 	morphologyEx(mask,mask1,2,element);
  131.     return mask;
  132. }
  133.  
  134. Mat Cornerdetect(Mat image, Mat corners,int numcorners)
  135. {
  136. 	Mat grayscale;
  137. 	#pragma omp parallel sections
  138. 	{	
  139. 		#pragma omp section 
  140. 		{ cvtColor(image,grayscale,CV_BGR2GRAY); }
  141. 		//goodfeatures(grayimage, output to store corners, quality factor, distance factor)
  142. 		#pragma omp section 
  143. 		{ goodFeaturesToTrack(grayscale,corners,numcorners,0.01,70); }   
  144. 		// Mark these corners on the original image
  145. 	}
  146.  
  147. 	if(corners.rows!=0)
  148. 	{
  149. 		cornerSubPix(grayscale, corners, Size(11, 11), Size(-1, -1), TermCriteria(CV_TERMCRIT_EPS | CV_TERMCRIT_ITER, 30, 0.1)); 
  150. 	}
  151. 	return corners; 
  152. }
  153.  
Compilation error #stdin compilation error #stdout 0s 0KB
comments (?)
stdin
copy 
Standard input is empty
compilation info 
prog.cpp:2:22: fatal error: opencv\cv.h: No such file or directory
 #include<opencv\cv.h>
                      ^
compilation terminated.
stdout
copy 
Standard output is empty
https://ideone.com/cyJSOo
language:
C++ (gcc 8.3)
created:
10 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!