OpenCV 4.13.0
开源计算机视觉库 (Open Source Computer Vision)
正在加载...
正在搜索...
未找到匹配项
线特征教程

在本教程中,将展示如何:

  • 使用 BinaryDescriptor 接口提取直线并将其存储在 KeyLine 对象中
  • 使用相同的接口为每个提取的直线计算描述符
  • 使用 BynaryDescriptorMatcher 确定来自不同图像的描述符之间的匹配项

直线提取和描述符计算

在以下代码片段中,展示了如何从图像中检测直线。LSD 提取器使用 LSD_REFINE_ADV 选项进行初始化;其余参数保留其默认值。使用全 1 掩码以接受所有提取的直线,最后,对第 0 层金字塔(octave 0)提取的直线使用随机颜色进行显示。

1/*M///////////////////////////////////////////////////////////////////////////////////////
2 //
3 // 重要提示:在下载、复制、安装或使用之前请阅读。
4 //
5 // 下载、复制、安装或使用本软件即表示您同意本许可协议。
6 // 如果您不同意本许可协议,请勿下载、安装、
7 // 复制或使用本软件。
8 //
9 //
10 // 许可协议
11 // 用于 OpenCV (Open Source Computer Vision Library)
12 //
13 // Copyright (C) 2014, Biagio Montesano, 保留所有权利。
14 // 第三方版权归其各自所有者所有。
15 //
16 // 在满足以下条件的情况下,允许以源码和二进制形式重新分发和使用,无论是否经过修改:
17 // 需满足以下条件:
18 //
19 // * 源码的再分发必须保留上述版权声明、
20 // 此条件列表以及以下免责声明。
21 //
22 // * 二进制形式的再分发必须在分发包提供的文档
23 // 和/或其他材料中保留上述版权声明、
24 // 此条件列表以及以下免责声明。
25 //
26 // * 未经事先具体书面许可,不得使用版权持有人的名称来认可或推广
27 // 衍生自本软件的产品。
28 //
29 // 本软件由版权持有人及贡献者“按原样”提供,不提供任何明示或暗示的保证,包括但不限于
30 // 对适销性和特定用途适用性的暗示保证。
31 // 在任何情况下,Intel 公司或贡献者均不对任何直接、
32 // 间接、附带、特殊、典型或后果性损害
33 // (包括但不限于替代商品或服务的采购;
34 // 使用、数据或利润损失;或业务中断)承担责任,无论其是如何造成的,
35 // 且基于任何责任理论,无论是合同、严格责任、
36 // 或侵权行为(包括疏忽或其他原因),即使已被告知
37 // 存在此类损害的可能性。
38 //M*/
39 //
40 //M*/
41
42#include <iostream>
43#include <opencv2/opencv_modules.hpp>
44
45#ifdef HAVE_OPENCV_FEATURES2D
46
49#include <opencv2/imgproc.hpp>
51#include <opencv2/highgui.hpp>
52
53using namespace cv;
54using namespace cv::line_descriptor;
55using namespace std;
56
57static const char* keys =
58{ "{@image_path | | 图像路径 }" };
59
60static void help()
61{
62 cout << "\n本示例展示了 BinaryDescriptor 类提供的直线提取功能\n"
63 << "请使用如下格式的命令运行本示例:\n" << "./example_line_descriptor_lines_extraction <输入图像路径>" << endl;
64}
65
66int main( int argc, char** argv )
67{
68 /* 从命令行获取参数 */
69 CommandLineParser parser( argc, argv, keys );
70 String image_path = parser.get<String>( 0 );
71
72 if( image_path.empty() )
73 {
74 help();
75 return -1;
76 }
77
78 /* 加载图像 */
79 cv::Mat imageMat = imread( image_path, 1 );
80 if( imageMat.data == NULL )
81 {
82 std::cout << "错误,无法加载图像。请检查其路径" << std::endl;
83 return -1;
84 }
85
86 /* 创建一个全 1 二进制掩码 */
87 cv::Mat mask = Mat::ones( imageMat.size(), CV_8UC1 );
88
89 /* 使用默认参数创建 LSDDetector 对象的指针 */
90 Ptr<LSDDetector> bd = LSDDetector::createLSDDetector();
91
92 /* 创建一个结构来存储提取的直线 */
93 vector<KeyLine> lines;
94
95 /* 提取直线 */
96 cv::Mat output = imageMat.clone();
97 bd->detect( imageMat, lines, 2, 1, mask );
98
99 /* 绘制从 octave 0 提取的直线 */
100 if( output.channels() == 1 )
101 cvtColor( output, output, COLOR_GRAY2BGR );
102 for ( size_t i = 0; i < lines.size(); i++ )
103 {
104 KeyLine kl = lines[i];
105 if( kl.octave == 0)
106 {
107 /* 获取随机颜色 */
108 int R = ( rand() % (int) ( 255 + 1 ) );
109 int G = ( rand() % (int) ( 255 + 1 ) );
110 int B = ( rand() % (int) ( 255 + 1 ) );
111
112 /* 获取直线端点 */
113 Point pt1 = Point2f( kl.startPointX, kl.startPointY );
114 Point pt2 = Point2f( kl.endPointX, kl.endPointY );
115
116 /* 绘制直线 */
117 line( output, pt1, pt2, Scalar( B, G, R ), 3 );
118 }
119
120 }
121
122 /* 在图像上显示直线 */
123 imshow( "LSD lines", output );
124 waitKey();
125}
126
127#else
128
129int main()
130{
131 std::cerr << "OpenCV 编译时未包含 features2d 模块" << std::endl;
132 return 0;
133}
134
135#endif // HAVE_OPENCV_FEATURES2D
专为命令行解析设计。
定义 utility.hpp:890
n 维密集数组类
定义于 mat.hpp:840
CV_NODISCARD_STD Mat clone() const
创建数组及其底层数据的完整副本。
MatSize size
定义 mat.hpp:2226
uchar * data
指向数据的指针
定义 mat.hpp:2206
int channels() const
返回矩阵通道数。
std::string String
定义 cvstd.hpp:151
std::shared_ptr< _Tp > Ptr
定义 cvstd_wrapper.hpp:23
#define CV_8UC1
定义 interface.h:88
GMat mask(const GMat &src, const GMat &mask)
将掩码应用于矩阵。
void imshow(const String &winname, InputArray mat)
在指定窗口中显示图像。
int waitKey(int delay=0)
等待按键操作。
Mat imread(const String &filename, int flags=IMREAD_COLOR_BGR)
从文件加载图像。
void cvtColor(InputArray src, OutputArray dst, int code, int dstCn=0, AlgorithmHint hint=cv::ALGO_HINT_DEFAULT)
将图像从一种颜色空间转换为另一种。
void line(InputOutputArray img, Point pt1, Point pt2, const Scalar &color, int thickness=1, int lineType=LINE_8, int shift=0)
绘制连接两个点的线段。
int main(int argc, char *argv[])
定义 highgui_qt.cpp:3
定义 descriptor.hpp:77
定义 core.hpp:107
STL 命名空间。
表示直线的类。
定义 descriptor.hpp:105
float endPointY
定义 descriptor.hpp:131
float startPointX
定义 descriptor.hpp:128
float endPointX
定义 descriptor.hpp:130
float startPointY
定义 descriptor.hpp:129
int octave
定义 descriptor.hpp:114

这是从著名的摄影师(cameraman)图像中获得的结果

替换文本

另一种提取直线的方法是使用 LSDDetector 类;该类使用 LSD 提取器来计算直线。要获得此结果,只需修改上述代码片段中的以下几行:

// 创建 LSDDetector 对象的指针
Ptr<LSDDetector> lsd = LSDDetector::createLSDDetector();
// 计算直线
std::vector<KeyLine> keylines;
lsd->detect( imageMat, keylines, mask );

这是 LSD 检测器在摄影师图片上再次返回的结果

替换文本

检测到关键线(keylines)后,可以按如下所示计算它们的描述符

1/*M///////////////////////////////////////////////////////////////////////////////////////
2 //
3 // 重要提示:在下载、复制、安装或使用之前请阅读。
4 //
5 // 下载、复制、安装或使用本软件即表示您同意本许可协议。
6 // 如果您不同意本许可协议,请勿下载、安装、
7 // 复制或使用本软件。
8 //
9 //
10 // 许可协议
11 // 用于 OpenCV (Open Source Computer Vision Library)
12 //
13 // Copyright (C) 2014, Biagio Montesano, 保留所有权利。
14 // 第三方版权归其各自所有者所有。
15 //
16 // 在满足以下条件的情况下,允许以源码和二进制形式重新分发和使用,无论是否经过修改:
17 // 需满足以下条件:
18 //
19 // * 源码的再分发必须保留上述版权声明、
20 // 此条件列表以及以下免责声明。
21 //
22 // * 二进制形式的再分发必须在分发包提供的文档
23 // 和/或其他材料中保留上述版权声明、
24 // 此条件列表以及以下免责声明。
25 //
26 // * 未经事先具体书面许可,不得使用版权持有人的名称来认可或推广
27 // 衍生自本软件的产品。
28 //
29 // 本软件由版权持有人及贡献者“按原样”提供,不提供任何明示或暗示的保证,包括但不限于
30 // 对适销性和特定用途适用性的暗示保证。
31 // 在任何情况下,Intel 公司或贡献者均不对任何直接、
32 // 间接、附带、特殊、典型或后果性损害
33 // (包括但不限于替代商品或服务的采购;
34 // 使用、数据或利润损失;或业务中断)承担责任,无论其是如何造成的,
35 // 且基于任何责任理论,无论是合同、严格责任、
36 // 或侵权行为(包括疏忽或其他原因),即使已被告知
37 // 存在此类损害的可能性。
38 //M*/
39 //
40 //M*/
41
42#include <iostream>
43#include <opencv2/opencv_modules.hpp>
44
45#ifdef HAVE_OPENCV_FEATURES2D
46
49#include <opencv2/imgproc.hpp>
51#include <opencv2/highgui.hpp>
52
53using namespace cv;
54using namespace cv::line_descriptor;
55
56
57static const char* keys =
58{ "{@image_path | | 图像路径 }" };
59
60static void help()
61{
62 std::cout << "\n本示例展示了 BinaryDescriptor 类提供的直线提取和描述符计算功能\n"
63 << "请使用如下格式的命令运行本示例:\n" << "./example_line_descriptor_compute_descriptors <输入图像路径>"
64 << std::endl;
65}
66
67int main( int argc, char** argv )
68{
69 /* 从命令行获取参数 */
70 CommandLineParser parser( argc, argv, keys );
71 String image_path = parser.get<String>( 0 );
72
73 if( image_path.empty() )
74 {
75 help();
76 return -1;
77 }
78
79 /* 加载图像 */
80 cv::Mat imageMat = imread( image_path, 1 );
81 if( imageMat.data == NULL )
82 {
83 std::cout << "错误,无法加载图像。请检查其路径" << std::endl;
84 }
85
86 /* 创建二进制掩码 */
87 cv::Mat mask = Mat::ones( imageMat.size(), CV_8UC1 );
88
89 /* 使用默认参数创建 BinaryDescriptor 对象的指针 */
90 Ptr<BinaryDescriptor> bd = BinaryDescriptor::createBinaryDescriptor();
91
92 /* 计算直线 */
93 std::vector<KeyLine> keylines;
94 bd->detect( imageMat, keylines, mask );
95
96 /* 计算描述符 */
97 cv::Mat descriptors;
98
99 bd->compute( imageMat, keylines, descriptors);
100
101}
102
103#else
104
105int main()
106{
107 std::cerr << "OpenCV 编译时未包含 features2d 模块" << std::endl;
108 return 0;
109}
110
111#endif // HAVE_OPENCV_FEATURES2D

描述符之间的匹配

如果我们从两幅不同的图像中提取了描述符,就可以在它们之间搜索匹配项。一种方法是将每个输入查询描述符与其最接近距离的描述符精确匹配:

1/*M///////////////////////////////////////////////////////////////////////////////////////
2 //
3 // 重要提示:在下载、复制、安装或使用之前请阅读。
4 //
5 // 下载、复制、安装或使用本软件即表示您同意本许可协议。
6 // 如果您不同意本许可协议,请勿下载、安装、
7 // 复制或使用本软件。
8 //
9 //
10 // 许可协议
11 // 用于 OpenCV (Open Source Computer Vision Library)
12 //
13 // Copyright (C) 2014, Biagio Montesano, 保留所有权利。
14 // 第三方版权归其各自所有者所有。
15 //
16 // 在满足以下条件的情况下,允许以源码和二进制形式重新分发和使用,无论是否经过修改:
17 // 需满足以下条件:
18 //
19 // * 源码的再分发必须保留上述版权声明、
20 // 此条件列表以及以下免责声明。
21 //
22 // * 二进制形式的再分发必须在分发包提供的文档
23 // 和/或其他材料中保留上述版权声明、
24 // 此条件列表以及以下免责声明。
25 //
26 // * 未经事先具体书面许可,不得使用版权持有人的名称来认可或推广
27 // 衍生自本软件的产品。
28 //
29 // 本软件由版权持有人及贡献者“按原样”提供,不提供任何明示或暗示的保证,包括但不限于
30 // 对适销性和特定用途适用性的暗示保证。
31 // 在任何情况下,Intel 公司或贡献者均不对任何直接、
32 // 间接、附带、特殊、典型或后果性损害
33 // (包括但不限于替代商品或服务的采购;
34 // 使用、数据或利润损失;或业务中断)承担责任,无论其是如何造成的,
35 // 且基于任何责任理论,无论是合同、严格责任、
36 // 或侵权行为(包括疏忽或其他原因),即使已被告知
37 // 存在此类损害的可能性。
38 //M*/
39 //
40 //M*/
41
42#include <iostream>
43#include <opencv2/opencv_modules.hpp>
44
45#ifdef HAVE_OPENCV_FEATURES2D
46
49#include <opencv2/imgproc.hpp>
51#include <opencv2/highgui.hpp>
52
53#define MATCHES_DIST_THRESHOLD 25
54
55using namespace cv;
56using namespace cv::line_descriptor;
57
58static const char* keys =
59{ "{@image_path1 | | 图像路径 1 }"
60 "{@image_path2 | | 图像路径 2 }" };
61
62static void help()
63{
64 std::cout << "\n本示例展示了 BinaryDescriptor 类提供的直线提取和描述符计算功能\n"
65 << "请使用如下格式的命令运行本示例:\n" << "./example_line_descriptor_compute_descriptors <输入图像 1 路径>"
66 << "<输入图像 2 路径>" << std::endl;
67
68}
69
70int main( int argc, char** argv )
71{
72 /* 从命令行获取参数 */
73 CommandLineParser parser( argc, argv, keys );
74 String image_path1 = parser.get<String>( 0 );
75 String image_path2 = parser.get<String>( 1 );
76
77 if( image_path1.empty() || image_path2.empty() )
78 {
79 help();
80 return -1;
81 }
82
83 /* 加载图像 */
84 cv::Mat imageMat1 = imread( image_path1, 1 );
85 cv::Mat imageMat2 = imread( image_path2, 1 );
86
87 if( imageMat1.data == NULL || imageMat2.data == NULL )
88 {
89 std::cout << "错误,无法加载图像。请检查其路径" << std::endl;
90 }
91
92 /* 创建二进制掩码 */
93 cv::Mat mask1 = Mat::ones( imageMat1.size(), CV_8UC1 );
94 cv::Mat mask2 = Mat::ones( imageMat2.size(), CV_8UC1 );
95
96 /* 使用默认参数创建 BinaryDescriptor 对象的指针 */
97 Ptr<BinaryDescriptor> bd = BinaryDescriptor::createBinaryDescriptor( );
98
99 /* 计算直线和描述符 */
100 std::vector<KeyLine> keylines1, keylines2;
101 cv::Mat descr1, descr2;
102
103 ( *bd )( imageMat1, mask1, keylines1, descr1, false, false );
104 ( *bd )( imageMat2, mask2, keylines2, descr2, false, false );
105
106 /* 从第一层金字塔中选择关键线及其描述符 */
107 std::vector<KeyLine> lbd_octave1, lbd_octave2;
108 Mat left_lbd, right_lbd;
109 for ( int i = 0; i < (int) keylines1.size(); i++ )
110 {
111 if( keylines1[i].octave == 0 )
112 {
113 lbd_octave1.push_back( keylines1[i] );
114 left_lbd.push_back( descr1.row( i ) );
115 }
116 }
117
118 for ( int j = 0; j < (int) keylines2.size(); j++ )
119 {
120 if( keylines2[j].octave == 0 )
121 {
122 lbd_octave2.push_back( keylines2[j] );
123 right_lbd.push_back( descr2.row( j ) );
124 }
125 }
126
127 /* 创建 BinaryDescriptorMatcher 对象 */
128 Ptr<BinaryDescriptorMatcher> bdm = BinaryDescriptorMatcher::createBinaryDescriptorMatcher();
129
130 /* 进行匹配 */
131 std::vector<DMatch> matches;
132 bdm->match( left_lbd, right_lbd, matches );
133
134 /* 筛选最佳匹配 */
135 std::vector<DMatch> good_matches;
136 for ( int i = 0; i < (int) matches.size(); i++ )
137 {
138 if( matches[i].distance < MATCHES_DIST_THRESHOLD )
139 good_matches.push_back( matches[i] );
140 }
141
142 /* 绘制匹配项 */
143 cv::Mat outImg;
144 cv::Mat scaled1, scaled2;
145 std::vector<char> mask( matches.size(), 1 );
146 drawLineMatches( imageMat1, lbd_octave1, imageMat2, lbd_octave2, good_matches, outImg, Scalar::all( -1 ), Scalar::all( -1 ), mask,
147 DrawLinesMatchesFlags::DEFAULT );
148
149 imshow( "Matches", outImg );
150 waitKey();
151 imwrite("/home/ubisum/Desktop/images/env_match/matches.jpg", outImg);
152 /* 创建一个 LSD 检测器 */
153 Ptr<LSDDetector> lsd = LSDDetector::createLSDDetector();
154
155 /* 检测直线 */
156 std::vector<KeyLine> klsd1, klsd2;
157 Mat lsd_descr1, lsd_descr2;
158 lsd->detect( imageMat1, klsd1, 2, 2, mask1 );
159 lsd->detect( imageMat2, klsd2, 2, 2, mask2 );
160
161 /* 为第一层金字塔的直线计算描述符 */
162 bd->compute( imageMat1, klsd1, lsd_descr1 );
163 bd->compute( imageMat2, klsd2, lsd_descr2 );
164
165 /* 从第一层金字塔中选择直线和描述符 */
166 std::vector<KeyLine> octave0_1, octave0_2;
167 Mat leftDEscr, rightDescr;
168 for ( int i = 0; i < (int) klsd1.size(); i++ )
169 {
170 if( klsd1[i].octave == 1 )
171 {
172 octave0_1.push_back( klsd1[i] );
173 leftDEscr.push_back( lsd_descr1.row( i ) );
174 }
175 }
176
177 for ( int j = 0; j < (int) klsd2.size(); j++ )
178 {
179 if( klsd2[j].octave == 1 )
180 {
181 octave0_2.push_back( klsd2[j] );
182 rightDescr.push_back( lsd_descr2.row( j ) );
183 }
184 }
185
186 /* 计算匹配项 */
187 std::vector<DMatch> lsd_matches;
188 bdm->match( leftDEscr, rightDescr, lsd_matches );
189
190 /* 筛选最佳匹配 */
191 good_matches.clear();
192 for ( int i = 0; i < (int) lsd_matches.size(); i++ )
193 {
194 if( lsd_matches[i].distance < MATCHES_DIST_THRESHOLD )
195 good_matches.push_back( lsd_matches[i] );
196 }
197
198 /* 绘制匹配项 */
199 cv::Mat lsd_outImg;
200 resize( imageMat1, imageMat1, Size( imageMat1.cols / 2, imageMat1.rows / 2 ), 0, 0, INTER_LINEAR_EXACT );
201 resize( imageMat2, imageMat2, Size( imageMat2.cols / 2, imageMat2.rows / 2 ), 0, 0, INTER_LINEAR_EXACT );
202 std::vector<char> lsd_mask( matches.size(), 1 );
203 drawLineMatches( imageMat1, octave0_1, imageMat2, octave0_2, good_matches, lsd_outImg, Scalar::all( -1 ), Scalar::all( -1 ), lsd_mask,
204 DrawLinesMatchesFlags::DEFAULT );
205
206 imshow( "LSD matches", lsd_outImg );
207 waitKey();
208
209
210}
211
212#else
213
214int main()
215{
216 std::cerr << "OpenCV 编译时未包含 features2d 模块" << std::endl;
217 return 0;
218}
219
220#endif // HAVE_OPENCV_FEATURES2D
Mat row(int y) const
为指定的矩阵行创建矩阵头。
void push_back(const _Tp &elem)
向矩阵底部添加元素。
用于指定图像或矩形大小的模板类。
定义 types.hpp:335
bool imwrite(const String &filename, InputArray img, const std::vector< int > &params=std::vector< int >())
将图像保存到指定文件。
void resize(InputArray src, OutputArray dst, Size dsize, double fx=0, double fy=0, int interpolation=INTER_LINEAR)
调整图像大小。
void drawLineMatches(const Mat &img1, const std::vector< KeyLine > &keylines1, const Mat &img2, const std::vector< KeyLine > &keylines2, const std::vector< DMatch > &matches1to2, Mat &outImg, const Scalar &matchColor=Scalar::all(-1), const Scalar &singleLineColor=Scalar::all(-1), const std::vector< char > &matchesMask=std::vector< char >(), int flags=DrawLinesMatchesFlags::DEFAULT)
绘制两幅图像中找到的关键线匹配项。

有时,我们可能对搜索给定输入描述符的最接近的 k 个描述符感兴趣。这需要稍微修改之前的代码:

// 准备用于承载匹配项的结构
std::vector<std::vector<DMatch> > matches;
// 进行 knn 匹配
bdm->knnMatch( descr1, descr2, matches, 6 );

在上述示例中,为每个查询返回最接近的 6 个描述符。在某些情况下,我们可能有一个搜索半径,并寻找与输入查询距离最多为 r 的所有描述符。之前的代码必须修改为:

// 准备用于承载匹配项的结构
std::vector<std::vector<DMatch> > matches;
// 计算匹配项
bdm->radiusMatch( queries, matches, 30 );

这是一个在从原始摄影师图像提取的描述符与其降采样(并模糊)版本之间进行匹配的示例

替换文本

查询内部数据库

BynaryDescriptorMatcher 类拥有一个内部数据库,可以用从不同图像中提取的描述符进行填充,并使用上一节描述的模式之一进行查询。可以使用 add 函数填充内部数据集;该函数不会直接将新数据添加到数据库中,只是将其局部存储。真正的更新发生在调用 train 函数或执行任何查询函数时,因为每个查询函数在查询前都会调用 train。查询时,内部数据库不仅返回所需的描述符,而且对于返回的每个匹配项,它都能指出该匹配描述符是从哪幅图像中提取的。以下代码描述了内部数据集用法的一个示例;在局部添加新描述符后,调用半径搜索。这会触发局部数据传输到数据集,进而对数据集进行查询。

1/*M///////////////////////////////////////////////////////////////////////////////////////
2 //
3 // 重要提示:在下载、复制、安装或使用之前请阅读。
4 //
5 // 下载、复制、安装或使用本软件即表示您同意本许可协议。
6 // 如果您不同意本许可协议,请勿下载、安装、
7 // 复制或使用本软件。
8 //
9 //
10 // 许可协议
11 // 用于 OpenCV (Open Source Computer Vision Library)
12 //
13 // Copyright (C) 2014, Biagio Montesano, 保留所有权利。
14 // 第三方版权归其各自所有者所有。
15 //
16 // 在满足以下条件的情况下,允许以源码和二进制形式重新分发和使用,无论是否经过修改:
17 // 需满足以下条件:
18 //
19 // * 源码的再分发必须保留上述版权声明、
20 // 此条件列表以及以下免责声明。
21 //
22 // * 二进制形式的再分发必须在分发包提供的文档
23 // 和/或其他材料中保留上述版权声明、
24 // 此条件列表以及以下免责声明。
25 //
26 // * 未经事先具体书面许可,不得使用版权持有人的名称来认可或推广
27 // 衍生自本软件的产品。
28 //
29 // 本软件由版权持有人及贡献者“按原样”提供,不提供任何明示或暗示的保证,包括但不限于
30 // 对适销性和特定用途适用性的暗示保证。
31 // 在任何情况下,Intel 公司或贡献者均不对任何直接、
32 // 间接、附带、特殊、典型或后果性损害
33 // (包括但不限于替代商品或服务的采购;
34 // 使用、数据或利润损失;或业务中断)承担责任,无论其是如何造成的,
35 // 且基于任何责任理论,无论是合同、严格责任、
36 // 或侵权行为(包括疏忽或其他原因),即使已被告知
37 // 存在此类损害的可能性。
38 //M*/
39 //
40 //M*/
41
42#include <iostream>
43#include <opencv2/opencv_modules.hpp>
44
45#ifdef HAVE_OPENCV_FEATURES2D
46
49#include <opencv2/imgproc.hpp>
51#include <opencv2/highgui.hpp>
52
53#include <vector>
54
55using namespace cv;
56using namespace cv::line_descriptor;
57
58static const std::string images[] =
59{ "cameraman.jpg", "church.jpg", "church2.png", "einstein.jpg", "stuff.jpg" };
60
61static const char* keys =
62{ "{@image_path | | 图像路径 }" };
63
64static void help()
65{
66 std::cout << "\n本示例展示了半径匹配(radius matching)的功能 " << "请使用如下格式的命令运行本示例:\n"
67 << "./example_line_descriptor_radius_matching <输入图像路径>/" << std::endl;
68}
69
70int main( int argc, char** argv )
71{
72 /* 从命令行获取参数 */
73 CommandLineParser parser( argc, argv, keys );
74 String pathToImages = parser.get < String > ( 0 );
75
76 /* 创建用于存放 KeyLines 和描述符的结构 */
77 int num_elements = sizeof ( images ) / sizeof ( images[0] );
78 std::vector < Mat > descriptorsMat;
79 std::vector < std::vector<KeyLine> > linesMat;
80
81 /* 创建 BinaryDescriptor 对象的指针 */
82 Ptr < BinaryDescriptor > bd = BinaryDescriptor::createBinaryDescriptor();
83
84 /* 计算直线和描述符 */
85 for ( int i = 0; i < num_elements; i++ )
86 {
87 /* 获取图像路径 */
88 std::stringstream image_path;
89 image_path << pathToImages << images[i];
90 std::cout << image_path.str().c_str() << std::endl;
91
92 /* 加载图像 */
93 Mat loadedImage = imread( image_path.str().c_str(), 1 );
94 if( loadedImage.data == NULL )
95 {
96 std::cout << "无法加载图像。" << std::endl;
97 help();
98 exit( -1 );
99 }
100
101 /* 计算直线和描述符 */
102 std::vector < KeyLine > lines;
103 Mat computedDescr;
104 bd->detect( loadedImage, lines );
105 bd->compute( loadedImage, lines, computedDescr );
106
107 descriptorsMat.push_back( computedDescr );
108 linesMat.push_back( lines );
109
110 }
111
112 /* 组成查询矩阵 */
113 Mat queries;
114 for ( size_t j = 0; j < descriptorsMat.size(); j++ )
115 {
116 if( descriptorsMat[j].rows >= 5 )
117 queries.push_back( descriptorsMat[j].rowRange( 0, 5 ) );
118
119 else if( descriptorsMat[j].rows > 0 && descriptorsMat[j].rows < 5 )
120 queries.push_back( descriptorsMat[j] );
121 }
122
123 std::cout << "已生成包含 " << queries.rows << " 个描述符的矩阵" << std::endl;
124
125 /* 创建 BinaryDescriptorMatcher 对象 */
126 Ptr < BinaryDescriptorMatcher > bdm = BinaryDescriptorMatcher::createBinaryDescriptorMatcher();
127
128 /* 填充匹配器 */
129 bdm->add( descriptorsMat );
130
131 /* 计算匹配项 */
132 std::vector < std::vector<DMatch> > matches;
133 bdm->radiusMatch( queries, matches, 30 );
134 std::cout << "匹配项样本大小 " << matches.size() << std::endl;
135
136 for ( int i = 0; i < (int) matches.size(); i++ )
137 {
138 for ( int j = 0; j < (int) matches[i].size(); j++ )
139 {
140 std::cout << "匹配: " << matches[i][j].queryIdx << " " << matches[i][j].trainIdx << " " << matches[i][j].distance << std::endl;
141 }
142
143 }
144
145}
146
147#else
148
149int main()
150{
151 std::cerr << "OpenCV 编译时未包含 features2d 模块" << std::endl;
152 return 0;
153}
154
155#endif // HAVE_OPENCV_FEATURES2D
int rows
行数和列数,当矩阵超过 2 维时为 (-1, -1)
定义 mat.hpp:2204