Halcon中的坐标系特点及XLD的镜像转换

我们知道，Halcon中的坐标系的原点在左上角，而一般二维平面坐标系的原点在左下角。那么Halcon中坐标系和一般的二维坐标系有什么区别呢？我通过下面这个例子来分析。

 gen_image_const (Image, 'byte', , )
 dev_set_draw ('margin')

 *点1
 gen_circle (Circle1, 10, 10, )
 disp_message (, '(10 ,10)', 'image', , , 'white', 'false')
 *点2
 gen_circle (Circle2, 100, 200, )
 disp_message (, '(100 ,200)', 'image', , , 'white', 'false')
 *点3
 gen_circle (Circle3, 400, 100, )
 disp_message (, '(400 ,100)', 'image', , , 'white', 'false')

 union2 (Circle1,Circle2, RegionUnion)
 union2 (RegionUnion,Circle2, RegionUnion)
 union2 (RegionUnion,Circle3, RegionUnion)
 gen_contour_region_xld (RegionUnion, Contours, 'border')
 write_contour_xld_dxf (Contours,'C:/Users/happy xia/Desktop/XLD的镜像变换/threePoint.dxf')

这三个点在Halcon中的位置如下图：

上面的程序最终生成了一个名为threePoint.dxf的文件。我用文本查看器（notepad++、记事本等）打开查看它的数据：

红框中“10”下面的数据代表点的x值，“20”下面的数据代表点的y值，即截图里有两个点(98.5, 396.5)、(99.5, 396.5)，这两个点都是XLD圆上的点，很明显这两个点所在圆的圆心坐标是(100, 400)。但是我们看最开始的那张图，发现只有点(400, 100)，并没有点(100, 400)！

根据上面的事实现象，我们不难发现Halcon中的坐标系实际是这样的：

根据导出的dxf文件中的坐标点反推，之前的那张图的点的坐标实际上是这样的：

这两个坐标系的区别如下图。不难看出，它们是上下镜像的关系（只需要用CAD看图软件打开threePoint.dxf文件一看便知）。

上文详细分析了Halcon坐标系的特点，由于有这样的特点，Halcon中的XLD和导出的dxf文件轮廓就存在上下镜像（上下颠倒）的关系，那么我们需要解决的一个问题就是：如何进行XLD的镜像转换？

镜像，对应的单词是“mirror”，我们在Halcon的帮助中搜索与之相关的算子，仅有mirror_image和mirror_region。

根据我前面的博客文章，我们知道，如果XLD是闭合的，那么可将XLD转为Region，然后再转回XLD，因此可将Region作为中介，来镜像变换XLD。

方法一：以Region为中介，通过mirror_region算子间接对闭合XLD进行镜像转换。

原XLD如下图所示。

 gen_image_const (Image, 'byte', , )
 read_contour_xld_dxf (Contours,'unionContour.dxf', [], [], DxfStatus)

 *方法1：先转换为Region，再镜像Region，最后转回到XLD
 gen_region_contour_xld (Contours, Region, 'margin')
 mirror_region (Region, RegionMirror, 'row', 6000)
 gen_contour_region_xld (RegionMirror, ContoursMirror, 'border')

效果如下：

方法二：由于XLD都是由点组成的，因此可以点对点镜像重绘

 gen_image_const (Image, 'byte', , )
 read_contour_xld_dxf (Contours,'unionContour.dxf', [], [], DxfStatus)

 *方法2：点对点镜像重绘
 count_obj (Contours, Num)
 Step := 5
 axisValue := 6000

 gen_empty_obj (Contour2)
 for i :=  to Num by
     select_obj (Contours, Contour, i)
     get_contour_xld (Contour, Row, Col)
     Row1 := []
     Col1 := []
     for j :=  to |Row|- by Step

         Row1:=[Row1,axisValue - Row[j]]
         Col1:=[Col1,Col[j]]

     endfor

     test_closed_xld (Contour, IsClosed)
     if (IsClosed = )
         Row1:=[Row1,axisValue - Row[]]
         Col1:=[Col1,Col[]]

     endif

     gen_contour_polygon_xld (Contour1, Row1, Col1)
     smooth_contours_xld (Contour1, SmoothedContours, )
     concat_obj (Contour2, SmoothedContours, Contour2)
 endfor

效果如下：

对于方法二程序有疑问的，可以参考我之前写的文章：http://www.cnblogs.com/xh6300/p/7414256.html

上面两种方法有一些小问题，不够完美。方法一主要的局限是只适用于闭合的XLD，方法二对点进行处理的时候，为提高运算速度进行了间隔采样，导致最终得到的XLD轮廓和初始XLD有细微差异。

其实Halcon还提供了一种最佳解决方案，即仿射变换：hom_mat2d_reflect —— 变换后的对象与初始对象相对于两点确定的一条直线彼此对称，其实就是镜像。

方法三：使用hom_mat2d_reflect 直接进行镜像的仿射变换

 gen_image_const (Image, 'byte', , )
 read_contour_xld_dxf (Contours,'unionContour.dxf', [], [], DxfStatus)

 *方法3：使用仿射变换镜像
 hom_mat2d_identity (HomMat2DIdentity)
 hom_mat2d_reflect (HomMat2DIdentity, 3000, 0, 3000, 100, HomMat2DReflect)
 affine_trans_contour_xld (Contours, ContoursAffinTrans, HomMat2DReflect)

效果如下：