如何在 openlayers 4.x 中获取多边形内的所有点？

Question

我想使用 Openlayers 集成一个功能 4.x 是我想获取地图上多边形内的所有点。目前我可以获得多边形本身的所有坐标。

但我想要多边形内的所有坐标或点。如果我解释得更多，那就意味着我希望地图上的所有点或坐标都被多边形区域包围。

Answer 1

将 vectorSource.forEachFeatureIntersectingExtent() 与 polygonGeom.getExtent() 结合使用。这将为您提供其中大多数的范围捷径。之后，您将需要自己实现 point-in-poly（网上有很多资源）或使用像 https://github.com/bjornharrtell/jsts 这样的库。 OpenLayers 仅提供与范围的几何交集。

Answer 2

我经常将 turf.js 库与 openlayers 放在一起，专门用于此类任务。我的大部分几何图形都是 geojson 中的 nativley，因此 turf.js 非常适合。如果你有一个 geojson FeatureCollection。您可以迭代 .features 数组（或者甚至只是 [x, y] 点的任何数组）并检查每一个以查看它是否在您的多边形内。如果有帮助，我可以做一个工作 fiddle。

// In this example I'm looking for all features
// that have AT LEAST ONE point within  
// the world extent (WGS84)
const polygon = turf.bboxPolygon([-180, -90, 180, 90])
myGeoJson.features.forEach(function(feature){
    const points = turf.explode(feature);
    const pointsWithin = turf.pointsWithinPolygon(points, polygon);
    if(pointsWithin.features && pointsWithin.features.length){
        // feature has AT LEAST ONE point inside the polygon
        // I can see what points by iterating the
        // pointsWithin.features array
    }else{
        // feature has ZERO points inside the polgyon
    }
});

http://turfjs.org/docs#pointsWithinPolygon

Answer 3

基于@willsters 的回答，如果已经确定了候选者并且您只是在寻找点（几何是范围），则可以在相反方向使用 forEachFeatureIntersectingExtent 来查看这些点是否与多边形的几何相交。

var candidates = [];
source.forEachFeatureIntersectingExtent(myPolygon.getGeometry().getExtent(),function(feature){
    if (feature.getGeometry().get('type') == 'Point') {
        candidates.push(feature);
    }
});

var selected = [];
candidates.forEach(function(candidate){
    source.forEachFeatureIntersectingExtent(candidate.getGeometry().getExtent(),function(feature){
        if (feature === myPolygon) {
            selected.push(candidate);
        }
    });
});

同样对于单个坐标点，我认为它可以一步完成：

var selected = [];
source.forEachFeatureIntersectingExtent(myPolygon.getGeometry().getExtent(),function(feature){
    if (feature.getGeometry().get('type') == 'Point' &&
        myPolygon.getGeometry().intersectsCoordinate(feature.getGeometry().get('coordinates')) {
            candidates.push(selected);
    }
});

关于划分为像这样的单元格会为包含多边形的 10x10 网格的每个单元格生成针推点。如果只有部分单元格与多边形相交，则单元格中心的针推可能在几何体之外。

var extent = myPolygon.getGeometry().getExtent();
for (var i=extent[0]; i<extent[2]; i+=(extent[2]-extent[0])/10) {
    for (var j=extent[1]; j<extent[3]; j+=(extent[3]-extent[1])/10) {
        var cellExtent = [i,j,i+(extent[2]-extent[0])/10),j+(extent[3]-extent[1])/10];
        source.forEachFeatureIntersectingExtent(cellExtent,function(feature){
            if (feature === myPolygon) {
               var pinPush = new ol.feature(new ol.geom.Point(ol.extent.getCenter(cellExtent)));
               source.addFeature(pinPush); 
            }
        });
    }
}

如何在 openlayers 4.x 中获取多边形内的所有点？

How to get all points inside a polygon in openlayers 4.x?

javascript

geometry

openlayers