基于COORD四参数坐标转换

1、原理

两个平面坐标系通过旋转、缩放、平移即可从一个坐标系转换到另外一个坐标系。

转换前坐标点（x2,y2）,转换后坐标点（x1,y1）；

转换公式如下：

求解方程一般用最小二乘法计算，这里提供计算工具，也可以自行代码。

2、计算方法

（1） 先记录两个坐标系的公共点坐标（至少两对）

（2） 打开COORD工具（见博客资源：坐标转换工具（四参数七参数计算）COORD软件.zip），选择计算四参数

（3） 输入公共点坐标信息

（4） 计算得到转换参数

DX表示x轴平移参数

DY表示y轴平移参数

T表示旋转角的弧度值

K表示缩放系数

（5） 代入公式进行坐标转换

这里提供的是python代码批处理坐标

def transfer(x,y):dx=391090.650080dy=2473060.653080T=0.017260K=0.999697x2=dx+K*(math.cos(T)*x-math.sin(T)*y)y2=dy+K*(math.sin(T)*x+math.cos(T)*y)return(x2,y2)

也可用LISP语言对CAD文件整体进行坐标转换

(defun c:sz2cgcs();定义一个函数 (setq pt1 (getpoint "First Pt:"))(setq pt2 (getpoint "Second Pt:"))(setq ss (ssget "c" pt1 pt2))(setq jd (* (/ 0.017260 pi) 180))(command "rotate" ss "" "0,0" jd)(command "scale" ss "" "0,0" 0.999697 )(command "move" ss "" "0,0" "391090.650080,2473060.653080" )(princ);清除 nil) ;函数定义结束