js版根据经纬度计算多边形面积（墨卡托投影）

左手的ㄟ右手 2022-05-28 10:12 623阅读 0赞

var earthRadiusMeters = 6371000.0;
    var metersPerDegree = 2.0 * Math.PI * earthRadiusMeters / 360.0;
    var radiansPerDegree = Math.PI / 180.0;
    var degreesPerRadian = 180.0 / Math.PI;
    var pointArr;
    
    $(document).ready(function() {
    	pointArr = new Array();
    	b();
    });
    
    
    function calculateArea(points) {
    	if (points.length > 2) {
    		var areaMeters2 = PlanarPolygonAreaMeters2(points);
    		if (areaMeters2 > 1000000.0) {
    			areaMeters2 = SphericalPolygonAreaMeters2(points);
    			alert("面积为" + areaMeters2 + "平方米");
    		}
    	}
    }
    
    /*球面多边形面积计算*/
    function SphericalPolygonAreaMeters2(points) {
    	var totalAngle = 0;
    	for (var i = 0; i < points.length; i++) {
    		var j = (i + 1) % points.length;
    		var k = (i + 2) % points.length;
    		totalAngle += Angle(points[i], points[j], points[k]);
    	}
    	var planarTotalAngle = (points.length - 2) * 180.0;
    	var sphericalExcess = totalAngle - planarTotalAngle;
    	if (sphericalExcess > 420.0) {
    		totalAngle = points.length * 360.0 - totalAngle;
    		sphericalExcess = totalAngle - planarTotalAngle;
    	} else if (sphericalExcess > 300.0 && sphericalExcess < 420.0) {
    		sphericalExcess = Math.abs(360.0 - sphericalExcess);
    	}
    	return sphericalExcess * radiansPerDegree * earthRadiusMeters * earthRadiusMeters;
    }
    
    /*角度*/
    function Angle(p1, p2, p3) {
    	var bearing21 = Bearing(p2, p1);
    	var bearing23 = Bearing(p2, p3);
    	var angle = bearing21 - bearing23;
    	if (angle < 0) {
    		angle += 360;
    	}
    	return angle;
    }
    
    
    /*方向*/
    function Bearing(from, to) {
    	var lat1 = from[1] * radiansPerDegree;
    	var lon1 = from[0] * radiansPerDegree;
    	var lat2 = to[1] * radiansPerDegree;
    	var lon2 = to[0] * radiansPerDegree;
    	var angle = -Math.atan2(Math.sin(lon1 - lon2) * Math.cos(lat2), Math.cos(lat1) * Math.sin(lat2) - Math.sin(lat1) * Math.cos(lat2) * Math.cos(lon1 - lon2));
    	if (angle < 0) {
    		angle += Math.PI * 2.0;
    	}
    	angle = angle * degreesPerRadian;
    	return angle;
    }
    
    /*平面多边形面积*/
    function PlanarPolygonAreaMeters2(points) {
    	var a = 0;
    	for (var i = 0; i < points.length; ++i) {
    		var j = (i + 1) % points.length;
    		var xi = points[i][0] * metersPerDegree * Math.cos(points[i][1] * radiansPerDegree);
    		var yi = points[i][1] * metersPerDegree;
    		var xj = points[j][0] * metersPerDegree * Math.cos(points[j][1] * radiansPerDegree);
    		var yj = points[j][1] * metersPerDegree;
    		a += xi * yj - xj * yi;
    	}
    	return Math.abs(a / 2);
    }
    
    function b() {
    	var s = "112.523197631836,37.868892669677734;112.5170669555664,37.8605842590332;112.52099609375,37.849857330322266;112.54137420654297,37.8512732521875;112.5351180302734,37.858699798583984";
    	var s1 = new Array()
    	s1 = s.split(";");
    	for (var i = 0; i < s1.length; i++) {
    		var ss = s1[i];
    		var temp = ss.split(",");
    		var point = new Array();
    		point.push(Number(temp[0]), Number(temp[1]));
    		pointArr.push(point);
    	}
    	calculateArea(pointArr);

\}

转载链接：https://www.cnblogs.com/grimm/p/5097383.html

发表评论取消回复

表情：

评论列表（有 0 条评论，623人围观）

还没有评论，来说两句吧...

相关阅读

相关墨卡托投影介绍

一、墨卡托投影墨卡托投影，又称正轴等角圆柱投影，由荷兰地图学家墨卡托(G.Mercator)于1569年创拟。假设地球被套在一个圆柱中，赤道与圆柱相切，然后在地球中心放

蔚落/ 2023年09月30日 13:43/ 0 赞/ 144 阅读

相关 Unity 经纬度转空间坐标及墨卡托坐标转经纬度等

注：此案例中有用到Mapbox相关库 1. 创建由于Vector2d 结构体原因：由于Unity中自带的Vector2是float类型，而我们需要用一个结构体来储存精

拼搏现实的明天。/ 2023年01月23日 13:50/ 0 赞/ 51 阅读

相关经纬度，墨卡托等坐标转换

[1.经纬度转墨卡托][1.] [1.1Java：][1.1Java] [1.2Scala：][1.2Scala] [2.墨卡

雨点打透心脏的1/2处/ 2023年01月13日 01:46/ 0 赞/ 507 阅读

相关横轴墨卡托投影参数_几种投影的主要参数【转】

Gauss Kruger(高斯-克吕格投影)：除中央经线和赤道为直线外，其他经线均为对称于中央经线的曲线。该投影没有角度变形，在长度和面积上变形也很小，中央经线无变形，自中央

冷不防/ 2023年01月09日 12:57/ 0 赞/ 365 阅读

相关墨卡托经纬度坐标转换

在GIS中，常会遇到不同投影坐标系之间的转换，学习了解下面的方法进行实现。 public static double[] Mercator2LonLat(double

水深无声/ 2022年10月17日 04:18/ 0 赞/ 330 阅读

相关理解墨卡托投影原理

http://blog.csdn.net/ryfdizuo/article/details/20793669 背景介绍墨卡托投影思想是GIS的基石，GPS设备采集的数

悠悠/ 2022年09月24日 12:20/ 0 赞/ 568 阅读

相关墨卡托投影、高斯-克吕格投影、UTM投影及我国分带方法

一、墨卡托投影、高斯-克吕格投影、UTM投影 1．墨卡托(Mercator)投影 ![12700d7706eg214.jpg][] 墨卡托(Mercator)投影，是一

本是古典何须时尚/ 2022年07月16日 08:14/ 0 赞/ 395 阅读

相关 js版根据经纬度计算多边形面积（墨卡托投影）

var earthRadiusMeters = 6371000.0; var metersPerDegree = 2.0 Math.PI earthRadi

左手的ㄟ右手/ 2022年05月28日 10:12/ 0 赞/ 624 阅读

相关坐标转换：墨卡托 & 经纬度

坐标转换：墨卡托 & 经纬度 var lonlat={x:114.514864,y:38.04232}; var mercator={x:12747736.3

快来打我*/ 2021年09月29日 16:14/ 0 赞/ 1718 阅读

相关墨卡托投影

微信搜索：“二十同学” 公众号，欢迎关注一条不一样的成长之路一、墨卡托投影墨卡托投影，又称正轴等角圆柱投影，由荷兰地图学家墨卡托(G.Merca

「爱情、让人受尽委屈。」/ 2021年09月29日 01:36/ 0 赞/ 647 阅读