首页 > 编程知识 正文

大地坐标系转换火星坐标系公式,空间坐标转换为大地坐标

时间:2023-05-04 04:14:09 阅读:274586 作者:2065

WGS84:国际坐标系,为一种大地坐标系,也是目前广泛使用的GPS全球卫星定位系统使用的坐标系。
GCJ02:火星坐标系,是由中国国家测绘局制订的地理信息系统的坐标系统。由WGS84坐标系经加密后的坐标系。
BD09:为百度坐标系,在GCJ02坐标系基础上再次加密。其中bd09ll表示百度经纬度坐标,bd09mc表示百度墨卡托米制坐标

上代码

const x_PI = 52.35987755982988; const PI = 3.1415926535897932384626; const a = 6378245.0; const ee = 0.00669342162296594323; 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 * 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); $dlng = ($dlng * 180.0) / (self::a / $sqrtmagic * cos($radlat) * self::PI); $mglat = $lat + $dlat; $mglng = $lng + $dlng; dump($mglng.','.$mglat);//转换后得经纬度 // } } private 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 * self::PI) + 20.0 * sin(2.0 * $lng * 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; } private 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 * self::PI) + 20.0 * sin(2.0 * $lng * self::PI)) * 2.0 / 3.0; $ret += (20.0 * sin($lng * self::PI) + 40.0 * sin($lng / 3.0 * self::PI)) * 2.0 / 3.0; $ret += (150.0 * sin($lng / 12.0 * self::PI) + 300.0 * sin($lng / 30.0 * self::PI)) * 2.0 / 3.0; return $ret; }

版权声明:该文观点仅代表作者本人。处理文章:请发送邮件至 三1五14八八95#扣扣.com 举报,一经查实,本站将立刻删除。