我要评论前言
在当今数字化时代，地理空间数据的应用愈发广泛，而海图作为重要的地理信息资源，在航海、海洋研究等领域扮演着关键角色。s57 海图文件作为一种标准化的海图数据格式，承载着丰富的海洋地理信息，如何高效、准确地读取和解析这些数据，成为了众多开发者和研究人员关注的焦点。近期有朋友与我交流，想了解s57海图文件的读取与发布。海图文件我其实没有接触过，所以对其也是一无所知。因此今天，就让我们一同走进这场基于 java 和 gdal 的 s57 海图文件解析实战之旅，探索其中的技术奥秘与实践要点。
在开始实战之前，我们需要对相关技术背景有所了解。java 作为一种成熟且广泛应用的编程语言，以其强大的跨平台特性、丰富的类库和良好的性能，为处理各种复杂任务提供了坚实基础。而 gdal（geospatial data abstraction library）则是一个开源的地理空间数据抽象库，它支持多种地理空间数据格式的读写操作，包括 s57 海图文件。将 java 与 gdal 结合，能够充分发挥两者的优势，实现对 s57 海图文件的高效解析。通过本次实战，我们将逐步深入，从搭建开发环境、配置相关依赖，到编写代码实现对 s57 海图文件的读取、解析以及数据提取等操作，全方位展示整个过程，让读者能够清晰地掌握每个环节的关键点。
这场实战不仅仅是对技术的实践应用，更是一次对地理空间数据处理领域的深入探索。在解析 s57 海图文件的过程中，我们会遇到各种数据结构、编码方式以及地理空间概念等问题，通过解决这些问题，我们能够更好地理解地理空间数据的内在逻辑和处理方法。此外，本次实战也将分享一些常见问题的解决方案和优化技巧，帮助读者在实际开发中避免遇到类似的困境。无论你是地理信息系统领域的开发者，还是对地理空间数据处理感兴趣的程序员，相信通过本文的介绍，你都将收获满满，为今后在相关领域的项目开发和个人学习奠定坚实的基础。
一、海图文件简介
为了让大家对海图文件有一个简单的了解，本节将对海图文件进行一个简单的介绍。由于之前也没有接触过海图文件，因此也是自己的学习过程。信息记录比较简单，在以后的学习中不断补充完善。首先给大家介绍依稀海图文件是什么？然后以s57为例具体介绍海图文件的文件类型。
1、什么是海图文件
海图文件是一种用于表示海洋地理信息的数字文件格式，它包含了海洋地形、水深、航行障碍物、海岸线、航标等重要地理信息。这些信息对于航海安全、海洋资源开发、海洋科学研究以及军事应用等领域至关重要。海图文件通常以矢量或栅格数据的形式存储，能够为船舶导航系统、地理信息系统（gis）以及相关海洋应用提供基础数据支持。
2、s57海图文件
s57 海图文件是一种国际标准化的矢量海图数据格式，由国际水道测量组织（iho）制定并广泛应用于全球海洋领域。s57 文件格式具有以下特点：
- 标准化：s57 格式遵循 iho 的 s-57 标准，确保了不同国家和地区之间的海图数据可以无缝交换和互操作。
- 矢量数据：s57 文件以矢量形式存储地理信息，包括点、线、面等几何对象，能够提供高精度的地理数据表示。
- 丰富的属性信息：除了几何数据外，s57 文件还包含了大量的属性信息，如水深、地形类型、航行警告等，这些属性信息为用户提供了更全面的海洋地理信息。
- 分层存储：s57 文件支持分层存储，用户可以根据需要选择显示或处理特定的图层，提高了数据处理的灵活性和效率。
- 兼容性：s57 文件格式被广泛应用于各种航海设备和地理信息系统软件中，具有良好的兼容性。
s57 海图文件的这些特点使其成为海洋地理信息领域的重要数据格式之一。通过解析 s57 文件，开发者可以提取出丰富的地理信息，用于开发各种海洋应用，如船舶导航系统、海洋环境监测、海洋资源管理等。
二、qgis读取s57实现
在了解了s57海图文件后，接下来以qgis为例，来看看海图文件展示出来是什么样子。从图层信息和属性信息查看两个方面进行介绍。
1、图层展示
为演示海图文件的加载，从c站获取网友提供的海图文件。下载后得到的s57海图信息如下：
然后在qgis中打开这些文件，与shp文件不一样。海图文件通常包含51个图层，以c1313100.000这个文件为例。使用qgis打开后可以看到有以下提示：
点击“ok”按钮将图层添加到图层组中，如下图所示：
在qgis的左边可以看到所有的图层信息，我们可以打开任意的一个图层文件查看文件信息。如下图所示：
2、属性查看
图层信息比较多，这里我们不一一的对其进行查看。跟shp文件一样，我们来看看其属性信息都包含哪些呢？在qgis中我们先来看看图层有哪些属性。在信息中可以看到以下信息：
可以看到这个图层大约有30个字段。再来看看属性表中这些信息度包含哪些内容。
这些数据跟我们常用的数据库类似，都是一条一条的记录。
三、gdal读取s57实现
除了可以使用qgis来打开读取s57文件，我们还可以使用gdal来进行海图文件的读取。目前很多对海图文件的读取都是c或者c++为主，还有python的。但是介绍java的读取的不多，因此本文也算是一个简单的补充。对广大的java gis 开发这有一个参考。
1、maven引用
这里使用一种比较简单的方式引入gdal的jar，需要注意的是，在java中要提前配置好gdal的环境才可以，关于gdal的java环境配置，可以翻阅之前的博客，有详细的介绍。这里贴出jar引用代码：
<project xmlns="http://maven.apache.org/pom/4.0.0"
xmlns:xsi="http://www.w3.org/2001/xmlschema-instance"
xsi:schemalocation="http://maven.apache.org/pom/4.0.0 http://maven.apache.org/xsd/maven-4.0.0.xsd">
<modelversion>4.0.0</modelversion>
<groupid>com.yelang</groupid>
<artifactid>gdal_demo1</artifactid>
<version>0.0.1-snapshot</version>
<name>gdal_demo1</name>
<description>gdal的第一次试验</description>
<properties>
<project.build.sourceencoding>utf-8</project.build.sourceencoding>
<maven.compiler.source>1.8</maven.compiler.source>
<maven.compiler.target>1.8</maven.compiler.target>
</properties>
<dependencies>
<dependency>
<groupid>org.gdal</groupid>
<artifactid>gdal</artifactid>
<version>3.4.3</version>
<scope>system</scope>
<systempath>${project.basedir}/lib/gdal.jar</systempath>
</dependency>
<dependency>
<groupid>net.sf.ucanaccess</groupid>
<artifactid>ucanaccess</artifactid>
<version>4.0.4</version>
</dependency>
<dependency>
<groupid>junit</groupid>
<artifactid>junit</artifactid>
<version>4.11</version>
<scope>test</scope>
</dependency>
<!-- https://mvnrepository.com/artifact/de.micromata.jak/javaapiforkml -->
<dependency>
<groupid>de.micromata.jak</groupid>
<artifactid>javaapiforkml</artifactid>
<version>2.2.1</version>
</dependency>
<dependency>
<groupid>cn.hutool</groupid>
<artifactid>hutool-all</artifactid>
<version>5.7.11</version>
</dependency>
<!-- https://mvnrepository.com/artifact/com.alibaba.fastjson2/fastjson2 -->
<dependency>
<groupid>com.alibaba.fastjson2</groupid>
<artifactid>fastjson2</artifactid>
<version>2.0.47</version>
</dependency>
<!-- https://mvnrepository.com/artifact/org.apache.commons/commons-compress -->
<dependency>
<groupid>org.apache.commons</groupid>
<artifactid>commons-compress</artifactid>
<version>1.21</version>
</dependency>
<dependency>
<groupid>com.vividsolutions</groupid>
<artifactid>jts</artifactid>
<version>1.13</version>
</dependency>
<dependency>
<groupid>com.google.code.gson</groupid>
<artifactid>gson</artifactid>
<version>2.6.2</version>
</dependency>
</dependencies>
</project>2、java读取
为了实现java调用gdal实现s57文件的读取，我们首先来看看应该调用什么驱动来进行读取。首先看看gdal官网有没有提供相关说明：
在官网上可以看到有对应的驱动，点开s57后，可以看到以下的描述：
在gdal中要实现s57文件的读取的驱动使用“s57”即可。下面给出具体的java读取代码：
package com.yelang.s57;
import java.util.arraylist;
import java.util.hashmap;
import java.util.list;
import java.util.map;
import org.gdal.gdal.gdal;
import org.gdal.ogr.datasource;
import org.gdal.ogr.feature;
import org.gdal.ogr.featuredefn;
import org.gdal.ogr.fielddefn;
import org.gdal.ogr.geometry;
import org.gdal.ogr.layer;
import org.gdal.ogr.ogr;
import org.gdal.osr.spatialreference;
import org.junit.test;
public class gdals57case {
@test
public void reads572() {
// 指定文件的名字和路径
string strvectorfile = "path/海图文件/testchart/testchart/c1313100.000";
// 注册所有的驱动
ogr.registerall();
// 为了支持中文路径，请添加下面这句代码
gdal.setconfigoption("gdal_filename_is_utf8", "yes");
// 为了使属性表字段支持中文，请添加下面这句
gdal.setconfigoption("shape_encoding", "cp936");
// 读取数据，这里以esri的shp文件为例
string strdrivername = "s57";
// 创建一个文件，根据strdrivername扩展名自动判断驱动类型
org.gdal.ogr.driver odriver = ogr.getdriverbyname(strdrivername);
if (odriver == null) {
system.out.println(strdrivername + " 驱动不可用！\n");
return;
}
datasource datasource = odriver.open(strvectorfile);
int totalcount = datasource.getlayercount();
for (int i = 0; i < datasource.getlayercount(); i++) {
layer layer = datasource.getlayer(i);
spatialreference srf = layer.getspatialref();
string srfstr = null == srf ? "空" : "空间参考坐标系：" + srf.getattrvalue("authority", 0)+ srf.getattrvalue("authority", 1);
string srid = null == srf ? "-" : srf.getattrvalue("authority", 1);
system.out.println("图层名称：<==>" + layer.getname() + "\t" + srfstr + "\t" + srid);
long featurecount = layer.getfeaturecount();
system.out.println("数据量大小==>" + featurecount);
// 获取所有的字段名
list<string> fieldnames = new arraylist<string>();
featuredefn featuredefn = layer.getlayerdefn();
for (int j = 0; j < featuredefn.getfieldcount(); j++) {
fielddefn fielddefn = featuredefn.getfielddefn(j);
string fields = fielddefn.getname();
fieldnames.add(fields);
}
//system.out.println("所有字段名：" + fieldnames);
// 遍历读取属性
feature feature = null;
while ((feature = layer.getnextfeature()) != null) {
list<string> fieldvalue = new arraylist<string>();
geometry geometry = feature.getgeometryref();
if(null != geometry) {
string geojson = feature.getgeometryref().exporttojson();
system.out.println("geojson ：" + geojson);
}
stringbuffer sbf = new stringbuffer(1024);
// 遍历字段名读取所有属性
for (string fieldname : fieldnames) {
string field = feature.getfieldasstring(fieldname);
fieldvalue.add(field);
//system.out.println(fieldname +"===>"+field);
sbf.append(fieldname + "=" + field + ";");
}
system.out.println(sbf.tostring());
}
system.out.println("---------------------------------------------------");
}
system.out.println("总图层数：" + totalcount);
datasource.delete();
gdal.gdaldestroydrivermanager();
}
}可以看到，对s57海图文件的解析与shp文件的解析基本一致。读取layer，然后就获取feature，获取属性解析输出即可。
3、读取成果输出
最后来看看使用java调用gdal对s57文件进行解析后的属性和空间信息，在控制台输出如下：
图层名称：<==>c_aggr 空 - 数据量大小==>2 rcid=2544;prim=255;grup=2;objl=400;rver=1;agen=71;fidn=-1375794613;fids=73;lnam=0047adff0a4b0049;lnam_refs=(34:0047583709140094,00470394b3e5002a,0047583709150094,004719a4545f0028,0047b20623220027,00475837090e0094,0047639bc6210029,00475837090f0094,0047b696f02a001d,0047583709170094,0047fd1733ef001a,0047583709160094,00477604421b001b,0047460172eb001f,0047f83ed80b0003,004700757e65004d,0047583709180094,0047cfea4be70020,0047583709190094,004792801e360025,00475837091b0094,0047722ab7040021,00475837091a0094,00471af81a8f0024,0047d5a9b55c001c,0047c35f10640048,0047566d45aa0091,0047ebd02f660047,0047befd2a3a0022,00475705dc260092,004700757e78004d,00470112a9790012,0047599563800001,00474bdbe3b40002);ffpt_rind=(34:3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3);nobjnm=;objnam=;inform=;ninfom=;ntxtds=;picrep=;scamax=;scamin=;txtdsc=;recdat=;recind=;sordat=;sorind=; rcid=2815;prim=255;grup=2;objl=400;rver=1;agen=71;fidn=950957764;fids=1;lnam=004738ae76c40001;lnam_refs=(7:00470dda54c0005a,0047c7096211005b,0047010366dd0011,00471e6602920033,00471e6602930033,00470112a8fe0012,00470112a90d0012);ffpt_rind=(7:3,3,3,3,3,3,3);nobjnm=;objnam=;inform=;ninfom=;ntxtds=;picrep=;scamax=;scamin=;txtdsc=;recdat=;recind=;sordat=;sorind=; ---------------------------------------------------
四、总结
以上就是本文的主要内容，本次实战，我们将逐步深入，从搭建开发环境、配置相关依赖，到编写代码实现对 s57 海图文件的读取、解析以及数据提取等操作，全方位展示整个过程，让读者能够清晰地掌握每个环节的关键点。这场实战不仅仅是对技术的实践应用，更是一次对地理空间数据处理领域的深入探索。海图文件繁多，可以做很多深入的研究，受限文章篇幅，这里不进行赘述。抛砖引玉，大家可以对功能进行扩展。行文仓促，定有不足之处，欢迎各位朋友在评论区批评指正，不胜感激。
以上就是java使用gdal解析s57海图文件全面详解的详细内容，更多关于java解析海图文件的资料请关注代码网其它相关文章！
发表评论