php 坐标转换_php wgs84转2000-程序员宅基地

技术标签: 坐标转换  高德坐标转换  火星坐标系转换  PHP  百度坐标转换  谷歌坐标转换  

1、首先,这下面是各个坐标系之间的转换

/**
     * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换
     * 即 百度 转 谷歌、高德
     * @param bd_lon
     * @param bd_lat
     * @returns array
     */
    public function bd09togcj02($bd_lon, $bd_lat)
    {
        $x_pi = 3.14159265358979324 * 3000.0 / 180.0;
        $x = $bd_lon - 0.0065;
        $y = $bd_lat - 0.006;
        $z = sqrt($x * $x + $y * $y) - 0.00002 * sin($y * $x_pi);
        $theta = atan2($y, $x) - 0.000003 * cos($x * $x_pi);
        $gg_lng = $z * cos($theta);
        $gg_lat = $z * sin($theta);
        return array($gg_lng, $gg_lat);
    }

    /**
     * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换
     * 即谷歌、高德 转 百度
     * @param lng
     * @param lat
     * @returns array
     */
    public function gcj02tobd09($lng, $lat)
    {
        $z = sqrt($lng * $lng + $lat * $lat) + 0.00002 * sin($lat * $this->x_PI);
        $theta = atan2($lat, $lng) + 0.000003 * cos($lng * $this->x_PI);
        $bd_lng = $z * cos($theta) + 0.0065;
        $bd_lat = $z * sin($theta) + 0.006;
        return array($bd_lng, $bd_lat);
    }

    /**
     * WGS84转GCj02
     * @param lng
     * @param lat
     * @returns array
     */
    public function wgs84togcj02($lng, $lat)
    {
        if ($this->out_of_china($lng, $lat)) {
            return array($lng, $lat);
        } else {
            $dlat = $this->transformlat($lng - 105.0, $lat - 35.0);
            $dlng = $this->transformlng($lng - 105.0, $lat - 35.0);
            $radlat = $lat / 180.0 * $this->PI;
            $magic = sin($radlat);
            $magic = 1 - $this->ee * $magic * $magic;
            $sqrtmagic = sqrt($magic);
            $dlat = ($dlat * 180.0) / (($this->a * (1 - ee)) / ($magic * $sqrtmagic) * $this->PI);
            $dlng = ($dlng * 180.0) / ($this->a / $sqrtmagic * cos($radlat) * $this->PI);
            $mglat = $lat + $dlat;
            $mglng = $lng + $dlng;
            return array($mglng, $mglat);
        }
    }

    /**
     * GCJ02 转换为 WGS84
     * @param lng
     * @param lat
     * @returns array
     */
    public function gcj02towgs84($lng, $lat)
    {
        if ($this->out_of_china($lng, $lat)) {
            return array($lng, $lat);
        } else {
            $dlat = $this->transformlat($lng - 105.0, $lat - 35.0);
            $dlng = $this->transformlng($lng - 105.0, $lat - 35.0);
            $radlat = $lat / 180.0 * $this->PI;
            $magic = sin($radlat);
            $magic = 1 - $this->ee * $magic * $magic;
            $sqrtmagic = sqrt($magic);
            $dlat = ($dlat * 180.0) / (($this->a * (1 - $this->ee)) / ($magic * $sqrtmagic) * $this->PI);
            $dlng = ($dlng * 180.0) / ($this->a / $sqrtmagic * cos($radlat) * $this->PI);
            $mglat = $lat + $dlat;
            $mglng = $lng + $dlng;
            return array($lng * 2 - $mglng, $lat * 2 - $mglat);
        }
    }

    public function transformlat($lng, $lat)
    {
        $ret = -100.0 + 2.0 * $lng + 3.0 * $lat + 0.2 * $lat * $lat + 0.1 * $lng * $lat + 0.2 * sqrt(abs($lng));
        $ret += (20.0 * sin(6.0 * $lng * $this->PI) + 20.0 * sin(2.0 * $lng * $this->PI)) * 2.0 / 3.0;
        $ret += (20.0 * sin($lat * $this->PI) + 40.0 * sin($lat / 3.0 * $this->PI)) * 2.0 / 3.0;
        $ret += (160.0 * sin($lat / 12.0 * $this->PI) + 320 * sin($lat * $this->PI / 30.0)) * 2.0 / 3.0;
        return $ret;
    }

    public function transformlng($lng, $lat)
    {
        $ret = 300.0 + $lng + 2.0 * $lat + 0.1 * $lng * $lng + 0.1 * $lng * $lat + 0.1 * sqrt(abs($lng));
        $ret += (20.0 * sin(6.0 * $lng * $this->PI) + 20.0 * sin(2.0 * $lng * $this->PI)) * 2.0 / 3.0;
        $ret += (20.0 * sin($lng * $this->PI) + 40.0 * sin($lng / 3.0 * $this->PI)) * 2.0 / 3.0;
        $ret += (150.0 * sin($lng / 12.0 * $this->PI) + 300.0 * sin($lng / 30.0 * $this->PI)) * 2.0 / 3.0;
        return $ret;
    }

    /**
     * 判断是否在国内,不在国内则不做偏移
     * @param lng
     * @param lat
     * @returns boolean
     */
    public function out_of_china($lng, $lat)
    {
        return ($lng < 72.004 || $lng > 137.8347) || (($lat < 0.8293 || $lat > 55.8271) || false);
    }




2、还需要设置变量

  private $x_PI = 0;
   private $PI = 3.1415926535897932384626;
   private $a = 6378245.0;
   private $ee = 0.00669342162296594323;

public function __construct()
{
       $this->x_PI = 3.14159265358979324 * 3000.0 / 180.0;
}

3、查出经纬度,经过转换就可以了


版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/YuYan_wang/article/details/72675460

智能推荐

ZZULIOJ 1188: 选票统计(一)(结构体专题)_c++题目描述 某单位进行选举,有5位候选人:zhang、wang、zhao、liu、miao。编写-程序员宅基地

文章浏览阅读1.2k次,点赞3次,收藏4次。zzulioj 1188: 选票统计(一)(结构体专题)题目描述某单位进行选举,有5位候选人:zhang、wang、zhao、liu、miao。编写程序,统计每人所得的票数。要求每人的信息里包括两部分:name和votes,分别描述姓名和所得票数。每个人的信息用一个结构体来表示,5个人的信息使用结构体数组。输入首先输入一个整数n,表示一张选票,接下来n行,每行是一个由小写英文字母组成的字符串,表示以上5个候选人之一。输出输出5行,按zhang、wang、zhao、liu、miao的顺序输出5个候_c++题目描述 某单位进行选举,有5位候选人:zhang、wang、zhao、liu、miao。编写程

使用vlc播放m3u8网络视频教程_vlc怎么播放m3u8视频-程序员宅基地

文章浏览阅读10w+次。例子视频网址https://20200407/1kDcVooU/hls/index.m3u8打开vlc视频播放器点击媒体选择打开网络串流将网址填入点击播放即可_vlc怎么播放m3u8视频

GICI-LIB源码阅读(一)程序简介、编译调试、配置文件、车载数据集-程序员宅基地

文章浏览阅读1.6k次,点赞12次,收藏23次。GICI-LIB 全称GNSS/INS/CameraIntegrated Navigation Library,是上海交大最新开源的一套基于图优化的 GNSS+INS+Camera 集成导航定位库。为了阐明 GNSS 的算法模型,加快在多源融合应用中针对 GNSS 的开发效率,我们开源了 GICI-LIB,并辅以详尽的文档和全面的数据集。GICI-LIB 以可扩展的设计理念,实现了GIC传感器之间多种形式的松紧组合。评估结果表明,GIC 系统能够在多种复杂环境下,提供分米到米级的高精度导航。_gici-lib

MySQL中的SET数据类型详解_mysql set-程序员宅基地

文章浏览阅读856次,点赞3次,收藏4次。创建一个SET类型的列时,你需要定义可能的值集合。这里可以存储任何这些值的组合,包括全部都不选。当向SET类型的列插入数据时,你可以插入定义的任何值的组合。_mysql set

matlab print dmeta,[转载]转:MATLAB导出高清晰图片-程序员宅基地

文章浏览阅读143次。一个写论文用的小技巧,MATLAB产生的图片如何导出来合适大小或者符合编辑部要求的图片:使用print命令,如下:用于LaTeX的话,可以使用eps图片:print -deps fig1用于word的话,使用分辨率比较高的tif或者png格式即可:print -dtiff -r500 fig1产生的图片fig1插入到word中就可以满足要求的哦~~-------------------------..._matlab中print-dbitmap是什么意思

高效能人士七个习惯及提升_高效人士的七个习惯 发展目标-程序员宅基地

文章浏览阅读1.1k次。1. 积极主动:他们知道要为自己的抉择负责,做出的选择总是基于原则和价值观,而不是基于情绪或受限于周围(天气环境,社会天气)。 自我感觉:可以自我选择是否痛苦,当你面对痛苦,囫囵吞枣般咽下这苦口良药,并承认是自己选择了痛苦之后,你的确依然可以选择不痛苦(就当没咽下痛苦这回事,从而更别去因自己的痛苦影响别人)。因此伤害我们的并非悲惨遭遇本身,而是我们对于悲惨遭遇的回应!尽管这些_高效人士的七个习惯 发展目标

随便推点

android 标题字体大小,android setTitle怎么实现字体的大小-程序员宅基地

文章浏览阅读279次。烧仙草VB主题里添加如下,可以修改title的字体的大小和样式比如没有ActionBar的Activity中,定义如下的样式false@style/WindowTitle#fff20spbold如果是在ActionBar存在的前提下,主题样式要如下@style/TextAppearance@color/main_color@style/ActionBar@style/TextAppearance..._安卓settitle 设置字体大小

[Unity基础]打包与读取AssetBundle_解包出来的prefab怎么读取-程序员宅基地

文章浏览阅读3.3k次。原文链接:http://www.xuanyusong.com/archives/24051.打包prefabusing UnityEngine;using System.Collections;using UnityEditor;//选择的物体放在不同的ABpublic class CreateAB_Respectively : MonoBehaviour {_解包出来的prefab怎么读取

关于JVM-三色标记算法剖析-程序员宅基地

文章浏览阅读424次。白色:表示对象未被垃圾收集器访问过,这是可达性分析开始的阶段,所有的对象都是白色,如果分析结束阶段,还是白色,即代表不可达。灰色:表示对象被垃圾收集器访问过,但是至少还有一个引用没有被扫描过。黑色:表示对象已被垃圾收集器访问过了,并且这个对象所有对象都被扫描过了,它是安全存活的,如果有其它对象引用指向了黑色对象,无须重新扫描一遍。黑色对象不可能直接指向某个白色对象。

如何快速调出软键盘_天生我材必有用 | 如何快速的计算和调用防火阀、调节阀、铝合金风口等材料价格...-程序员宅基地

文章浏览阅读291次。应用场景:在项目中,通常会遇到防火阀、调节阀、铝合金风口等材料,已知其规格尺寸,需要通过计算公式,得出其价格,如下表所示。在品茗胜算软件里,在“工料汇总”界面调用我材助手,找到对应材料,点击“风阀风口”,软件内置了计算公式,自动读取规格尺寸,自动算出价格,点击“应用价格”,价格就直接应用到项目中了,省去了翻找公式、敲计算器的时间。当然也可以点击“扫描件”,调出当期这些材料的计算公式。如果..._品茗计价中风口信息价换算在哪里

设计模式系列-第二章(设计原则-SOLID)_设计模式第二章完整代码-程序员宅基地

文章浏览阅读108次。前言:SOLID 原则并非单纯的 1 个原则,而是由 5 个设计原则组成的,它们分别是:单一职责原则、开闭原则、里式替换原则、接口隔离原则和依赖反转原则,依次对应 SOLID 中的 S、O、L、I、D 这 5 个英文字母。一、单一职责原则(SRP) 单一职责原则的英文是 Single Responsibility Principle,缩写为 SRP。一个类或者模块只负责完成一个职责(或者功能)。单一职责原则是为了实现代码高内聚、低耦合,提高代码的复用性、可读性、可维护性。单一职责原则通过..._设计模式第二章完整代码

mysql-workbench双击后无法启动解决方案_workbench双击无反应-程序员宅基地

文章浏览阅读3.3k次,点赞2次,收藏8次。好久没更新了 这个问题我遇到了两次,一次是腾讯云服务器上,一次是我的电脑上,都是安装后点mysql-workbench没有任何反应。网上别的解决方案我就不说了,要是好使你们也不会找到这里。**很简单!你们首先去mysql-workbench得安装目的地哪里,找到这个东西**双击它!看看是不是这个结果然后去下载这个.dll文件之后把这个文件放到C:\Windows\System32 ..._workbench双击无反应

推荐文章

热门文章

相关标签