1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768 |
- import turfbbox from "@turf/bbox";
- import booleanPointInPolygon from "@turf/boolean-point-in-polygon";
- import rbush from "rbush";
- /**
- * Merges a specified property from a FeatureCollection of points into a
- * FeatureCollection of polygons. Given an `inProperty` on points and an `outProperty`
- * for polygons, this finds every point that lies within each polygon, collects the
- * `inProperty` values from those points, and adds them as an array to `outProperty`
- * on the polygon.
- *
- * @name collect
- * @param {FeatureCollection<Polygon>} polygons polygons with values on which to aggregate
- * @param {FeatureCollection<Point>} points points to be aggregated
- * @param {string} inProperty property to be nested from
- * @param {string} outProperty property to be nested into
- * @returns {FeatureCollection<Polygon>} polygons with properties listed based on `outField`
- * @example
- * var poly1 = turf.polygon([[[0,0],[10,0],[10,10],[0,10],[0,0]]]);
- * var poly2 = turf.polygon([[[10,0],[20,10],[20,20],[20,0],[10,0]]]);
- * var polyFC = turf.featureCollection([poly1, poly2]);
- * var pt1 = turf.point([5,5], {population: 200});
- * var pt2 = turf.point([1,3], {population: 600});
- * var pt3 = turf.point([14,2], {population: 100});
- * var pt4 = turf.point([13,1], {population: 200});
- * var pt5 = turf.point([19,7], {population: 300});
- * var pointFC = turf.featureCollection([pt1, pt2, pt3, pt4, pt5]);
- * var collected = turf.collect(polyFC, pointFC, 'population', 'values');
- * var values = collected.features[0].properties.values
- * //=values => [200, 600]
- *
- * //addToMap
- * var addToMap = [pointFC, collected]
- */
- function collect(polygons, points, inProperty, outProperty) {
- var rtree = rbush(6);
- var treeItems = points.features.map(function (item) {
- var _a;
- return {
- minX: item.geometry.coordinates[0],
- minY: item.geometry.coordinates[1],
- maxX: item.geometry.coordinates[0],
- maxY: item.geometry.coordinates[1],
- property: (_a = item.properties) === null || _a === void 0 ? void 0 : _a[inProperty],
- };
- });
- rtree.load(treeItems);
- polygons.features.forEach(function (poly) {
- if (!poly.properties) {
- poly.properties = {};
- }
- var bbox = turfbbox(poly);
- var potentialPoints = rtree.search({
- minX: bbox[0],
- minY: bbox[1],
- maxX: bbox[2],
- maxY: bbox[3],
- });
- var values = [];
- potentialPoints.forEach(function (pt) {
- if (booleanPointInPolygon([pt.minX, pt.minY], poly)) {
- values.push(pt.property);
- }
- });
- poly.properties[outProperty] = values;
- });
- return polygons;
- }
- export default collect;
|