1.一种球面三角形离散格网编码向地理经纬度坐标的快速转换方法,其特征在于:依次包括以下步骤:(1)输入格网编码;
(2)根据编码的剖分层次L确定控制剖分层次LC;
(3)在低于控制剖分层次的阶段,采用递归逼近方法确定格网编码在控制剖分层次上所对应的控制三角形单元V0V1V2及目标编码格元在控制三角形V0V1V2内的行列号(X,Y);
(4)当转换到控制剖分层次LC时,采用直接映射法将控制三角形V0V1V2采用等分方式直接划分成格网,而后根据控制三角形内的行列号(X,Y)找到目标三角形格元;
(5)计算目标三角形格元中心点坐标,将其从直角坐标转换到经纬度坐标,输出完成转换;
步骤(3)中通过递归逼近方法获得控制剖分层次上的控制三角形单元V0V1V2后,根据剩余的剖分次数N=L-LC,采用直接映射方法得到目标剖分层次上的三角形单元,参照平面栅格,将控制三角形V0V1V2采用等分方式直接划分成格网,然后根据控制三角形内的行列号(X,Y)找到目标三角形格元,具体流程如下:(4.1)将V0V1、V0V2等分成2N段,在V0V1上分别取第Y和Y+1个等分点e1、e2,在V0V2上分别取第Y,Y+1个等分点e3、e4;
(4.2)将e1e3等分成Y段,e2e4等分成Y+1段,而后在e1e3上分别取第 个等分点d1、d2,在e2e4上分别取第 个等分点d3、d4;
(4.3)若X为偶数,则目标层次上的三角形三个顶点为d1、d3、d4;若为奇数,则三个顶点为d4、d1、d2。
2.根据权利要求1所述的球面三角形离散格网编码向地理经纬度坐标的快速转换方法,其特征在于:步骤(3)中在低于控制剖分层次的阶段时,采用递归逼近法确定格网编码在控制剖分层次上所对应的控制三角形单元V0V1V2目标编码格元在控制三角形V0V1V2内的行列号(X,Y),具体过程为:(3.1)根据格网编码获取对应初始三角形单元,用其顶点初始化控制三角形V0V1V2;
(3.2)初始化格网剖分计数l=0,根据格网编码与行列号的对应关系,将格网编码转换为行列号(X,Y);
(3.3)若当前控制三角形的层次l等于控制剖分层次LC,输出控制三角形单元V0V1V2,及剩余(X,Y),递归方法结束;否则,就进入步骤(3.4);此处剩余(X,Y)是指控制三角形内的行号和列号;
(3.4)计算控制三角形单元三条弧边的中点M0、M1、M2,通过X、Y判断下一层次控制三角形在当前控制三角形中的位置,若Y<2L-l-1,则为当前控制三角形的上部子三角形,进入步骤(3.5);若X<2Y-2L-l,则为当前控制三角形的左部子三角形,令Y=Y-2L-l-1,进入步骤(3.5);若X<2L-l-1,则为当前控制三角形的中部子三角形,令X=X-2×Y+2L-l-1,Y=2L-l-1-Y-1,进入步骤(3.5);若前三个判断条件均不满足,则为当前控制三角形的上部子三角形,令X=X-2L-l,Y=Y-2L-l-1,进入步骤(3.5);
(3.5)更新控制三角形单元的顶点坐标V0、V1、V2,剖分计数l=l+1;返回步骤(3.3)。
3.根据权利要求1所述的球面三角形离散格网编码向地理经纬度坐标的快速转换方法,其特征在于:基于步骤(4.3)所得的三个顶点计算三角形中心点坐标,并将其从直角坐标转换到经纬度坐标,完成转换。