[英]How to split a large JTS polygon into multiple 1 km by 1 km geometries?
[英]How to split a JTS polygon
我有一個大多邊形,我想找到與多邊形相交的特征,但由於多邊形太大,我遇到超時異常。
我試圖研究 JTS 方法,但不知道如何使用它。
final List<Coordinate> coordinates = List.of(new Coordinate(0, 0), new Coordinate(-1, 1),
new Coordinate(1, 3), new Coordinate(2, 3), new Coordinate(3, 1), new Coordinate(0, 0));
final GeometryFactory factory = new GeometryFactory();
final Polygon polygon = factory.createPolygon(coordinates.toArray(new Coordinate[0]));
final Geometry envelope = polygon.getEnvelope();
有人可以指導我如何拆分多邊形 object 嗎?
有許多(無數種)方法可以分割多邊形,但我會這樣做,通過計算包絡線的中線並將新包絡線與多邊形相交。
public List<Geometry> split(Polygon p) {
List<Geometry> ret = new ArrayList<>();
final Envelope envelope = p.getEnvelopeInternal();
double minX = envelope.getMinX();
double maxX = envelope.getMaxX();
double midX = minX + (maxX - minX) / 2.0;
double minY = envelope.getMinY();
double maxY = envelope.getMaxY();
double midY = minY + (maxY - minY) / 2.0;
Envelope llEnv = new Envelope(minX, midX, minY, midY);
Envelope lrEnv = new Envelope(midX, maxX, minY, midY);
Envelope ulEnv = new Envelope(minX, midX, midY, maxY);
Envelope urEnv = new Envelope(midX, maxX, midY, maxY);
Geometry ll = JTS.toGeometry(llEnv).intersection(p);
Geometry lr = JTS.toGeometry(lrEnv).intersection(p);
Geometry ul = JTS.toGeometry(ulEnv).intersection(p);
Geometry ur = JTS.toGeometry(urEnv).intersection(p);
ret.add(ll);
ret.add(lr);
ret.add(ul);
ret.add(ur);
return ret;
}
這為您的多邊形提供了這個:
如果您在 output 上再次調用它:
在生產設置中,您需要進行一些錯誤檢查以確保您可以處理生成的點。
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.