2f:[{"id":"8a259fc1bd96ec206db5ec13f6dd8773c0af4ef879a5f51328b4e9988006346db328ad1e59065ee1687f4acceff0bb28879cbb90feb7db00fc854e0963c3a158","href":"https://static.observableusercontent.com/files/8a259fc1bd96ec206db5ec13f6dd8773c0af4ef879a5f51328b4e9988006346db328ad1e59065ee1687f4acceff0bb28879cbb90feb7db00fc854e0963c3a158","name":"tl_2024_us_county@1.json","type":"application/json","lastModified":1730927424117.737,"size":2290470},{"id":"431e0e64a134bec795262e437f8a2a661c7b92f73af547b2fb07d375547ea30028354a77b9c4ee3d5ef6511f49b2ae6ccda41aaf0b1dfc9bd05cd456105e3d8a","href":"https://static.observableusercontent.com/files/431e0e64a134bec795262e437f8a2a661c7b92f73af547b2fb07d375547ea30028354a77b9c4ee3d5ef6511f49b2ae6ccda41aaf0b1dfc9bd05cd456105e3d8a","name":"ne_50m_admin_0_countries_lakes.json","type":"application/json","lastModified":1730929617055.776,"size":63067}] 2d:{"title":"Clip counties file","cells":[{"id":0,"value":"# Clip counties file","pinned":false,"hidden":false,"mode":"md"},{"id":76,"value":"htl.html`\n

\n ${pct.toFixed(1)}%\n \n ${d3.format(\",\")(countyGeoClipped.features.length)} of ${d3.format(\",\")(countyGeo.features.length)} counties\n

`","pinned":false,"hidden":false,"mode":"ojs"},{"id":19,"value":"{\n const svg = d3.create(\"svg\")\n .attr(\"width\", width)\n .attr(\"height\", height);\n\n svg.selectAll(\".county\")\n .data(countyGeo.features)\n .join(\"path\")\n .attr(\"class\", \"county\")\n .attr(\"d\", path)\n .attr(\"fill\", \"#f5f5f5\")\n .attr(\"stroke\", \"#c0c0c0\");\n\n svg.selectAll(\".outline\")\n .data(usGeo.features)\n .join(\"path\")\n .attr(\"class\", \"outline\")\n .attr(\"d\", path)\n .attr(\"fill\", \"none\")\n .attr(\"stroke\", \"black\");\n\n svg.selectAll(\".clipped\")\n .data(countyGeoClipped.features)\n .join(\"path\")\n .attr(\"class\", \"clipped\")\n .attr(\"d\", path)\n .attr(\"fill\", \"lightblue\")\n .attr(\"stroke\", \"steelblue\");\n \n return svg.node();\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":17,"value":"path = d3.geoPath(projection)","pinned":true,"hidden":false,"mode":"ojs"},{"id":15,"value":"projection = d3.geoAlbersUsa().fitSize([width, height], countyGeo)","pinned":true,"hidden":false,"mode":"ojs"},{"id":12,"value":"height = width * 0.63","pinned":true,"hidden":false,"mode":"ojs"},{"id":7,"value":"countyTopo = FileAttachment(\"tl_2024_us_county@1.json\").json()","pinned":true,"hidden":false,"mode":"ojs"},{"id":9,"value":"countyGeo = ({\n type: \"FeatureCollection\",\n features: topojson\n .feature(countyTopo, countyTopo.objects.tl_2024_us_county)\n .features\n .sort((a, b) => d3.ascending(+a.properties.INTPTLON, +b.properties.INTPTLON))\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":53,"value":"usTopo = FileAttachment(\"ne_50m_admin_0_countries_lakes.json\").json()","pinned":true,"hidden":false,"mode":"ojs"},{"id":55,"value":"usGeo = topojson.feature(usTopo, usTopo.objects.ne_50m_admin_0_countries_lakes)","pinned":true,"hidden":false,"mode":"ojs"},{"id":39,"value":"countyGeoClipped = {\n const out = {\n type: \"FeatureCollection\",\n features: []\n }\n\n for (let i = 0, l = countyGeo.features.length; i < l; i++) {\n const feature = countyGeo.features[i];\n\n out.features.push(geoIntersection(usGeo.features[0], feature))\n \n if (i % 10 === 0) {\n yield out;\n }\n }\n\n yield out;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":81,"value":"pct = countyGeoClipped.features.length / countyGeo.features.length * 100","pinned":true,"hidden":false,"mode":"ojs"},{"id":33,"value":"## Utilities","pinned":false,"hidden":false,"mode":"md"},{"id":35,"value":"function geoIntersection(clipper, feature){\n const union = clip.intersection(\n feature.geometry.coordinates,\n clipper.geometry.coordinates\n );\n\n const clipped = rewind({\n type: \"Feature\",\n properties: feature.properties,\n geometry: {\n type: union.length === 1 ? \"Polygon\" : \"MultiPolygon\",\n coordinates: union.length === 1 ? union[0] : union\n }\n }, 1);\n\n clipped.properties.pct = d3.geoArea(clipped) / d3.geoArea(feature)\n\n return clipped\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":29,"value":"## Imports","pinned":false,"hidden":false,"mode":"md"},{"id":31,"value":"clip = require(\"polygon-clipping@0.15.7\")","pinned":true,"hidden":false,"mode":"ojs"},{"id":37,"value":"rewind = require(\"https://bundle.run/geojson-rewind@0.3.1\");","pinned":true,"hidden":false,"mode":"ojs"}],"files":"$W2f","stdlib":"1","titling":"auto"}