美文网首页
常用地图经纬度转换,以及遇到的问题和解决方式

常用地图经纬度转换,以及遇到的问题和解决方式

作者: 你可以叫我老白 | 来源:发表于2022-04-27 17:28 被阅读0次

    坐标系分类

    第一种分类

    1、

    GPS,WGS-84,原始坐标体系。一般用国际标准的GPS记录仪记录下来的坐标,都是GPS的坐标。很可惜,在中国,任何一个地图产品都不允许使用GPS坐标,据说是为了保密。

    2、

    GCJ-02,国测局02年发布的坐标体系。又称“火星坐标”。在中国,必须至少使用GCJ-02的坐标体系。比如谷歌,腾讯,高德都在用这个坐标体系。GCJ-02也是国内最广泛使用的坐标体系。

    3、其他坐标体系。一般都是由GCJ-02进过偏移算法得到的。这种体系就根据每个公司的不同,坐标体系都不一样了。比如,百度和搜狗就使用自己的坐标体系,与其他坐标体系不兼容。

    第二种分类

    首先明白,所有坐标体系的原点,都是非洲。

    1、经纬度。这个是球面坐标,对于北京来说,就是(116.38817139.935961)这样的坐标。比如腾讯、高德、百度都是这样的经纬度坐标。谷歌是经纬度顺序写反的经纬度坐标。

    如果是度分秒坐标,需要进行转换,才能得到这样的经纬度坐标。详见坐标转换。

    2、墨卡托坐标。平面坐标,相当于是直线距离,数字一般都比较大,像这样的。(215362.00021333335

    99526.00034912192)

    墨卡托坐标,主要用于程序的后台计算。直线距离嘛,加加减减几乎计算方便。

    搜狗地图API就是直接使用的墨卡托坐标。

    经纬度转换工具类

    static double pi = 3.14159265358979324;

    static double a = 6378245.0;

    static double ee = 0.00669342162296594323;

    public final static double x_pi = pi * 3000.0 / 180.0;

    /**

    * 高德经纬度转百度经纬度

    */

    public static double[] gcj2bd(double lat, double lon) {

    double z = Math.sqrt(lon * lon + lat * lat) + 0.00002 * Math.sin(lat * x_pi);

    double theta = Math.atan2(lat, lon) + 0.000003 * Math.cos(lon * x_pi);

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

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

    return new double[] { bd_lat, bd_lon };

    }

    /**

    * 百度经纬度转高德

    */

    public static double[] bd2gcj(double lat, double lon) {

    double x = lon - 0.0065, y = lat - 0.006;

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

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

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

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

    return new double[] { gg_lat, gg_lon };

    }

    /**

    * 国标经纬度转高德经纬度

    */

    public static double[] wgs2gcj(double lat, double lon) {

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

    double dLon = transformLon(lon - 105.0, lat - 35.0);

    double radLat = lat / 180.0 * pi;

    double magic = Math.sin(radLat);

    magic = 1 - ee * magic * magic;

    double sqrtMagic = Math.sqrt(magic);

    dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * pi);

    dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * pi);

    double mgLat = lat + dLat;

    double mgLon = lon + dLon;

    return new double[]{ mgLat, mgLon };

    }

    /**

    * 高德转到GPS

    */

    public static double[] gcj2wgc(double lat, double lon) {

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

    double dLon = transformLon(lon - 105.0, lat - 35.0);

    double radLat = lat / 180.0 * pi;

    double magic = Math.sin(radLat);

    magic = 1 - ee * magic * magic;

    double sqrtMagic = Math.sqrt(magic);

    dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * pi);

    dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * pi);

    double mgLat = lat + dLat;

    double mgLon = lon + dLon;

    double lontitude = lon * 2 - mgLon;

    double latitude = lat * 2 - mgLat;

    return new double[]{latitude, lontitude};

    }

    /**

    * 国标转百度

    */

    public static double[] wgs2bd(double lat, double lon) {

    double[] wgs2gcj = wgs2gcj(lat, lon);

    return gcj2bd(wgs2gcj[0], wgs2gcj[1]);

    }

    /**

    * 经纬度转换

    */

    private static double transformLat(double lat, double lon) {

    double ret = -100.0 + 2.0 * lat + 3.0 * lon + 0.2 * lon * lon + 0.1 * lat * lon + 0.2 * Math.sqrt(Math.abs(lat));

    ret = getRet(lat, lon, ret);

    ret += (160.0 * Math.sin(lon / 12.0 * pi) + 320 * Math.sin(lon * pi / 30.0)) * 2.0 / 3.0;

    return ret;

    }

    private static double getRet(double lat, double lon, double ret) {

    ret += (20.0 * Math.sin(6.0 * lat * pi) + 20.0 * Math.sin(2.0 * lat * pi)) * 2.0 / 3.0;

    ret += (20.0 * Math.sin(lon * pi) + 40.0 * Math.sin(lon / 3.0 * pi)) * 2.0 / 3.0;

    return ret;

    }

    private static double transformLon(double lat, double lon) {

    double ret = 300.0 + lat + 2.0 * lon + 0.1 * lat * lat + 0.1 * lat * lon + 0.1 * Math.sqrt(Math.abs(lat));

    ret = getRet(lat, lat, ret);

    ret += (150.0 * Math.sin(lat / 12.0 * pi) + 300.0 * Math.sin(lat / 30.0 * pi)) * 2.0 / 3.0;

    return ret;

    }

    PS:因为经纬度的获取是依靠卫星的,所以屋内无法获取经纬度。并且哪怕你站在一个点移动,经纬度的获取也是有飘忽的,允许误差大概在10米以内即可。

    软硬件结合时要提前沟通好使用的经纬度分类。比如软件使用的是高德坐标,而硬件获取的则是其他坐标系或者使用的安卓原生坐标,尤其是后者转换后仍旧不属于任何一种坐标,所以尽量避免硬件传输安卓原生坐标,可以通过相关地图,比如高德的sdk进行转换。

    创作不易,如果这篇文章对你有用,请点个赞谢谢♪(・ω・)ノ!

    相关文章

      网友评论

          本文标题:常用地图经纬度转换,以及遇到的问题和解决方式

          本文链接:https://www.haomeiwen.com/subject/zsudyrtx.html