一种基于GeoSOT网格的星地数据交互方法_许明魁.pdf

一种基于GeoSOT网格的星地数据交互方法摘要：由于卫星数据具有较大的数据量，数据分发传输过程中，现有的影像数据压缩技术和数据分层分块技术存在较多数据冗余，并不适合一些窄带通信，导致数据请求分发效率低下。鉴于此，提出了一种基于GeoSOT剖分网格的星地数据交互方法，以剖分网格为参考框架，利用网格编码对空间区位进行多尺度表达，计算包含关系，离散化的卫星数据基于网格编码实现多尺度按需分发；研究了分发服务模式，以满足不同用户需要；最后通过实验验证了该方法的可行性。关键词：剖分网格；网格编码；数据请求；数据分发中图分类号：P228文献标志码：B文章编号：1672-4623（2023）03-0036-05Interactive Method of Satellite and Earth Data Based on GeoSOT GridXU Mingkui1,YANG Jianpeng2,CHEN Chunhua3,LIU Guanghui3(1.Wuhan Jiwei Spatial Information Technology Research Institute Co.,Ltd.,Wuhan 433070,China;2.Shenzhen Tanglang Technology Co.,Ltd.,Shenzhen 518038,China;3.Hubei Surveying and Mapping Engineering Institute,Wuhan 430074,China)Abstract:Because of the large amount of satellite data,the existing image data compression technology and data layered block technology havemore data redundancy in the process of data distribution and transmission.It s not suitable for narrowband communication,it will cause the prob-lem of inefficient in data request and distribution.In view of this,we put forward a interactive method of satellite and earth data based on Geo-SOT global subdivision grid.Taking GeoSOT global subdivision grid as the reference frame,we used grid coding to express spatial location atmultiple scales and calculate inclusion relation.We distributed the discrete satellite data on demand at multiple scales based on grid coding.Then,we studied the distribution service model,to meet the needs of different users.Finally,we verified the feasibility of this method by experiments.Key words:subdivision grid,grid coding,data request,data distribution收稿日期：2021-11-10；修回日期：2021-12-21。项目来源：精密工程与工业测量国家测绘地理信息局重点实验室开放基金资助项目（PF2017-4）。第一作者简介：许明魁（1985），主要研究方向为卫星遥感数据处理与应用，E-mail：。引文格式：许明魁,杨建朋,陈春华,等.一种基于GeoSOT网格的星地数据交互方法J.地理空间信息,2023,21(3):36-40.doi:10.3969/j.issn.1672-4623.2023.03.008Mar.,2023Vol.21,No.3地 理 空 间 信 息GEOSPATIAL INFORMATION2023 年 3 月第21卷第 3 期（1.武汉吉威空间信息技术研究院有限公司，湖北 武汉 433070；2.深圳市塘朗科技有限公司，广东 深圳 518038；3.湖北省测绘工程院，湖北 武汉 430074）许明魁1，杨建朋2，陈春华3，刘光辉3近年来，随着各行各业对卫星数据需要的不断增长，现有的以“轨道条带”、“景”、“图幅”为单位进行的数据组织与分发对通信条件提出了很高的要求，尤其是对于计算资源和通信链路资源相对较低的手持终端用户来说增加了数据分发的难度1；而用户往往只对影像内的局部区域感兴趣2或只需更新一定时间间隔之后的影像变化信息，目前现有的数据分发模式效率低下3。鉴于此，本文提出了一种基于剖分网格的星地数据交互方法，通过GeoSOT剖分网格组织与管理卫星数据，基于网格编码请求与分发数据，通过不同的服务模式满足各类用户对卫星数据交互的需求。1基于GeoSOT网格编码的卫星数据请求与分发根据不同用户的实际需要，以不同大小的剖分网格单元为信息服务单元，随时随地、并行地为不同任务不同需求的用户提供大到整个地球、小到某个具体目标实体的空间信息灵活调度与分发服务；同时，基于剖分标识编码全球统一、编码唯一的特性，用户能以剖分标识申领空间数据。1.1GeoSOT区位标识的多尺度表达GeoSOT区位标识的多尺度表达特性主要体现在GeoSOT区位标识对空间信息对象区位结构的分层或分支面片标识方式以及区位面片标识编码上4，见图1。1）基于分层面片标识表达的多尺度特性。在GeoSOT区位标识的分层面片标识方式下，一方面可采用不同的剖分层级对空间信息对象进行多层次多尺度区位标识与表达，反映同一空间信息对象在不同剖分层级上的分层区位结构信息；另一方面位于最低剖分层级的分层面片可逐层向上聚合，反映空间信息对第21卷第3期许明魁等：一种基于GeoSOT网格的星地数据交互方法象区位的多尺度信息。2）基于分支面片标识表达的多尺度特性。在GeoSOT区位标识的分支面片标识方式下，一方面体现了同一空间信息对象在不同空间尺度上的多层次区位结构信息；另一方面位于最底层的最小“尺度面片”也可形成对空间信息对象的区域覆盖，逐层向上聚合。另外，GeoSOT区位标识编码也具有层次继承关系，不仅可隐含区位标识的剖分层级，还可反映当前区位面片标识的尺度信息以及由此剖分层级向上的各尺度层信息5。分层面片标识表达分支面片标识表达图1GeoSOT区位标识的多尺度表达特性示意图1.2基于剖分编码的区域包含关系计算在全球剖分框架下，利用面片集合表达地理对象并对其进行剖分编码，可解决对象编码不统一的问题，有利于实现空间信息区域包含关系的计算，从而进一步实现不同层次空间信息之间的关联6。1.2.1基本原理以剖分面片为空间关系计算的基本单位，计算对象空间关系，可分析不同面片集之间的包含关系。其分析对象主要是处于同一级别的剖分面片集合，针对区域范围与尺度差异较大的不同层次空间对象，若将其统一到同一剖分层级进行分析，会产生计算量过大或分析精度损失等问题7。本文以对象表达多尺度剖分编码模型为基础，选用网格划分与编码更加合理的GeoSOT剖分框架，进一步扩展剖分空间关系计算方法，可分析不同级别剖分面片集合之间的空间包含关系，并设计了具体的算法流程。以此为基础，建立跨尺度层次间的空间信息区域包含关系。1.2.2算法与流程为分析不同尺度间对象剖分编码的包含关系，首先需要获取特定地理对象的最小外接矩形（MBR）。其具体方法步骤为：1）根据对象的空间位置，判断编码的层级，并确定编码code：计算实体最小外接矩形R，通过该地理实体的最小/最大经纬度坐标获得，令R的长边为r；确定网格层级，若r小于等于k级网格边长ak并大于k+1级网格边长ak+1，则将k作为基本的剖分层级；得到R所占据的经向和纬向网格数，计算R左下和右上顶点所在的k+1级网格编码，再通过二进制位运算得到R占据的经向和纬向网格数m和n，进而得到覆盖R的网格集合C；聚合为上一级网格，当经向和纬向跨度均为3个网格以上时，网格集合C中至少有 4 个k+1级网格可聚合为k级网格；确定网格层级码L和定位网格码O，若k+1级网格集合C无法聚合为k级网格，则L=k+1，定位网格选取居中并更接近GeoSOT定位原点的k+1级网格，将其编码作为定位网格码O，若网格集合C可聚合为一个或多个k级网格，则L=k，定位网格按首先居中、其次更接近GeoSOT定位原点的原则在k级网格中选取，将其编码作为定位网格码O，定位网格编码的计算可通过C左下角网格编码的二进制位运算得到。2）判断编码长短，按照“异或”操作进行判断。若k层网格的编码位数为n，k层网格的n位和k+1层网格的对应n位一致，则k层包含k+1层，如图2所示。包含CellBCellACellA编码：1001010001110001CellB编码：100101000111000000000000“异或”位运算实现包含关系计算图2异或操作流程示意图3）精确化判断包含关系。网格层级的判断较粗糙，为了更加精细化判断，首先选取精度更高一级的层级网格进行掩码计算，然后确定对象的形状，再利用掩码进行“异或”计算或操作，判断上下层级对象的精确包含关系。首先确定最小外接多边形的区位掩码，然后在区位编码所确定的GeoSOT网格最小外包矩阵内，在自定义的网格尺度下，由单元的具体位置形成的二维数据矩阵即为单元区位掩码（图3），占据的网格数据为1，不占据的网格数据为0。单元区位掩码生成步骤为：在区位编码所确定的网格外包矩形（尺度为L1）中，根据自定义的网格尺度L2，生成该网格外包矩形所对应的L2尺度下的新的网格外包矩形（原来2518的网格，可能变为5036）；在新的网格外包矩形中，将多边形对象坐标串中的每一个坐标都转化为L2尺度下的GeoSOT网格码，并将相应的网格值设37地理空间信息第21卷第3期为1（图4）；依照线填充算法，可找出不动产单元的边界网格，并将这些边界网格值设为1。图3单元区位掩码多边形对象坐标串：(B1,L1),(B2,L2),(B3,L3),图4多边形对象坐标串转化4）在层级包含关系的基础上，利用区位外接矩形的掩码实现对象精确包含关系的分析判断。1.3基于剖分编码的多尺度卫星信息分发地球剖分GeoSOT网格划分了32个层级，具备优良的多尺度特性8。在数据分发中可充分利用GeoSOT剖分网格的多尺度性，实现数据多尺度按需分发，更好地满足不同类型用户多尺度请领、多带宽分发的应用需要。1.4基于剖分编码的小面片卫星信息自动分发用户可根据终端设备的导航定位模块自动获得所在区域的导航定位编码，并定位到剖分网格。以该剖分网格为条件，查询获得相关数据后，用户选择需要传输的数据，对数据进行切分，只将用户所在网格或感兴趣的区域（查询）数据进行网络传输。如果用户连续关注某区域，可将最新影像与用户历史影像进行变化检测，只对变化面片的数据进行传输。由于小面片数据量非常小，受带宽限制小，可实现数据的自动分发。小面片自动分发是在卫星数据地球剖分技术体系下进行的，后端服务数据均按GeoSOT网格进行逻辑剖分和索引9；在通信资源和导航定位资源的支持下，后端服务平台根据终端当前所在位置，随时主动推送用户所在位置的最新小面片数据，实现后端平台对用户终端的随时随地服务支持。传统方法直接传输原始数据，冗余量超过50%，有时甚至多达数倍；多尺度分发方法可根据用户关心区域大小，按网格进行剖分，只传输关心网格数据，冗余量可低于20%，减少带宽浪费，若连续关注相同区域，则可只分发变化的小网格数据，冗余量