2b:null 2c:["$","title",null,{"children":"Voronoi Incircles Trippy Edition | Christopher Pietsch | Observable"}] 31:T578,canvas = { const context = DOM.context2d(width, height); while (true) { context.fillStyle = "#00000002"; 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(); context.beginPath(); //voronoi.render(context); //voronoi.renderBounds(context); context.strokeStyle = "#bbb"; context.stroke(); context.beginPath(); //voronoi.delaunay.renderPoints(context); context.fillStyle = "#000"; 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(); } }32:[] 2f:{"id":"e0b33a28a750de84","slug":"voronoi-incircles-trippy-edition","title":"Voronoi Incircles Trippy Edition","thumbnail":"2d29a885b9c8d5714961d81ac46cdb13c0c91ce2ac76ee576649938a38e8955a","access_level":"public","create_time":"$D2021-09-15T12:46:20.942Z","update_time":"$D2021-09-15T13:02:50.597Z","view_time":null,"first_public_time":"$D2021-09-15T12:46:47.155Z","fork":{"id":"ba0d23af0509620f","slug":"voronoi-incircles","title":"Voronoi Incircles","owner":{"id":"8c2e47252fd4995c","login":"mbostock","name":"Mike Bostock"}},"stars":1,"starred":false,"owner":{"type":"team","id":"dbaac3ef74ea8f0c","login":"cpietsch","name":"Christopher Pietsch","avatar_url":"https://avatars-next.observableusercontent.com/5c893db7cf5c2699bda61c80eaa28fc9ebd560ec948b9bec441e6eefe34c25f8","tier":"starter_2024"},"version":352,"publish_version":352,"role":null,"authored":false,"authors":[{"type":"individual","id":"42a00b3d53c425cb","login":"cpietsch","name":"Christopher Pietsch","avatar_url":"https://avatars-next.observableusercontent.com/5c893db7cf5c2699bda61c80eaa28fc9ebd560ec948b9bec441e6eefe34c25f8","tier":null}],"members":[{"type":"individual","id":"42a00b3d53c425cb","login":"cpietsch","name":"Christopher Pietsch","avatar_url":"https://avatars-next.observableusercontent.com/5c893db7cf5c2699bda61c80eaa28fc9ebd560ec948b9bec441e6eefe34c25f8","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"subscribed":null,"forks":1,"collection":null,"type":"notebook","body":{"title":"Voronoi Incircles Trippy Edition","cells":[{"id":0,"value":"md`# Voronoi Incircles Trippy Edition\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":"$31","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":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":"$W32","stdlib":"1","titling":"auto"}} 30:"$2f:body"