坐标系转换及判断点是否在多边形内

最新推荐文章于 2024-04-02 09:48:51 发布

zhou_bo88

最新推荐文章于 2024-04-02 09:48:51 发布

阅读量173

点赞数

文章标签： javascript

本文链接：https://blog.csdn.net/zhou_bo88/article/details/123358094

版权

本文介绍了地理坐标系统之间的转换，包括WGS84、GCJ02和BD09坐标系。提供了从WGS84到GCJ02，以及GCJ02到BD09的转换函数，并详细解释了坐标偏移的原理。此外，还包含了一个判断点是否在多边形内的辅助函数。

摘要由CSDN通过智能技术生成

/*
坐标转换程序
WGS84坐标系：即地球坐标系，国际上通用的坐标系。Earth
GCJ02坐标系：即火星坐标系，WGS84坐标系经加密后的坐标系。Mars
BD09坐标系：即百度坐标系，GCJ02坐标系经加密后的坐标系。 Bd09
搜狗坐标系、图吧坐标系等，估计也是在GCJ02基础上加密而成的。
百度地图API 百度坐标
腾讯搜搜地图API 火星坐标
搜狐搜狗地图API 搜狗坐标*
阿里云地图API 火星坐标
图吧MapBar地图API 图吧坐标
高德MapABC地图API 火星坐标
灵图51ditu地图API 火星坐标
*/

var x_earth_pi=3.14159265358979324 * 3000.0 / 180.0;

//火星坐标（GCJ02）转百度坐标（bd09II）
function marsTobaidu(mars_point){
   var baidu_point={lon:0,lat:0};
   var x=mars_point.lon;
   var y=mars_point.lat;
   var z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * x_earth_pi);
   var theta = Math.atan2(y, x) + 0.000003 * Math.cos(x * x_earth_pi);
   baidu_point.lon = z * Math.cos(theta) + 0.0065;
   baidu_point.lat = z * Math.sin(theta) + 0.006;
   return baidu_point;
}

//百度坐标（bd09II）转火星坐标（GCJ02）
function baiduTomars(baidu_point){
var mars_point={lon:0,lat:0};
var x=baidu_point.lon-0.0065;
var y=baidu_point.lat-0.006;
var z=Math.sqrt(x*x+y*y)- 0.00002 * Math.sin(y * x_earth_pi);
var theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * x_earth_pi);
mars_point.lon=z * Math.cos(theta);
mars_point.lat=z * Math.sin(theta);
return mars_point;
}

// 地球坐标系（WGS84）转火星坐标系（GCJ02）
var earth_pi = 3.14159265358979324;
/*判断是否在国内，不在国内则不做偏移*/
function outOfChina(lon, lat)
{
if ((lon < 72.004 || lon > 137.8347)&&(lat < 0.8293 || lat > 55.8271)){
return true;
}else {
return false;
}
}
function transformLat(x,y)
{
var ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * Math.sqrt(Math.abs(x));
ret += (20.0 * Math.sin(6.0 * x * earth_pi) + 20.0 * Math.sin(2.0 * x * earth_pi)) * 2.0 / 3.0;
ret += (20.0 * Math.sin(y * earth_pi) + 40.0 * Math.sin(y / 3.0 * earth_pi)) * 2.0 / 3.0;
ret += (160.0 * Math.sin(y / 12.0 * earth_pi) + 320 * Math.sin(y * earth_pi / 30.0)) * 2.0 / 3.0;
return ret;
}

//火星坐标系GCJ02转地球坐标系WGS84
function transformGCJ2WGS(gcjLat, gcjLon) {
   let d = delta(gcjLat, gcjLon)
   return {
   'lat': gcjLat - d.lat,
   'lon': gcjLon - d.lon
   }
}
function delta(lat, lon) {
let a = 6378245.0 // a: 卫星椭球坐标投影到平面地图坐标系的投影因子。
let ee = 0.00669342162296594323 // ee: 椭球的偏心率。
let dLat = transformLat(lon - 105.0, lat - 35.0)
let dLon = transformLon(lon - 105.0, lat - 35.0)
let radLat = lat / 180.0 * earth_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) * earth_pi)
dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * earth_pi)
return {
'lat': dLat,
'lon': dLon
}
}
function transformLat(x, y) {
let ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * Math.sqrt(Math.abs(x))
ret += (20.0 * Math.sin(6.0 * x * earth_pi) + 20.0 * Math.sin(2.0 * x * earth_pi)) * 2.0 / 3.0
ret += (20.0 * Math.sin(y * earth_pi) + 40.0 * Math.sin(y / 3.0 * earth_pi)) * 2.0 / 3.0
ret += (160.0 * Math.sin(y / 12.0 * earth_pi) + 320 * Math.sin(y * earth_pi / 30.0)) * 2.0 / 3.0
return ret
}
function transformLon(x, y) {
let ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * Math.sqrt(Math.abs(x))
ret += (20.0 * Math.sin(6.0 * x * earth_pi) + 20.0 * Math.sin(2.0 * x * earth_pi)) * 2.0 / 3.0
ret += (20.0 * Math.sin(x * earth_pi) + 40.0 * Math.sin(x / 3.0 * earth_pi)) * 2.0 / 3.0
ret += (150.0 * Math.sin(x / 12.0 * earth_pi) + 300.0 * Math.sin(x / 30.0 * earth_pi)) * 2.0 / 3.0
return ret
}

//判断点是否在多边形内
function isInPolygon(checkPoint, polygonPoints) {
var counter = 0;
var i;
var xinters;
var p1, p2;
var pointCount = polygonPoints.length;
p1 = polygonPoints[0];

for (i = 1; i <= pointCount; i++) {
p2 = polygonPoints[i % pointCount];
if (
checkPoint[0] > Math.min(p1[0], p2[0]) &&
checkPoint[0] <= Math.max(p1[0], p2[0])
) {
if (checkPoint[1] <= Math.max(p1[1], p2[1])) {
if (p1[0] != p2[0]) {
xinters =
(checkPoint[0] - p1[0]) *
(p2[1] - p1[1]) /
(p2[0] - p1[0]) +
p1[1];
if (p1[1] == p2[1] || checkPoint[1] <= xinters) {
counter++;
}
}
}
}
p1 = p2;
}
if (counter % 2 == 0) {
return false;
} else {
return true;
}
}