计算矩形方位角（两点连线计算方位角）

一、说在前面

方位角的定义是正北方向与两点连线间的夹角，本质是需要两个点，计算两点的xy坐标差值进行atan即可，但地球是个圆，两点间的连线是个弧线不是直线，也导致其连线与正北方向的夹角是始终变化的，因此短距离方位角的计算比较适合使用在平面坐标下计算。

二、计算开始

1、获得两点（GM、Arcgis）

（1）获取中线

一个矩形有四条边五个点（起点和终点在同一点，但是是两点），因此使用GM基于面图层获得方向中线。

（2）获取一点

通过GM，将上一步获得的中线转成两点。

通过Arcgis属性挑选，mod（字段，2)=1，筛选出奇数或偶数，删除。

（3）获取中点

通过Arcgis，基于第一步获取的中线获得中点，数据管理工具-要素-要素折点转点。

（4）获取两个图层的经纬度坐标。

通过Arcgis，数据管理工具-要素-添加XY坐标。

2、计算方位角

（1）将两个图层经纬度数据导出

中点作为原点，边点作为角点。通过dbf得到两点坐标。

（2）通过excel计算方位角

首先将经纬度转换为墨卡托投影坐标
其次通过excel表格计算atan2(注意xy位置）
然后判断象限，增补角度
处理异常角度