【GDAL/OGR】利用GDAL/OGR读取shp文件并转换为json文件(Java版)

前言:

对于GIS开发者来说,GDAL/OGR是最熟悉不过的开源GIS库了,GDAL(Geospatial Data Abstraction
Library)是一个在X/MIT许可协议下的开源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式。它还有一系列命令行工具来进行数据转换和处理。OGR是GDAL项目的一个分支,功能与GDAL类似,只不过它提供对矢量数据的支持。

GDAL官网

编译的dll文件下载以及各个语言的相应开发包,这个地址里面有各个不同版本的编译文件,里面会详细介绍里面的各个文件目录的文件都是什么文件。例如:这个版本的编译好的文件目录详情

在开发前,需要安装GDAL(其实就是配置相应的环境变量)

GDAL的安装:
在上面第三个网址中,里面有相应的安装介绍,读者可以自己相应的安装自己需要的组件和开发环境。
按照官网的说法如下:

In order to run the utility programs the compiled dll-s should be available to load during the execution. Therefore the PATH environment parameter should contain the full path to /bin directory of the package or the dlls in this directory should be copied into the directory from which the application (executable) is running

这里有两种安装方法:
(1)把自己需要用到的dll文件路径配置到环境变量的Path变量中(主要是bin目录下的dll文件,这部分文件是common dll,以及bin/gdal/java目录下的dll文件),至于jar包就需要引入到项目中。
(2)将需要的dll文件放到项目的根目录下

这里需要使用java调用dll文件,至于其原理,可以参考:

Java 本地调用(java调用window DLL里的函数)

主要代码:

/**
     * 将shp文件转换为json数据文件输出,这个json后面缺少反括号,要自己加上去,是gdal的bug还是?
     * @param shpfilePath
     * @param outGeojsonPath
     */
    public static void parseShpfile2GeojsonFile(String shpfilePath,String outGeojsonPath) {
        // 注册所有的驱动
        ogr.RegisterAll();
        // 为了支持中文路径,请添加下面这句代码
        gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "YES");
        // 为了使属性表字段支持中文,请添加下面这句
        gdal.SetConfigOption("SHAPE_ENCODING", "");

        // 打开文件
        DataSource ds = ogr.Open(shpfilePath, 0);
        if (ds == null) {
            System.out.println("打开文件失败!");
            return;
        }
        System.out.println("打开文件成功!");
        Driver dv = ogr.GetDriverByName("GeoJSON");
        if (dv == null) {
            System.out.println("打开驱动失败!");
            return;
        }
        System.out.println("打开驱动成功!");
        dv.CopyDataSource(ds, outGeojsonPath);
        System.out.println("转换成功!");
    }

注:

就笔者的实验结果,转换后的json文件的json格式有一点问题,就是json格式的最后面缺少反大括号,不知道这是笔者使用的这个gdal版本的bug还是其他情况。读者根据自己的实际情况做处理,同时也希望有了解的同学,指导交流一下。


另外,笔者看到另外一种不使用代码方式转换文件的方法,就是使用ogr中的命令行进行格式转换。
具体网址参考:

ogr2ogr命令行使用

评论中提到的用法

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
以下是使用JavaGDAL/OGR打开矢量文件并输出每个面外界矩形范围内的POINT数据的代码示例: ```java import org.gdal.ogr.DataSource; import org.gdal.ogr.Driver; import org.gdal.ogr.Layer; import org.gdal.ogr.Feature; import org.gdal.ogr.Geometry; import org.gdal.ogr.ogr; public class GDALOGRExample { public static void main(String[] args) { // 初始化GDAL/OGR ogr.RegisterAll(); // 打开矢量文件 String filename = "path/to/vector/file.shp"; DataSource dataSource = ogr.Open(filename); // 获取第一个图层 Layer layer = dataSource.GetLayer(0); // 遍历每个面 Feature feature; while ((feature = layer.GetNextFeature()) != null) { Geometry geometry = feature.GetGeometryRef(); // 获取外界矩形 Geometry envelope = geometry.GetEnvelope(); double minX = envelope.GetX(0); double maxX = envelope.GetX(1); double minY = envelope.GetY(0); double maxY = envelope.GetY(1); // 输出每个面外界矩形范围内的POINT数据 Layer pointLayer = dataSource.ExecuteSQL("SELECT * FROM " + layer.GetName() + " WHERE " + "X > " + minX + " AND X < " + maxX + " AND Y > " + minY + " AND Y < " + maxY, null, "SQLITE"); Feature pointFeature; while ((pointFeature = pointLayer.GetNextFeature()) != null) { Geometry pointGeometry = pointFeature.GetGeometryRef(); System.out.println(pointGeometry.GetX() + "," + pointGeometry.GetY()); } pointLayer.delete(); feature.delete(); geometry.delete(); } layer.delete(); dataSource.delete(); } } ``` 代码解释: 1. 首先需要初始化GDAL/OGR,可以调用`ogr.RegisterAll()`方法进行初始化; 2. 打开矢量文件,可以使用`ogr.Open()`方法; 3. 获取第一个图层,可以使用`dataSource.GetLayer()`方法; 4. 遍历每个面,可以使用`layer.GetNextFeature()`方法获取每个面的Feature对象; 5. 获取每个面的外界矩形范围,可以使用`geometry.GetEnvelope()`方法获取外界矩形,然后获取最小X、最大X、最小Y、最大Y; 6. 输出每个面外界矩形范围内的POINT数据,可以使用SQL语句过滤出符合条件的数据,然后使用`dataSource.ExecuteSQL()`方法执行SQL语句并返回一个新的Layer对象,遍历该Layer对象的每个Feature对象,获取每个Feature对象的Geometry对象,并输出其坐标; 7. 最后需要释放资源,可以调用各个对象的`delete()`方法进行释放。 注意事项: 1. 需要引入GDAL/OGRJava库,可以使用Maven或手动下载并添加到classpath中; 2. 矢量文件的路径需要修改为实际的路径; 3. 输出结果可以根据需求进行修改。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值