SuperMap iObjects Java 三维中量算方位角

 翔翔

目录

一、背景介绍

二、方位角

三、代码剖析

四、源码

---

一、背景介绍

在三维中量算长度、高度、面积等功能已经不能满足实际项目中的基本需求了，有的可能想要量算一条直线与正北方向的夹角也叫方位角，有的可能想要量算两条线段的夹角等等。

本文从实际需求出发，实现了在三维场景中量算直线的方位角的功能。先看下量算效果：

二、方位角

从某点的指北经线起，依顺时针方向到目标方向线之间的水平夹角（如图所示θ，可以将其看成是指南针所指示的角度），也即是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代码实现绘制时，实时显示方位角即可。详细代码可参考第四节源码

四、源码

https://github.com/songxiang/SuperMap-MeasureHeading

https://gitee.com/bazingaga/SuperMap-MeasureHeading