当前位置:首页 > 技术分享 > 正文内容

地图经纬度转换工具(天地图,高德地图,百度地图,腾讯地图)

admin2年前 (2023-02-01)技术分享4314

uniapp 中不同地图坐标系经纬度转换(天地图,高德地图,百度地图,腾讯地图)

 WGS-84:是国际标准,GPS坐标(Google Earth使用、或者GPS模块、天地图)
 GCJ-02:中国坐标偏移标准,Google Map、高德、腾讯使用
 BD-09:百度坐标偏移标准,Baidu Map使用


function mapTool() {  
    this.x_PI = 3.14159265358979324 * 3000.0 / 180.0,  
    this.PI = 3.1415926535897932384626,  
    this.ee = 0.00669342162296594323,  
    this.a = 6378245.0;  
}  
/**  
 * GCJ02 转换为 WGS84  
 * @param lng  
 * @param lat  
 * @returns {*[]}  
 */  mapTool.prototype.gcj02towgs84 = function(lng, lat) {  
        if (this.out_of_china(lng, lat)) {  
            return [lng, lat]  
        } else {  
            let dlat = this.transformlat(lng - 105.0, lat - 35.0);  
            let dlng = this.transformlng(lng - 105.0, lat - 35.0);  
            let radlat = lat / 180.0 * this.PI;  
            let magic = Math.sin(radlat);  
            magic = 1 - this.ee * magic * magic;  
            let sqrtmagic = Math.sqrt(magic);  
            dlat = (dlat * 180.0) / ((this.a * (1 - this.ee)) / (magic * sqrtmagic) * this.PI);  
            dlng = (dlng * 180.0) / (this.a / sqrtmagic * Math.cos(radlat) * this.PI);  
            let mglat = lat + dlat;  
            let mglng = lng + dlng;  
            return [lng * 2 - mglng, lat * 2 - mglat]  
        }  
    },  
    /**  
     * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换  
     * 即谷歌、高德 转 百度  
     * @param lng  
     * @param lat  
     * @returns {*[]}  
     */  
    mapTool.prototype.gcj02tobd09 = function(lng, lat) {  
        let z = Math.sqrt(lng * lng + lat * lat) + 0.00002 * Math.sin(lat * this.x_PI);  
        let theta = Math.atan2(lat, lng) + 0.000003 * Math.cos(lng * this.x_PI);  
        let bd_lng = z * Math.cos(theta) + 0.0065;  
        let bd_lat = z * Math.sin(theta) + 0.006;  
        return [bd_lng, bd_lat]  
    },  
    /**  
     * WGS84转GCj02  
     * @param lng  
     * @param lat  
     * @returns {*[]}  
     */  
    mapTool.prototype.wgs84togcj02 = function(lng, lat) {  
        if (this.out_of_china(lng, lat)) {  
            return [lng, lat]  
        } else {  
            let dlat = this.transformlat(lng - 105.0, lat - 35.0);  
            let dlng = this.transformlng(lng - 105.0, lat - 35.0);  
            let radlat = lat / 180.0 * this.PI;  
            let magic = Math.sin(radlat);  
            magic = 1 - this.ee * magic * magic;  
            let sqrtmagic = Math.sqrt(magic);  
            dlat = (dlat * 180.0) / ((this.a * (1 - this.ee)) / (magic * sqrtmagic) * this.PI);  
            dlng = (dlng * 180.0) / (this.a / sqrtmagic * Math.cos(radlat) * this.PI);  
            let mglat = lat + dlat;  
            let mglng = lng + dlng;  
            return [mglng, mglat]  
        }  
    },  
    /**  
     * GCJ02 转换为 WGS84  
     * @param lng  
     * @param lat  
     * @returns {*[]}  
     */  
    mapTool.prototype.gcj02towgs84 = function(lng, lat) {  
        if (this.out_of_china(lng, lat)) {  
            return [lng, lat]  
        } else {  
            let dlat = this.transformlat(lng - 105.0, lat - 35.0);  
            let dlng = this.transformlng(lng - 105.0, lat - 35.0);  
            let radlat = lat / 180.0 * this.PI;  
            let magic = Math.sin(radlat);  
            magic = 1 - this.ee * magic * magic;  
            let sqrtmagic = Math.sqrt(magic);  
            dlat = (dlat * 180.0) / ((this.a * (1 - this.ee)) / (magic * sqrtmagic) * this.PI);  
            dlng = (dlng * 180.0) / (this.a / sqrtmagic * Math.cos(radlat) * this.PI);  
            let mglat = lat + dlat;  
            let mglng = lng + dlng;  
            return [lng * 2 - mglng, lat * 2 - mglat]  
        }  
    },  
    mapTool.prototype.transformlat = function(lng, lat) {  
        let ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + 0.1 * lng * lat + 0.2 * Math.sqrt(Math.abs(lng));  
        ret += (20.0 * Math.sin(6.0 * lng * this.PI) + 20.0 * Math.sin(2.0 * lng * this.PI)) * 2.0 / 3.0;  
        ret += (20.0 * Math.sin(lat * this.PI) + 40.0 * Math.sin(lat / 3.0 * this.PI)) * 2.0 / 3.0;  
        ret += (160.0 * Math.sin(lat / 12.0 * this.PI) + 320 * Math.sin(lat * this.PI / 30.0)) * 2.0 / 3.0;  
        return ret  
    },  
    mapTool.prototype.transformlng = function(lng, lat) {  
        let ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + 0.1 * lng * lat + 0.1 * Math.sqrt(Math.abs(lng));  
        ret += (20.0 * Math.sin(6.0 * lng * this.PI) + 20.0 * Math.sin(2.0 * lng * this.PI)) * 2.0 / 3.0;  
        ret += (20.0 * Math.sin(lng * this.PI) + 40.0 * Math.sin(lng / 3.0 * this.PI)) * 2.0 / 3.0;  
        ret += (150.0 * Math.sin(lng / 12.0 * this.PI) + 300.0 * Math.sin(lng / 30.0 * this.PI)) * 2.0 / 3.0;  
        return ret  
    },  
    /**  
     * 判断是否在国内,不在国内则不做偏移  
     * @param lng  
     * @param lat  
     * @returns {boolean}  
     */  
    mapTool.prototype.out_of_china = function(lng, lat) {  
        return (lng < 72.004 || lng > 137.8347) || ((lat < 0.8293 || lat > 55.8271) || false);  
    },  
    module.exports.mapTool = mapTool;

php类

<?php
/**
 * @name Coordinate
 * @desc 坐标转换
 */
class Coordinate
{
    const x_PI = 52.35987755982988;
    const PI = 3.1415926535897932384626;
    const a = 6378245.0;
    const ee = 0.00669342162296594323;
 
    /**
     * 百度坐标系(BD-09) 转 火星坐标系(GCJ-02)
     * @param 
     * @return 
     **/
    public static function bd09ToGcj02($bd_lon, $bd_lat)
    {
        $x = $bd_lon - 0.0065;
        $y = $bd_lat - 0.006;
        $z = sqrt($x * $x + $y * $y) - 0.00002 * sin($y * self::x_PI);
        $theta = atan2($y, $x) - 0.000003 * cos($x * self::x_PI);
        $g_lon = $z * cos($theta);
        $g_lat = $z * sin($theta);
        return array('lon' => $g_lon, 'lat' => $g_lat);
    }
 
    /**
     * 火星坐标系(GCJ-02) 转 百度坐标系(BD-09)
     * 即谷歌、高德 转 百度
     * @param 
     * @return 
     **/
    public static function gcj02Tobd09($g_lon, $g_lat)
    {
        $z = sqrt($g_lon * $g_lon + $g_lat * $g_lat) + 0.00002 * sin($g_lat * self::x_PI);
        $theta = atan2($g_lat, $g_lon) + 0.000003 * cos($g_lon * self::x_PI);
        $bd_lon = $z * cos($theta) + 0.0065;
        $bd_lat = $z * sin($theta) + 0.006;
        return array('lon' => $bd_lon, 'lat' => $bd_lat);
    }
 
    /**
     * WGS84 转 GCj02
     * @param 
     * @return 
     **/
    public static function wgs84ToGcj02($w_lon, $w_lat)
    {
        $dlat = self::transFormLat($w_lon - 105.0, $w_lat - 35.0);
        $dlon = self::transFormLon($w_lon - 105.0, $w_lat - 35.0);
        $radlat = $w_lat / 180.0 * self::PI;
        $magic = sin($radlat);
        $magic = 1 - self::ee * $magic * $magic;
        $sqrtmagic = sqrt($magic);
        $dlat = ($dlat * 180.0) / ((self::a * (1 - self::ee)) / ($magic * $sqrtmagic) * self::PI);
        $dlon = ($dlon * 180.0) / (self::a / $sqrtmagic * cos($radlat) * self::PI);
        $g_lat = $w_lat + $dlat;
        $g_lon = $w_lon + $dlon;
        return array('lon' => $g_lon, 'lat' => $g_lat);
    }
 
    /**
     * GCJ02 转换为 WGS84
     * @param 
     * @return 
     **/
    public static function gcj02ToWgs84($g_lon, $g_lat)
    {
        $dlat = self::transFormLat($g_lon - 105.0, $g_lat - 35.0);
        $dlon = self::transFormLon($g_lon - 105.0, $g_lat - 35.0);
        $radlat = $g_lat / 180.0 * self::PI;
        $magic = sin($radlat);
        $magic = 1 - self::ee * $magic * $magic;
        $sqrtmagic = sqrt($magic);
        $dlat = ($dlat * 180.0) / ((self::a * (1 - self::ee)) / ($magic * $sqrtmagic) * self::PI);
        $dlon = ($dlon * 180.0) / (self::a / $sqrtmagic * cos($radlat) * self::PI);
        $w_lat = $g_lat + $dlat;
        $w_lon = $g_lon + $dlon;
        return array('lon' => $g_lon * 2 -$w_lon, 'lat' => $g_lat * 2 -$w_lat);
    }
 
    /**
     * BD09 转换为 WGS84
     * @param 
     * @return 
     **/
    public static function bd09ToWgs84($bd_lon, $bd_lat)
    {
        $gcj02 = self::bd09ToGcj02($bd_lon, $bd_lat);
        $g_lon = $gcj02['lon'];
        $g_lat = $gcj02['lat'];
        $wgs84 = self::gcj02ToWgs84($g_lon, $g_lat);
        return $wgs84;
    }
 
    /**
     * WGS84 转换为 BD09
     * @param 
     * @return 
     **/
    public static function wgs84ToBd09($w_lon,$w_lat){
        $gcj02 = self::wgs84ToGcj02($w_lon,$w_lat);
        $g_lon = $gcj02['lon'];
        $g_lat = $gcj02['lat'];
        $bd09 = self::gcj02Tobd09($g_lon,$g_lat);
        return $bd09;
    } 
 
    /**
     * 转换纬度
     * @param 
     * @return 
     **/
    protected static function transFormLat($lon, $lat)
    {
        $ret = -100.0 + 2.0 * $lon + 3.0 * $lat + 0.2 * $lat * $lat + 0.1 * $lon * $lat + 0.2 * sqrt(abs($lon));
        $ret += (20.0 * sin(6.0 * $lon * self::PI) + 20.0 * sin(2.0 * $lon * self::PI)) * 2.0 / 3.0;
        $ret += (20.0 * sin($lat * self::PI) + 40.0 * sin($lat / 3.0 * self::PI)) * 2.0 / 3.0;
        $ret += (160.0 * sin($lat / 12.0 * self::PI) + 320 * sin($lat * self::PI / 30.0)) * 2.0 / 3.0;
        return $ret;
    }
 
    /**
     * 转换经度
     * @param 
     * @return 
     **/
    protected static function transFormLon($lon, $lat)
    {
        $ret = 300.0 + $lon + 2.0 * $lat + 0.1 * $lon * $lon + 0.1 * $lon * $lat + 0.1 * sqrt(abs($lon));
        $ret += (20.0 * sin(6.0 * $lon * self::PI) + 20.0 * sin(2.0 * $lon * self::PI)) * 2.0 / 3.0;
        $ret += (20.0 * sin($lon * self::PI) + 40.0 * sin($lon / 3.0 * self::PI)) * 2.0 / 3.0;
        $ret += (150.0 * sin($lon / 12.0 * self::PI) + 300.0 * sin($lon / 30.0 * self::PI)) * 2.0 / 3.0;
        return $ret;
    }
}
?>


扫描二维码推送至手机访问。

版权声明:本文由小刚刚技术博客发布,如需转载请注明出处。

本文链接:https://blog.bitefu.net/post/461.html

分享给朋友:

相关文章

Chrome 开启自带多线程下载支持 提升下载速度

Chrome 开启自带多线程下载支持 提升下载速度

Chrome 是当前性能最高的浏览器之一,但在下载速度方面,Chrome 可一点都不快。当前 Chrome 的下载仍是单线程机制,速度非常慢。对此,大家不得不使用第三方的下载工具,例如迅雷、IDM 等支持多线程下载的下载软件,来提高下载速度...

linux centos svn搭建及网站同步更新

以前上线的项目都是用ftp来更新的.后来技术人员多了,ftp的弊端就出来了.所以在这里给大家介绍一下我们现在正在使用的svn同步修改更新的方案.原理很简单主要是用到了svn的hooks功能.即本地提交到服务器的svn服务器.然后svn服务器...

360浏览器如何禁止切换到兼容模式[本地终极办法]

360浏览器如何禁止切换到兼容模式[本地终极办法]

浏览器内核控制标签meta说明背景介绍由于众所周知的原因,国内的主流浏览器都是双核浏览器:基于Webkit的内核用于常用网站的高速浏览,基于IE的内核主要用于部分网银、政府、办公系统等网站的正常使用。以360浏览器为例,我们优先通过Webk...

安装Windows 10X 教你如何安装Win10X正式版 及下载地址

安装Windows 10X 教你如何安装Win10X正式版 及下载地址

安装Windows 10X 教你如何安装Win10X正式版:Windows 10X是Windows 10操作系统的新版本,主要针对双屏电脑。由于即将运行Windows 10X的双屏电脑(例如即将面世的Surface Neo)的开发遇到挫折,...

VirtualXposed 不支持32位应用 32位无法安装问题解决办法

VirtualXposed 不支持32位应用 32位无法安装问题解决办法

VirtualXposed 是基于VirtualApp 和 epic 在非ROOT环境下运行Xposed模块的实现(支持5.0~10.0)。与 Xposed 相比,目前 VirtualXposed 有...

抢先体验太阳谷!20H1、20H2、21H1替换“Dev版新图标”

尽量不要替换shell32.dll.mun和zipfldr.dll.mun,这两老哥可能会带来无法预测的风险!不解除被替换文件的硬链接,已确定会导致无法安装质量更新!(详见H大测评)单替换imageres.dll.mun新图标覆盖率基本可达...

发表评论

访客

看不清,换一张

◎欢迎参与讨论,请在这里发表您的看法和观点。