PHP收藏本版 +发表新主题

php各个坐标系转换

  1. <?php

  2. namespace Common\Controller;

  3. /**
  4. * Description of Map
  5. *
  6. */
  7. class MapController
  8. {
  9.     const x_PI = 52.35987755982988;
  10.     const PI   = 3.1415926535897932384626;
  11.     const a    = 6378245.0;
  12.     const ee   = 0.00669342162296594323;

  13.         /**
  14.          * WGS84转GCj02(北斗转高德)
  15.          * @param lng
  16.          * @param lat
  17.          * @returns {*[]}
  18.          */
  19.     public function wgs84togcj02($lng, $lat)
  20.     {
  21.         if ($this->out_of_china($lng, $lat)) {
  22.             return array($lng, $lat);
  23.         } else {
  24.             $dlat      = $this->transformlat($lng - 105.0, $lat - 35.0);
  25.             $dlng      = $this->transformlng($lng - 105.0, $lat - 35.0);
  26.             $radlat    = $lat / 180.0 * self::PI;
  27.             $magic     = sin($radlat);
  28.             $magic     = 1 - self::ee * $magic * $magic;
  29.             $sqrtmagic = sqrt($magic);
  30.             $dlat      = ($dlat * 180.0) / ((self::a * (1 - self::ee)) / ($magic * $sqrtmagic) * self::PI);
  31.             $dlng      = ($dlng * 180.0) / (self::a / $sqrtmagic * cos($radlat) * self::PI);
  32.             $mglat     = $lat + $dlat;
  33.             $mglng     = $lng + $dlng;
  34.             return array($mglng, $mglat);
  35.         }
  36.     }
  37.         /**
  38.          * GCJ02 转换为 WGS84 (高德转北斗)
  39.          * @param lng
  40.          * @param lat
  41.          * @return array(lng, lat);
  42.          */
  43.     public function gcj02towgs84($lng, $lat)
  44.     {
  45.         if ($this->out_of_china($lng, $lat)) {
  46.             return array($lng, $lat);
  47.         } else {
  48.             $dlat      = $this->transformlat($lng - 105.0, $lat - 35.0);
  49.             $dlng      = $this->transformlng($lng - 105.0, $lat - 35.0);
  50.             $radlat    = $lat / 180.0 * self::PI;
  51.             $magic     = sin($radlat);
  52.             $magic     = 1 - self::ee * $magic * $magic;
  53.             $sqrtmagic = sqrt($magic);
  54.             $dlat      = ($dlat * 180.0) / ((self::a * (1 - self::ee)) / ($magic * $sqrtmagic) * self::PI);
  55.             $dlng      = ($dlng * 180.0) / (self::a / $sqrtmagic * cos($radlat) * self::PI);
  56.             $mglat     = $lat + $dlat;
  57.             $mglng     = $lng + $dlng;
  58.             return array($lng * 2 - $mglng, $lat * 2 - $mglat);
  59.         }
  60.     }

  61.         /**
  62.          * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换
  63.          * 即 百度 转 谷歌、高德
  64.          * @param bd_lon
  65.          * @param bd_lat
  66.          * @returns
  67.          */
  68.     public function bd09togcj02($bd_lon, $bd_lat)
  69.     {
  70.         $x_pi   = 3.14159265358979324 * 3000.0 / 180.0;
  71.         $x      = $bd_lon - 0.0065;
  72.         $y      = $bd_lat - 0.006;
  73.         $z      = sqrt($x * $x + $y * $y) - 0.00002 * sin($y * $x_pi);
  74.         $theta  = atan2($y, $x) - 0.000003 * cos($x * $x_pi);
  75.         $gg_lng = $z * cos($theta);
  76.         $gg_lat = $z * sin($theta);

  77.         return array($gg_lng, $gg_lat);
  78.     }

  79.         /**
  80.          * GCJ-02 转换为 BD-09 (火星坐标系 转百度即谷歌、高德 转 百度)
  81.          * @param $lng
  82.          * @param $lat
  83.          * @returns array(bd_lng, bd_lat)
  84.          */
  85.     public function gcj02tobd09($lng, $lat)
  86.     {
  87.         $z      = sqrt($lng * $lng + $lat * $lat) + 0.00002 * Math . sin($lat * x_PI);
  88.         $theta  = Math . atan2($lat, $lng) + 0.000003 * Math . cos($lng * x_PI);
  89.         $bd_lng = $z * cos($theta) + 0.0065;
  90.         $bd_lat = z * sin($theta) + 0.006;
  91.         return array($bd_lng, $bd_lat);
  92.     }

  93.     private function transformlat($lng, $lat)
  94.     {
  95.         $ret = -100.0 + 2.0 * $lng + 3.0 * $lat + 0.2 * $lat * $lat + 0.1 * $lng * $lat + 0.2 * sqrt(abs($lng));
  96.         $ret += (20.0 * sin(6.0 * $lng * self::PI) + 20.0 * sin(2.0 * $lng * self::PI)) * 2.0 / 3.0;
  97.         $ret += (20.0 * sin($lat * self::PI) + 40.0 * sin($lat / 3.0 * self::PI)) * 2.0 / 3.0;
  98.         $ret += (160.0 * sin($lat / 12.0 * self::PI) + 320 * sin($lat * self::PI / 30.0)) * 2.0 / 3.0;
  99.         return $ret;
  100.     }
  101.     private function transformlng($lng, $lat)
  102.     {
  103.         $ret = 300.0 + $lng + 2.0 * $lat + 0.1 * $lng * $lng + 0.1 * $lng * $lat + 0.1 * sqrt(abs($lng));
  104.         $ret += (20.0 * sin(6.0 * $lng * self::PI) + 20.0 * sin(2.0 * $lng * self::PI)) * 2.0 / 3.0;
  105.         $ret += (20.0 * sin($lng * self::PI) + 40.0 * sin($lng / 3.0 * self::PI)) * 2.0 / 3.0;
  106.         $ret += (150.0 * sin($lng / 12.0 * self::PI) + 300.0 * sin($lng / 30.0 * self::PI)) * 2.0 / 3.0;
  107.         return $ret;
  108.     }

  109.     private function rad($param)
  110.     {
  111.         return $param * self::PI / 180.0;
  112.     }
  113.         /**
  114.          * 判断是否在国内,不在国内则不做偏移
  115.          * @param $lng
  116.          * @param $lat
  117.          * @returns {boolean}
  118.          */
  119.     private function out_of_china($lng, $lat)
  120.     {
  121.         return ($lng < 72.004 || $lng > 137.8347) || (($lat < 0.8293 || $lat > 55.8271) || false);
  122.     }

  123. }

求验证
扫描二维码,手机查看
声明:本文来源于互联网,观点仅代表作者本人,不代表欢乐你我,真实性请妥善甄别。
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则