CN108052629A - 一种基于高精度dem数据的快速海陆判断方法 - Google Patents
一种基于高精度dem数据的快速海陆判断方法 Download PDFInfo
- Publication number
- CN108052629A CN108052629A CN201711377020.0A CN201711377020A CN108052629A CN 108052629 A CN108052629 A CN 108052629A CN 201711377020 A CN201711377020 A CN 201711377020A CN 108052629 A CN108052629 A CN 108052629A
- Authority
- CN
- China
- Prior art keywords
- block
- extra large
- index
- land
- file
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 12
- 241000132092 Aster Species 0.000 claims abstract description 9
- 235000013399 edible fruits Nutrition 0.000 claims description 15
- 230000006835 compression Effects 0.000 claims description 14
- 238000007906 compression Methods 0.000 claims description 14
- 230000006837 decompression Effects 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000006399 behavior Effects 0.000 description 3
- 238000000205 computational method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000013441 quality evaluation Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003653 coastal water Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T1/00—General purpose image data processing
- G06T1/60—Memory management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- Remote Sensing (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Instructional Devices (AREA)
Abstract
本发明提供了一种基于全球高精度DEM数据的快速海陆判断方法，属于地理信息系统领域，先将全球ASTER GDEM数据文件按照一定的空间索引方案进行插值重构，在此基础上对重构后的DEM数据文件进行二值化和压缩并合并成一个文件，最后利用合并后的数据文件进行实时海陆判断。该方法可对海量地理坐标进行快速的海陆判断，可广泛应用于航船舶交通信息坐标异常判断、船舶交通行为可视化与分析等领域。
Description
技术领域
本发明涉及地理信息系统领域，特别是一种基于高精度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、对全球ASTER GDEM数据文件进行插值重构获得重构后的全球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、对全球ASTER GDEM数据文件进行插值重构
(1)全球空间索引方案
ASTER GDEM数据覆盖范围为北纬83°到南纬83°之间的所有陆地区域，比以往任何地形图都要广得多，达到了地球陆地表面的99％。其空间分辨率为1″，垂直精度达20米，水平精度达30米，是迄今为止可为用户提供的最完整的、精度最高的全球数字高程模型数据。ASTER GDEM在全球范围内按照1°×1°进行分片，每个分片包含两个压缩文件，一个数字高程模型(DEM)文件和一个质量评估(QA)文件。所有DEM文件均以GeoTIFF格式存储，像素个数为3601×3601，海平面高程值为0，无效像素值为-9999。
按照ASTER GDEM数据的分片方式，其全球DEM数据文件共计64800个，每个文件大小为24.7M字节，数据总量达1.5T字节。数据总量大但同时每个文件又很小，非常不便于进行频繁的IO操作，而且不利于快速的海陆判断，这就需要对这些DEM数据文件按照一定的空间索引方案进行插值重构。
本发明设计了如下的空间索引方案：将全球平均分成4096×2048块，每100×100块合并为一个新的DEM文件，块内顶点个数为256×256，每个顶点的高程值用一个float型变量存储。按照这样的方案，数据水平分辨率为(360°÷4096)÷(256-1)≈1.2″，与原有水平分辨率1″基本相当；每个文件大小为2.44G字节，覆盖地理范围约为8.79°，可以大大提高频繁IO操作的效率。
(2)插值重构
首先，在每个块内对256×256个顶点的高程值进行双线性插值。需要注意的是，原始ASTER GDEM数据中可能会因为云遮挡等原因出现无效像素值-9999，用该值参与双线性插值，是没有任何意义的。我们的做法是在用四个原始顶点进行双线性插值过程中，一旦有两个或两个以上顶点为无效像素，则认为待插值点为无效像素点，直接将其高程值记录为-9999，以做标识。
接下来，是要将插值后的每100×100块存储为一个新的DEM数据文件，后缀名为demx，全球共计861个这样的文件。每个demx文件分为三部分：文件头、块索引、块数据(如图1所示)。其中，文件头存储数据的版本号和文件内的分块数目；每个块索引存储该块数据的字节总数和对应的块数据在文件的偏移地址；每个块数据存储该块的元信息和块内所有顶点(256×256个)的高程值。
2、对重构后的全球DEM数据文件进行二值化、压缩、合并
(1)海洋和陆地块的二值化
重构后的全球DEM数据分块进行存储，接下来需要根据块数据中顶点的高程值来判断该块的地形特征，分为四种情况：海洋、陆地、海陆接边、未知。
如果所有顶点的高程值全部小于等于0，则该块为海洋，块标识记为0；
如果所有顶点的高程值全部大于0，则该块为陆地，块标识记为1；
否则，该块为海陆接边，块标识记为2；
如果块中没有任何顶点，即字节数为0，则该块为未知(这些块位于最后一个demx文件的最后若干列，即该文件不满100×100块)，块标识记为3。
如果用n表示某顶点在块内的序号，用k表示其在8个相邻顶点中的序号，那么该点的海陆二值化结果Vbit为：
其所在字节Vbyte的计算方法为：Vbyte＝Vbyte|Vbit。
如果块为海洋、陆地或未知，情况简单，只需要在原有块索引区将块标识记为一个整数即可，不再需要字节数约为256K字节的块数据，再加上全球海洋面积占地球表面积的71.8％，这种二值化做法就可以大大减少全球DEM数据的字节数。
如果块为海陆接边，情况就要复杂一些，仅仅记录一个块标识还不够，还需要对块内每个顶点继续进行海陆二值化，以便明确每个顶点到底位于海洋还是陆地。
(2)海陆接边块的压缩
对于海陆接边块中的每个顶点，如果高程值小于等于0，则认为该顶点位于海洋，标识为0；否则，认为其位于陆地，标识为1。为减少数据量，每个顶点的标识占用1个bit，每8个相邻的顶点占用1个字节。这样，海陆接边块的数据量就从256K字节减少到了8K字节。为进一步减少数据量，还需要对二值化后的海陆接边块数据(8K字节)进行Zip压缩，可使用开源库zlib完成该操作，将compress2()函数中的参数level设置为9，表示使用最高的压缩率。
实验结果表明，经过二值化和Zip压缩之后，每个demx文件的大小从2.44G锐减到1M字节以下，这就大大减少了数据量，提高了IO操作的效率，便于接下来的快速海陆判断。
(3)合并为一个数据文件
经过上述步骤处理之后的每个demx文件，不超过1M字节，为进一步提高IO操作效率，可将全球861个demx文件合并为一个全球DEM数据文件，全球DEM数据文件格式如图2所示，将全球DEM数据文件栅格化后的结果如图3所示。图2中的全球DEM数据文件格式与图1中的文件格式很相似，主要区别在于块数据区不再需要记录海洋、陆地、未知这样三类块的数据了，只需要记录海陆接边块的Zip压缩数据即可。注意：在头文件中需要记录四类块(海洋、陆地、海陆接边、未知)的总数，以简化后续的统计工作。
3、利用合并后的全球DEM数据文件进行实时海陆判断
(1)在计算机内存中构建一级索引和二级索引
ASTER GDEM数据经过重构、二值化、压缩、合并之后，就可以存储到计算机硬盘，准备用于实时的海陆判断了。但要想快速判断出一个点位于海洋还是陆地，不可能频繁地到硬盘中去访问该文件，还需要将文件中的相关信息读取到计算机内存，因为计算机内存访问速度毕竟比IO访问速度要快得多。
本发明在计算机内存中构建了两级索引，分别如下：
一级索引：用于记录从文件中读取的所有块的类型标识，每块占据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位用于记录该块在二级索引中的序号，即该块在所有海陆接边中的自然顺序。
如果用k表示某个海陆接边块在所有海陆接边块中的序号，那么该块的一级索引值Index的3个字节依次为：
V＝(k＜＜1)|0x01
Index[0]＝(BYTE)(V＞＞16)
Index[1]＝(BYTE)((V＜＜8)＞＞16)
Index[2]＝(BYTE)((V＜＜16)＞＞16)。
(2)基于两级索引的快速海陆判断
在内存中构建完两级索引之后，就可以进行快速海陆判断了。如果是海洋块或陆地块，可直接认定待判断的点位于海洋或陆地，算法结束；如果是海陆接边块，还需要解析一级索引的3个字节，计算其在二级索引中的序号。如果分别用Index[0]、Index[1]、Index[2]表示解析后的一级索引的3个字节，那么该块在二级索引中的序号IndexSecond为：
IndexSecond＝[((int)Index[0]＜＜16)|((int)Index[1]＜＜8)|(Index[2])]＞＞1
根据该二级索引序号，即可到缓存或文件中获取该块的海陆二值化数据，进一步确定该点位于海洋还是陆地。如图4所示，为整个实时判断算法的流程图，主要包括以下几个步骤：
①根据输入点的经纬度坐标，计算其所在的地形块，并从内存中获取该块的一级索引信息；
②解析一级索引的3个字节，如果从高位到低位依次为0、0、0，则说明该块为海洋，从而知该点位于海洋，本次判断结束；如果从高位到低位依次为0、0、2，则说明该块为陆地，从而知该点位于陆地，本次判断结束；如果从高位到低位依次为0、0、4，则说明该块为未知，从而知该点地形未知，本次判断结束；否则，该块为海陆接边，该点位于海洋或者陆地，继续进行步骤③；
③根据海陆接边块的全球编号，到缓存中查找是否有该块解压后的海陆二值化数据，如果有，直接转向步骤④；否则，根据一级索引中前面23个bit位的值，反算其在二级索引中的序号，并根据该序号提取其在内存中的二级索引信息，根据二级索引信息中记录的文件偏移地址，到文件对应位置读取该块数据并解压为海陆二值化数据，然后存储到缓存中；
④根据该点的经纬度坐标，计算其在该海陆接边块中的行列号，根据行列号计算该点位于二值化数据的哪个字节，以及位于该字节的哪个bit位，如果该bit位为0，表示该点位于海洋，否则，表示该点位于陆地，本次判断结束。
本发明的技术方案不限于上述具体实施例的限制，凡是根据本发明的技术方案做出的技术变形，均落入本发明的保护范围之内。
Claims (2)
1.一种基于高精度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，表示该点位于海洋，否则，表示该点位于陆地，本次判断结束。
2.根据权利要求1所述的基于高精度DEM数据的快速海陆判断方法，其特征在于：所述步骤S12中用每个块内的四个原始顶点进行双线性插值过程中，一旦有两个或两个以上原始顶点为无效像素，则认为待插值点为无效像素点，直接将其高程值记录为-9999，以做标识。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711377020.0A CN108052629B (zh) | 2017-12-19 | 2017-12-19 | 一种基于高精度dem数据的快速海陆判断方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711377020.0A CN108052629B (zh) | 2017-12-19 | 2017-12-19 | 一种基于高精度dem数据的快速海陆判断方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108052629A true CN108052629A (zh) | 2018-05-18 |
| CN108052629B CN108052629B (zh) | 2021-07-06 |
Family
ID=62130099
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711377020.0A Active CN108052629B (zh) | 2017-12-19 | 2017-12-19 | 一种基于高精度dem数据的快速海陆判断方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108052629B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109886134A (zh) * | 2019-01-25 | 2019-06-14 | 北京市遥感信息研究所 | 高分辨率全球数字高程模型图像的海陆信息索引编码方法 |
| CN111913179A (zh) * | 2020-08-20 | 2020-11-10 | 上海无线电设备研究所 | 一种提高星载雷达高度计近海岸观测能力的方法 |
| CN117312471A (zh) * | 2023-09-26 | 2023-12-29 | 中国人民解放军91977 部队 | 海量位置点的海陆属性判断方法及装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103678658A (zh) * | 2013-12-24 | 2014-03-26 | 北京中科大洋科技发展股份有限公司 | 一种地形高程数据存储和读取的装置 |
| CN104063870A (zh) * | 2014-07-04 | 2014-09-24 | 中国科学院大学 | 一种基于扫描线检测的海陆模板自动分割方法及其应用 |
| GB2529271A (en) * | 2014-08-15 | 2016-02-17 | Satavia Ltd | Meteorological hazard identification apparatus, moving craft and methods |
| CN105512622A (zh) * | 2015-12-01 | 2016-04-20 | 北京航空航天大学 | 一种基于图分割和监督学习的可见光遥感图像海陆分割方法 |
| CN105513041A (zh) * | 2015-10-28 | 2016-04-20 | 深圳大学 | 一种大幅面遥感图像海陆分割的方法及系统 |
| CN107340503A (zh) * | 2017-07-02 | 2017-11-10 | 中国航空工业集团公司雷华电子技术研究所 | 一种基于数字高程地图抑制虚假海面目标的方法 |
-
2017
- 2017-12-19 CN CN201711377020.0A patent/CN108052629B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103678658A (zh) * | 2013-12-24 | 2014-03-26 | 北京中科大洋科技发展股份有限公司 | 一种地形高程数据存储和读取的装置 |
| CN104063870A (zh) * | 2014-07-04 | 2014-09-24 | 中国科学院大学 | 一种基于扫描线检测的海陆模板自动分割方法及其应用 |
| GB2529271A (en) * | 2014-08-15 | 2016-02-17 | Satavia Ltd | Meteorological hazard identification apparatus, moving craft and methods |
| CN105513041A (zh) * | 2015-10-28 | 2016-04-20 | 深圳大学 | 一种大幅面遥感图像海陆分割的方法及系统 |
| CN105512622A (zh) * | 2015-12-01 | 2016-04-20 | 北京航空航天大学 | 一种基于图分割和监督学习的可见光遥感图像海陆分割方法 |
| CN107340503A (zh) * | 2017-07-02 | 2017-11-10 | 中国航空工业集团公司雷华电子技术研究所 | 一种基于数字高程地图抑制虚假海面目标的方法 |
Non-Patent Citations (5)
| Title |
|---|
| N. OTHMAN ET AL: "Using ASTER GDEM and SRTM digital elevation models to generate contour lines over rugged terrain of Makkah", 《IEEE》 * |
| WENJIAN NI ET AL: "Investigation of forest height retrieval using SRTM-DEM and ASTER-GDEM", 《IEEE》 * |
| 杜莹 等: "全球多分辨率虚拟地形环境的金字塔模型研究", 《系统方正学报》 * |
| 杜莹 等: "设计模式在三维地形环境仿真平台中的应用", 《测绘工程》 * |
| 武玉国 等: "大规模地形TIN模型的LOD算法设计与实现", 《系统仿真学报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109886134A (zh) * | 2019-01-25 | 2019-06-14 | 北京市遥感信息研究所 | 高分辨率全球数字高程模型图像的海陆信息索引编码方法 |
| CN111913179A (zh) * | 2020-08-20 | 2020-11-10 | 上海无线电设备研究所 | 一种提高星载雷达高度计近海岸观测能力的方法 |
| CN111913179B (zh) * | 2020-08-20 | 2023-04-11 | 上海无线电设备研究所 | 一种提高星载雷达高度计近海岸观测能力的方法 |
| CN117312471A (zh) * | 2023-09-26 | 2023-12-29 | 中国人民解放军91977 部队 | 海量位置点的海陆属性判断方法及装置 |
| CN117312471B (zh) * | 2023-09-26 | 2024-05-28 | 中国人民解放军91977部队 | 海量位置点的海陆属性判断方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108052629B (zh) | 2021-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhao et al. | A method for simplifying ship trajectory based on improved Douglas–Peucker algorithm | |
| CN113553682B (zh) | 一种数据驱动的多层次船舶航路网络构建方法 | |
| CN108052629A (zh) | 一种基于高精度dem数据的快速海陆判断方法 | |
| CN112562415A (zh) | 船舶航线推荐方法 | |
| CN109949197B (zh) | 一种面向海河联运的船岸一体化系统 | |
| CN112347218B (zh) | 无人船艇的环境地图生成方法及无人船艇感知系统 | |
| CN115577324B (zh) | 雷达和船舶自动识别系统数据融合方法及系统、存储介质 | |
| CN111428916B (zh) | 一种海上救援船舶巡航路径规划方法 | |
| CN114936332A (zh) | 一种风电场水域船舶交通流分析的方法、装置和存储介质 | |
| CN117037547B (zh) | 一种航道关键节点的拥堵程度评价方法、系统及存储介质 | |
| KR101214996B1 (ko) | 해양지리 정보 구축 및 서비스 방법 | |
| CN104392113B (zh) | 一种近海海面冷空气大风风速的估算方法 | |
| CN109886134A (zh) | 高分辨率全球数字高程模型图像的海陆信息索引编码方法 | |
| CN117082474A (zh) | 一种科考船实时获取海洋环境预报数据的系统 | |
| Weixi et al. | A Case Research on Vulnerability of logistics system in the Tianjin port | |
| Ribeiro et al. | Integration of high-resolution metocean forecast and observing systems at Port of Santos | |
| CN114529834A (zh) | 小目标对象的检测方法和检测装置 | |
| CN111028482B (zh) | 远洋航行的提醒方法、装置及系统 | |
| CN103678658A (zh) | 一种地形高程数据存储和读取的装置 | |
| CN109241207A (zh) | 一种在地图上显示数据的方法及装置 | |
| CN105279069A (zh) | 一种基于索引矩阵的目标快速告警判断方法 | |
| CN108597013B (zh) | 一种特定区域气象水文资料抽点填图方法 | |
| Lin et al. | Shoreline data based sea-land segmentation method for on-orbit ship detection from panchromatic images | |
| Nugraha et al. | Spatiotemporal Analysis of Urban Land Cover: Case Study-Pekanbaru City, Indonesia | |
| CN115376013B (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 |