翔翔
目录
一、背景介绍
在三维中量算长度、高度、面积等功能已经不能满足实际项目中的基本需求了,有的可能想要量算一条直线与正北方向的夹角也叫方位角,有的可能想要量算两条线段的夹角等等。
本文从实际需求出发,实现了在三维场景中量算直线的方位角的功能。先看下量算效果:
二、方位角
从某点的指北经线起,依顺时针方向到目标方向线之间的水平夹角(如图所示θ,可以将其看成是指南针所指示的角度),也即是OPN平面与OPQ平面的所构成的二面角大小。
这里默认地球是一个圆球,通过数学方法推导出θ的计算公式:
三、代码剖析
根据上诉推导出来的公式,编译成Java代码如下:
/**
* 计算方向角
* @param start 起始点坐标
* @param end 终止点坐标
* @return 方位角,以度为单位
*/
private double Complex(Point3D start, Point3D end) {
double numerator = Math.sin(Math.toRadians(end.x - start.x)) * Math.cos(end.y);
double denominator = Math.cos(Math.toRadians(start.y)) * Math.sin(Math.toRadians(end.y))
- Math.sin(Math.toRadians(start.y)) * Math.cos(Math.toRadians(end.y)) * Math.cos(Math.toRadians(end.x - start.x));
double x = Math.atan2(Math.abs(numerator), Math.abs(denominator));
double result = x;
/**
* 判断线的走向,确定方位角
*/
if (end.x > start.x) {
if (end.y > start.y) result = x;
else if (end.y < start.y) result = Math.PI - x;
else result = Math.PI;
} else if (end.x < start.x) {
if (end.y > start.y) result = 2 * Math.PI - x;
else if (end.y < start.y) result = Math.PI + x;
else result = Math.PI * 3 / 2;
} else {
if (end.y > start.y) result = 0;
else if (end.y < start.y) result = Math.PI;
else result = 0;
}
return result * 180 / Math.PI;
}
上诉代码已经能够计算两个经纬度点连线的方向角了,再结合SuperMap iObjects Java代码实现绘制时,实时显示方位角即可。详细代码可参考第四节源码