29:null 2a:["$","title",null,{"children":"Voronoi Incircles | Aleksander | Observable"}] 2f:T462,canvas = { const context = DOM.context2d(width, height); while (true) { context.fillStyle = "#eee"; context.fillRect(0, 0, width, height); context.beginPath(); for (const cell of voronoi.cellPolygons()) { const [x, y, r] = polygonIncircle(cell); if (r < 1.5) continue; context.moveTo(x + r - 1.5, y); context.arc(x, y, r - 1.5, 0, 2 * Math.PI); } context.fillStyle = "#fff"; context.fill(); yield context.canvas; for (let i = 0; i < n; ++i) { const x = i << 1; const y = x + 1; positions[x] += velocities[x]; positions[y] += velocities[y]; if (positions[x] < -margin) positions[x] += width + margin * 2; else if (positions[x] > width + margin) positions[x] -= width + margin * 2; if (positions[y] < -margin) positions[y] += height + margin * 2; else if (positions[y] > height + margin) positions[y] -= height + margin * 2; velocities[x] += 0.2 * (Math.random() - 0.5) - 0.01 * velocities[x]; velocities[y] += 0.2 * (Math.random() - 0.5) - 0.01 * velocities[y]; } voronoi.update(); } }30:[{"id":"2d747586c41e42c49f2d57b9d87da83a409e9cbf1a68be3e4f9be39a54ada1f982d68b0c8bc21413e852d706fe7a9cb3ea9ef9edcf364b3cde7a1a4e03c30143","href":"https://static.observableusercontent.com/files/2d747586c41e42c49f2d57b9d87da83a409e9cbf1a68be3e4f9be39a54ada1f982d68b0c8bc21413e852d706fe7a9cb3ea9ef9edcf364b3cde7a1a4e03c30143","name":"Meisje_met_de_parel.jpg","type":"image/jpeg","lastModified":1676968369847.541,"size":114245}] 2d:{"id":"c7879a54e79dc540","slug":null,"title":"Voronoi Incircles","thumbnail":"bb053e5d9aea0bd6342c49a4ce4e835a3f9f5b57c6c0be64cb20795b11a81772","access_level":"public","create_time":"$D2023-02-21T08:08:07.301Z","update_time":"$D2023-02-21T08:38:55.795Z","view_time":null,"first_public_time":"$D2023-02-21T08:08:07.301Z","fork":{"id":"ba0d23af0509620f","slug":"voronoi-incircles","title":"Voronoi Incircles","owner":{"id":"8c2e47252fd4995c","login":"mbostock","name":"Mike Bostock"}},"stars":0,"starred":false,"owner":{"type":"team","id":"f8bd23fb65b0c2a8","login":"fallenartist","name":"Aleksander","avatar_url":"https://avatars-next.observableusercontent.com/bfafc2ba125265418f3da4af2d2aed36cdd227f1b2db87e9da76a4a3add844d6","tier":"starter_2024"},"version":347,"publish_version":null,"role":null,"authored":false,"authors":[{"type":"individual","id":"9b54e3c43bd8e2da","login":"fallenartist","name":"Aleksander","avatar_url":"https://avatars-next.observableusercontent.com/bfafc2ba125265418f3da4af2d2aed36cdd227f1b2db87e9da76a4a3add844d6","tier":null}],"members":[{"type":"individual","id":"9b54e3c43bd8e2da","login":"fallenartist","name":"Aleksander","avatar_url":"https://avatars-next.observableusercontent.com/bfafc2ba125265418f3da4af2d2aed36cdd227f1b2db87e9da76a4a3add844d6","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"subscribed":null,"forks":0,"collection":null,"type":"notebook","body":{"title":"Voronoi Incircles","cells":[{"id":0,"value":"md`# Voronoi Incircles\n\nFor each Voronoi cell, the [largest contained circle](/@mbostock/convex-polygon-incircle) is shown in white.`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":6,"value":"$2f","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":314,"value":"positions = Float64Array.from({length: n * 2}, (_, i) => Math.random() * (i & 1 ? height : width))","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":315,"value":"velocities = new Float64Array(n * 2)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":317,"value":"voronoi = new d3.Delaunay(positions).voronoi([0.5, 0.5, width - 0.5, height - 0.5])","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":24,"value":"function polygonIncircle(polygon) {\n let circle = [NaN, NaN, 0];\n for (let i = 0, n = polygon.length - 1; i < n; ++i) {\n const pi0 = polygon[i], pi1 = polygon[i + 1];\n for (let j = i + 1; j < n; ++j) {\n const pj0 = polygon[j], pj1 = polygon[j + 1];\n search: for (let k = j + 1; k < n; ++k) {\n const pk0 = polygon[k], pk1 = polygon[k + 1];\n const c = circleTangent(pi0, pi1, pj0, pj1, pk0, pk1);\n if (!(c[2] > circle[2])) continue;\n for (let l = 0; l < n; ++l) {\n if (l === i || l === j || l === k) continue;\n const d = pointLineDistance(c, polygon[l], polygon[l + 1]);\n if (d + 1e-6 < c[2]) continue search;\n }\n circle = c;\n }\n }\n }\n return circle;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":27,"value":"function circleTangent(p0, p1, p2, p3, p4, p5) {\n const b0 = lineLineBisect(p0, p1, p3, p2);\n const b1 = lineLineBisect(p2, p3, p5, p4);\n const i = lineLineIntersect(...b0, ...b1);\n return [...i, pointLineDistance(i, p0, p1)];\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":62,"value":"function pointLineDistance([x0, y0], [x2, y2], [x1, y1]) {\n const x21 = x2 - x1, y21 = y2 - y1;\n return (y21 * x0 - x21 * y0 + x2 * y1 - y2 * x1) / Math.sqrt(y21 * y21 + x21 * x21);\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":255,"value":"function lineLineBisect([x0, y0], [x1, y1], [x2, y2], [x3, y3]) {\n const x02 = x0 - x2, y02 = y0 - y2;\n const x10 = x1 - x0, y10 = y1 - y0, l10 = Math.sqrt(x10 ** 2 + y10 ** 2);\n const x32 = x3 - x2, y32 = y3 - y2, l32 = Math.sqrt(x32 ** 2 + y32 ** 2);\n const ti = (x32 * y02 - y32 * x02) / (y32 * x10 - x32 * y10);\n const xi = x0 + ti * x10, yi = y0 + ti * y10;\n return [[xi, yi], [xi + x10 / l10 + x32 / l32, yi + y10 / l10 + y32 / l32]];\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":29,"value":"function lineLineIntersect([x0, y0], [x1, y1], [x2, y2], [x3, y3]) {\n const x02 = x0 - x2, y02 = y0 - y2;\n const x10 = x1 - x0, y10 = y1 - y0;\n const x32 = x3 - x2, y32 = y3 - y2;\n const t = (x32 * y02 - y32 * x02) / (y32 * x10 - x32 * y10);\n return [x0 + t * x10, y0 + t * y10];\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":336,"value":"samplePoint = (p, ctx) => {\n let sampleSize = 10;\n let imageData = ctx.getImageData(p[0], p[1], sampleSize, sampleSize);\n let rgb = { r: 0, g: 0, b: 0 };\n let len = imageData.data.length;\n let s = 0;\n let i = 0;\n while (i < len) {\n rgb.r += imageData.data[i + 0];\n rgb.g += imageData.data[i + 1];\n rgb.b += imageData.data[i + 2];\n i += 4;\n s += 1;\n }\n rgb.r = rgb.r / s;\n rgb.g = rgb.g / s;\n rgb.b = rgb.b / s;\n\n return `rgb(${rgb.r}, ${rgb.g}, ${rgb.b})`;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":338,"value":"image = FileAttachment(\"Meisje_met_de_parel.jpg\").image();","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":344,"value":"cols = d3.scaleOrdinal(d3.schemeAccent)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":13,"value":"height = 600","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":320,"value":"margin = 60","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":150,"value":"n = 100","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":3,"value":"d3 = require(\"d3-delaunay@5\")","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"}],"files":"$W30","stdlib":"1","titling":"auto"}} 2e:"$2d:body"