java 利用 gdal 对shpfile 进行裁剪(一)

简要说明

在java,使用gdal ,使用裁剪范围的geojson 对 目标shpfile 进行裁剪,裁剪结果为一份新的shpfile;前置条件是geojson的坐标系和 shpfile 的坐标系一致,不然裁剪结果可能为空。

maven依赖

<!--需要安装完gdal后,本地install gdal包才能使用 --><!--gdal安装可参考 https://blog.csdn.net/qq_41613913/article/details/135743562 --><dependency><groupId>org.gdal</groupId><artifactId>gdal</artifactId><version>3.7.3</version></dependency>

样例代码

public static void main(String[] args) {String geojsonStr = "your geojson";if(StringUtils.isEmpty(geojsonStr )) {System.out.println("geojson 字符串为空");return;}JSONObject jsonObject = JSONUtil.parseObj(overlayRange);String features = jsonObject.getStr("features");//将geoJson转成Geometry集合对象Geometry geometry = ogr.CreateGeometryFromJson(JsonUtil.toJsonString(JsonUtil.jsonStringToList(features).get(0).get("geometry")));//临时的裁剪结果目录String clippingParentDirFileStr = "D:/tmp/clipping" + File.separator ;DataSource originalDataSource = null;try {File clippingParentDir = new File(clippingParentDirFileStr);if (!clippingParentDir.exists()) {clippingParentDir.mkdirs();}//打开数据org.gdal.ogr.Driver shpDriver = ogr.GetDriverByName("ESRI Shapefile");//被裁剪shp文件String originalFilePath = "D:/data/shp/xxxx.shp";originalDataSource = shpDriver.Open(originalFilePath, 0);//通过hutools工具获取文件名String originalName = FileUtil.getName(originalFilePath);String layerName = originalName.replace(".shp", "");//拿到图层Layer originalLayer = originalDataSource.GetLayer(layerName);//结果图层文件String outputFilePath = clippingParentDirFileStr + "output.shp";DataSource outDataSource = shpDriver.CreateDataSource(outputFileStr);//设置坐标系String originalSrid = originalLayer.GetSpatialRef().GetAttrValue("AUTHORITY", 1);if("0".equals(originalSrid)) {System.out.println("源文件不存在坐标系");return;}SpatialReference sr = new SpatialReference();sr.ImportFromEPSG(Integer.valueOf(originalSrid));Layer outlayer = outDataSource.CreateLayer("output", sr);//设置字段FeatureDefn originalfeatureDefn = originalLayer.GetLayerDefn();int fieldCount = originalfeatureDefn.GetFieldCount();for (int i = 0; i < fieldCount; i++) {FieldDefn defn = originalfeatureDefn.GetFieldDefn(i);outlayer.CreateField(defn);}//根据相交结果进行裁剪long featureCount = originalLayer.GetFeatureCount();Feature feature2 = null;for (int i = 0; i < featureCount ;i++) {feature2 = originalLayer.GetFeature(i);Geometry geometry2 = feature2.GetGeometryRef();if (geometry.Intersect(geometry2)) {Geometry intersection = geometry.Intersection(geometry2);feature2.SetGeometry(intersection);outlayer.CreateFeature(feature2);}}outlayer.delete();outDataSource.delete();shpDriver.delete();originalDataSource.delete();} catch (Exception e) {e.printStackTrace();}}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/pingmian/39120.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

STL的六大组件

一.总体概念 STL&#xff08;Standard Template Library&#xff0c;标准模板库&#xff09;是C标准库的一部分&#xff0c;提供了丰富且高效的数据结构和算法。STL主要由六大组件组成&#xff0c;它们是&#xff1a; 容器&#xff08;Containers&#xff09;&#xff1a;STL提…

ai智能电话机器人的语音识别究竟是如何实现的,机器人部署,识别次数惠

什么是智能语音识别系统&#xff1f;语音识别实际就是将人类说话的内容和意思转化为计算机可读的输入&#xff0c;例如按键、二进制编码或者字符序列等。与说话人的识别不同&#xff0c;后者主要是识别和确认发出语音的人并非其中所包含的内容。语音识别的目的就是让机器人听懂…

RPA自动化实践:发票处理流程自动化

在日常的企业运营中&#xff0c;发票处理是一项繁琐且重复性高的工作。本文将通过一个具体的RPA项目——自动化发票处理流程&#xff0c;展示RPA技术如何提高效率和准确性。 引言 RPA自动化不仅仅是技术的应用&#xff0c;更是一种工作方式的革新。通过自动化发票处理流程&am…

1.(vue3.x+vite)实现卷帘效果

前端技术社区总目录(订阅之前请先查看该博客) 1:效果预览 2:代码编写 <template><div style="width

如何使用Maxscript访问C#类库?

本教程帮助你一步一步集成maxscript文档和C#类库&#xff0c;以便你能够在maxscript提供的相当有限的环境中访问dotnet框架的强大功能。这让你可以使用maxscript做一些功能更强大的事情&#xff0c;比如访问数据库、获取web部署的内容等等。 还是直接上教程实惠&#xff1a; …

昇思25天学习打卡营第6天|关于函数与神经网络梯度相关技术探讨

目录 Python 库及 MindSpore 相关模块和类的导入 函数与计算图 微分函数与梯度计算 Stop Gradient Auxiliary data 神经网络梯度计算 Python 库及 MindSpore 相关模块和类的导入 Python 中的 numpy 库被成功导入&#xff0c;并简称为 np。numpy 在科学计算领域应用广泛&#x…

SQLyog脚本无限试用重置脚本

文章目录 引言脚本(win)必要操作、说明 引言 SQLyog 需要po jie&#xff0c;但是网上的没看到很好使的&#xff0c;直接下的官方。能处理14天试用也是很ok的。 脚本(win) echo offREM SQLyog注册表key&#xff0c;可能跟你的不一样&#xff0c;如果不一样&#xff0c;请替换…

ai绘画一条作品变现1400+,怎么做一个赚钱的AI绘画账号?

大家都知道现在AI很火&#xff0c;变现的玩法也多种多样&#xff0c;但一说到AI&#xff0c;大家就下意识认为这东西离我太远&#xff0c;自己没有那么高学历&#xff0c;不会英文&#xff0c;不会用AI模型等。 其实Ai没有大家想象得那么难&#xff0c;尤其在AI绘画这块&#…

kvm 虚拟机增加内存、cpu

目录 kvm 虚拟机增加内存、cpu kvm 虚拟机增加内存、cpu #关闭虚拟机 virsh shutdown <虚拟机名称或ID>#如果虚拟机未能正常关闭&#xff0c;可以使用强制关闭命令&#xff1a; virsh destroy <虚拟机名称或ID>#编辑配置 virsh edit <虚拟机名称或ID><vcp…

Phoenix Omid Timestamp Oracle 组件实现原理

Omid Timestamp Oracle 组件实现原理 作用 生成全局单调递增的时间戳&#xff0c;支持获取操作和崩溃恢复。 功能 1.生成全局单调递增的时间戳(支持崩溃恢复)apinext返回下一个时间戳getLast返回最后一个分配的时间戳(当前时间戳)实现方式TimestampOracleImpl单调递增的时间…

DNS访问百度

DNS&#xff0c;英文全称是 domain name system&#xff0c;域名解析系统&#xff0c;它的作用也很明确&#xff0c;就是域名和 IP 相互映射。 假设你要查询 baidu.com 的 IP 地址: 首先会查找浏览器的缓存,看看是否能找到 baidu.com 对应的IP地址&#xff0c;找到就直接返回&…

【第七节】C/C++排序算法

目录 前言 一、冒泡排序 二、选择排序 三、插入排序 四、希尔排序 五、归并排序 六、快速排序 七、 堆排序 八、计数排序 九、桶排序 十、基数排序 前言 排序算法可以大致分为两大类&#xff1a;比较类排序和非比较类排序。以下是这两大类中一些常见的排序算法示例&…

ChatGPT-4o医学应用、论文撰写、数据分析与可视化、机器学习建模、病例自动化处理、病情分析与诊断支持

2022年11月30日&#xff0c;可能将成为一个改变人类历史的日子——美国人工智能开发机构OpenAI推出了聊天机器人ChatGPT-3.5&#xff0c;将人工智能的发展推向了一个新的高度。2023年11月7日&#xff0c;OpenAI首届开发者大会被称为“科技界的春晚”&#xff0c;吸引了全球广大…

WPF布局控件

目录 Grid StackPanel WrapPanel DockPanel UniformGrid Canvas&InkCanvas Canvas InkCanvas Border Grid 属性 ShowGridLines&#xff1a;显示边线 ColumnDefinitions 列集合 表示有几列下面就写几个ColumnDefinition Width 宽&#xff1a;如果写具体数字则表…

Spring-cloud-openfeign解码器Decoder接口(后置拦截器)

使用feign调用第三方的http服务&#xff0c;对方返回response&#xff0c;之后这个Decoder接口会将对方的返回值&#xff0c;序列化成我们的返回值&#xff0c;例如下面的代码中&#xff0c;为什么我们能拿到User类型&#xff0c;而不是一个String类型&#xff0c;这就是Decode…

科普文:一文搞懂jvm实战(一)Runtime实时监控jvm

概叙 Java Runtime 类是 Java 标准库中的关键类之一。它提供了对当前Java虚拟机&#xff08;JVM&#xff09;实例的访问和控制&#xff0c;允许程序动态地修改和管理运行时环境。 Java Runtime 是Java虚拟机&#xff08;JVM&#xff09;的一个实例&#xff0c;代表了正在执行Ja…

Kafka 进阶指南

Kafka 进阶指南 引言 在掌握了 Kafka 的基本概念和操作后&#xff0c;我们可以进一步探索 Kafka 的高级特性和使用技巧&#xff0c;以提高其性能、可扩展性和可靠性。本指南将介绍 Kafka 的进阶主题&#xff0c;包括性能调优、扩展策略、数据复制、日志压缩、流处理和安全性。…

nrm工具——npm源管理器

nrm 是一个 npm 源管理器&#xff0c;允许你快速地在 npm源间切换。 一开始觉得这个工具没用&#xff0c;一般使用npm config set registry XXXXX命令就可以切换源地址 但是突然发现&#xff01;&#xff01;&#xff01;&#xff01; 淘宝镜像地址改了&#xff01;&#xf…

Spring Boot 实现 AOP 动态热插拔功能并附DEMO源码

&#x1f604; 19年之后由于某些原因断更了三年&#xff0c;23年重新扬帆起航&#xff0c;推出更多优质博文&#xff0c;希望大家多多支持&#xff5e; &#x1f337; 古之立大事者&#xff0c;不惟有超世之才&#xff0c;亦必有坚忍不拔之志 &#x1f390; 个人CSND主页——Mi…

python中的包和模块

目录 一、包与模块 二、第三方包的安装 2.1 pip install 2.2使用 curl 管道 2.3其他安装方法 三、导入单元的构成 3.1pip的使用 四、模块的缓存 一、包与模块 Python 中除了函数库以外&#xff0c;还有非常多且优秀的第三方库、包、模块。 模块Module&#xff1a;以…