* This example program shows how to find scaled and rotated shape models.

 dev_update_pc ('off')

 dev_update_window ('off')

 dev_update_var ('off')

 read_image (Image, 'green-dot')

 get_image_size (Image, Width, Height) 获取了图像大小

 dev_close_window ()

 dev_open_window (, , Width, Height, 'black', WindowHandle)

 dev_set_color ('red')

 dev_display (Image)

 threshold (Image, Region, , ) 对图像进行二值化

 connection (Region, ConnectedRegions)  区域生长得到连通域

 select_shape (ConnectedRegions, SelectedRegions, 'area', 'and', , ) 通过面积进行筛选，得到里面的圆

 fill_up (SelectedRegions, RegionFillUp)   对圆进行填充

 dilation_circle (RegionFillUp, RegionDilation, 5.5)  对填充区域进行膨胀

 reduce_domain (Image, RegionDilation, ImageReduced) ROI操作得到imagereduced

 create_scaled_shape_model (ImageReduced, , rad(-), rad(), 'auto', 0.8, 1.0, 'auto', 'none', 'ignore_global_polarity', , , ModelID)
//基于区域创建匹配模型，得到模型的ID modelID

 get_shape_model_contours (Model, ModelID, ) 基于模型ID 得到模型的轮廓 model

 area_center (RegionFillUp, Area, RowRef, ColumnRef) 获得里面圆中心位置，相对于全图来说

 vector_angle_to_rigid (, , , RowRef, ColumnRef, , HomMat2D) //vector_angle_to_rigid只需要一个点对及一个角度对即可计算刚性变换矩阵，所以可利用find_shape_model的结果
//HomMat2D 通过顶点得到其变换矩阵 

 affine_trans_contour_xld (Model, ModelTrans, HomMat2D) //对XLD轮廓(contour)进行一个任意二维仿射变换。 将model轮廓转换为相对于全图的，得到ModelTrans

 dev_display (Image)

 dev_display (ModelTrans) //在全局中显示模型区域 ，注意modeltrans 为xld 轮廓！！！

 read_image (ImageSearch, 'green-dots')

 dev_display (ImageSearch)
//将之前模型ID 作为模板，

 find_scaled_shape_model (ImageSearch, ModelID, rad(-), rad(), 0.8, 1.0, 0.5, , 0.5, 'least_squares', , 0.8, Row, Column, Angle, Scale, Score)
//得到图像中模板的点 位置 和 角度 尺度 和匹配值

 for I :=  to |Score|- by

     hom_mat2d_identity (HomMat2DIdentity)

     hom_mat2d_translate (HomMat2DIdentity, Row[I], Column[I], HomMat2DTranslate)

     hom_mat2d_rotate (HomMat2DTranslate, Angle[I], Row[I], Column[I], HomMat2DRotate)

     hom_mat2d_scale (HomMat2DRotate, Scale[I], Scale[I], Row[I], Column[I], HomMat2DScale)
2*3 刚性变换矩阵

     affine_trans_contour_xld (Model, ModelTrans, HomMat2DScale)  把model 转换为modeltrans 

     gen_region_contour_xld(ModelTrans,Region2,'filled')      填充区域

     intensity ( Region2, ImageSearch, Mean, Deviation )      获得区域的灰度平均值和方差值

     dev_display (ModelTrans)

 endfor

 clear_shape_model (ModelID)

 HTuple HalHwndView1;

   void DispImage( const Hobject& hoImage);

 HalHwndView1=NULL;

 void CSpayTestDlg::DispImage(const Hobject& hoImage)

 {

 HTuple tmpWidth;

 HTuple tmpHeight;

 get_image_size(hoImage,&tmpWidth,&tmpHeight);

 Halcon::set_part(HalHwndView1, , , tmpHeight - , tmpWidth - );//

 Halcon::clear_window(HalHwndView1);

 Halcon::disp_obj(hoImage, HalHwndView1);

 }

 CRect rect;

 GetDlgItem(IDC_SHOW_WINDOW)->GetClientRect(rect);

 HWND HwndView1=GetDlgItem(IDC_SHOW_WINDOW)->GetSafeHwnd();

 Halcon::open_window(,,rect.Width(),rect.Height(),Hlong(HwndView1),"visible","",&HalHwndView1);

Hobject  Image;

HTuple  Width, Height, WindowHandle, ModelID;

Hobject  Region, ConnectedRegions, SelectedRegions;

Hobject  RegionFillUp, RegionDilation, ImageReduced, Model;

HTuple  Angle, Scale, Score, I, HomMat2DIdentity, HomMat2DTranslate;

HTuple  HomMat2DRotate, HomMat2DScale;

  HTuple  Area, RowRef, ColumnRef, HomMat2D, Row, Column;

    Hobject  ModelTrans;

void CSpayTestDlg::OnBnClickedButtonDetection()

{

    // TODO: 在此添加控件通知处理程序代码

   Image=IplImageToHImage(m_workImg);

    DispImage(Image);

   Halcon::set_color(HalHwndView1,"red");

   Halcon::set_draw(HalHwndView1,"margin");

   get_image_size(Image, &Width, &Height);

   
通过鼠标选择区域

   //HTuple  Row1,Column1,Row2,Column2,h_Area,h_Row,h_Column;

   //draw_rectangle1(HalHwndView1,&Row1,&Column1,&Row2,&Column2);

   //gen_rectangle1(&RegionDilation,Row1,Column1,Row2,Column2);

   //area_center(RegionDilation,&h_Area,&h_Row,&h_Column);

   //reduce_domain(Image,RegionDilation,&ImageReduced);

   threshold(Image, &Region, , );

   connection(Region, &ConnectedRegions);

   select_shape(ConnectedRegions, &SelectedRegions, "area", "and", , );

   fill_up(SelectedRegions, &RegionFillUp);

   dilation_circle(RegionFillUp, &RegionDilation, 5.5);

   reduce_domain(Image, RegionDilation, &ImageReduced);

   //roi 区域 ImageReduced

   create_scaled_shape_model(ImageReduced, , HTuple(-).Rad(), HTuple().Rad(),

       "auto", 0.8, 1.0, "auto", "none", "ignore_global_polarity", , , &ModelID);//模型是个数据 ModelID

   get_shape_model_contours(&Model, ModelID, );//获取模型的轮廓 Model中

   Halcon::write_shape_model(ModelID,"Model.shm");

   area_center(RegionFillUp, &Area, &RowRef, &ColumnRef);

   vector_angle_to_rigid(, , , RowRef, ColumnRef, , &HomMat2D);

   affine_trans_contour_xld(Model, &ModelTrans, HomMat2D);//roi 转换到原图的过程 基于2d仿射 将 model模型 轮廓 变化到trans上

   Halcon::disp_obj(ModelTrans,HalHwndView1);

}

     Image=IplImageToHImage(m_workImg);

     DispImage(Image);

     Halcon::set_color(HalHwndView1,"green");

     HTuple  ModelIDadd;

      Hobject  ModelTransadd;

     Halcon::read_shape_model( "Model.shm",&ModelIDadd);

     find_scaled_shape_model(Image, ModelIDadd, HTuple(-).Rad(), HTuple().Rad(),

         0.8, 1.0, 0.5, , 0.5, "least_squares", , 0.8, &Row, &Column, &Angle, &Scale,

         &Score);

     Hobject  Region2;

     HTuple  Mean, Deviation;

     for (I=; I<=(Score.Num())-; I+=)

     {

         hom_mat2d_identity(&HomMat2DIdentity);

         hom_mat2d_translate(HomMat2DIdentity, HTuple(Row[I]), HTuple(Column[I]), &HomMat2DTranslate);

         hom_mat2d_rotate(HomMat2DTranslate, HTuple(Angle[I]), HTuple(Row[I]), HTuple(Column[I]),

             &HomMat2DRotate);

         hom_mat2d_scale(HomMat2DRotate, HTuple(Scale[I]), HTuple(Scale[I]), HTuple(Row[I]),

             HTuple(Column[I]), &HomMat2DScale);

         affine_trans_contour_xld(Model, &ModelTransadd, HomMat2DScale);

         gen_region_contour_xld(ModelTrans,&Region2,"filled");

         intensity ( Region2, Image, &Mean, &Deviation );

         double mean=Mean[].D();

         double deviation=Deviation[].D();

         disp_obj(ModelTransadd,HalHwndView1);

     }

     clear_shape_model(ModelIDadd);