坐标系转换基础知识

地图应用、定位服务、GIS系统开发中，经常会遇到不同坐标体系之间的数据不兼容问题。例如，GPS设备采集的数据通常采用WGS84坐标系，而国内互联网地图服务大多使用GCJ-02或BD09坐标系。如果直接将不同坐标系的数据叠加显示，会出现几十米甚至数百米的位置偏移。World Airport Maps

Spring Boot作为企业级Java开发框架，常用于构建地图服务接口、空间数据处理平台以及位置服务后台。通过在Spring Boot项目中集成坐标转换算法，可以实现GCJ-02、BD09、CGCS2000与WGS84之间的数据转换，为地图应用提供统一的数据基准。

常见坐标系特点如下：

WGS84（World Geodetic System 1984）

全球通用地理坐标系，GPS设备默认输出通常采用该坐标系，也是国际GIS系统广泛使用的标准。
GCJ-02（火星坐标系）

国内互联网地图常用坐标体系，在WGS84基础上经过偏移处理，高德地图、腾讯地图等服务通常采用该坐标系。GIS工具箱
BD09（百度坐标系）

百度地图在GCJ-02基础上再次加密偏移形成的坐标体系，主要应用于百度地图生态。百度开发者中心
CGCS2000（中国大地坐标系2000）

我国国家大地坐标参考系统，与WGS84非常接近，主要用于测绘、国土、规划等专业GIS场景。二者之间通常认为存在厘米级差异。Hylab

Spring Boot实现坐标转换的整体设计

在Spring Boot项目中实现坐标转换，一般采用服务层封装方式，将坐标转换逻辑独立出来。

典型架构如下：


Controller
    |
    |
CoordinateService
    |
    |
CoordinateConvertUtil
    |
    |
转换算法

例如：

用户请求：


POST /coordinate/convert

{
    "type":"GCJ02_TO_WGS84",
    "lng":116.404,
    "lat":39.915
}

服务返回：


{
    "lng":116.397428,
    "lat":39.90923
}

通过统一接口，可以方便支持：

GCJ02 → WGS84
WGS84 → GCJ02
BD09 → GCJ02
GCJ02 → BD09
BD09 → WGS84
CGCS2000 → WGS84

创建Spring Boot坐标对象

首先定义坐标实体类：

Java
@Data
public class Coordinate {

    /**
     * 经度
     */
    private double longitude;

    /**
     * 纬度
     */
    private double latitude;

    public Coordinate(double longitude, double latitude) {
        this.longitude = longitude;
        this.latitude = latitude;
    }
}

该对象用于统一封装各种坐标系的数据。

实现GCJ-02转WGS84

GCJ-02转换到WGS84属于逆向计算。

由于GCJ-02算法属于非线性偏移，无法简单通过减去固定值完成，需要采用迭代逼近方式获取结果。

核心流程：

判断坐标是否位于中国区域。
根据GCJ-02正向算法计算预测值。
不断修正误差。
获取接近真实WGS84坐标。

实现代码：

Java
public class CoordinateConvertUtil {


    private static final double EARTH_RADIUS = 6378245.0;

    private static final double EE = 0.00669342162296594323;


    public static Coordinate gcj02ToWgs84(
            double lng,
            double lat) {


        if (outOfChina(lng, lat)) {
            return new Coordinate(lng, lat);
        }


        double threshold = 0.000000001;

        double mglng = lng;
        double mglat = lat;

        double dLng;
        double dLat;


        while (true) {

            Coordinate temp =
                    wgs84ToGcj02(
                            mglng,
                            mglat);


            dLng = temp.getLongitude() - lng;
            dLat = temp.getLatitude() - lat;


            if (Math.abs(dLng) < threshold
                    && Math.abs(dLat) < threshold) {

                break;
            }


            mglng -= dLng;
            mglat -= dLat;
        }


        return new Coordinate(
                mglng,
                mglat);
    }
}

这种方式比简单偏移修正更加准确，适用于地图定位、轨迹回放等业务场景。

实现WGS84转GCJ-02

WGS84转GCJ02属于正向转换。

核心公式：

Java
public static Coordinate wgs84ToGcj02(
        double lng,
        double lat) {


    if (outOfChina(lng, lat)) {
        return new Coordinate(lng, lat);
    }


    double dLat =
            transformLat(
                    lng - 105.0,
                    lat - 35.0);


    double dLng =
            transformLng(
                    lng - 105.0,
                    lat - 35.0);


    double radLat =
            lat / 180.0 * Math.PI;


    double magic =
            Math.sin(radLat);


    magic =
            1 - EE * magic * magic;


    double sqrtMagic =
            Math.sqrt(magic);


    dLat =
            (dLat * 180.0)
            /
            ((EARTH_RADIUS
            * (1 - EE))
            /
            (magic * sqrtMagic)
            * Math.PI);


    dLng =
            (dLng * 180.0)
            /
            (EARTH_RADIUS
            /
            sqrtMagic
            * Math.cos(radLat)
            * Math.PI);


    double mgLat = lat + dLat;
    double mgLng = lng + dLng;


    return new Coordinate(
            mgLng,
            mgLat);
}

实现BD09与GCJ02互转

百度坐标BD09是在GCJ02基础上的再次偏移，因此转换关系如下：


WGS84
 |
 |
GCJ02
 |
 |
BD09

GCJ02转BD09：

Java
public static Coordinate gcj02ToBd09(
        double lng,
        double lat) {


    double z =
            Math.sqrt(
                    lng * lng
                    +
                    lat * lat)
            +
            0.00002 *
            Math.sin(lat * Math.PI);


    double theta =
            Math.atan2(lat, lng)
            +
            0.000003 *
            Math.cos(lng * Math.PI);


    double bdLng =
            z * Math.cos(theta)
            + 0.0065;


    double bdLat =
            z * Math.sin(theta)
            + 0.006;


    return new Coordinate(
            bdLng,
            bdLat);
}

BD09转GCJ02：

Java
public static Coordinate bd09ToGcj02(
        double lng,
        double lat) {


    double x =
            lng - 0.0065;


    double y =
            lat - 0.006;


    double z =
            Math.sqrt(
                    x*x+y*y)
            -
            0.00002 *
            Math.sin(y*Math.PI);


    double theta =
            Math.atan2(y,x)
            -
            0.000003 *
            Math.cos(x*Math.PI);


    double gcjLng =
            z*Math.cos(theta);


    double gcjLat =
            z*Math.sin(theta);


    return new Coordinate(
            gcjLng,
            gcjLat);
}

CGCS2000转换WGS84实现方式

CGCS2000和WGS84属于不同的大地坐标基准，但二者非常接近。

在普通互联网地图应用中，通常可以直接认为：


CGCS2000 ≈ WGS84

简单转换：

Java
public static Coordinate cgcs2000ToWgs84(
        double lng,
        double lat){

    return new Coordinate(
            lng,
            lat);
}

不过在高精度测绘项目中，需要使用七参数转换模型：


X2 = DX + (1+M)X1 - RzY1 + RyZ1
Y2 = DY + RzX1 + (1+M)Y1 - RxZ1
Z2 = DZ - RyX1 + RxY1 + (1+M)Z1

实际工程需要结合：

椭球参数
平移参数
旋转参数
比例因子

进行三维坐标转换。

Spring Boot封装转换接口

创建Service：

Java
@Service
public class CoordinateService {


    public Coordinate convert(
            String type,
            Coordinate coordinate){


        switch(type){

            case "GCJ02_TO_WGS84":

                return CoordinateConvertUtil
                    .gcj02ToWgs84(
                    coordinate.getLongitude(),
                    coordinate.getLatitude());


            case "WGS84_TO_GCJ02":

                return CoordinateConvertUtil
                    .wgs84ToGcj02(
                    coordinate.getLongitude(),
                    coordinate.getLatitude());


            case "GCJ02_TO_BD09":

                return CoordinateConvertUtil
                    .gcj02ToBd09(
                    coordinate.getLongitude(),
                    coordinate.getLatitude());


            case "BD09_TO_GCJ02":

                return CoordinateConvertUtil
                    .bd09ToGcj02(
                    coordinate.getLongitude(),
                    coordinate.getLatitude());


            default:

                throw new RuntimeException(
                        "不支持的坐标转换类型");
        }
    }
}

Controller：

Java
@RestController
@RequestMapping("/coordinate")
public class CoordinateController {


    @Autowired
    private CoordinateService service;


    @PostMapping("/convert")
    public Coordinate convert(
            @RequestBody CoordinateRequest request){

        return service.convert(
                request.getType(),
                new Coordinate(
                request.getLng(),
                request.getLat()));
    }
}

这样即可提供REST接口供前端地图系统调用。

坐标转换中的常见问题

1. 为什么转换后仍然存在偏差？

坐标转换误差来源主要包括：

GPS设备自身误差
地图服务数据误差
算法实现差异
坐标采集时间不同

GCJ-02和BD09属于偏移坐标体系，不同实现可能存在一定误差。开发工具

2. 是否所有地区都需要GCJ02转换？

不需要。

GCJ02主要针对中国大陆区域。

代码中通常增加判断：

Java
private static boolean outOfChina(
        double lng,
        double lat){

    return lng < 72.004
            || lng > 137.8347
            || lat < 0.8293
            || lat > 55.8271;
}