@dra2020/baseclient 1.0.28 → 1.0.31

package/dist/geo/geo.d.ts

@@ -4,6 +4,12 @@ export declare type GeoProperties = geojson.GeoJsonProperties;
 export declare type GeoFeature = geojson.Feature;
 export declare type GeoFeatureArray = GeoFeature[];
 export declare type GeoFeatureCollection = geojson.FeatureCollection;
+export declare type GeoCentroidMap = {
+    [geoid: string]: {
+        x: number;
+        y: number;
+    };
+};
 export declare function geoEnsureID(col: GeoFeatureCollection): void;
 export declare function geoCollectionToMap(col: GeoFeatureCollection): GeoFeatureMap;
 export declare function geoMapToCollection(map: GeoFeatureMap): GeoFeatureCollection;

package/dist/poly/all.d.ts

@@ -14,3 +14,4 @@ export * from './topo';
 export * from './selfintersect';
 export * from './shamos';
 export * from './pointinpoly';
+export * from './mapto';

package/dist/poly/mapto.d.ts

@@ -0,0 +1,3 @@
+import * as G from '../geo/all';
+export declare function polyMapToByCentroid(districts: G.GeoFeatureCollection, centroids: G.GeoCentroidMap): any;
+export declare function polyMapTo(districts: G.GeoFeatureCollection, blocks: G.GeoFeatureCollection): any;

package/lib/geo/geo.ts

@@ -6,6 +6,7 @@ export type GeoProperties = geojson.GeoJsonProperties;
 export type GeoFeature = geojson.Feature;
 export type GeoFeatureArray = GeoFeature[];
 export type GeoFeatureCollection = geojson.FeatureCollection;
+export type GeoCentroidMap = { [geoid: string]: { x: number, y: number } };
 
 export function geoEnsureID(col: GeoFeatureCollection): void
 {
@@ -19,6 +20,11 @@ export function geoEnsureID(col: GeoFeatureCollection): void
     props.forEach(p => { if (prop === undefined && f.properties[p] !== undefined) prop = p; });
     if (prop)
       col.features.forEach(f => { f.properties.id = f.properties[prop] });
+    else
+    {
+      let n = 1;
+      col.features.forEach(f => { f.properties.id = String(n++) });
+    }
   }
 }
 

package/lib/poly/all.ts

@@ -14,3 +14,4 @@ export * from './topo';
 export * from './selfintersect';
 export * from './shamos';
 export * from './pointinpoly';
+export * from './mapto';

package/lib/poly/mapto.ts

@@ -0,0 +1,73 @@
+import * as G from '../geo/all';
+import * as P from './poly';
+import * as PP  from './polypack';
+import * as PL from './polylabel';
+import * as BB from './boundbox';
+import { polyContainsPoint } from './pointinpoly';
+
+function setLabels(c: G.GeoFeatureCollection): void
+{
+  c.features.forEach(f => {
+      if (f.properties.labelx === undefined)
+      {
+        let {x,y} = PL.polyLabel(f);
+        f.properties.labelx = x;
+        f.properties.labely = y;
+      }
+    });
+}
+
+// polyMapTo:
+//
+//  Given a set of underlying blocks (or precincts), map them to the district (or any other feature)
+//  they are contained in.
+//
+//  If a block maps to multiple districts or no district, it is left out of the result.
+//
+//  Note that the algorithm is to see where the centroid of the block overlaps. So a block
+//  could potentially overlap multiple districts and this would return only the one containing
+//  the centroid.
+//
+//  The thinking is that if you want fine granularity, use real blocks as the base collection.
+//  If you provide precincts, you will get a result on precinct granularity.
+//
+//  The return value is an object that maps the block feature ids to the district feature id.
+//
+export function polyMapToByCentroid(districts: G.GeoFeatureCollection, centroids: G.GeoCentroidMap): any
+{
+  let map: any = {};
+
+  // Cache district boundboxes for quick containment exclusion
+  let bbDistricts: BB.BoundBox[] = districts.features.map(f => BB.boundbox(f));
+
+  // Walk over blocks, mapping centroid to district
+  Object.keys(centroids).forEach(blockid => {
+      let x = centroids[blockid].x;
+      let y = centroids[blockid].y;
+
+      let fIn: G.GeoFeature[] = [];
+      districts.features.forEach((fDistrict: G.GeoFeature, i: number) => {
+          if (BB.boundboxContains(bbDistricts[i], x, y))
+            if (polyContainsPoint(fDistrict, x, y))
+              fIn.push(fDistrict);
+        });
+
+      if (fIn.length == 1)
+        map[blockid] = fIn[0].properties.id;
+    });
+
+  return map;
+}
+
+export function polyMapTo(districts: G.GeoFeatureCollection, blocks: G.GeoFeatureCollection): any
+{
+  // Cache labelx, labely if necessary
+  setLabels(blocks);
+
+  let centroids: G.GeoCentroidMap = {};
+  blocks.features.forEach(f => {
+    centroids[f.properties.id] = { x: f.properties.labelx, y: f.properties.labely }
+    });
+
+  return polyMapToByCentroid(districts, centroids);
+}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@dra2020/baseclient",
-  "version": "1.0.28",
+  "version": "1.0.31",
   "description": "Utility functions for Javascript projects.",
   "main": "dist/baseclient.js",
   "types": "./dist/all/all.d.ts",