kikkimo / WgsToGcj

This class WGStoGCJ implements geodetic coordinate transform between geodetic coordinate in WGS84 coordinate system and geodetic coordinate in GCJ-02 coordinate system.

This class WGStoGCJ implements geodetic coordinate transform between geodetic coordinate in WGS84 coordinate system and geodetic coordinate in GCJ-02 coordinate system with high precision.

Supports 4 ways to convert coordinates form GCJ-02 coordinate system to WGS84 coordinate system:

• Simple linear iteration
• Analytical derivation
• Numerical derivation
• Automatic derivation(needed Ceres library)

WGS84 to GCJ-02: $$\eqalign{ & \Delta lo{n_0} = lo{n_{WGS84}} - 105 \cr & \Delta la{t_0} = la{t_{WGS84}} - 35 \cr & \Delta lo{n_1} = 300 + \Delta lo{n_0} + 2\Delta la{t_0} + 0.1\Delta lo{n_0}^2 + 0.1\Delta lo{n_0} \cdot \Delta la{t_0} + 0.1\sqrt {\left| {\Delta lo{n_0}} \right|} \cr & + \frac{2}{3}\left( {20\sin \left( {6\pi \cdot \Delta lo{n_0}} \right) + 20\sin \left( {2\pi \cdot \Delta lo{n_0}} \right)} \right) \cr & + \frac{2}{3}\left( {20\sin \left( {\pi \cdot \Delta lo{n_0}} \right) + 40\sin \left( {\frac{{\pi \cdot \Delta lo{n_0}}}{3}} \right)} \right) \cr & + \frac{2}{3}\left( {150\sin \left( {\frac{{\pi \cdot \Delta lo{n_0}}}{{12}}} \right) + 300\sin \left( {\frac{{\pi \cdot \Delta lo{n_0}}}{{30}}} \right)} \right) \cr & \Delta la{t_1} = - 100 + 2\Delta lo{n_0} + 3\Delta la{t_0} + 0.2\Delta la{t_0}^2 + 0.1\Delta lo{n_0} \cdot \Delta la{t_0} + 0.2\sqrt {\left| {\Delta lo{n_0}} \right|} \cr & + \frac{2}{3}\left( {20\sin \left( {6\pi \cdot \Delta lo{n_0}} \right) + 20\sin \left( {2\pi \cdot \Delta lo{n_0}} \right)} \right) \cr & + \frac{2}{3}\left( {20\sin \left( {\pi \cdot \Delta la{t_0}} \right) + 40\sin \left( {\frac{{\pi \cdot \Delta la{t_0}}}{3}} \right)} \right) \cr & + \frac{2}{3}\left( {160\sin \left( {\frac{{\pi \cdot \Delta la{t_0}}}{{12}}} \right) + 320\sin \left( {\frac{{\pi \cdot \Delta la{t_0}}}{{30}}} \right)} \right) \cr & {B_{WGS84}} = Deg2Rad(la{t_{WGS84}}) \cr & W = \sqrt {1 - {e^2}\sin {{\left( {{B_{WGS84}}} \right)}^2}} \cr & N = \frac{a}{W} \cr & \Delta lo{n_2} = \frac{{180\Delta lo{n_1}}}{{\pi N \cdot \cos \left( {{B_{WGS84}}} \right)}} \cr & \Delta la{t_2} = \frac{{180\Delta la{t_1} \cdot {W^2}}}{{\pi N \cdot \left( {1 - {e^2}} \right)}} \cr & lo{n_{GCJ - 02}} = lo{n_{WGS84}} + \Delta lo{n_2} \cr & la{t_{GCJ - 02}} = la{t_{WGS84}} + \Delta la{t_2} \cr} $$

Notice:

1. Use C++ standard 17.
2. Ceres is not a mandatory option for this project, and we suggest to set 'USE_CERES = OFF'. If set check 'USE_CERES', you should set the following librarys path: Ceres Glog GFlags Eigen

• Author: kikkimo
• Email: [e-mail]
• Details: https://blog.csdn.net/gudufuyun/article/details/106738942
• School of Remote Sensing and Information Engineering,WuHan University,Wuhan, Hubei, P.R. China. 430079
• Copyright (c) 2020. All rights reserved.

This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version.

这个类实现了WGS84大地坐标与GCJ-02大地坐标的相互转换，具有较高的转换精度，可通过迭代实现误差1mm以内的逆转换。 支持4种方式将GCJ-02坐标转换为WGS84坐标:

• 简单线性迭代
• 解析求导
• 数值求导
• 自动求导(需借助Ceres库)

注意事项： 1、需使用C++17标准 2、Ceres并非本工程的必选项，使用CMake创建工程时，建议不要勾选USE_CERES。如果勾选了USE_CERES，需要将以下库引入至本工程： Ceres Glog GFlags Eigen

• 作者：kikkimo
• 邮箱：[e-mail]
• 原理：https://blog.csdn.net/gudufuyun/article/details/106738942
• 地址：武汉大学遥感信息工程学院，**湖北省武汉市，430079
• Copyright (c) 2020. All rights reserved.