本发明专利技术提供了一种基于全球高精度DEM数据的快速海陆判断方法，属于地理信息系统领域，先将全球ASTER GDEM数据文件按照一定的空间索引方案进行插值重构，在此基础上对重构后的DEM数据文件进行二值化和压缩并合并成一个文件，最后利用合并后的数据文件进行实时海陆判断。该方法可对海量地理坐标进行快速的海陆判断，可广泛应用于航船舶交通信息坐标异常判断、船舶交通行为可视化与分析等领域。
A fast sea land judgment method based on high precision DEM data
The invention provides a fast sea land judgment method based on global high precision DEM data. It belongs to the field of geographic information system. First, the global ASTER GDEM data file is interpolated and reconstructed according to a certain spatial index scheme. On this basis, the reconstructed DEM data file is two valued and compressed and merged into a file. Finally, we use the merged data files to conduct real-time sea land judgment. This method can be used to judge the sea and land quickly, and can be widely used in the field of navigation information anomaly detection, visualization and analysis of ship traffic behavior.
【技术实现步骤摘要】
一种基于高精度DEM数据的快速海陆判断方法
本专利技术涉及地理信息系统领域，特别是一种基于高精度DEM数据的快速海陆判断方法。
技术介绍
海洋运输由于运量大、成本低、能耗小、污染少等优势，在我国乃至全球的综合交通运输体系中都占有非常重要的地位。航运管理信息化和智能化系统的相关研究一直是国际航运界关注的热点问题，而实时获得航船舶交通信息是构建智能航运系统的关键。船舶自动识别系统(AIS)是目前获取船舶交通信息的一种重要手段，它主要由岸基设施和船载设备共同组成，是一种广播式的船－船、船－岸信息自动交互工具，是重要的船舶助航导航设备。海事监管部门通过岸基AIS基站接收船台发送的AIS报文，实现对船舶航行动态的实时跟踪与监控。AIS船台发送的报文通常包括：MMSI编号、船舶类型、船舶在航状态、外形尺寸等静态信息和船舶速度、航向角、经纬度坐标等动态信息。由于AIS报文包含的信息容量非常大，所以近年来成为分析水上交通行为与风险的重要数据来源。然而，AIS缺乏完整的信息验证机制，在实际应用中经常会由于AIS传感器接收GPS信号过程中产生的数据错误、AIS信号传输过程中产生的数据错误、AIS设备自身原因产生的数据错误、AIS网络通信阻塞导致AIS信号发射延迟等问题，导致其动态信息中含有某些错误信息。例如，船舶在正常航行状态下的经纬度坐标有时会异常出现在陆地部分，就是一个比较常见的错误。这种错误的坐标信息会干扰海事监管，导致对水面交通态势的误判，同时也会对以AIS数据为基础的船舶交通行为可视化与分析产生困扰。因此，需要对实时接收的AIS数据中的经纬度坐标进行海陆判断。此外，AIS报文的更新频率非常高，正常情况下，船舶每隔2秒-6分钟就会广播1次AIS信息，大量船舶不停广播，导致海量AIS数据产生。仅我国海事部门搜集近海船舶动态数据，单日数据量峰值就高达千万条，3个月信息总量超过亿条。这就迫切需要提供一种准确高效的计算方法，有效解决大数据情况下的快速海陆判断问题。
技术实现思路
本专利技术的目的是要克服现有技术中利用矢量地图数据进行海陆判断效率低下的问题，提供一种基于高精度DEM数据的快速海陆判断方法，可有效解决大数据情况下的快速海陆判断问题。为达到上述目的，本专利技术是按照以下技术方案实施的：一种基于高精度DEM数据的快速海陆判断方法，包括以下步骤：S1、对全球ASTERGDEM数据文件进行插值重构获得重构后的全球DEM数据文件：S11、将全球平均分成4096×2048块，每100×100块合并为一个新的DEM文件，块内顶点个数为256×256；S12、在每个块内对256×256个顶点的高程值进行双线性插值，在用每个块内的四个原始顶点进行双线性插值的该过程中，一旦有两个或两个以上原始顶点为无效像素，则认为待插值点为无效像素点，直接将其高程值记录为-9999，以做标识；S13、将双线性插值后的每100×100块存储为一个新的DEM数据文件，后缀名为demx，全球共计861个这样的文件，每个demx文件分为三部分：文件头、块索引、块数据；其中，文件头存储数据的版本号和文件内的分块数目；每个块索引存储该块数据的字节总数和对应的块数据在文件的偏移地址；每个块数据存储该块的元信息和块内所有顶点的高程值；S2、对插值重构后的全球DEM数据文件进行二值化、压缩、合并：S21、二值化：根据步骤S1中每个demx文件的块数据中顶点的高程值来判断该块的地形特征，分为四种情况：海洋、陆地、海陆接边、未知，依次用整数标识记为0、1、2、3；对于海陆接边块，对海陆接边块内每个顶点继续进行海陆二值化，如果高程值小于等于0，则认为该顶点位于海洋，标识为0；否则，认为其位于陆地，标识为1；S22、压缩：对于海陆接边块中的每个顶点，每个顶点的标识占用1个bit，每8个相邻的顶点占用1个字节，然后再对二值化后的海陆接边块数据进行Zip压缩；S23、合并：将经过二值化和压缩后的861个demx文件合并为一个全球DEM数据文件；S3、利用合并后的全球DEM数据文件进行实时海陆判断：S31、将步骤S2中经过二值化、压缩、合并的全球DEM数据文件存储到计算机的硬盘，在计算机内存中构建一级索引和二级索引，其中一级索引用于记录从文件中读取的所有块的类型标识，每块占据3个字节，即24个bit位；二级索引用于记录所有海陆接边块的原始索引信息，即块的字节数和块在文件中的偏移地址；根据块的地形类型，一级索引情况如下：如果块为海洋，一级索引的24个bit位中最后3位记为000，其它位均为0，即3个字节从高位到低位依次为0、0、0；如果块为陆地，一级索引的24个bit位中最后3位记为010，其它位均为0，即3个字节从高位到低位依次为0、0、2；如果块为未知，一级索引的24个bit位中最后3位记为100，其它位均为0，即3个字节从高位到低位依次为0、0、4；如果块为海陆接边，一级索引的24个bit位中最后1位记为1，其它23个bit位用于记录该块在二级索引中的序号，即该块在所有海陆接边中的自然顺序；S32、在内存中构建完两级索引之后，进行快速海陆判断：S321、根据输入点的经纬度坐标，计算其所在的地形块，并从内存中获取该块的一级索引信息；S322、解析一级索引的3个字节，如果从高位到低位依次为0、0、0，则说明该块为海洋，从而知该点位于海洋，本次判断结束；如果从高位到低位依次为0、0、2，则说明该块为陆地，从而知该点位于陆地，本次判断结束；如果从高位到低位依次为0、0、4，则说明该块为未知，从而知该点地形未知，本次判断结束；否则，该块为海陆接边，该点位于海洋或者陆地，继续进行步骤③；S323、根据海陆接边块的全球编号，到缓存中查找是否有该块解压后的海陆二值化数据，如果有，直接转向步骤④；否则，根据一级索引中前面23个bit位的值，反算其在二级索引中的序号，并根据该序号提取其在内存中的二级索引信息，根据二级索引信息中记录的文件偏移地址，到文件对应位置读取该块数据并解压为海陆二值化数据，然后存储到缓存中；S324、根据该点的经纬度坐标，计算其在该海陆接边块中的行列号，根据行列号计算该点位于二值化数据的哪个字节，以及位于该字节的哪个bit位，如果该bit位为0，表示该点位于海洋，否则，表示该点位于陆地，本次判断结束。与现有技术相比，本专利技术可对海量地理坐标进行快速的海陆判断，可广泛应用于航船舶交通信息坐标异常判断、船舶交通行为可视化与分析等领域。附图说明图1为本专利技术按照全球空间索引方案重构的demx文件格式。图2为本专利技术合并后的全球demx文件格式。图3为本专利技术合并后的全球demx文件栅格化结果。图4为本专利技术基于两级索引的快速海陆判断算法流程图。具体实施方式下面结合具体实施例对本专利技术作进一步描述，在此专利技术的示意性实施例以及说明用来解释本专利技术，但并不作为对本专利技术的限定。1、对全球ASTERGDEM数据文件进行插值重构(1)全球空间索引方案ASTERGDEM数据覆盖范围为北纬83°到南纬83°之间的所有陆地区域，比以往任何地形图都要广得多，达到了地球陆地表面的99％。其空间分辨率为1″，垂直精度达20米，水平精度达30米，是迄今为止可为用户提供的最完整的、精度最高的全球数字高程模型数据。AST本文档来自技高网...
【技术保护点】
一种基于高精度DEM数据的快速海陆判断方法，其特征在于，包括以下步骤：S1、对全球ASTER GDEM数据文件进行插值重构获得重构后的全球DEM数据文件：S11、将全球平均分成4096×2048块，每100×100块合并为一个新的DEM文件，块内顶点个数为256×256；S12、在每个块内对256×256个顶点的高程值进行双线性插值；S13、将双线性插值后的每100×100块存储为一个新的DEM数据文件，后缀名为demx，全球共计861个这样的文件，每个demx文件分为三部分：文件头、块索引、块数据；其中，文件头存储数据的版本号和文件内的分块数目；每个块索引存储该块数据的字节总数和对应的块数据在文件的偏移地址；每个块数据存储该块的元信息和块内所有顶点的高程值；S2、对插值重构后的全球DEM数据文件进行二值化、压缩、合并：S21、二值化：根据步骤S1中每个demx文件的块数据中顶点的高程值来判断该块的地形特征，分为四种情况：海洋、陆地、海陆接边、未知，依次用整数标识记为0、1、2、3；对于海陆接边块，对海陆接边块内每个顶点继续进行海陆二值化，如果高程值小于等于0，则认为该顶点位于海洋，标识为0；否则，认为其位于陆地，标识为1；S22、压缩：对于海陆接边块中的每个顶点，每个顶点的标识占用1个bit，每8个相邻的顶点占用1个字节，然后再对二值化后的海陆接边块数据进行Zip压缩；S23、合并：将经过二值化和压缩后的861个demx文件合并为一个全球DEM数据文件；S3、利用合并后的全球DEM数据文件进行实时海陆判断：S31、将步骤S2中经过二值化、压缩、合并的全球DEM数据文件存储到计算机的硬盘，在计算机内存中构建一级索引和二级索引，其中一级索引用于记录从文件中读取的所有块的类型标识，每块占据3个字节，即24个bit位；二级索引用于记录所有海陆接边块的原始索引信息，即块的字节数和块在文件中的偏移地址；根据块的地形类型，一级索引情况如下：如果块为海洋，一级索引的24个bit位中最后3位记为000，其它位均为0，即3个字节从高位到低位依次为0、0、0；如果块为陆地，一级索引的24个bit位中最后3位记为010，其它位均为0，即3个字节从高位到低位依次为0、0、2；如果块为未知，一级索引的24个bit位中最后3位记为100，其它位均为0，即3个字节从高位到低位依次为0、0、4；如果块为海陆接边，一级索引的24个bit位中最后1位记为1，其它23个bit位用于记录该块在二级索引中的序号，即该块在所有海陆接边中的自然顺序；S32、在内存中构建完两级索引之后，进行快速海陆判断：S321、根据输入点的经纬度坐标，计算其所在的地形块，并从内存中获取该块的一级索引信息；S322、解析一级索引的3个字节，如果从高位到低位依次为0、0、0，则说明该块为海洋，从而知该点位于海洋，本次判断结束；如果从高位到低位依次为0、0、2，则说明该块为陆地，从而知该点位于陆地，本次判断结束；如果从高位到低位依次为0、0、4，则说明该块为未知，从而知该点地形未知，本次判断结束；否则，该块为海陆接边，该点位于海洋或者陆地，继续进行步骤③；S323、根据海陆接边块的全球编号，到缓存中查找是否有该块解压后的海陆二值化数据，如果有，直接转向步骤④；否则，根据一级索引中前面23个bit位的值，反算其在二级索引中的序号，并根据该序号提取其在内存中的二级索引信息，根据二级索引信息中记录的文件偏移地址，到文件对应位置读取该块数据并解压为海陆二值化数据，然后存储到缓存中；S324、根据该点的经纬度坐标，计算其在该海陆接边块中的行列号，根据行列号计算该点位于二值化数据的哪个字节，以及位于该字节的哪个bit位，如果该bit位为0，表示该点位于海洋，否则，表示该点位于陆地，本次判断结束。...
【技术特征摘要】
1.一种基于高精度DEM数据的快速海陆判断方法，其特征在于，包括以下步骤：S1、对全球ASTERGDEM数据文件进行插值重构获得重构后的全球DEM数据文件：S11、将全球平均分成4096×2048块，每100×100块合并为一个新的DEM文件，块内顶点个数为256×256；S12、在每个块内对256×256个顶点的高程值进行双线性插值；S13、将双线性插值后的每100×100块存储为一个新的DEM数据文件，后缀名为demx，全球共计861个这样的文件，每个demx文件分为三部分：文件头、块索引、块数据；其中，文件头存储数据的版本号和文件内的分块数目；每个块索引存储该块数据的字节总数和对应的块数据在文件的偏移地址；每个块数据存储该块的元信息和块内所有顶点的高程值；S2、对插值重构后的全球DEM数据文件进行二值化、压缩、合并：S21、二值化：根据步骤S1中每个demx文件的块数据中顶点的高程值来判断该块的地形特征，分为四种情况：海洋、陆地、海陆接边、未知，依次用整数标识记为0、1、2、3；对于海陆接边块，对海陆接边块内每个顶点继续进行海陆二值化，如果高程值小于等于0，则认为该顶点位于海洋，标识为0；否则，认为其位于陆地，标识为1；S22、压缩：对于海陆接边块中的每个顶点，每个顶点的标识占用1个bit，每8个相邻的顶点占用1个字节，然后再对二值化后的海陆接边块数据进行Zip压缩；S23、合并：将经过二值化和压缩后的861个demx文件合并为一个全球DEM数据文件；S3、利用合并后的全球DEM数据文件进行实时海陆判断：S31、将步骤S2中经过二值化、压缩、合并的全球DEM数据文件存储到计算机的硬盘，在计算机内存中构建一级索引和二级索引，其中一级索引用于记录从文件中读取的所有块的类型标识，每块占据3个字节，即24个bit位；二级索引用于记录所有海陆接边块的原始索引信息，即块的字节数和块在文件中的偏移地址；根据块的地形类型，一级索引情况如下：如果块为海洋，一级索引的24个bit位中最后3位记为0...
【专利技术属性】
技术研发人员：杜莹，武玉国，蒲欢欢，马亚明，孙汀，韦原原，
申请(专利权)人：郑州师范学院，
类型：发明
国别省市：河南,41
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。