geotools外接矩形
2023-04-18 15:39:40 时间
geom.getEnvelope() 得到外接矩形,不一定是面积最小;可以对多边形的每一条边求外接矩形,然后比较得到最小外接矩形
参考:https://blog.csdn.net/qq_40985985/article/details/127034254
ArrayList<Polygon> polys = new ArrayList<Polygon>(); Polygon bbox = (Polygon) geom.getEnvelope(); polys.add(bbox); geoms.add(geom.getEnvelope()); List<Geometry> geoms = new ArrayList<Geometry>(); //Polygon poly = null; Geometry polyUnion = null; int i=0; while(itr.hasNext()){ SimpleFeature sf = itr.next(); //MultiPolygon geom = (MultiPolygon)sf.getDefaultGeometry(); Geometry geom = (Geometry) sf.getDefaultGeometry(); Polygon bbox = (Polygon) geom.getEnvelope(); //polys.add(bbox); geoms.add(geom.getEnvelope()); if(i==0){ polyUnion = bbox }else { polyUnion = polyUnion.union(bbox); } i++; } itr.close();
//shpReadWrite.CreateNewShapeFile
求并集。。
为什么要这么迂回呢?为什么不直接对原始多边形求并集,然后求外接矩形。。
List<Geometry> geoms = new ArrayList<Geometry>(); Geometry polyUnion = null; int i=0; while(itr.hasNext()){ SimpleFeature sf = itr.next(); Geometry geom = (Geometry) sf.getDefaultGeometry(); if(i==0){ Geometry geom2 = geom.intersection(geom);//避免自交 polyUnion = geom2; }else { Geometry geom2 = geom.intersection(geom);//避免自交 polyUnion = polyUnion.union(geom2); } i++; } itr.close();
s.geom.TopologyException: found non-noded intersection between LINESTRING........geotools空间分析存在自相交情况时的错误解决:先和自己自交,再求并集。。
>>更简单的办法:
SimpleFeatureCollection fc = contentFeatureSource.getFeatures();
double minX = fc.getBounds().getMinX();
...
相关文章
- 【技术种草】cdn+轻量服务器+hugo=让博客“云原生”一下
- CLB运维&运营最佳实践 ---访问日志大洞察
- vnc方式登陆服务器
- 轻松学排序算法:眼睛直观感受几种常用排序算法
- 十二个经典的大数据项目
- 为什么使用 CDN 内容分发网络?
- 大数据——大数据默认端口号列表
- Weld 1.1.5.Final,JSR-299 的框架
- JavaFX 2012:彻底开源
- 提升as3程序性能的十大要点
- 通过凸面几何学进行独立于边际的在线多类学习
- 利用行动影响的规律性和部分已知的模型进行离线强化学习
- ModelLight:基于模型的交通信号控制的元强化学习
- 浅谈Visual Source Safe项目分支
- 基于先验知识的递归卡尔曼滤波的代理人联合状态和输入估计
- 结合网络结构和非线性恢复来提高声誉评估的性能
- 最佳实践丨云开发CloudBase多环境管理实践
- TimeVAE:用于生成多变量时间序列的变异自动编码器
- 具有线性阈值激活的神经网络:结构和算法
- 内网渗透之横向移动 -- 从域外向域内进行密码喷洒攻击