本博客算法及代碼參考自賈志剛老師的《OpenCV圖像處理-小案例實戰》,若涉及侵權問題,望通知,會第一時間刪除。
算法功能:
1.圖像角度傾斜矯正 (基於仿射變換)
2.去掉多餘的邊(輪廓查找+ROI提取)
原始圖像如下:
算法思路:
一、進行圖像角度糾正
二、取出ROI區域,去掉多餘的白邊
代碼實現:
/*
=======圖像旋轉+切邊=======
*/
#include <opencv2/opencv.hpp>
#include <iostream>
#include <math.h>
using namespace cv;
using namespace std;
Mat Check_Skew(Mat&);
void FindROI(Mat&);
int threshold_value = 100;
int max_level = 255;
const char* output_win = "Contours Result";
const char* roi_win = "Final Result";
int main(int argc, char** argv) {
Mat src = imread("D:/VS2015_Projects/opencv_workspace/img/img_skew.jpg");
if (src.empty()) {
printf("could not load image...\n");
return -1;
}
namedWindow("input image", CV_WINDOW_AUTOSIZE);
imshow("input image", src);
namedWindow(output_win, CV_WINDOW_AUTOSIZE);
Mat img_skew = Check_Skew(src);
// namedWindow(roi_win, CV_WINDOW_AUTOSIZE);
//createTrackbar("Threshold:", output_win, &threshold_value, max_level, FindROI);
FindROI(img_skew);
waitKey(0);
return 0;
}
//角度矯正
Mat Check_Skew(Mat& src) {
Mat gray_src,canny_output;
cvtColor(src, gray_src, COLOR_BGR2GRAY);
//邊緣檢測
Canny(gray_src, canny_output, threshold_value, threshold_value * 2, 3, false);
//輪廓查找
vector<vector<Point>> contours;
vector<Vec4i> hireachy;
findContours(canny_output, contours, hireachy, RETR_TREE, CHAIN_APPROX_SIMPLE, Point(0, 0));
Mat drawImg = Mat::zeros(src.size(), CV_8UC3);
float maxw = 0;
float maxh = 0;
double degree = 0;
//角度獲取
for (size_t t = 0; t < contours.size(); t++) {
RotatedRect minRect = minAreaRect(contours[t]);
degree = abs(minRect.angle);
if (degree > 0) {
maxw = max(maxw, minRect.size.width);
maxh = max(maxh, minRect.size.height);
}
}
//輪廓繪製
RNG rng(12345);
for (size_t t = 0; t < contours.size(); t++) {
RotatedRect minRect = minAreaRect(contours[t]);
if (maxw == minRect.size.width && maxh == minRect.size.height) {
degree = minRect.angle;
Point2f pts[4];
minRect.points(pts);
Scalar color = Scalar(rng.uniform(0, 255), rng.uniform(0, 255), rng.uniform(0, 255));
for (int i = 0; i < 4; i++) {
line(drawImg, pts[i], pts[(i + 1) % 4], color, 2, 8, 0);
}
}
}
printf("max contours width : %f\n", maxw);
printf("max contours height : %f\n", maxh);
printf("max contours angle : %f\n", degree);
imshow(output_win, drawImg);
//獲得旋轉矩陣
Point2f center(src.cols / 2, src.rows / 2);
Mat rotm = getRotationMatrix2D(center, degree, 1.0);
//旋轉圖像
Mat dst;
warpAffine(src, dst, rotm, src.size(), INTER_LINEAR, 0, Scalar(255, 255, 255));
//顯示結果
imshow("Correct Image", dst);
return dst;
}
//去邊
void FindROI(Mat& img) {
Mat gray_src;
cvtColor(img, gray_src, COLOR_BGR2GRAY);
//邊緣檢測
Mat canny_output;
Canny(gray_src, canny_output, threshold_value, threshold_value * 2, 3, false);
//輪廓查找
vector<vector<Point>> contours;
vector<Vec4i> hireachy;
findContours(canny_output, contours, hireachy, RETR_TREE, CHAIN_APPROX_SIMPLE, Point(0, 0));
//繪製輪廓
int minw = img.cols*0.75;
int minh = img.rows*0.75;
RNG rng(12345);
Mat drawImage = Mat::zeros(img.size(), CV_8UC3);
Rect bbox;
for (size_t t = 0; t < contours.size(); t++) {
//查找可傾斜的最小外接矩
RotatedRect minRect = minAreaRect(contours[t]);
//獲得傾斜角度
float degree = abs(minRect.angle);
if (minRect.size.width > minw && minRect.size.height > minh && minRect.size.width < (img.cols - 5)) {
printf("current angle : %f\n", degree);
Point2f pts[4];
minRect.points(pts);
bbox = minRect.boundingRect();
Scalar color = Scalar(rng.uniform(0, 255), rng.uniform(0, 255), rng.uniform(0, 255));
for (int i = 0; i < 4; i++) {
line(drawImage, pts[i], pts[(i + 1) % 4], color, 2, 8, 0);
}
}
}
imshow(output_win, drawImage);
//提取ROI區域
if (bbox.width > 0 && bbox.height > 0) {
Mat roiImg = img(bbox);
imshow(roi_win, roiImg);
}
return;
}
實現效果