专利名称：球面菱形网格递归剖分方法
技术领域：
本文提出一种结合经纬线并且用正多面体逐级递归的思想在球面直接剖分的方 法。
背景技术：
目前，GIS正面临着海量全球空间数据库、全球性问题研究以及位置相关信息社 会化服务等方面的巨大挑战。传统的基于地图的空间信息表达、组织、管理和发布方式已 不能满足全球空间信息管理的需要。随着航天遥感应用研究与理论探索的深人，尤其是在 “数字地球(Digital Earth，DE)”提出后，地理信息系统(GIS)所处理的不只是某一局部区 域的信息，而是全球范围内连续的、多层次的、动态的环境资源和社会信息，传统平面数据 模型逐渐暴露出投影复杂、缺乏多尺度数据集成管理等局限性，不能满足大范围甚至全球 多分辨率海量数据管理的需要。全球空间数据剖分系统(Global Geodata Partitioning System,G2PS)即是在这一背景下提出。其研究如何将地球(或球面)剖分为等面积和等形 状的层次状面片，并实现高效空间数据的表达和管理。目前关于地球的剖分方法，根据剖分单元的形状特征将球面剖分系统分为规则剖 分系统和不规则剖分系统，其中规则剖分系统包括基于地理坐标系的球面格网系统和基于 正多面体剖分的球面格网系统。1. 1基于地理坐标系的球面网格人们最早使用的球面网格为经纬网格(如图1所示)。基于大地坐标系统，将地球 表面按相等经纬度增量划分成网格单元，所有的数据都与这些网格单元有关。它是适合人 们认识地球习惯的一种模型，是现有大量空间数据集、处理算法以及软件的基础，也是目前 GIS软件中处理算法的基础。但是这种格网存在面积变形、形状变形且由赤道向南北极递增，在南北两级网格 退化成三角形而不是矩形，格网单元大小分布不均勻，不利于多分辨率、不同比例尺数据之 间的转换，复杂的邻接关系不利于模拟应用。1.2正多面体的球面网格剖分20世纪80年代末以来，国内外学者对基于正多面体的球面网格剖分研究较多， 其基本方法是以球体的内接正多面体在球面的投影为基础，对球面空间进行连续的递归 剖分，包括正四面体、正六面体、正八面体、正十二面体和正二十面体以及14个半规则立 体(如图2所示)。以Dutton的基于八面体的四分三角形格网(Octahedral-Quaternary Triangular Mesh, 0-QTM)和 Fekete 的基于正二十面体的球面四叉树(Sphere QuadTree, SQT)为代表，其主要优点是对地球表面进行无缝、多级的格网划分，使全球空间数据能忽略 投影的影响，在管理连续的、多层次和多分辨率的动态全球海量数据方面有着传统平面格 网不具备的优势。但是由球面的几何特性可知，区域剖分后产生的各个球面三角形在大小和形状上 都不可避免地存在一定变形，各个球面三角形的上下方向不一致，边长和面积不全相等，即没有一种剖分方法能使得球面各剖分单元具有完全相等的几何特征(如面积、边长、 形状)，只能达到近似相等。特别是基于正多面体的剖分与地理坐标之间的转换比较复杂， 不符合人们的一贯思维。一些学者先后总结了评价网格系统优劣的指标，主要有以下几点(1)网格能以任意分辨率无缝无叠地覆盖全球；(2)不同分辨率的网格能够形成一个高度一致的层次结构；(3)网格系统对应一套有效的编码方案；(4)网格系统与地理坐标之间的转换关系较简单；(5)单元形状和结构完全一致。设计出各项指标都很优秀的网格系统一直是地学领域研究人员孜孜以求的目标。 但是，由于球面具有特殊的几何属性，因此任何一个网格系统都不能同时满足上述条件。而 在实际应用中，对各种指标的要求也不一样，合格的网格系统应该根据需要在各种标准之 间取得平衡。基于地理坐标剖分产生的传统网格系统在理论和技术上起步较早，几乎所有的 GIS软件都采用了这种方式；基于多面体剖分产生的新型网格系统种类繁多，每一种都针 对不同的应用背景，绝大多数成果还仅限于学术研究，鲜有工程化。主要是因为地理坐标符 合人们的思维习惯，计算简单直观且与现有各类数据转换比较方便。基于多面体剖分的全 球网格是一种新型数据结合，尽管在气候模拟、制图综合等分散领域有部分积累，但系统性 和实用性不强，理论体系不够完善
发明内容
基于以上两种方法的优缺点，并权衡以上五项评价指标，本发明提出一种新的剖 分思路——基于地理坐标的球面菱形网格(Sphere Rhombus Grid, SRG)递归剖分方法。1、基于地理坐标的球面菱形网格SRG剖分思想SRG结合以上两种格网系统的剖分方法，充分利用了经纬线，无需内接正多面体， 无需投影，直接在球面用弧线连接经线和纬线，进行多级递归菱形剖分，实现多分辨率层次 性无缝无叠的剖分。2、SRG剖分过程(1)0 级剖分首先用0°、90°、180°、270°经线把地球分成四等份，从0°开始自东向西分分 别编码为0、1、2、3(如错误！未找到引用源。)(2) 1 级剖分以0号球面为例(如错误！未找到引用源。所示)，赤道将其等为上下两部分，赤道 以北的0°、90°经线是其上半部分的两条边，赤道以南的0°、90°经线是其下半部分的 两条边。分别取各边的中点北极点与赤道中间的纬线即北纬45°与上半部分的两条边交 于两点Al、A2，恰好是上半部分两条边的中点；同样南极点与赤道中间的纬线即南纬45° 与下半部分交于两点Bi、B2，恰是下半部分两条边的中点。0°、90°经线中间的45°经线 与赤道的交点01就是0号球面的中心点。用平滑弧线连接点Al和01，使AlOl这条弧线与赤道以北的90°经线平行；用平滑弧线连接点A2和01，使A201这条弧线与赤道以北的0°经线平行；用平滑弧线连接点Bl 和01，使BlOl这条弧线与赤道以南的90°经线平行；用平滑弧线连接点B2和01，使B201 这条弧线与赤道以南的0°经线平行。这样0号球面就被分为四个近似菱形球面，按照上下 左右的顺序分别编码为00、01、02、03。其它球面类似，这样经过一级剖分全球被分为4X4 =16个菱形球面。(3) 2 级剖分对于00菱形球面，左右两个端点A1、A2都在北纬45°纬线上，北纬45°纬线与上 顶点所在纬线(即北极点)的中间纬线——北纬67. 5°与端点Al、A2上方两条边交于两 点，这两点恰是端点Al、A2上方两条边的中点；同样北纬45°纬线与下顶点所在纬线(即 赤道)的中间纬线——北纬22. 5°与端点A1、A2下方两条边交于两点，恰是端点A1、A2下 方两条边的中点。00菱形球面的中心点取左右两个端点Al、A2的连线(即北纬45°线的 一部分)的中点即北纬45°，东经45°的点(错误！未找到引用源。中点P)。00菱形球面各条边的中点与00菱形球面的中心点分别用弧线连接，四条弧线按 照连接的方向或者平行于AlOl或者平行于A201。这样00菱形球面就被近似等分为四个小 菱形，按照上下左右的顺序编码为000、001、002、003(如错误！未找到引用源。)。其它四 边形也可以按照上述方法进行二级剖分，01四边形就被分为010、011、012、013四个小的菱 形球面，这样经过2级剖分，全球可分为4 X 42 = 64个近似菱形球面。(4)多级剖分以此类推，逐层细分，每个菱形球面各边的中点都取其各边两个端点所在两条纬 线的中间纬线与这条边的交点，每个菱形球面的中心点都取该菱形球面左右两个端点的连 线的中点，该连线属于某条纬线的一部分；菱形球面的中心点与各边中点分别连线，形成4 条弧线，这4条弧线有的平行于0°经线的一部分，有的平行于90°经线的一部分，有的平 行于180°经线的一部分，有的平行于270°经线的一部分。这4条弧线将菱形球面剖分成 4个小的菱形球面。经过η级剖分，全球可分为4X4n个菱形球面(三级剖分见错误！未找 到引用源。)。以上提供了剖分的具体方法，但对于全球来说不同地区地形差异很大，而且不同 地区的军事关注度差异也很大，因此，针对不同区域分辨率需求的不同，可采用不同级别的 剖分。3、SRG剖分编码以上剖分过程中已进行了编码，具体编码如下η级剖分产生的菱形球面的编码长度为n+1，编码都是由0、1、2、3中的几位数字 组合而成，每个菱形球面编码从左到右按照从低级剖分到高级剖分顺序排列，不同级别之 间有层次性。每个菱形球面都可分成四个小的菱形球面，四个小的菱形球面按照上下左右 的顺序分别对应0、1、2、3。编码的首位是几就代表该菱形球面面片在几号球面上，以后的 每一位都是在上一级四边形编码的基础上多编一位。假设第k层某个剖分面片的编码为 a0ala2. . .ak，其中al_ak是k层编码，取值为0，1，2，3 ;a0由0级剖分产生，取值是0，1，2， 3。每个剖分面片都有唯一的编码与之对应，编码的长短反映了剖分的层次，剖分单元编码 具有空间位置相关性。例如1号菱形球面经过二级剖分后产生的编码如错误！未找到引用 源。所示。整个编码的过程与DUTTON的QTM编码相似。
5
4、SRG 的优点4. 1菱形剖分单元更合理剖分单元采用菱形，因为菱形结构类似于正方形格网，具有一致的方向性、径向对 称和平移相和性，并不依赖于从多面体表面到球体或椭球体表面的映射方法，可以直接利 用基于平面四叉树的许多算法，因而在空间操作特别是邻近搜索方面更容易实现。另外，菱 形格网还能和一个等面积格网或一个基于其他标准的、最优化的性能格网一起用。4. 2编码与地理坐标之间的转换较简单首先SRG方法借鉴了传统基于地理坐标系的平面网格剖分方法，每个剖分单元的 端点和中心点都对应特定的地理坐标点，剖分过程与经纬线进行密切结合；其次菱形单元 方向一致，中心点直接取左右两个端点连线的中点，该连线为某条纬线的一部分。计算方 便；还有编码按照一定的规则，每个菱形球面面片都与唯一的编码与之对应。4. 3剖分单元变形小本方法结合了地理坐标但是汲取了正多面体的逐级递归的剖分思路，每个面片都 可以以任意分辨率剖分，采用了特定的弧线连接，剖分单元从赤道到两极都是相似的四边 形，变形很小。4. 4无需投影操作简单本方法直接在球面操作，无需内接正多面体，无需投影，操作简单。能以任意分别 率，基本实现了无缝无叠的全球剖分。每个剖分单元的顶点和中心点对应的地理坐标都容易求得，从而使剖分单元的编 码与地理坐标之间的转换相对比较简单
图1是传统球面网格；
图2是正多面体球面网格剖分
图3是SRGO级剖分；
图4是SRGl级剖分；
图5是SRG2级剖分；
图6是SRG3级剖分；
图7是SRG2级剖分编。
具体实施例方式本发明提供了一种新的基于地理坐标的球面菱形网格递归剖分方法，利用经纬 线，直接在球面上用弧线连接经线和纬线，进行0-n级的多级递归的菱形剖分。1、首先进行0级剖分用0°、90°、180°、270°经线把地球分成四等份，从0°经线开始自东向西分别 给这四等份编码0、1、2、3。2、进行1-n级的多级递归菱形剖分3、给每个剖分的菱形球面编码剖分单元边界与经纬度进行了紧密结合，使剖分单元的编码与地理坐标之间的转换相对容易；同时汲取了正多面体逐级递归剖分的思路，使剖分单元从赤道到两级的面积 和形状变形减小，而且可以实现任意分辨率的剖分；另外SRG方法还具有操作简单的特点。
权利要求
球面菱形网格递归剖分方法，其特征在于利用经纬线，直接在球面用弧线连接经线和纬线，进行多级菱形剖分，首先进行0级剖分，然后进行1 n级的多级菱形剖分，并给每个剖分的菱形球面编码。
2.根据权利要求1所述的球面菱形网格递归剖分方法，其特征在于0级剖分的过程为 用0°、90°、180°、270°经线把地球分成四等份，从0°经线开始自东向西分别给这四等 份编码0、1、2、3。
3.根据权利要求1所述的球面菱形网格递归剖分方法，其特征在于多级菱形剖分的过 程为每个菱形球面各边取中点，每个菱形球面取中心点，菱形球面的中心点与各边中点分 别连线，形成4条弧线，这4条弧线将菱形球面剖分成4个小的菱形球面，经过η级剖分，全 球可分为4 X 4η个菱形球面。
4.根据权利要求1所述的球面菱形网格递归剖分方法，其特征在于编码过程如下-M级 剖分产生的菱形球面的编码长度为η+1，编码都是由0、1、2、3中的几位数字组合而成，每个 菱形球面编码从左到右按照从低级剖分到高级剖分顺序排列；每个菱形球面都可分成四个 小的菱形球面，四个小的菱形球面按照上下左右的顺序分别对应0、1、2、3 ；编码的首位是 几就代表该菱形球面面片在几号球面上，以后的每一位都是在上一级四边形编码的基础上 多编一位。
全文摘要
本发明提出一种结合经纬线并且用正多面体逐级递归的思想在球面直接剖分的方法。该方法利用了经纬线，直接在球面用弧线连接经线和纬线，进行多级递归菱形剖分。首先进行0级剖分，用0°、90°、180°、270°经线把地球分成四等份，从0°经线开始自东向西分分别给这四等份编码为0、1、2、3。然后进行1-n级的多级递归菱形剖分。并给每个剖分的菱形球面编码。本方法直接在球面操作，无需内接正多面体，无需投影，操作简单。能以任意分别率，基本实现了无缝无叠的全球剖分。
文档编号G09B29/00GK101964157SQ20101029304
公开日2011年2月2日 申请日期2010年9月27日 优先权日2010年9月27日
发明者孙兆林, 张玉梅, 徐晖, 徐欣, 聂洪山 申请人:中国人民解放军国防科学技术大学