* ------------------------------------------------------------------------------------------------
* This example program uses shape-based matching to align ROIs for the measure
* tool, which then inspects individual razor blades.
* The program can be run in two modes: (1) with the full affine transformation
* (2) using translate_measure
* Modify the next line to switch between the modes.
USING_TRANSLATE_MEASURE := 1
* ------------------------------------------------------------------------------------------------
* general configuration of HDevelop
dev_update_window ('on')//设置窗口不自动更新
* image acquisition and window size
read_image (ModelImage, 'razors1')//读取一张图像
get_image_pointer1 (ModelImage, Pointer, Type, Width, Height)//获得图像的指针和宽高
dev_close_window ()//关闭窗口
dev_open_window (0, 0, Width, Height, 'white', WindowHandle)//打开窗口
dev_set_part (0, 0, Height - 1, Width - 1)//设置图像在窗口显示的位置
dev_display (ModelImage)//在窗口显示Object变量
* colors and other settings for the visualization
dev_set_color ('cyan')//设置窗口显示内容的颜色
dev_set_draw ('margin')//设置区域是填充还是显示轮廓
dev_set_line_width (2)//设置线宽
stop ()
* ------------------- start of the application ----------------
* -> select the model object
Row1 := 46
Column1 := 57
Row2 := 79
Column2 := 94
gen_rectangle1 (ROIPart1, Row1, Column1, Row2, Column2)//生成矩形区域
gen_rectangle1 (ROIPart2, Row1 + 364, Column1 + 13, Row2 + 364, Column2 + 13)//生成矩形区域
union2 (ROIPart1, ROIPart2, ModelROI)//将两个区域联合起来变成一个区域
area_center (ModelROI, Area, CenterROIRow, CenterROIColumn)//计算区域面积和中心点
dev_display (ModelImage)//显示图像,显示Object变量
dev_display (ModelROI)//显示区域
stop ()
* -> create the model
reduce_domain (ModelImage, ModelROI, ImageROI)//抠图
create_shape_model (ImageROI, 4, 0, 0, 'auto', 'none', 'use_polarity', 30, 10, ModelID)//创建图形模板
inspect_shape_model (ImageROI, ShapeModelImage, ShapeModelRegion, 1, 30)//
get_shape_model_contours (ShapeModel, ModelID, 1)//获得模板轮廓
dev_clear_window ()//清空窗口
dev_set_color ('blue')//设置窗口显示内容的颜色
dev_display (ShapeModelRegion)//显示区域
stop ()
* step 1: create variables describing the measurement ROIs and display them
Rect1Row := 244
Rect1Col := 73
DistColRect1Rect2 := 17
Rect2Row := Rect1Row
Rect2Col := Rect1Col + DistColRect1Rect2
RectPhi := rad(90)
RectLength1 := 122
RectLength2 := 2
gen_rectangle2 (MeasureROI1, Rect1Row, Rect1Col, RectPhi, RectLength1, RectLength2)//生成矩形区域
gen_rectangle2 (MeasureROI2, Rect2Row, Rect2Col, RectPhi, RectLength1, RectLength2)//生成矩形区域
dev_display (ModelImage)
dev_set_color ('yellow')
dev_display (MeasureROI1)
dev_display (MeasureROI2)
* translate measurement ROIs to lie on XLD model (without clipping!)
get_system ('clip_region', OriginalClipRegion)
set_system ('clip_region', 'false')
move_region (MeasureROI1, MeasureROI1Ref, -CenterROIRow, -CenterROIColumn)//平移区域,平移的距离,而不是平移的目标点
move_region (MeasureROI2, MeasureROI2Ref, -CenterROIRow, -CenterROIColumn)
set_system ('clip_region', OriginalClipRegion)
DistRect1CenterRow := Rect1Row - CenterROIRow
DistRect1CenterCol := Rect1Col - CenterROIColumn
DistRect2CenterRow := Rect2Row - CenterROIRow
DistRect2CenterCol := Rect2Col - CenterROIColumn
if (USING_TRANSLATE_MEASURE != 0)
* -> measure objects are created only once in advance and then translated later
gen_measure_rectangle2 (Rect1Row, Rect1Col, RectPhi, RectLength1, RectLength2, Width, Height, 'bilinear', MeasureHandle1)//生成测量矩形
gen_measure_rectangle2 (Rect2Row, Rect2Col, RectPhi, RectLength1, RectLength2, Width, Height, 'bilinear', MeasureHandle2)
endif
stop ()
* step 2: find the objects in another image
read_image (SearchImage, 'razors2')
dev_display (SearchImage)
find_shape_model (SearchImage, ModelID, 0, 0, 0.8, 0, 0.5, 'least_squares', 0, 0.7, RowCheck, ColumnCheck, AngleCheck, Score)
if (|Score| > 0)
for i := 0 to |Score| - 1 by 1
* step 3: determine the affine transformation
vector_angle_to_rigid (0, 0, 0, RowCheck[i], ColumnCheck[i], AngleCheck[i], MovementOfObject)//计算刚性放射变换矩阵,从(0,0,0)变换到目标点。
affine_trans_contour_xld (ShapeModel, ModelAtNewPosition, MovementOfObject)//执行轮廓仿射变换
dev_display (ModelAtNewPosition)
* step 4: measure width and distance of the teeth
* -> display the moved ROIs
affine_trans_region (MeasureROI1Ref, MeasureROI1AtNewPosition, MovementOfObject, 'constant')//执行区域仿射变换
affine_trans_region (MeasureROI2Ref, MeasureROI2AtNewPosition, MovementOfObject, 'constant')//执行区域仿射变换
dev_display (MeasureROI1AtNewPosition)
dev_display (MeasureROI2AtNewPosition)
affine_trans_pixel (MovementOfObject, DistRect1CenterRow, DistRect1CenterCol, Rect1RowCheck, Rect1ColCheck)//执行点仿射变换
affine_trans_pixel (MovementOfObject, DistRect2CenterRow, DistRect2CenterCol, Rect2RowCheck, Rect2ColCheck)//执行点仿射变换
if (USING_TRANSLATE_MEASURE != 0)
* -> translate the already created measure objects
translate_measure (MeasureHandle1, Rect1RowCheck, Rect1ColCheck)//平移测量矩形到目标点
translate_measure (MeasureHandle2, Rect2RowCheck, Rect2ColCheck)
measure_pairs (SearchImage, MeasureHandle1, 2, 25, 'negative', 'all', RowEdge11, ColEdge11, Amp11, RowEdge21, ColEdge21, Amp21, Width1, Distance1)//测量一对
measure_pairs (SearchImage, MeasureHandle2, 2, 25, 'negative', 'all', RowEdge12, ColEdge12, Amp12, RowEdge22, ColEdge22, Amp22, Width2, Distance2)
else
* -> create new measure objects and destroy them after the measurement
RectPhiCheck := RectPhi + AngleCheck[i]
gen_measure_rectangle2 (Rect1RowCheck, Rect1ColCheck, RectPhiCheck, RectLength1, RectLength2, Width, Height, 'bilinear', MeasureHandle1)//生成测量矩形
gen_measure_rectangle2 (Rect2RowCheck, Rect2ColCheck, RectPhiCheck, RectLength1, RectLength2, Width, Height, 'bilinear', MeasureHandle2)
* step 5: perform the measurement
measure_pairs (SearchImage, MeasureHandle1, 2, 25, 'negative', 'all', RowEdge11, ColEdge11, Amp11, RowEdge21, ColEdge21, Amp21, Width1, Distance1)//测量一对
measure_pairs (SearchImage, MeasureHandle2, 2, 25, 'negative', 'all', RowEdge12, ColEdge12, Amp12, RowEdge22, ColEdge22, Amp22, Width2, Distance2)
close_measure (MeasureHandle1)//删除测量矩形
close_measure (MeasureHandle2)
endif
* step 6: check for too short or missing teeth
NumberTeeth1 := |Width1|
NumberTeeth2 := |Width2|
dev_set_color ('red')
if (NumberTeeth1 < 37)
for j := 0 to NumberTeeth1 - 2 by 1
if (Distance1[j] > 4.0)
RowFault := round(0.5 * (RowEdge11[j + 1] + RowEdge21[j]))
ColFault := round(0.5 * (ColEdge11[j + 1] + ColEdge21[j]))
disp_rectangle2 (WindowHandle, RowFault, ColFault, 0, 4, 4)//显示一个矩形
dev_open_window (0, Width + 20, 80, 80, 'black', WindowHandleZoom)//打开一个窗口
dev_set_part (RowFault - 10, ColFault - 10, RowFault + 10, ColFault + 10)
dev_display (SearchImage)
disp_rectangle2 (WindowHandleZoom, RowFault, ColFault, 0, 4, 4)
stop ()
dev_close_window ()
dev_set_part (0, 0, Height - 1, Width - 1)
endif
endfor
endif
if (NumberTeeth2 < 37)
for j := 0 to NumberTeeth2 - 2 by 1
if (Distance2[j] > 4.0)
RowFault := round(0.5 * (RowEdge12[j + 1] + RowEdge22[j]))
ColFault := round(0.5 * (ColEdge12[j + 1] + ColEdge22[j]))
disp_rectangle2 (WindowHandle, RowFault, ColFault, 0, 4, 4)
dev_open_window (0, Width + 20, 80, 80, 'black', WindowHandleZoom)
dev_set_part (RowFault - 10, ColFault - 10, RowFault + 10, ColFault + 10)
dev_display (SearchImage)
disp_rectangle2 (WindowHandleZoom, RowFault, ColFault, 0, 4, 4)
stop ()
dev_close_window ()
dev_set_part (0, 0, Height - 1, Width - 1)
endif
endfor
endif
dev_set_color ('yellow')
stop ()
endfor
endif
* ------------------- end of the application -----------------
* clean up
if (USING_TRANSLATE_MEASURE != 0)
close_measure (MeasureHandle1)//删除测量矩形
close_measure (MeasureHandle2)
endif
dev_update_window ('on')
clear_shape_model (ModelID)//清除模板
该例程实现了测量剃须刀片的检测,
例程中使用了两种方式,第一种是全部使用放射变换,第二种是平移测量矩形。
该例程中使用了inspect_shape_model算子,该算子可以查看图像的每层金字塔的状态,用来在创建模板时确定NumLevels(金字塔层数)和Contrast(对比度)这两个参数,方便地用于create_shape_model,create_scaled_shape_model或create_aniso_shape_model中。模型的表示是在多个图像金字塔级别上创建的,其中级别的数量由NumLevels确定。
具体用法可以参考该文章