1.一种基于GIS的空间距离校核三维仿真方法,其特征在于,所述方法包括如下步骤:步骤1.预处理GIS中的散乱点,确定点序列;步骤2.根据所述点序列,初始化凸壳链表、边链表及三角链表;步骤3.将未处理的所述散乱点加入三角网;步骤4.优化所述三角网,得到Delaunay三角网;步骤5.根据所述Delaunay三角网,在GIS场景中进行空间距离计算并进行校验;所述步骤1包括:1-1.在散乱点的点集中,筛选出y坐标值最小的点,并判断其数量;若所述y坐标值最小的点只有1个,则该点为选点;若所述y坐标值最小的点大于1个,则选取其中x坐标最小的点为选点;1-2.连接所述选点与所述点集中其他各点,得到有向线段;1-3.分别计算各所述有向线段与正半轴的夹角值和长度;1-4.按所述夹角值的递增以及长度递减进行词典式排序,得到点序列;所述步骤2包括:2-1.在所述点序列中,依次查找不在有向边上的点,将第一个所述有向边添加到凸壳链表和边链表中,令所述凸壳链表中凸壳边数及所述边链表中的有向边数均为i-1;其中,i为点序列总数;2-2.将第二个所述有向边添加到所述边链表中,令所述边链表中的有向边数为i-2;2-3.将新产生的三角形添加到三角形链表中,并设置相关边的相邻三角形信息;令所述三角形链表中的三角形的数目为i-2;2-4.将所述选点添加到所述凸壳链表中;2-5.输出所述凸壳链表、边链表及三角链表,结束初始化;所述步骤3包括:依次判断未处理的所述散乱点是否已经联结成三角网;若是,并继续处理下一点;反之,则在所述凸壳链表中沿最后添加进去的所述有向边位置往前搜索直至所述凸壳链表结束,并修改相关边的相邻三角形信息,令凸壳链表中的凸壳边数+l,并继续处理下一点,直至所有所述散乱点都加入至三角网中。
2.如权利要求1所述的方法,其特征在于,将所述散乱点的点集及预处理后的所述点序列保存在有序点表中。
3.如权利要求1所述的方法,其特征在于,所述步骤4包括:在三角网的联结过程中不断对子三角网进行优化处理、或将点集联结成三角网后再统一对边链表的每条边进行一次局部优化处理;形成所述Delaunay三角网。
4.如权利要求1所述的方法,其特征在于,所述步骤5包括:5-1.按下式(1)计算线路中任一点与地形地物间广义距离d ij (q):其中,q为线路中任一点;x il 为第l个点的第i维坐标;x jl 为l个点的j维坐标;n为点的总数;l为1到n个点中的任意一点;当q=2时,该值为欧氏距离;当q=1时,该值为绝对距离、即城市街坊距离d ij (1):当q→∞时,该值为契比雪夫距离d ij (∞):d ij (∞)=max{|x il -x jl |}l=1,2,...,n (3)根据所述契比雪夫距离求出固定点之间距离,并根据极限逼近的方法求解导线与地形地物间的最小距离;5-2.按下式(4),根据点到直线段距离的计算来确定点到线的距离d PL :其中,L为线状物体;P为点状物体;d Px 为点状物体P到线状物体L上任一点的距离;x为线状物体L上任一点;5-3.根据两个线状物体L1,L2间的距离可以定义为L1中点P1(x1,y1)与L2中点P2(x2,y2)之间距离的极小值,且L1,L2均表现为折线;计算出线-线物体距离d:求得极大值,再计算方差曲线;5-4.设 为角P1的度数; 为角P2的度数;λ1为P1点的极角;λ2为P2点的极角;P为极点,球面三角形中点P1及P2中的以弧度计的两条边和一夹角已知:∠P=λ 1 -λ 2 (8)其中,PP 1 为球面三角形中的一边;PP 2 为球面三角形中的另一边;∠P为球面三角形中的两边夹角;则根据球面三角形的余弦定理,有:其中,P 1 P 2 为弧度,且P 1 P 2 实际长度为R×cos -1 (P 1 P 2 ),R为球面三角形半径;根据球面三角形的余弦定理,有:cosPP 3 =cosPP 1 cosP 1 P 3 +sinPP 1 sinP 1 P 3 cosP 1 (10)其中,P 3 为由P、P 1 、P 3 所组成的球面三角形中的一点, 为P3点的极角;若先给定 则求得 再由球面三角正弦定理求出∠P 3 ;其中,∠P 3 为由P、P 1 、P 3 所组成的球面三角形中的P 3 点的度数;再由耐普尔公式求出P 1 P 3 :5-5.根据peano规则,对于平面上规则分布的正方形格网点,采用旋转函数进行排列次序。
5.如权利要求4所述的方法,其特征在于,所述步骤5-5中的所述排列次序的方式包括:a.按Morton次序和Pi次序将初始空间划分为用四叉树来表示的子空间;b.在Morton次序和Row次序中由空间坐标求得任何一个空间物体的顺序号和地址;c.在Pi次序和Row-Prime次序中,两者权值相等。