《方位角算法改进.doc》由会员分享，可在线阅读，更多相关《方位角算法改进.doc（3页珍藏版）》请在金锄头文库上搜索。
1、方位角又称地平经度(Azimuth (angle)缩写Az)，是在平面上量度物体之间的角度差的方法之一。是从某点的指北方向线起，依顺时针方向到目标方向线之间的水平夹角。 方位角是在地震工作中是测量地震震中位置的一个重要数据，本台用的msdp软件，测出的方位角会有偏差，为此，在实际工作中，我们会参考已经给出的震中位置的经纬度来反算出精确的方位角，作为工作之用。如图一所示，设A点本台位置，B点为震中位置，那么CAB 图一 图二就是B点的方位角，可以记作，设A（A，A）， B(B, B)， D(D，D)，表示纬度和经度值，均为弧度值，如果是角度值，先转换成弧度值。根据球面三角形余弦公式有：COSAO
2、B=SINASINB + COSACOSBCOS(B-A)A点到B点的距离: AOB=arcCOSAOB 那么 AOB=AOB180/PI()B点方位角： SINCAB = COSBSIN(B-A)/ 1-COS2AOB 求反函数 CAB=arcSINCAB 那么 CAB=CAB180/PI()如图二在一象限 方位角为 a , 在第二象限 方位角为 180-|a| ,在第三象限 方位角为 180+|a|,在第四象限 方位角为 360-|a|。如果是在平面上，B点的纬度大于A点的纬度，可以判断方位角可能在一或四象限，B点纬度小于A点纬度，可以判断方位角可能在二或三象限。这个结论在平面成立，但是球
3、面上却不成立，因为不同经度上的一个和A点相同纬度值的点对应的方位角不是90度，如本台A点的纬度为26.506，某一条经度线上的一个点纬度是26.506，这个点和A点的连线构成的方位角，不等于90度。不等于90度，就说明这个点和A点的连线不垂直于A点，所以某一条经线上的待测点的纬度值大于或小于26.506，是不能判断他所在的象限的。如果能在任意一条经线上找到一个点，这个点和A点的连线构成的方位角是90度，那么这个点的纬度值大于或小于26.506，就可以判断方位角处在哪个象限了。如图一，设D点和B点同在一条经线上，这个D点和A点的连线构成的方位角CAD是90度，求D点的纬度值。D点的D=B，公式如
4、下：SIND =COS(B-A) SINA /1-SIN2A SIN2(B-A)上式右边的经纬度值都是已知的，那么 D=arcSIND 就求出了 D点的纬度值。结论：BD 方位角 在一或四象限 BD 并且 -78.251BD 并且 （ 101.749B180 或B-78.251 ） 方位角 在一象限 BD 并且 -78.251B101.749 方位角 在三象限BD 并且 （ 101.749B180 或B-78.251 ） 方位角 在二象限举例：一， 时间 震级 纬度 经度 深度 地点 2012-11-08 00:35:49 7.3 14.1 -92.0 30 危地马拉附近海域 2012-10-24 08:45:34 6.6 10.0 -85.4 30 哥斯达黎加这是简单比较两个纬度值算出的方位角：是用本章所述算出的方位角： 上述两个地震是发生在中美洲的极远震，他们的纬度都小于26.506，但方位角应该在一象限，而不应在第二象限。二、纬度 经度 深度 地点 25.5 121.5 30 台湾岛附近这是简单比较两个纬度值算出的方位角：是用本章所述算出的方位角：比较可以看出，方位角差距很大，台湾不可能在相对于本台方位角为178.8的位置。这个位置几乎都在攀枝花的正下方了。
