本文实例为大家分享了Opencv透视变换综合实例的具体代码，供大家参考，具体内容如下
案例背景：对下面发生畸变的图像进行校正
方案思路：灰度二值化分割，闭操作，寻找轮廓，霍夫直线检测，直线排序，直线方程，直线交点，透视矩阵，透视变换。
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
|
#include<opencv2\opencv.hpp>using namespace cv;using namespace std;int main(int arc, char** argv) { Mat src = imread("1.jpg"); namedWindow("input", CV_WINDOW_AUTOSIZE); imshow("input", src); //灰度化 Mat grayImg; cvtColor(src, grayImg, CV_BGR2GRAY); //二值化 Mat binaryImg; threshold(grayImg, binaryImg, 0, 255, THRESH_BINARY_INV | THRESH_OTSU); //闭操作 Mat kernel = getStructuringElement(MORPH_RECT,Size(3,3)); morphologyEx(binaryImg, binaryImg, MORPH_CLOSE,kernel,Point(-1,-1) ,3); imshow("output", binaryImg); //寻找轮廓 Mat draw = Mat::zeros(src.size(), CV_8UC3); vector<vector<Point>>contours; findContours(binaryImg, contours, RETR_EXTERNAL, CHAIN_APPROX_SIMPLE, Point()); for (int i = 0; i < contours.size(); i++) { Rect rect = boundingRect(contours[i]); if (rect.width < src.cols / 2 && rect.height < src.rows / 2)continue; drawContours(draw, contours, i, Scalar(0, 0, 255), 2); } imshow("output2", draw); //霍夫直线检测 vector<Vec4i> lines; cvtColor(draw, draw, CV_BGR2GRAY); HoughLinesP(draw, lines, 1, CV_PI / 180, src.rows/2,src.rows/2,0); Mat draw2 = Mat::zeros(src.size(), CV_8UC3); for (int j = 0; j < lines.size(); j++) { Vec4i ln = lines[j]; line(draw2, Point(ln[0], ln[1]), Point(ln[2], ln[3]), Scalar(0, 255, 0), 2); } printf("number of line:%d\n", lines.size()); imshow("output3", draw2); //寻找与定位直线 Vec4i topLine,bottomLine,leftLine,rightLine; for (int j = 0; j < lines.size(); j++) { Vec4i ln = lines[j]; if (ln[1] < src.rows / 2 && ln[3] < src.rows / 2) { topLine = ln; } if (ln[1] > src.rows / 2 && ln[3] > src.rows / 2) { bottomLine = ln; } if (ln[0] < src.cols / 2 && ln[2] < src.cols / 2) { leftLine = ln; } if (ln[0] > src.cols / 2 && ln[2] > src.cols / 2) { rightLine = ln; } } cout << "topLine:p1(x,y)=" << topLine[0] << "," << topLine[1]<<" " << "p2(x,y)=" << topLine[2] << "," << topLine[3] << endl; cout << "bottomLine:p1(x,y)=" << bottomLine[0] << "," << bottomLine[1] << " " << "p2(x,y)=" << bottomLine[2] << "," << bottomLine[3] << endl; cout << "leftLine:p1(x,y)=" << leftLine[0] << "," << leftLine[1] << " " << "p2(x,y)=" << leftLine[2] << "," << leftLine[3] << endl; cout << "rightLine:p1(x,y)=" << rightLine[0] << "," << rightLine[1] << " " << "p2(x,y)=" << rightLine[2] << "," << rightLine[3] << endl; //求解直线方程 float k1, c1; k1 = float((topLine[3] - topLine[1])) / float(topLine[2] - topLine[0]); c1 = topLine[1] - k1*topLine[0]; float k2, c2; k2 = float((bottomLine[3] - bottomLine[1])) / float(bottomLine[2] - bottomLine[0]); c2 = bottomLine[1] - k2*bottomLine[0]; float k3, c3; k3 = float((leftLine[3] - leftLine[1])) / float(leftLine[2] - leftLine[0]); c3 = leftLine[1] - k3*leftLine[0]; float k4, c4; k4 = float((rightLine[3] - rightLine[1])) / float(rightLine[2] - rightLine[0]); c4 = rightLine[1] - k4*rightLine[0]; //求解直线交点 Point p1; p1.x = (int)((c1 - c3) / (k3 - k1)); p1.y = (int)(k1*p1.x + c1); Point p2; p2.x = (int)((c1 - c4) / (k4 - k1)); p2.y = (int)(k1*p2.x + c1); Point p3; p3.x = (int)((c2 - c4) / (k4 - k2)); p3.y = (int)(k2*p3.x + c2); Point p4; p4.x = (int)((c2 - c3) / (k3 - k2)); p4.y = (int)(k2*p4.x + c2); cout << "左上角：" << p1.x << "," << p1.y << endl; cout << "右上角：" << p2.x << "," << p2.y << endl; cout << "右下角：" << p3.x << "," << p3.y << endl; cout << "左下角：" << p4.x << "," << p4.y << endl; //画出交点 circle(draw2, p1, 2, Scalar(0, 0, 255)); circle(draw2, p2, 2, Scalar(0, 0, 255)); circle(draw2, p3, 2, Scalar(0, 0, 255)); circle(draw2, p4, 2, Scalar(0, 0, 255)); imshow("output4", draw2); //透视变换 vector<Point2f> srcCorners(4); srcCorners[0] = p1; srcCorners[1] = p2; srcCorners[2] = p3; srcCorners[3] = p4; vector<Point2f> dstCorners(4); dstCorners[0] = Point(0, 0); dstCorners[1] = Point(src.cols, 0); dstCorners[2] = Point(src.cols, src.rows); dstCorners[3] = Point(0, src.rows); Mat warpMartrix = getPerspectiveTransform(srcCorners, dstCorners);//Mat warpMartrix = findHomography(srcCorners, dstCorners); Mat result = Mat::zeros(src.size(), -1); warpPerspective(src, result, warpMartrix, result.size(),INTER_LINEAR); imshow("output5", result); waitKey(0); return 0; } |
原图像
二值化闭操作
寻找轮廓
霍夫直线
直线及其交点
效果图
以上就是本文的全部内容，希望对大家的学习有所帮助，也希望大家多多支持服务器之家。
原文链接：https://blog.csdn.net/qq_24946843/article/details/82772654
