怎么快速求图像的旋转角度

如何快速求图像的旋转角度 在线等
在平面上有很多相同的图(不重叠),但是旋转角度都不同,如何以其中一个图作为标准,算出其他凸显的旋转角度呢?
如图:怎么快速求图像的旋转角度
以左下角的图为参照,求其他图的旋转角度,图上标出的角度是我用软件刻意旋转的,目标就是要求出这些角度,允许误差可以较大。
------解决方案--------------------
求最小外接矩形,再求矩形的旋转角度
------解决方案--------------------
怎么快速求图像的旋转角度
代码有的地方,还是稍微有点问题:::




// FindRotation-angle.cpp : 定义控制台应用程序的入口点。
//

// findContours.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"



#include <iostream>
#include <vector>
#include <opencv2/opencv.hpp> 
#include <opencv2/core/core.hpp>
#include <opencv2/imgproc/imgproc.hpp>
#include <opencv2/highgui/highgui.hpp>
//#include "highlight"
//#include "highgui.h"



#define PI 3.1415926

int main()
{
// Read input binary image
cv::Mat image = cv::imread("test.jpg",0);
if (!image.data)
return 0; 

cv::namedWindow("Binary Image");
cv::imshow("Binary Image",image);



// 从文件中加载原图  
   IplImage *pSrcImage = cvLoadImage("test.jpg", CV_LOAD_IMAGE_UNCHANGED);  
  
   // 转为2值图

 cvThreshold(pSrcImage,pSrcImage,200,255,cv::THRESH_BINARY_INV);
   

   image = cv::Mat(pSrcImage,true);

   cv::imwrite("binary.jpg",image);

// Get the contours of the connected components
std::vector<std::vector<cv::Point>> contours;
cv::findContours(image, 
contours, // a vector of contours 
CV_RETR_EXTERNAL, // retrieve the external contours
CV_CHAIN_APPROX_NONE); // retrieve all pixels of each contours

// Print contours' length
std::cout << "Contours: " << contours.size() << std::endl;
std::vector<std::vector<cv::Point>>::const_iterator itContours= contours.begin();
for ( ; itContours!=contours.end(); ++itContours) 
{

std::cout << "Size: " << itContours->size() << std::endl;
}

// draw black contours on white image
cv::Mat result(image.size(),CV_8U,cv::Scalar(255));
cv::drawContours(result,contours,
-1, // draw all contours
cv::Scalar(0), // in black
2); // with a thickness of 2

cv::namedWindow("Contours");
cv::imshow("Contours",result);









// Eliminate too short or too long contours
int cmin= 100;  // minimum contour length
int cmax= 1000; // maximum contour length
std::vector<std::vector<cv::Point>>::const_iterator itc= contours.begin();
while (itc!=contours.end()) {

if (itc->size() < cmin 
------解决方案--------------------
 itc->size() > cmax)
itc= contours.erase(itc);
else 
++itc;
}

// draw contours on the original image
cv::Mat original= cv::imread("test.jpg");
cv::drawContours(original,contours,
-1, // draw all contours
cv::Scalar(255,255,0), // in white
2); // with a thickness of 2

cv::namedWindow("Contours on Animals");
cv::imshow("Contours on Animals",original);



// Let's now draw black contours on white image
result.setTo(cv::Scalar(255));
cv::drawContours(result,contours,
-1, // draw all contours
cv::Scalar(0), // in black
1); // with a thickness of 1
image= cv::imread("binary.jpg",0);

// testing the bounding box 
cv::Rect r0= cv::boundingRect(cv::Mat(contours[0]));
cv::rectangle(result,r0,cv::Scalar(0),2);



// testing the enclosing circle 
float radius;
cv::Point2f center;
cv::minEnclosingCircle(cv::Mat(contours[1]),center,radius);
cv::circle(result,cv::Point(center),static_cast<int>(radius),cv::Scalar(0),2);

// cv::RotatedRect rrect= cv::fitEllipse(cv::Mat(contours[1]));
// cv::ellipse(result,rrect,cv::Scalar(0),2);

// testing the approximate polygon
std::vector<cv::Point> poly;
cv::approxPolyDP(cv::Mat(contours[2]),poly,5,true);

std::cout << "Polygon size: " << poly.size() << std::endl;

// Iterate over each segment and draw it
std::vector<cv::Point>::const_iterator itp= poly.begin();
while (itp!=(poly.end()-1)) {
cv::line(result,*itp,*(itp+1),cv::Scalar(0),2);