2f:T449,{ const outside = hexagons.filter(h3Index => { const bounds = h3.h3ToGeoBoundary(h3Index); return bounds.every(([lat, lng]) => h3.geoToH3(lat, lng, parentRes) !== parent) }).length; const parentArea = h3.cellArea(parent, h3.UNITS.km2); const fmtKm2 = (num) => num.toFixed(2) + ' km²'; const details = { 'Parent area': fmtKm2(parentArea), 'Total child area': fmtKm2(hexagons.reduce( (sum, h3Index) => sum + h3.cellArea(h3Index, h3.UNITS.km2), 0 )), 'Child area outside parent': fmtKm2(outsideArea), 'Child % outside parent': (outsideArea / parentArea * 100).toFixed(2) + '%', 'Children inside parent': hexagons.length - outside, 'Childen outside parent': outside }; return html` ${Object.keys(details).map(k => html` `)}

${k}	${details[k]}

` }30:[] 2d:{"title":"H3 Hierarchical (Non)Containment","cells":[{"id":0,"value":"md`# H3 Hierarchical (Non)Containment\n\nWhile H3 has a hierarchical relationship between cells at different resolutions, the hexagon cell shape doesn't allow for exact containment of a cell's descendants. Between each level, a fixed percentage of the children's area lies outside of the parent, and the same percentage of the parent's area is not included within its children. The margin of error is relatively stable across multiple resolutions. At > 2 resolutions finer than the parent, entire cells lie outside the parent's boundary.\n\nThis issue is an unavoidable consequence of using a hexagonal grid rather than a square one, a trade-off for the other advantages of using hexagons. You can find [more discussion of these tradeoffs in the docs](https://h3geo.org/docs/usecases).\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":59,"value":"viewof childDepth = slider({\n min: 0, \n max: 5, \n step: 1, \n value: 1, \n title: \"Child Depth\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":5,"value":"viewof map = {\n const container = html`

`;\n yield container; // Give the container dimensions.\n \n const map = container.value = new mapboxgl.Map({\n container,\n center: [config.lng, config.lat],\n zoom: config.zoom,\n style: 'mapbox://styles/mapbox/light-v9'\n });\n \n map.on('load', () => {\n container.value = map;\n container.dispatchEvent(new CustomEvent(\"input\"));\n });\n \n invalidation.then(() => map.remove());\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":121,"value":"$2f","pinned":false,"hidden":false,"mode":"ojs"},{"id":204,"value":"md`---`","pinned":false,"hidden":false,"mode":"ojs"},{"id":73,"value":"parentRes = h3.h3GetResolution(parent)","pinned":false,"hidden":false,"mode":"ojs"},{"id":221,"value":"h3Resolution = parentRes + childDepth","pinned":false,"hidden":false,"mode":"ojs"},{"id":93,"value":"parent = \"851fa093fffffff\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":32,"value":"hexagons = {\n return h3.h3ToChildren(parent, h3Resolution);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":233,"value":"parentLatLng = h3.h3ToGeo(parent)","pinned":true,"hidden":false,"mode":"ojs"},{"id":46,"value":"config = ({\n lat: parentLatLng[0],\n lng: parentLatLng[1],\n zoom: 6\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":25,"value":"{\n renderH3IndexesToMap(map, \"hexes\", hexagons, {\n \"fill-color\": \"steelblue\",\n \"fill-opacity\": 0.2\n });\n renderH3IndexesToMap(map, \"parent\", [parent], {\n \"fill-color\": \"red\",\n \"fill-opacity\": 0.2\n });\n\n // Uncomment to double-check the difference\n // renderMapLayer({\n // map,\n // id: 'outside',\n // data: outsideDiff,\n // type: 'fill',\n // paint: {\n // 'fill-color': 'black',\n // 'fill-opacity': 0.4\n // }\n // });\n\n return md`*{map rendering logic}*`;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":169,"value":"features = ({\n children: geojson2h3.h3SetToMultiPolygonFeature(hexagons),\n parent: geojson2h3.h3ToFeature(parent)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":176,"value":"outsideDiff = {\n const parentFeature = geojson2h3.h3ToFeature(parent);\n return {\n type: 'FeatureCollection',\n features: hexagons.map(h3Index => difference(geojson2h3.h3ToFeature(h3Index), parentFeature))\n .filter(Boolean)\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":179,"value":"outsideArea = getArea(outsideDiff) / 1e6;","pinned":false,"hidden":false,"mode":"ojs"},{"id":241,"value":"{\n console.time(\"perf\");\n for (let i = 0; i < 10000000; i++) {\n h3.geoToH3(parentLatLng, 11);\n }\n console.timeEnd(\"perf\");\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":6,"value":"md`## Dependencies`","pinned":false,"hidden":false,"mode":"ojs"},{"id":10,"value":"h3 = require('h3-js@3.7.2')","pinned":false,"hidden":false,"mode":"ojs"},{"id":16,"value":"import {mapboxApiToken} from '@nrabinowitz/mapbox-token'","pinned":false,"hidden":false,"mode":"ojs"},{"id":28,"value":"import {renderH3IndexesToMap, renderMapLayer, h3ToGeoJSON} from '@nrabinowitz/mapbox-utils'","pinned":false,"hidden":false,"mode":"ojs"},{"id":13,"value":"mapboxgl = {\n let mapboxgl = await require('mapbox-gl');\n mapboxgl.accessToken = mapboxApiToken;\n return mapboxgl;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":56,"value":"import {slider} from \"@jashkenas/inputs\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":150,"value":"geojson2h3 = require('https://bundle.run/geojson2h3@1.2.0')","pinned":false,"hidden":false,"mode":"ojs"},{"id":140,"value":"difference = require('https://bundle.run/@turf/difference@6.0.1')","pinned":false,"hidden":false,"mode":"ojs"},{"id":142,"value":"getArea = {\n const module = await require('https://bundle.run/@turf/area@6.0.1');\n return module.default;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":18,"value":"// CSS import required for the Mapbox map\nhtml``","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W30","stdlib":"1","titling":"auto"}