[发明专利]大地线拱高限差约束的大地线内插方法

摘要

大地线拱高限差约束的大地线内插算法，包括以下步骤：读入大地线起点、终点的经纬度坐标；设定大地线拱高限差；计算未予标记的各折线段大地线长度、起始方位角及经纬度坐标中点；计算该点处的地球椭球面切点动径；计算未予标记的各折线段起点、终点的空间直角坐标；计算未予标记的各折线段大地线拱高；顺次判断未予标记的折线段大地线拱高与大地线拱高限差的关系，若折线段大地线拱高大于大地线拱高限差，则对折线段进行对分加密内插处理；若折线段大地线拱高小于等于大地线拱高限差，则该折线段不予内插处理。本发明通过大地线展绘精度与内插距离的自适应调整与匹配，实现了在指定精度阈值条件下任意长度大地线的精密展绘。

主权项

大地线拱高限差约束的大地线内插方法，其特征在于：包括以下步骤：a、输入大地线起点(P)的经纬度坐标(B1,L1)和终点(Q)的经纬度坐标(B2,L2)；b、确定起点(P)、终点(Q)在地球椭球面上构成的大地线(w)拱高限差(hΩ)；c、计算折线w'中未予标记的各折线段w′i(i∈[0,n])(w′0＝PT1、w′i＝TiTi+1、w′n＝TnQ)的大地线长度S(w′i)、起始方位角(Ai)及经纬度坐标中点起点(P)、终点(Q)及内插点Ti(i∈[1,n])相连构成折线(w')，n表示大地线w上的内插点数，n＝0,1,2...；d、依据折线(w')中未予标记的各折线段w′i(i∈[0,n])的经纬度坐标中点，计算折线(w')中未予标记的折线段(w′i)经纬度坐标中点处的地球椭球面切点动径(ri)；e、以球面代替地球椭球面，在替代球面上计算折线(w')中未予标记的折线段(w')大地线拱高(hi)；f、顺次判断折线(w')中未予标记的折线段w′i(i∈[0,n])大地线拱高(hi)与大地线拱高限差(hΩ)的关系，若hi＞hΩ，则内插出该折线段的大地线中点，遍历所有折线段后转至步骤c；若hi≤hΩ，则该折线段不予内插处理。