当前位置:网站首页>天地图坐标系转高德坐标系 WGS84转GCJ02
天地图坐标系转高德坐标系 WGS84转GCJ02
2022-08-04 01:31:00 【sr_www】
let [lng, lat] = [116.391358342, 39.9049662090001]
console.log([lng, lat].toString())
console.log(wgs84Togcj02(lng, lat).toString())
function wgs84Togcj02(lng, lat) {
if (out_of_china(lng, lat)) {
console.log('国外')
return [lng, lat]
}
//定义一些常量
//GCJ02 转换为 WGS84
const x_PI = 3.14159265358979324 * 3000.0 / 180.0
const PI = 3.1415926535897932384626
const a = 6378245.0
const ee = 0.00669342162296594323
let dlat = transformlat(lng - 105.0, lat - 35.0)
let dlng = transformlng(lng - 105.0, lat - 35.0)
let radlat = lat / 180.0 * PI
let magic = Math.sin(radlat)
magic = 1 - ee * magic * magic
let sqrtmagic = Math.sqrt(magic)
dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI)
dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * PI)
var mglat = lat + dlat
var mglng = lng + dlng
return [mglng, mglat]
//return [lng * 2 - mglng, lat * 2 - mglat]
}
function out_of_china(lng, lat) {
return (lng < 72.004 || lng > 137.8347) || ((lat < 0.8293 || lat > 55.8271) || false)
}
function transformlat(lng, lat) {
const x_PI = 3.14159265358979324 * 3000.0 / 180.0
const PI = 3.1415926535897932384626
const a = 6378245.0
const ee = 0.00669342162296594323
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 * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0
ret += (20.0 * Math.sin(lat * PI) + 40.0 * Math.sin(lat / 3.0 * PI)) * 2.0 / 3.0
ret += (160.0 * Math.sin(lat / 12.0 * PI) + 320 * Math.sin(lat * PI / 30.0)) * 2.0 / 3.0
return ret
}
function transformlng(lng, lat) {
const x_PI = 3.14159265358979324 * 3000.0 / 180.0
const PI = 3.1415926535897932384626
const a = 6378245.0
const ee = 0.00669342162296594323
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 * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0
ret += (20.0 * Math.sin(lng * PI) + 40.0 * Math.sin(lng / 3.0 * PI)) * 2.0 / 3.0
ret += (150.0 * Math.sin(lng / 12.0 * PI) + 300.0 * Math.sin(lng / 30.0 * PI)) * 2.0 / 3.0
return ret
}
边栏推荐
- VR全景拍摄线上展馆,3D全景带你沉浸体验
- Getting started with MATLAB 3D drawing command plot3
- C 学生管理系统_添加学生
- FeatureNotFound( bs4.FeatureNotFound: Couldn‘t find a tree builder with the features you requested:
- Modulo operation (MOD)
- Analysis: What makes the Nomad Bridge hack unique
- LeetCode第三题(Longest Substring Without Repeating Characters)三部曲之三:两次优化
- Observability:你所需要知道的关于 Syslog 的一些知识
- Array_Sliding window | leecode brushing notes
- 114. How to find the cause of Fiori Launchpad routing error by single-step debugging
猜你喜欢
随机推荐
Web3 安全风险令人生畏?应该如何应对?
How to find the cause of Fiori Launchpad routing errors by single-step debugging
字符串的排列
在Activity中获取另一个XML文件的控件
【虚拟户生态平台】虚拟化平台安装时遇到的坑
Google Earth Engine ——利用公开的河流数据计算河流的有效宽度
redis中常见的问题(缓存穿透,缓存雪崩,缓存击穿,redis淘汰策略)
互斥锁、读写锁、自旋锁,以及原子操作指令xaddl、cmpxchg的使用场景剖析
Vant3—— 点击对应的name名称跳转到下一页对应的tab栏的name的位置
typescript55 - generic constraints
ThreadLocal
企业虚拟偶像产生了实质性的价值效益
appium软件自动化测试框架
工程制图复习题(带答案)
即席查询——Presto
C 学生管理系统_添加学生
通用的测试用例编写大全(登录测试/web测试等)
typescript52-简化泛型函数调用
如何通过API接口从淘宝(或天猫店)复制宝贝到拼多多接口代码对接教程
敏捷交付的工程效能治理