【halcon每日一练】

测量仿射跟随

程序功能：测量不同视角下，芯片引脚长度，输出相应测量信息。
思路：

1. 创建ROI模板，计算仿射变换矩阵
2. 利用点仿射创建测量矩形
3. 输出计算结果

halcon代码

*创建形状模板，实现轮廓对象显示

read_image (Image, 'C:/Users/njtl003/Desktop/board-01.png')
*创建模板ROI区域
gen_rectangle2 (Rectangle, 243, 297, 0, 100, 40)
area_center (Rectangle, Area, Row, Column)
*创建测量ROI区域-设定偏移值
RowOffset:=100
ColOffset:=5
Rec2Phi:=0
Rec2Len1:=170
Rec2Len2:=5
gen_rectangle2 (Rectangle2, Row-RowOffset, Column+ColOffset, Rec2Phi, Rec2Len1, Rec2Len2)
gen_rectangle2 (Rectangle3, Row+RowOffset, Column+ColOffset, Rec2Phi, Rec2Len1, Rec2Len2)
*抠图创建模板匹配模板
reduce_domain (Image, Rectangle, ImageReduced)
create_shape_model (ImageReduced, 'auto', 0, 6.28, 'auto', 'auto', 'use_polarity', \
                    'auto', 'auto', ModelID)
get_shape_model_contours (shapeModelContours, ModelID, 1)
find_shape_model(Image,ModelID,0,6.28,0.5,1,0.5,'least_squares',\
                 0,0.9,Row1, Column1, Angle, Score)
*创建仿射变换矩阵
hom_mat2d_identity (HomMat2DIdentity)
hom_mat2d_scale (HomMat2DIdentity, 1, 1, 0, 0, HomMat2DScale)
hom_mat2d_rotate (HomMat2DScale, Angle, 0, 0, HomMat2DRotate)
hom_mat2d_translate (HomMat2DRotate, Row1, Column1, HomMat2DTranslate)

affine_trans_contour_xld (shapeModelContours, ContoursAffineTrans, HomMat2DTranslate)

*显示结果
dev_display (Image)
dev_get_window (WindowHandle)
dev_set_color ('green')
dev_display (shapeModelContours)
dev_set_color ('red')
dev_set_draw ('margin')
dev_set_line_width (3)
dev_display (Rectangle2)
dev_display (Rectangle3)
dev_set_color ('yellow')
disp_message (WindowHandle, 'Next Test the Picture', 'window', 12, 12, 'black', 'true')
disp_continue_message (WindowHandle, 'black', 'true')
stop()

*载入测试图片；实现形状匹配；

* Image Acquisition 01: Code generated by Image Acquisition 01
list_files ('C:/Users/njtl003/Desktop/11111', ['files','follow_links'], ImageFiles)
tuple_regexp_select (ImageFiles, ['\\.(tif|tiff|gif|bmp|jpg|jpeg|jp2|png|pcx|pgm|ppm|pbm|xwd|ima|hobj)$','ignore_case'], ImageFiles)
for Index := 0 to |ImageFiles| - 1 by 1
    dev_set_window (WindowHandle)
    read_image (Image, ImageFiles[Index])
    get_image_size (Image, Width, Height)
    dev_set_part (0, 0, Height-1, Width-1)          
    dev_display (Image)
    find_shape_model (Image, ModelID, 0, 6.28, 0.5, 1, 0.5, 'least_squares', 0, 0.9, shapeCenterRow, shapeCenterColumn, shapeAngle, shapeScore)

    *如果找到模板则显示模板结果和仿射测量矩形
    if(|shapeScore|>0)
        dev_set_color ('red')
        vector_angle_to_rigid (0, 0, 0, shapeCenterRow, shapeCenterColumn, shapeAngle, HomMat2D)
        affine_trans_contour_xld (shapeModelContours, ContoursAffineTrans1, HomMat2D)
        dev_display (ContoursAffineTrans1)
     *对测量矩形中心点进行仿射变换
     affine_trans_pixel (HomMat2D, -RowOffset, ColOffset, RowTrans1, ColTrans1)
     affine_trans_pixel (HomMat2D, RowOffset, ColOffset, RowTrans2, ColTrans2)
     *显示仿射变换后的rec2矩形
     gen_rectangle2 (Rectangle1Affine, RowTrans1, ColTrans1, shapeAngle, Rec2Len1, Rec2Len2)
     gen_rectangle2 (Rectangle2Affine, RowTrans2, ColTrans2, shapeAngle, Rec2Len1, Rec2Len2)

     dev_set_color ('blue')
     dev_display (Rectangle1Affine)
     dev_display (Rectangle2Affine)

*1D测量实现；结果显示

     *创建测量矩形
     
     gen_measure_rectangle2 (RowTrans1, ColTrans1, shapeAngle, Rec2Len1, Rec2Len2, Width, Height, 'nearest_neighbor', MeasureHandle1)
     gen_measure_rectangle2 (RowTrans2, ColTrans2, shapeAngle, Rec2Len1, Rec2Len2, Width, Height, 'nearest_neighbor', MeasureHandle2)

     *测量矩形1应用-可通过测量助手调试后填入相关参数
     measure_pairs (Image, MeasureHandle1, 2, 90, 'positive', \
                    'all', RowEdgeFirst1, ColumnEdgeFirst1, AmplitudeFirst1,\
                    RowEdgeSecond1, ColumnEdgeSecond1, AmplitudeSecond1, \
                    IntraDistance1, InterDistance1)
     
     *测量矩形2应用
     measure_pairs (Image, MeasureHandle2, 2, 90, 'positive', \
                    'all', RowEdgeFirst2, ColumnEdgeFirst2, AmplitudeFirst2,\
                    RowEdgeSecond2, ColumnEdgeSecond2, AmplitudeSecond2, \
                    IntraDistance2, InterDistance2)
     
     close_measure (MeasureHandle1)
     close_measure (MeasureHandle2)

     *显示第一个测量矩形左边缘
     disp_line(WindowHandle, RowEdgeFirst1-Rec2Len2*cos(shapeAngle), ColumnEdgeFirst1-Rec2Len2*sin(shapeAngle), RowEdgeFirst1+Rec2Len2*cos(shapeAngle), ColumnEdgeFirst1+Rec2Len2*sin(shapeAngle))
     *显示第一个测量矩形右边缘
     disp_line(WindowHandle, RowEdgeSecond1-Rec2Len2*cos(shapeAngle), ColumnEdgeSecond1-Rec2Len2*sin(shapeAngle), RowEdgeSecond1+Rec2Len2*cos(shapeAngle), ColumnEdgeSecond1+Rec2Len2*sin(shapeAngle))
     *显示第二个测量矩形左边缘
     disp_line(WindowHandle, RowEdgeFirst2-Rec2Len2*cos(shapeAngle), ColumnEdgeFirst2-Rec2Len2*sin(shapeAngle), RowEdgeFirst2+Rec2Len2*cos(shapeAngle), ColumnEdgeFirst2+Rec2Len2*sin(shapeAngle))
     *显示第二个测量矩形右边缘
     disp_line(WindowHandle, RowEdgeSecond2-Rec2Len2*cos(shapeAngle), ColumnEdgeSecond2-Rec2Len2*sin(shapeAngle), RowEdgeSecond2+Rec2Len2*cos(shapeAngle), ColumnEdgeSecond2+Rec2Len2*sin(shapeAngle))
     NumLead:=|IntraDistance1|+|IntraDistance2|
     MeanDistance:=mean([InterDistance1,InterDistance2])
     MinDistance:=min([InterDistance1,InterDistance2])
     dev_clear_window ()
     *显示共测得多少边缘对
     *每个边缘对的距离
     *芯片引脚平均距离
     *芯片引脚最小距离
     dev_display (Image)
     disp_message(WindowHandle,'共测得边缘对:'+NumLead$'2d','image',20,20,'green','true')
     disp_message(WindowHandle,'平均距离:'+MeanDistance$'3.2f','image',50,20,'red','true')
     disp_message(WindowHandle,'最小距离:'+MinDistance$'3.2f','image',80,20,'black','true')
    endif    
    stop ()
endfor
    *关闭形状匹配句柄
    clear_shape_model (ModelID)

结果显示