1c:{"title":"D3.js Potato Plan","cells":[{"id":0,"value":"# D3.js Potato Plan","pinned":false,"hidden":false,"mode":"md"},{"id":5,"value":"Testing how to create [potato plan-style maps](https://www.joostgrootens.nl/archive/1908_the_potato_plan_collection) in D3.js. Example below is based on [USGS NCLD](https://www.usgs.gov/centers/eros/science/national-land-cover-database) clipped to 100km around the Minneapolis-St.Paul area.\n\n

‼️ _This notebook is very much an experiment, not optimized for performance and not production-ready. Feedback and ideas for improvement are welcome._

🤖 _It was also created with the help of Claude Opus 4.7, but all mistakes are my own._

","pinned":false,"hidden":false,"mode":"md"},{"id":118,"value":"viewof targetRes = Inputs.range([180, 900], {step: 120, value: 540, label: html`Resolution (m)`})","pinned":false,"hidden":false,"mode":"ojs"},{"id":240,"value":"legend = {\n const codes = unique_values.filter(c => nlcdCodetoColor(c) !== \"none\");\n const perRow = width>800 ? 5 : 2;\n const colW = 200, rowH = 22, swatch = 16;\n const rows = Math.ceil(codes.length / perRow);\n\n const svg = d3.create(\"svg\")\n .attr(\"width\", perRow * colW)\n .attr(\"height\", rows * rowH + 6)\n .attr(\"style\", \"display:block; font:12px sans-serif;\");\n\n const items = svg.selectAll(\"g\")\n .data(codes)\n .join(\"g\")\n .attr(\"transform\", (d, i) =>\n `translate(${(i % perRow) * colW},${6 + Math.floor(i / perRow) * rowH})`);\n\n items.append(\"rect\")\n .attr(\"width\", swatch)\n .attr(\"height\", swatch)\n .attr(\"rx\", 3)\n .attr(\"fill\", d => nlcdCodetoColor(d))\n .attr(\"stroke\", \"#fff\").attr(\"stroke-width\", 0.5);\n\n items.append(\"text\")\n .attr(\"x\", swatch + 6)\n .attr(\"y\", swatch - 3)\n .attr(\"fill\", \"#333\")\n .text(d => nlcdCodetoClass(d));\n\n return svg.node();\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":116,"value":"$1e","pinned":false,"hidden":false,"mode":"ojs"},{"id":23,"value":"#### key variables","pinned":false,"hidden":false,"mode":"md"},{"id":41,"value":"unique_values = [...new Set(values)].sort((a, b) => a - b)","pinned":false,"hidden":false,"mode":"ojs"},{"id":38,"value":"values = (await image.readRasters())[0]","pinned":false,"hidden":false,"mode":"ojs"},{"id":60,"value":"n = image.getWidth()","pinned":false,"hidden":false,"mode":"ojs"},{"id":63,"value":"m = image.getHeight()","pinned":false,"hidden":false,"mode":"ojs"},{"id":131,"value":"scale_factor = Math.max(1, Math.round(targetRes / 30))","pinned":false,"hidden":false,"mode":"ojs"},{"id":29,"value":"#### helper functions","pinned":false,"hidden":false,"mode":"md"},{"id":68,"value":"contours = d3.contours()\n .size([image_resampled.width, image_resampled.height])\n .smooth(false)\n .thresholds([0.5])","pinned":false,"hidden":false,"mode":"ojs"},{"id":86,"value":"nlcdCodetoColor = function(nlcd_code){\n const map = {\n \"11\": \"#2A2B45\",\n \"21\": \"#E8CFC8\",\n \"22\": \"#E6b7A8\",\n \"23\": \"#E08F79\",\n \"24\": \"#DE7054\",\n \"31\": \"#F4F1EA\",\n \"41\": \"#A4B48C\",\n \"42\": \"#1C6330\",\n \"43\": \"#979B8E\",\n \"52\": \"#CCBA7C\",\n \"71\": \"#DFDFC2\",\n \"81\": \"#F5D37A\",\n \"82\": \"#E0D2CB\",\n \"90\": \"#A1BCCD\",\n \"95\": \"#7CAEAB\"\n };\n return map[nlcd_code] || nlcd_code;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":46,"value":"nlcdCodetoClass = function(nlcd_code){\n const map = {\n \"11\": \"Open Water\",\n \"21\": \"Developed, Open Space\",\n \"22\": \"Developed, Low Intensity\",\n \"23\": \"Developed, Medium Intensity\",\n \"24\": \"Developed, High Intensity\",\n \"31\": \"Barren Land\",\n \"41\": \"Deciduous Forest\",\n \"42\": \"Evergreen Forest\",\n \"43\": \"Mixed Forest\",\n \"52\": \"Shrub/Scrub\",\n \"71\": \"Grassland/Herbaceous\",\n \"81\": \"Pasture/Hay\",\n \"82\": \"Cultivated Crops\",\n \"90\": \"Woody Wetlands\",\n \"95\": \"Herbaceous Wetlands\"\n };\n return map[nlcd_code] || nlcd_code;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":173,"value":"blobPath = {\n const line = d3.line().curve(d3.curveBasisClosed);\n const step = 2; // higher = blobbier\n const minArea = 4; // higher = drop bigger stragglers\n return (poly) =>\n poly.coordinates\n .map(polygon => polygon\n .filter(ring => Math.abs(d3.polygonArea(ring)) >= minArea)\n .map(ring => {\n const pts = ring.filter((_, i) => i % step === 0);\n return pts.length > 2 ? line(pts) : line(ring);\n }).join(\"\"))\n .join(\"\");\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":325,"value":"roadProjection = {\n const [ox, oy] = image.getOrigin();\n const [px, py] = image.getResolution();\n const factor = targetRes / 30;\n return ([X, Y]) => [\n (X - ox) / px / factor,\n (Y - oy) / py / factor\n ];\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":27,"value":"#### data structures","pinned":false,"hidden":false,"mode":"md"},{"id":169,"value":"image_cleaned = {\n const { data, width: w, height: h } = image_resampled;\n const out = new Float64Array(w * h);\n for (let y = 0; y < h; y++) {\n for (let x = 0; x < w; x++) {\n const counts = new Map();\n for (let dy = -1; dy <= 1; dy++) {\n for (let dx = -1; dx <= 1; dx++) {\n const nx = x + dx, ny = y + dy;\n if (nx < 0 || ny < 0 || nx >= w || ny >= h) continue;\n const v = data[ny * w + nx];\n counts.set(v, (counts.get(v) || 0) + 1);\n }\n }\n let best = data[y * w + x], bestN = -1;\n for (const [v, c] of counts) if (c > bestN) { best = v; bestN = c; }\n out[y * w + x] = best;\n }\n }\n return { data: out, width: w, height: h };\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":144,"value":"image_resampled = {\n const factor = targetRes / 30;\n const dw = Math.ceil(n / factor);\n const dh = Math.ceil(m / factor);\n const out = new Float64Array(dw * dh);\n for (let y = 0; y < dh; y++) {\n for (let x = 0; x < dw; x++) {\n const counts = new Map();\n for (let dy = 0; dy < factor; dy++) {\n for (let dx = 0; dx < factor; dx++) {\n const sx = x * factor + dx, sy = y * factor + dy;\n if (sx >= n || sy >= m) continue;\n const v = values[sy * n + sx];\n counts.set(v, (counts.get(v) || 0) + 1);\n }\n }\n let best = 0, bestN = -1;\n for (const [v, c] of counts) if (c > bestN) { best = v; bestN = c; }\n out[y * dw + x] = best;\n }\n }\n return { data: out, width: dw, height: dh };\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":35,"value":"image = tiff.getImage()","pinned":false,"hidden":false,"mode":"ojs"},{"id":31,"value":"tiff = FileAttachment(\"Minneapolis_NLCD.tif\").arrayBuffer()\n .then(buffer => GeoTIFF.fromArrayBuffer(buffer))","pinned":false,"hidden":false,"mode":"ojs"},{"id":318,"value":"roads_clipped = {\n const b = raster_bbox;\n const inBox = ([x, y]) => x >= b.xmin && x <= b.xmax && y >= b.ymin && y <= b.ymax;\n const features = [];\n for (const f of roads.features) {\n const kept = f.geometry.coordinates.filter(seg => seg.some(inBox));\n if (kept.length) {\n features.push({\n type: \"Feature\",\n properties: f.properties,\n geometry: { type: \"MultiLineString\", coordinates: kept }\n });\n }\n }\n return { type: \"FeatureCollection\", features };\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":315,"value":"raster_bbox = {\n const [ox, oy] = image.getOrigin();\n const [px, py] = image.getResolution();\n return {\n xmin: ox, xmax: ox + px * n,\n ymax: oy, ymin: oy + py * m // py negative, so this is the bottom\n };\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":309,"value":"roads = FileAttachment(\"MN_roads.geojson\").json()","pinned":false,"hidden":false,"mode":"ojs"},{"id":21,"value":"### Dependencies","pinned":false,"hidden":false,"mode":"md"},{"id":18,"value":"GeoTIFF = require(\"geotiff@2\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":304,"value":"### Sources:\n[National Land Cover Database](https://www.usgs.gov/centers/eros/science/national-land-cover-database) from the United States Geological Service. [TIGER/Line Roads](https://www.census.gov/geographies/mapping-files/time-series/geo/tiger-line-file.html) from the United States Census Bureau.\n\n*All data retrieved 05/28/2026*","pinned":false,"hidden":false,"mode":"md"}],"files":"$W1f","stdlib":"1"}