基于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.650080
    dy=2473060.653080
    T=0.017260
    K=0.999697
    x2=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
)   ;函数定义结束