CN107229725B - 一种地理坐标向球面三角形离散格网编码的快速转换方法 - Google Patents
一种地理坐标向球面三角形离散格网编码的快速转换方法 Download PDFInfo
- Publication number
- CN107229725B CN107229725B CN201710416259.8A CN201710416259A CN107229725B CN 107229725 B CN107229725 B CN 107229725B CN 201710416259 A CN201710416259 A CN 201710416259A CN 107229725 B CN107229725 B CN 107229725B
- Authority
- CN
- China
- Prior art keywords
- triangle
- control
- point
- grid
- tempx
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
Landscapes
- Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开一种地理坐标向球面三角形离散格网编码的快速转换方法,包括以下步骤:输入经纬度坐标点P以及点P所在地图的比例尺;根据地图比例尺确定格网编码的目标剖分层次L与控制剖分层次LC;控制剖分层次前的递归逼近和控制剖分层次后的直接映射;将最终行列号转换为格网编码,并作为结果输出。本发明利用球面三角形离散格网的几何性质的变化规律,在保证转换精度的前提下,能够明显提高球面三角形离散格网编码与地理坐标的转换效率,且算法时间复杂度随着剖分层次的提高,并不呈线性增长,相较于传统算法具有明显优势。
Description
技术领域
本发明涉及地理信息技术,具体涉及一种地理坐标向球面三角形离散格网编码的快速转换方法。
背景技术
球面三角形离散格网是一种经典的全球离散格网系统。它以三角形为基本格网单元,基于正八面体或正二十面体对球面空间进行同构离散化处理,形成具有层次递归性、单元形态均匀性,无缝无叠的多尺度全球格网系统(Tong,2013)。经过众多学者长期的研究,球面三角形离散格网不仅成为了全球多分辨率空间数据建模的有效工具,还被广泛应用于空间数据索引与可视化(Goodchild,1992;Feket,1990)、地图综合模型(Dutton,1996)、全球导航模型(Lee,2000)、全球环境监测(Dutton,1999)、全球气候模拟(Sadourny,1968;Heikes,1995;1996;Thuburn,1997)等领域。
编码作为球面三角形离散格网单元的唯一数字标识,既表达了空间位置信息,又蕴含了空间尺度信息,为处理全球多分辨率海量空间数据提供了可能,是格网系统的重要组成部分(Dutton,1989,2000;Fekete,1990;Goodchild,1992;Otoo,1993;Lee,2000;Bartholdi,2001)。然而,现有的大部分空间数据是以地理坐标为主,许多传统空间信息系统也是以地理坐标为基础。因此,如何准确、高效的实现球面三角形离散格网编码与地理坐标的双向转换是球面三角形离散格网研究中的核心问题之一(Dutton,1999)。
由于球面是一个各向异性的二维紧致流形空间,在拓扑性质上与平面空间不同胚,属于不可展曲面(Kimerling,1999)。因此,球面空间无法与平面空间一样用相同的格网连续铺砌,难以借助格网原点与目标点的距离和角度关系直接推算出目标点对应的格网编码,通常采用递归逼近的方式逐层建立地理坐标和格网编码的对应关系。在此方面,相关学者基于正八面体球面三角形格网,先后提出了ZOT投影(Dutton,1991)、ETP投影(Goodchild,1992)、SQC(Otoo,1993)、RCA(Zhao,2006)、三向互化(Tong,2006)、EATP投影(Sun,2007)等方法。
这些算法充分利用了基于正八面体的球面三角形格网与地理坐标具有简易转换的特性(Sahr,2003),从不同角度简化了递归逼近时的算法实现难度,提高了双向转换的效率。然而,这些算法的本质仍是基于递归逼近的方法逐层建立地理坐标与格网编码的映射关系,算法时间复杂度与格网剖分层次呈线性关系。随着格网剖分层次的提高,算法效率明显下降,不利于高分辨率空间数据的集成。同时,这些改进方法仅适用于格网结构简单但几何性质较差的正八面体球面三角形格网,难以应用在格网结构复杂但几何性质良好的正二十面体球面三角形格网上,不具有通用性。
通过分析发现,尽管现有算法充分考虑到球面空间的特殊性对球面三角形离散格网系统中格网单元编码与地理坐标双向转换的影响,但却忽略了这一特殊性质在格网剖分层次提高过程中的变化规律。球面三角形离散格网的几何性质并不是固定不变的。球面三角形离散格网在最初的剖分层次中格网单元几何性质确实差异较大,但随着剖分层次的提高,则几何性质趋于相同(White,1998),局部区域的几何性质已近似于平面。若能根据格网几何性质的变化规律,寻找合适的控制剖分层次作为阈值,在低于该阈值的剖分层次中采用可以确保精度但效率较低的递归逼近转换方法,在高于该阈值的剖分层次采用类似平面格网的直接映射方法,则有可能构建一种兼顾精度与效率的新型混合式转换方法。
发明内容
发明目的:本发明的目的在于解决现有技术中存在的不足,提供一种地理坐标向球面三角形离散格网编码的快速转换方法,其转换过程主要包括两个重要部分,控制剖分层次前的递归逼近和控制剖分层次后的直接映射。
技术方案:本发明的一种地理坐标向球面三角形离散格网编码的快速转换方法,包括以下步骤:
(1)输入经纬度坐标点P以及点P所在地图的比例尺;
(2)根据地图比例尺确定格网编码的目标剖分层次L与控制剖分层次LC;
(3)在低于控制剖分层次LC时,为保证转换精度,采用递归逼近法得到点P在控制剖分层次LC上所对应的格网行号X列号Y和对应的控制三角形单元V0V1V2;
(4)当转换到控制剖分层次LC时,球面格网的曲率已趋于平面,故采用类似平面格网的直接映射方法在控制三角形V0V1V2内得到目标剖分层次上对应的最终格网行列号X和Y;
(5)将最终行列号转换为格网编码,并作为结果输出。
进一步的,所述步骤(3)的具体流程如下:
(3.1)将P转换成空间直角坐标,初始化格网编码的行号X、列号Y和剖分层次l:X=0,Y=0,l=0;
(3.2)通过点P与各初始三角形单元边界的关系确定其所在的初始三角形单元(例如:正八面体球面三角形离散格网的初始八个球面三角形,或正二十面体球面三角形离散格网的球面三角形),并用初始三角形单元的顶点与方向初始化控制三角形顶点V0、V1、V2和方向d;
(3.3)若当前控制三角形的剖分层次l等于控制剖分层次LC,则当前行列号即为点P在控制剖分层次上对应的格网行列号,同时得到对应的控制三角形V0V1V2,递归算法结束,进入直接映射;否则,进入步骤(3.4);
(3.4)以当前控制三角形边界V0V1的中点M0、V0V2的中点M2以及球心O构成平面OM0M2,求其三维空间中的平面方程;由于平面OM0M2经过球心O,则可得平面OM0M2的法向量根据法方程直接给出平面OM0M2的点法式方程xN×x+yN×y+zN×z=0;将点P转换为空间直角坐标后,带入平面方程,若结果大于0,则在平面上方,下一层次的控制三角形为当前控制三角形的上部子三角形,进入步骤(3.7);若结果小于0,则进入步骤(3.5);
(3.5)用同样方法判断P是否在平面OM1M0之上,若是,下一层次的控制三角形为当前控制三角形的左部子三角形,进入步骤(3.7),否则进入步骤(3.6);
(3.6)再用同样方法判断P是否在平面OM2M1之上,若是,下一层次的控制三角形为当前控制三角形的右部子三角形;否则,下一层次的控制三角形为当前控制三角形的中部子三角形,进入步骤(3.7)。
(3.7)按照下表中的行列号更新规则,更新行号X和列号Y,同时更新当前控制三角形块的顶点V0、V1、V2和方向d,l=l+1,返回步骤(3.3);:
进一步的,所述步骤(4)中,通过步骤(3)确定控制剖分层次三角形在初始三角形单元内对应的格网行列号X、Y后,在余下的剖分层次N=L-LC间采用直接映射方法获得待求点P所对应的三角形格网行列号:在控制剖分层次中确定的控制三角形单元V0V1V2内,待求点P所处的行列号tempX、tempY可由P与该控制三角形V0V1V2的相对位置求得,具体流程如下:
(4.1)在控制三角形单元V0V1V2内,过V0作V1V2的垂线交V1V2于M1,过待求点P且与V1V2平行的直线交V0M1于N1,计算列方向上的三角形格元高度进而获得控制三角形内的列号
(4.2)计算M2V2方向上的三角形格元高度从而待求点P相对于直线V0V1的偏移数结合点P所在列号tempY,点P所在的三角形格元就定位到两个格元之内,此时tempX=2×L1(左边三角形),或者tempX=2×L1+1(右边三角形);
(4.3)计算M3V1方向上的三角形格元高度从而待求点P相对于点V1的偏移数结合点P所在列号tempY,点P所在的三角形格元就定位到两个格元之内,此时tempX=2×(L2-2N+tempY)+1(左边三角形),或者tempX=2×(L2-2N+tempY)+2(右边三角形);
(4.4)步骤(4.2)和步骤(4.3)各确定两个三角形格元,它们的交集即为待求点P所在的三角形格元;先令tempX=2×L1,若tempX≠2×(L2-2N+tempY)+1且tempX≠2×(L2-2N+tempY)+2,则令tempX=2×L1+1;
(4.5)若控制三角形V0V1V2的方向与初始三角形方向相反,则令tempY=2N-tempY-1,tempX=2×tempY+tempX+1;
(4.6)finY=Y×2N+tempY,行列号finX,finY就为最终确定待求点P在初始格网中所对应的格网行列号。
有益效果:本发明依据三角形格网的几何特征,建立了在控制三角形V0V1V2内待求点P的相对位置与其对应格网编码的直接映射关系,该映射关系通过点P与控制三角形V0V1V2的相对位置关系直接定位到点P所属的目标层次三角形格元,并借助三个不同方向上的偏移量求得该格元的编码。
本发明利用球面三角形离散格网的几何性质的变化规律,在保证转换精度的前提下,能够明显提高球面三角形离散格网编码与地理坐标的转换效率,且算法时间复杂度随着剖分层次的提高,并不呈线性增长,相较于传统算法具有明显优势。
附图说明
图1为本发明的整体流程图;
图2为本发明中格网行列号的转换流程示意图;
图3为实施例中下一层次控制三角形位置判断示意图;
图4为实施例中子三角形位置分布;
图5为实施例中步骤(4.1)示意图;
图6为实施例中步骤(4.2)示意图;
图7为实施例中步骤(4.3)示意图;
图8为本发明中三角形格网的SQC编码与行列号对应关系示意图。
具体实施方式
下面对本发明技术方案进行详细说明,但是本发明的保护范围不局限于所述实施例。
如图1所示,本发明的一种地理坐标向球面三角形离散格网编码的快速转换方法,包括以下步骤:
(1)输入经纬度坐标点P以及点P所在地图的比例尺;
(2)根据地图比例尺确定格网编码的目标剖分层次L与控制剖分层次LC;
(3)在低于控制剖分层次LC时,为保证转换精度,采用递归逼近法得到点P在控制剖分层次LC上所对应的格网行号X列号Y和对应的控制三角形单元V0V1V2;
(4)当转换到控制剖分层次LC时,球面格网的曲率已趋于平面,故采用类似平面格网的直接映射方法在控制三角形V0V1V2内得到目标剖分层次上对应的最终格网行列号X和Y;
(5)将最终行列号转换为格网编码,并作为结果输出。
如图2所示,所述步骤(3)的具体流程如下:
(3.1)将P转换成空间直角坐标,初始化格网编码的行号X、列号Y和剖分层次l:X=0,Y=0,l=0;
(3.2)通过点P与各初始三角形单元边界的关系确定其所在的初始三角形单元(例如:正八面体球面三角形离散格网的初始八个球面三角形,或正二十面体球面三角形离散格网的球面三角形),并用初始三角形单元的顶点与方向初始化控制三角形顶点V0、V1、V2和方向d;
(3.3)若当前控制三角形的剖分层次l等于控制剖分层次LC,则当前行列号即为点P在控制剖分层次上对应的格网行列号,同时得到对应的控制三角形V0V1V2,递归算法结束,进入直接映射;否则,进入步骤(3.4);
(3.4)如图3所示,以当前控制三角形边界V0V1的中点M0、V0V2的中点M2以及球心O构成平面OM0M2,求其三维空间中的平面方程;由于平面OM0M2经过球心O,则可得平面OM0M2的法向量根据法方程直接给出平面OM0M2的点法式方程xN×x+yN×y+zN×z=0;将点P转换为空间直角坐标后,带入平面方程,若结果大于0,则在平面上方,下一层次的控制三角形为当前控制三角形的上部子三角形(如图4所示),进入步骤(3.7);若结果小于0,则进入步骤(3.5);
(3.5)用同样方法判断P是否在平面OM1M0之上,若是,下一层次的控制三角形为当前控制三角形的左部子三角形,进入步骤(3.7),否则进入步骤(3.6);
(3.6)再用同样方法判断P是否在平面OM2M1之上,若是,下一层次的控制三角形为当前控制三角形的右部子三角形;否则,下一层次的控制三角形为当前控制三角形的中部子三角形,进入步骤(3.7)。
(3.7)按照下表中的行列号更新规则,更新行号X和列号Y,同时更新当前控制三角形块的顶点V0、V1、V2和方向d,l=l+1,返回步骤(3.3);:
所述步骤(4)中,通过步骤(3)确定控制剖分层次三角形在初始三角形单元内对应的格网行列号X、Y后,在余下的剖分层次N=L-LC间采用直接映射方法获得待求点P所对应的三角形格网行列号:在控制剖分层次中确定的控制三角形单元V0V1V2内,待求点P所处的行列号tempX、tempY可由P与该控制三角形V0V1V2的相对位置求得,具体流程如下:
(4.1)如图5所示,在控制三角形单元V0V1V2内,过V0作V1V2的垂线交V1V2于M1,过待求点P且与V1V2平行的直线交V0M1于N1,计算列方向上的三角形格元高度进而获得控制三角形内的列号
(4.2)如图6所示,同理,计算M2V2方向上的三角形格元高度 从而待求点P相对于直线V0V1的偏移数结合点P所在列号tempY,点P所在的三角形格元就为两个粗线格元之一,此时tempX=2×L1(左边三角形),或者tempX=2×L1+1(右边三角形);
(4.3)如图7所示,同理,计算M3V1方向上的三角形格元高度 从而待求点P相对于点V1的偏移数结合点P所在列号tempY,点P所在的三角形格元就为两个粗线格元之一,此时tempX=2×(L2-2N+tempY)+1(左边三角形),或者tempX=2×(L2-2N+tempY)+2(右边三角形);
(4.4)步骤(4.2)和步骤(4.3)各确定两个三角形格元,它们的交集即为待求点P所在的三角形格元;先令tempX=2×L1,若tempX≠2×(L2-2N+tempY)+1且tempX≠2×(L2-2N+tempY)+2,则令tempX=2×L1+1;
(4.5)若控制三角形V0V1V2的方向与初始三角形方向相反,则令tempY=2N-tempY-1,tempX=2×tempY+tempX+1;
(4·6)finX=LX/2」×2N+1+tempX,finY=Y×2N+tempY,行列号finX,finY就为最终确定待求点P在初始格网中所对应的格网行列号。
最后,根据8所示的SQC编码与行列号的对应关系,将格网行列号转换为SQ编码,输出结果。
实施例1:
本实施例中,以南京师范大学仙林校区地理科学学院的经纬度坐标P(118.9139003°E,32.1145815°N)转换成1:1万比例尺地图(目标剖分层次为22,控制剖分层次为7)对应的正二十面体三角形离散格网SQC编码为例。
本实施例的转换过程分为两部分,控制剖分层次前的递归逼近和控制剖分层次后的直接映射。
一、递归逼近过程步骤如下:
步骤1:将P(118.9139003°E,32.1145815°N)转换成空间直角坐标P(-2620.8871131827618,4745.0076361872807,3402.3305648447690),初始化格网编码的行、列号和剖分层次:X=0,Y=0,l=0;
步骤2:通过点P与各初始三角形单元边界的关系确定其所在的初始三角形单元编号为2,并用编号为2的初始三角形单元的顶点与方向初始化控制三角形顶点V0、V1、V2和方向d;
步骤3:若当前控制三角形的剖分层次l等于控制剖分层次LC,则当前行列号即为点P在控制剖分层次上对应的格网行列号,同时得到对应的控制三角形V0V1V2,递归算法结束,进入直接映射;否则,进入下一步;
步骤4:如图2所示,以当前控制三角形边界V0V1的中点M0、V0V2的中点M2以及球心O构成平面OM0M2,求其三维空间中的平面方程。由于平面OM0M2经过球心O,则可得平面OM0M2的法向量根据法方程可以直接给出平面OM0M2的点法式方程xN×x+yN×y+zN×z=0。将点P转换为空间直角坐标后,带入平面方程。若结果大于0,则在平面上方,下一层次的控制三角形为当前控制三角形的上部子三角形如图3所示,进入步骤7;若结果小于0,则进入下一步;
步骤5:用同样方法判断P是否在平面OM1M0之上,若是,下一层次的控制三角形为当前控制三角形的左部子三角形,进入步骤7,否则进入下一步;
步骤6:再用同样方法判断P是否在平面OM2M1之上,若是,下一层次的控制三角形为当前控制三角形的右部子三角形;否则,下一层次的控制三角形为当前控制三角形的中部子三角形。进入下一步。
步骤7:按照下表中的规则,更新X和Y,同时更新当前控制三角形块的顶点V0、V1、V2和方向d,l=l+1,返回步骤3。
表1行列号更新规则
步骤3到步骤7是一个递归逼近过程,每次逼近行列号(X,Y)的值如下:(0,0)、(2,1)、(4,3)、(10,7)、(20,15)、(40,31)、(81,63)、(163,126)。
同时得到最终的控制三角形单元V0V1V2,其三个顶点的空间直角坐标分别为V0(-2627.6295552955821,4766.8453899067490,3366.4147024485728)、V1(-2599.5515195261419,4748.9227890212578,3413.2190438417356)、V2(-2649.1708091961154,4724.8164347521324,3408.5192799219058),并且该三角形方向为下。
二、直接映射过程步骤如下:
步骤1:如图4所示,在控制三角形单元V0V1V2内,过V0作V1V2的垂线交V1V2于M1,过待求点P且与V1V2平行的直线交V0M1于N1。计算列方向上的三角形格元高度从而获得控制三角形内的列号
步骤2:同理,如图5所示,计算M2V2方向上的三角形格元高度从而待求点P相对于直线V0V1的偏移数
步骤3:同理,如图6所示,计算M3V1方向上的三角形格元高度从而待求点P相对于点V1的偏移数
步骤4:先令tempX=2×L1=10646,因为tempX≠2×(L2-2N+tempY)+1且tempX≠2×(L2-2N+tempY)+2,则令tempX=2×L1+1=20647。
步骤5:因为控制三角形V0V1V2的方向与初始三角形方向相反,令tempY=2N-tempY-1=6616,tempX=2×tempY+tempX+1=33880。
步骤finY=Y×2N+tempY=4135384,行列号(5342296,4135384)就为最终确定待求点P在初始格网中所对应的格网行列号。
最后,如图8所示的SQC编码与行列号的对应关系,将格网行列号(5342296,4135384)转换为二进制的SQC格网编码(110111010101101000010100100101011001111000000),输出结果。
本发明利用了球面三角形离散格网的几何性质的变化规律,这一规律表明球面三角形离散格网的几何性质是不断变化的,随着剖分层次的提高不断趋于近似。在低剖分层次中,由于格网单元几何性质差异较大,且格网曲率与平面差异明显,因此采用现有的层次递归逼近转换方法,以确保编码与解码的正确性。在高剖分层次中,由于格网单元几何性质近似相同,格网曲率近乎于平面,因此采用这种类似平面格网的直接映射方法,在确保转换精度的前提下,提高转换的效率。
Claims (2)
1.一种地理坐标向球面三角形离散格网编码的快速转换方法,其特征在于:包括以下步骤:
(1)输入经纬度坐标点P以及点P所在地图的比例尺;
(2)根据地图比例尺确定格网编码的目标剖分层次L与控制剖分层次LC;
(3)在当前控制三角形的剖分层次低于控制剖分层次LC时,采用递归逼近法得到点P在控制剖分层次LC上所对应的格网行号X列号Y和对应的控制三角形单元V0V1V2;
(4)当行列号转换到控制剖分层次LC时,球面格网的曲率已趋于平面,采用直接映射方法在控制三角形V0V1V2内得到目标剖分层次上对应的最终格网行列号X和Y;
(5)将最终行列号转换为格网编码,并作为结果输出;
其中,步骤(4)通过步骤(3)确定控制剖分层次三角形在初始三角形单元内对应的格网行列号X、Y后,在余下的剖分层次N=L-LC间采用直接映射方法获得待求点P所对应的三角形格网行列号:在控制剖分层次中确定的控制三角形单元V0V1V2内,待求点P所处的行列号tempX、tempY可由P与该控制三角形V0V1V2的相对位置求得,具体流程如下:
(4.1)在控制三角形单元V0V1V2内,过V0作V1V2的垂线交V1V2于M1,过待求点P且与V1V2平行的直线交V0M1于N1,计算列方向上的三角形格元高度进而获得控制三角形内的列号
(4.2)计算M2V2方向上的三角形格元高度从而待求点P相对于直线V0V1的偏移数结合点P所在列号tempY,点P所在的三角形格元就定位到两个格元之内,此时tempX=2×L1,或者tempX=2×L1+1;
(4.3)计算M3V1方向上的三角形格元高度从而待求点P相对于点V1的偏移数结合点P所在列号tempY,点P所在的三角形格元就定位到两个格元之内,此时tempX=2×(L2-2N+tempY)+1,或者tempX=2×(L2-2N+tempY)+2;
(4.4)步骤(4.2)和步骤(4.3)各确定两个三角形格元,它们的交集即为待求点P所在的三角形格元;先令tempX=2×L1,若tempX≠2×(L2-2N+tempY)+1且tempX≠2×(L2-2N+tempY)+2,则令tempX=2×L1+1;
(4.5)若控制三角形V0V1V2的方向与初始三角形方向相反,则令tempY=2N-tempY-1,tempX=2×tempY+tempX+1;
(4.6)finY=Y×2N+tempY,行列号finX,finY就为最终确定待求点P在初始格网中所对应的格网行列号;
上述过程中,N1为过P点的V1V2平行线与V0M1的交点;M2为V0V1与V2M2的垂足,N2为过P点的V0V1平行线与V2M2的交点;M3为V0V2与V1M3的垂足,N3为过P点的V0V2平行线与V1M3的交点;V0V1垂直于M2V2,M2V2方向是指三角形以V0V1为底时的高方向;V0V2与M3V1垂直,M3V1方向是指三角形以V0V2为底时的高方向;|M2N2|、|V1N3|分别表示M2N2和V1N3的长度。
2.根据权利要求1所述的地理坐标向球面三角形离散格网编码的快速转换方法,其特征在于:所述步骤(3)的具体流程如下:
(3.1)将P转换成空间直角坐标,初始化格网编码的行号X、列号Y和剖分层次l:X=0,Y=0,l=0;
(3.2)通过点P与各初始三角形单元边界的关系确定其所在的初始三角形单元,并用初始三角形单元的顶点与方向初始化控制三角形顶点V0、V1、V2和方向d;
(3.3)若当前控制三角形的剖分层次l等于控制剖分层次LC,则当前行列号即为点P在控制剖分层次上对应的格网行列号,同时得到对应的控制三角形V0V1V2,递归算法结束,进入直接映射;否则,进入步骤(3.4);
(3.4)以当前控制三角形边界V0V1的中点M0、V0V2的中点M2以及球心O构成平面OM0M2,求其三维空间中的平面方程;由于平面OM0M2经过球心O,则可得平面OM0M2的法向量根据法方程直接给出平面OM0M2的点法式方程xN×x+yN×y+zN×z=0;将点P转换为空间直角坐标后,代入平面方程,若结果大于0,则在平面上方,下一层次的控制三角形为当前控制三角形的上部子三角形,进入步骤(3.7);若结果小于0,则进入步骤(3.5);其中,(xN,yN,zN)为法向量的坐标;
(3.5)用同样方法判断P是否在平面OM1M0之上,若是,下一层次的控制三角形为当前控制三角形的左部子三角形,进入步骤(3.7),否则进入步骤(3.6);其中,O为球心,M1为V1V2中点,M0为V0V1中点;
(3.6)再用同样方法判断P是否在平面OM2M1之上,若是,下一层次的控制三角形为当前控制三角形的右部子三角形;否则,下一层次的控制三角形为当前控制三角形的中部子三角形,进入步骤(3.7);
(3.7)按照下表中的行列号更新规则,更新行号X和列号Y,同时更新当前控制三角形的顶点V0、V1、V2和方向d,l=l+1,返回步骤(3.3);
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710416259.8A CN107229725B (zh) | 2017-06-06 | 2017-06-06 | 一种地理坐标向球面三角形离散格网编码的快速转换方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710416259.8A CN107229725B (zh) | 2017-06-06 | 2017-06-06 | 一种地理坐标向球面三角形离散格网编码的快速转换方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107229725A CN107229725A (zh) | 2017-10-03 |
CN107229725B true CN107229725B (zh) | 2019-07-09 |
Family
ID=59934901
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710416259.8A Active CN107229725B (zh) | 2017-06-06 | 2017-06-06 | 一种地理坐标向球面三角形离散格网编码的快速转换方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107229725B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109104611A (zh) * | 2018-07-11 | 2018-12-28 | 深圳超多维科技有限公司 | 一种用于球面与平面间映射的方法和装置 |
CN110084740B (zh) * | 2019-04-25 | 2023-05-23 | 北京工业大学 | 一种基于规格化结构的球面图像生成及转换方法 |
CN112634448B (zh) * | 2020-12-22 | 2022-12-30 | 郑州大学 | 全球离散点云系统驱动的空间格网通用构建方法 |
CN113114272B (zh) * | 2021-04-12 | 2023-02-17 | 中国人民解放军战略支援部队信息工程大学 | 全球瓦片一致的六边形格网数据结构编码方法与装置 |
CN115801024B (zh) * | 2022-11-24 | 2023-09-12 | 西南林业大学 | 局部等距性优化球面格网的编码方法、系统、装置及介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0587285B1 (en) * | 1992-09-11 | 1996-04-17 | Karsten Manufacturing Corporation | Golf ball with novel dimple pattern |
CN101599183A (zh) * | 2008-06-06 | 2009-12-09 | 吴立新 | 一种实现地球空间三维网格剖分的方法 |
CN102567439A (zh) * | 2011-10-09 | 2012-07-11 | 中国人民解放军国防科学技术大学 | Srg剖分编码与地理坐标的转换算法 |
-
2017
- 2017-06-06 CN CN201710416259.8A patent/CN107229725B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0587285B1 (en) * | 1992-09-11 | 1996-04-17 | Karsten Manufacturing Corporation | Golf ball with novel dimple pattern |
CN101599183A (zh) * | 2008-06-06 | 2009-12-09 | 吴立新 | 一种实现地球空间三维网格剖分的方法 |
CN102567439A (zh) * | 2011-10-09 | 2012-07-11 | 中国人民解放军国防科学技术大学 | Srg剖分编码与地理坐标的转换算法 |
Also Published As
Publication number | Publication date |
---|---|
CN107229725A (zh) | 2017-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107229725B (zh) | 一种地理坐标向球面三角形离散格网编码的快速转换方法 | |
CN107273466B (zh) | 球面三角形离散格网编码向地理经纬度坐标的快速转换方法 | |
WO2022047959A1 (zh) | 一种提升风环境的人工智能城市设计形态布局方法 | |
CN104318617B (zh) | 一种面向虚拟应急演练的三维地理场景仿真方法 | |
CN109682381A (zh) | 基于全向视觉的大视场场景感知方法、系统、介质及设备 | |
CN104121903B (zh) | 一种基于边界值问题的滚动航路规划方法 | |
CN113436319B (zh) | 城市室内三维语义模型的异形排列矩阵构建方法及其系统 | |
CN107665242A (zh) | 一种区域空间多尺度网格编码方法及装置 | |
CN109725327A (zh) | 一种多机构建地图的方法及系统 | |
CN108282743A (zh) | 室内定位方法、装置及系统 | |
CN102054166B (zh) | 一种新的用于户外增强现实系统的场景识别方法 | |
CN104469676A (zh) | 一种移动终端的定位方法及系统 | |
CN105844104A (zh) | 基于Morphing变形技术的中尺度涡旋可视化算法 | |
CN108830926A (zh) | 三维虚拟地球软件中空间动态粒子模型建立与可视化方法 | |
CN101582215A (zh) | 一种空间信息的多级九宫格网定位方法 | |
CN108665117A (zh) | 一种室内空间最短路径的计算方法、装置、终端设备以及存储介质 | |
CN105307264A (zh) | 一种无线传感器网络的移动节点定位方法 | |
CN105761310A (zh) | 一种天空可视域数字地图的模拟分析及图像显示方法 | |
CN113505189B (zh) | 一种基于多层级全球格网的面状目标栅格化方法 | |
CN107179525A (zh) | 一种基于泰森多边形的克里金插值的位置指纹构建方法 | |
CN111931934A (zh) | 一种基于改进遗传算法的海量控制点下仿射变换求解方法 | |
CN109741209A (zh) | 台风灾害下配电网多源数据融合方法、系统及存储介质 | |
CN105868393A (zh) | 一种城市中心微热岛数字地图的制作和显示方法 | |
Zhou et al. | Efficient encoding and decoding algorithm for triangular discrete global grid based on Hybrid Transformation Strategy | |
CN205068800U (zh) | 船舶适航预警系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |