地理坐标转换及校准方法

声明

出于工作需要写这么一篇文章。本人不是专业搞地理和数学的，文章中有不对或不准确的地方请指正，勿喷。

个人认为，在几个明确定义的坐标系之间，坐标应该有明确的转换方法（可能为近似），而不是使用这种拟合的方法。但有时一些临时使用的，或定义不明确的，或小范围的坐标系，这种拟合也不失为一种好的方法。

问题

有两个不同的地理坐标系，世界上某一点α在坐标系1上的坐标值为A，在坐标系2上的坐标值为B，如果知道A的坐标如何求B的坐标（或反之）呢，这就需要我们去找出坐标系1和坐标系2之间的映射关系。

原理

取一定数量的样本，在样本上使用拟合的方法，找出（方程）映射关系。

求解过程

我们先用数学语言表达整理一下问题和条件，有集合A{(x,y)|x属于坐标系1的经度,y属于坐标系1的纬度}，集合B{(m,n)|m属于坐标系2的经度,n属于坐标系2的纬度}，求从集合A到集合B的映射关系f，使f((x,y))=(m,n)，(x,y)属于A。

根据经验，A=>B有如下关系

m = A * x + B * y + C * x * x + D * y * y + E * x * y + F
n = G * x + H * y + I * x * x + J * y * y + K * x * y + L

而且，我们已经有了一些点的映射

(x1,y1) => (m1,n1)
(x2,y2) => (m2,n2)
...

分割线——————————————————————

这样我们就可以建立一个矩阵

  |  x1  y1  x1*x1  y1*y1  x1*y1  1  0   0   0      0      0        0  |     |  m1  |
  |  0   0   0      0      0      0  x1  y1  x1*x1  y1*y1  x1*y1    1  |  =  |  n1  |
  |  x2  y2  x2*x2  y2*y2  x2*y2  1  0   0   0      0      0        0  |     |  m2  |
  |  0   0   0      0      0      0  x2  y2  x2*x2  y2*y2  x2*y2    1  |     |  n2  |
  ...

然后最小二乘拟合求解即可。

上面的经验关系和拟合方法都可以根据不同情况进行调整。

打字困难，数学表达式难敲，凑合看。