华测坐标转换工具中的TGO水平平差法计算公式

先吐槽一下，首先这个TGO水平平差法我在网上真的没找到相关的计算公式，由于任务需要，导师要求转换后的结果要和该软件的结果保持一致，这就相当的困难了，公式没给，参数没给，我也不是学测绘专业的，导致这个事情相当的麻烦，更何况根据我这一段时间的学习以及编写，我个人认为不应该把某个企业做的东西叫做标杆，没有所谓的标杆。比如这个TGO水平平差法，并不是一个大众熟知、随便一搜就有定义的公认方法，我个人认为这个方法是华测自己下的一个定义，当然如果有解释说明可以留言@我，不胜感激。

说实话，什么叫TGO水平平差法？我到现在也没弄懂，以下公式由本人通过结果反推公式得到，经过了一系列数据检验，可以保证该公式的正确性。

TGO水平平差法一共7个参数，分别为北原点xo、东原点yo、北平移量r_x、东平移量r_y、旋转角度θ、比例因子S,输入值为x，y坐标，输出值为X,Y。

公式如下：

X = S * (cosθ * x - sinθ * (y - yo)) + (0.5 - cosθ) * xo * S + xo + r_x

Y = S * (sinθ * (x - xo) + cosθ * y) + (0.5 - cosθ) * yo * S + yo + r_y

结果图对照图：

此外笔者还复现了很多华测坐标转换相关的公式代码，均可有偿咨询。