Halcon学习之模板匹配

Halcon提供了许多的匹配方式：

匹配算法流程:
1.抠图，获取模板区域；
2.创建模板；
3.查找与模板相似区域，获得区域的坐标，角度，比例大小等结果；
4.清除模板

在实际应用中根据选择相应发匹配方式即可；大多数情况下使用基于相关性和基于形状这两种方式。

一、基于相关性
根据灰度做归一化的互相关匹配。可用于线性光照变换、边缘轮廓模糊的情况下。
基本流程：
*1采集图像
read_image (Image, ImageFiles[4])

*2准备模板
gen_circle (ROI_0, 249.543, 341.002, 165.223)
reduce_domain (Image, ROI_0, ImageReduced)
*3创建模板
create_ncc_model (ImageReduced, 4, 0,rad(180), ‘auto’, ‘use_polarity’, ModelID)
*获得模板区域
get_ncc_model_region (ModelRegion, ModelID)//该区域会被仿射变换到原点
get_ncc_model_params (ModelID, NumLevels, AngleStart, AngleExtent, AngleStep, Metric)
*4查找模板
find_ncc_model (Image, ModelID, 0, rad(270), 0.5, 1, 0.5, ‘true’, 0, Row, Column, Angle, Score)
*5显示
vector_angle_to_rigid (0, 0, 0, Row, Column, Angle, HomMat2D)
affine_trans_region (ModelRegion, RegionAffineTrans, HomMat2D, ‘nearest_neighbor’)
dev_display (Image)
dev_display (RegionAffineTrans)
*6清除模板
clear_ncc_model (ModelID)

重要算子参数说明：

create_ncc_model(Template : : NumLevels,AngleStart, AngleExtent, AngleStep, Metric : ModelID)
参数1：Template ，模板图像
参数2：NumLevels，金字塔层数（0，或’auto’都为自动确认）
参数3：AngleStart，起始角度（单位rad）
参数4：Template ，角度范围（单位rad）
参数5：Template ，角度步长（单位rad）
参数6：Metric ，极性；'use_polarity’为使用极性； 'ignore_global_polarity’为忽略极性
参数7：ModelID，模板句柄

find_ncc_model(Image : : ModelID, AngleStart, AngleExtent, MinScore, NumMatches, MaxOverlap, SubPixel, NumLevels : Row, Column, Angle, Score)
create_ncc_model(Template : : NumLevels,AngleStart, AngleExtent, AngleStep, Metric : ModelID)
参数1：Image ，图像
参数2：ModelID，句柄
参数3：AngleStart，起始角度（单位rad）
参数4：Template ，角度范围（单位rad）
参数5：MinScore，最小分数（0-1）
参数6：NumMatches，最大匹配数（0为All）
参数7：MaxOverlap，最大重叠率
参数8：SubPixel，亚像素精度
参数9：NumLevels ，金字塔层数
输出：
参数10：Row，匹配到区域中心的行
参数11：Column，匹配到区域中心的列
参数12：Angle，模板旋转到匹配到区域的角度（rad）
参数13：Score，得分（0-1）

二、基于形状
根据边缘方向梯度进行匹配，最常用的方式，可应用与边缘轮廓清晰的情况下
基本流程
*1采集图像
read_image (Image, ImageFiles[0])
*2准备模板
gen_rectangle1 (ROI_0, 189.637, 164.533, 309.563, 435.045)
reduce_domain (Image, ROI_0, ImageReduced)
*3创建模板
create_shape_model (ImageReduced, ‘auto’, 0, rad(360), 0.01, ‘auto’, ‘use_polarity’, ‘auto’, ‘auto’, ModelID)
*4查找
find_shape_model (Image, ModelID, 0, rad(360), 0.7, 1, 0.5, ‘none’, 0, 0.5, Row, Column, Angle, Score)
*5显示
if(|Row| > 0)
vector_angle_to_rigid (0, 0, 0, Row, Column, Angle, HomMat2D)
* affine_trans_region (ROI_0, RegionAffineTrans, HomMat2D, ‘nearest_neighbor’)
affine_trans_contour_xld (ModelContours, ContoursAffineTrans, HomMat2D)
dev_display (Image)
dev_display (ContoursAffineTrans)
disp_message (WindowHandle, Score, ‘window’, 12, 12, ‘black’, ‘true’)
Endif
*6清除
clear_shape_model (ModelID)
算法详解：
创建模板算子：
create_shape_model(Template : : NumLevels, AngleStart, AngleExtent, AngleStep, Optimization, Metric, Contrast, MinContrast : ModelID)
参数1：模板图像
参数2：金字塔最大层数 ‘0’默认，相当与‘auto’
参数3：模板的起始角度 rad
参数4：模板的角度范围 rad
参数5：角度步长 rad
参数6：优化方式，与速度有关
‘none’：不做处理；
‘no_pregeneration’：不预先存储到内存，会慢
‘pregeneration’：预先存储到内存，会快
‘point_reduction_high’, ‘point_reduction_low’, ‘point_reduction_medium’,；减少点数
参数7：极性
‘use_polarity’：使用极性；‘ignore_global_polarity’：忽略全局极性；
‘ignore_local_polarity’：忽略局部极性；'ignore_color_polarity’忽略颜色极性
参数8：模板图像的前景与背景的对比度
参数9：搜索到图像的最小对比度
参数10：模板句柄

寻找模板算子：
find_shape_model(Image : : ModelID, AngleStart, AngleExtent, MinScore, NumMatches, MaxOverlap, SubPixel, NumLevels, Greediness : Row, Column, Angle, Score)
参数1：搜索图像（可以多通道）
参数2：模型句柄
参数3：模板的起始角度
参数4：模板的角度范围；起始角度和角度范围需比创建模板时的小，不然则采用创建时的参数。
参数5：最小分数（模板与实例的相似度）
参数6：最大的匹配实例的数量；‘0’为全找
参数7：最大重叠度
参数8：亚像素精度
与插值算法相关：’
none’：不采用亚像素； ‘interpolation’,：简单插值
采用最小二乘法： ‘least_squares’, ‘least_squares_high’, ‘least_squares_very_high’
与变形相关：
‘max_deformation 1’, ‘max_deformation 2’, ‘max_deformation 3’,最大变形度为32；
使用该参数时，必须要亚像素；与插值搭配使用：【‘least_squares’,‘max_deformation 5’】
参数9：过程实际采用的金字塔层数
5：从第5层从上往下搜
【5，2】从第5层一直到第2层搜，找到最佳的匹配层数
【5，-2】从第5层搜索到第2层，返回第2层；如果最低层没有结果；采用容差模式（图像离焦时）。写‘0’表示与创建模板的参数一致。
参数10：贪婪度（0-1），0：慢但安全；1，快但可能找不到
参数11：被找到实例中心的行，
参数12：被找到实例中心的列，
参数13：模板旋转到实例的角度
参数14：得分

相关算子：
*检查模板
inspect_shape_model (ImageReduced, ModelImages, ModelRegions, 4, 30)
*获取模板的轮廓
get_shape_model_contours (ModelContours, ModelID, 1)
获取模板的参数
get_shape_model_params (ModelID, NumLevels, AngleStart, AngleExtent, AngleStep, ScaleMin, ScaleMax, ScaleStep, Metric, MinContrast)
*设置边界是否查找
set_system (‘border_shape_models’, ‘true’)