快捷導航

Opencv實現(xiàn)最小外接矩形和圓

更新時間：2020年07月21日 09:22:21 作者：東城青年

這篇文章主要為大家詳細介紹了Opencv實現(xiàn)最小外接矩形和圓，文中示例代碼介紹的非常詳細，具有一定的參考價值，感興趣的小伙伴們可以參考一下

本文實例為大家分享了Opencv實現(xiàn)最小外接矩形和圓的具體代碼，供大家參考，具體內(nèi)容如下

步驟：將一幅圖像先轉(zhuǎn)灰度，再canny邊緣檢測得到二值化邊緣圖像，再尋找輪廓，輪廓是由一系列點構(gòu)成的，要想獲得輪廓的最小外接矩形，首先需要得到輪廓的近似多邊形，用道格拉斯-普克抽稀（DP）算法，道格拉斯-普克抽稀算法，是將曲線近似表示為一系列點，并減少點的數(shù)量的一種算法。

該算法實現(xiàn)抽稀的過程是：

1）對曲線的首末點虛連一條直線，求曲線上所有點與直線的距離，并找出最大距離值dmax，用dmax與事先給定的閾值D相比：
2）若dmax<D，則將這條曲線上的中間點全部舍去；則該直線段作為曲線的近似，該段曲線處理完畢。

若dmax≥D，保留dmax對應(yīng)的坐標點，并以該點為界，把曲線分為兩部分，對這兩部分重復使用該方法，即重復1），2）步，直到所有dmax均<D，即完成對曲線的抽稀。

#include<opencv2/opencv.hpp>
using namespace cv;
using namespace std;
 
int value = 60;
RNG rng(1);
Mat src,gray_img,canny_img,dst;
void callback(int, void*);
int main(int arc, char** argv){ 
 src = imread("2.jpg"); 
 namedWindow("input",CV_WINDOW_AUTOSIZE);
 imshow("input", src);
 cvtColor(src, gray_img, CV_BGR2GRAY);
 namedWindow("output", CV_WINDOW_AUTOSIZE);
 createTrackbar("threshold", "output", &value, 255, callback);
 callback(0, 0);
 waitKey(0);
 return 0;
}
void callback(int, void*) { 
 Canny(gray_img, canny_img, value, 2 * value);
 
 vector<vector<Point>>contours;
 vector<Vec4i> hierarchy;
 findContours(canny_img, contours, hierarchy, RETR_EXTERNAL, CHAIN_APPROX_SIMPLE, Point(0, 0));
 
 vector<vector<Point>> contours_poly(contours.size());
 vector<Rect>poly_rects(contours.size());
 vector<Point2f>ccs(contours.size());
 vector<float>radius(contours.size());
 
 vector<RotatedRect> minRects(contours.size());
 vector<RotatedRect> myellipse(contours.size());
 
 for (int i = 0; i < contours.size(); i++) {
 approxPolyDP(contours[i], contours_poly[i], 20, true);//獲得點數(shù)比較少的近似多邊形
 poly_rects[i] = boundingRect(contours_poly[i]);//從近似多邊形獲得最小外接矩形
 minEnclosingCircle(contours_poly[i], ccs[i], radius[i]);//從近似多邊形獲得最小外接圓
 //多邊形點數(shù)大于5才能繪制帶方向的最小矩形和橢圓
 if (contours_poly[i].size() > 5) {
 minRects[i] = minAreaRect(contours_poly[i]);//從近似多邊形獲得帶方向的最小外接矩形
 myellipse[i] = fitEllipse(contours_poly[i]);//從近似多邊形獲得帶方向的最小外接橢圓
 } 
 }
 
 //繪制
 src.copyTo(dst);
 Point2f pts[4];
 for (int j = 0; j < contours.size(); j++) {
 Scalar color = Scalar(rng.uniform(0, 255), rng.uniform(0, 255), rng.uniform(0, 255));
 rectangle(dst, poly_rects[j], color, 2,8);
 circle(dst, ccs[j], (int)radius[j], color, 2,8);
 
 //繪制帶方向的最小外接矩形和橢圓
 if (contours_poly[j].size() > 5) {
 ellipse(dst, myellipse[j], color, 2);
 minRects[j].points(pts);
 for (int k = 0; k < 4; k++) {
 line(dst, pts[k], pts[(k + 1)%4], color, 2);
 }
 }
 }
 imshow("output", dst);
 
}

運行結(jié)果如下：