Cell.js 4.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126
  1. import {createBorderEdge} from "./Edge";
  2. import {cells, edges, epsilon} from "./Diagram";
  3. export function createCell(site) {
  4. return cells[site.index] = {
  5. site: site,
  6. halfedges: []
  7. };
  8. }
  9. function cellHalfedgeAngle(cell, edge) {
  10. var site = cell.site,
  11. va = edge.left,
  12. vb = edge.right;
  13. if (site === vb) vb = va, va = site;
  14. if (vb) return Math.atan2(vb[1] - va[1], vb[0] - va[0]);
  15. if (site === va) va = edge[1], vb = edge[0];
  16. else va = edge[0], vb = edge[1];
  17. return Math.atan2(va[0] - vb[0], vb[1] - va[1]);
  18. }
  19. export function cellHalfedgeStart(cell, edge) {
  20. return edge[+(edge.left !== cell.site)];
  21. }
  22. export function cellHalfedgeEnd(cell, edge) {
  23. return edge[+(edge.left === cell.site)];
  24. }
  25. export function sortCellHalfedges() {
  26. for (var i = 0, n = cells.length, cell, halfedges, j, m; i < n; ++i) {
  27. if ((cell = cells[i]) && (m = (halfedges = cell.halfedges).length)) {
  28. var index = new Array(m),
  29. array = new Array(m);
  30. for (j = 0; j < m; ++j) index[j] = j, array[j] = cellHalfedgeAngle(cell, edges[halfedges[j]]);
  31. index.sort(function(i, j) { return array[j] - array[i]; });
  32. for (j = 0; j < m; ++j) array[j] = halfedges[index[j]];
  33. for (j = 0; j < m; ++j) halfedges[j] = array[j];
  34. }
  35. }
  36. }
  37. export function clipCells(x0, y0, x1, y1) {
  38. var nCells = cells.length,
  39. iCell,
  40. cell,
  41. site,
  42. iHalfedge,
  43. halfedges,
  44. nHalfedges,
  45. start,
  46. startX,
  47. startY,
  48. end,
  49. endX,
  50. endY,
  51. cover = true;
  52. for (iCell = 0; iCell < nCells; ++iCell) {
  53. if (cell = cells[iCell]) {
  54. site = cell.site;
  55. halfedges = cell.halfedges;
  56. iHalfedge = halfedges.length;
  57. // Remove any dangling clipped edges.
  58. while (iHalfedge--) {
  59. if (!edges[halfedges[iHalfedge]]) {
  60. halfedges.splice(iHalfedge, 1);
  61. }
  62. }
  63. // Insert any border edges as necessary.
  64. iHalfedge = 0, nHalfedges = halfedges.length;
  65. while (iHalfedge < nHalfedges) {
  66. end = cellHalfedgeEnd(cell, edges[halfedges[iHalfedge]]), endX = end[0], endY = end[1];
  67. start = cellHalfedgeStart(cell, edges[halfedges[++iHalfedge % nHalfedges]]), startX = start[0], startY = start[1];
  68. if (Math.abs(endX - startX) > epsilon || Math.abs(endY - startY) > epsilon) {
  69. halfedges.splice(iHalfedge, 0, edges.push(createBorderEdge(site, end,
  70. Math.abs(endX - x0) < epsilon && y1 - endY > epsilon ? [x0, Math.abs(startX - x0) < epsilon ? startY : y1]
  71. : Math.abs(endY - y1) < epsilon && x1 - endX > epsilon ? [Math.abs(startY - y1) < epsilon ? startX : x1, y1]
  72. : Math.abs(endX - x1) < epsilon && endY - y0 > epsilon ? [x1, Math.abs(startX - x1) < epsilon ? startY : y0]
  73. : Math.abs(endY - y0) < epsilon && endX - x0 > epsilon ? [Math.abs(startY - y0) < epsilon ? startX : x0, y0]
  74. : null)) - 1);
  75. ++nHalfedges;
  76. }
  77. }
  78. if (nHalfedges) cover = false;
  79. }
  80. }
  81. // If there weren’t any edges, have the closest site cover the extent.
  82. // It doesn’t matter which corner of the extent we measure!
  83. if (cover) {
  84. var dx, dy, d2, dc = Infinity;
  85. for (iCell = 0, cover = null; iCell < nCells; ++iCell) {
  86. if (cell = cells[iCell]) {
  87. site = cell.site;
  88. dx = site[0] - x0;
  89. dy = site[1] - y0;
  90. d2 = dx * dx + dy * dy;
  91. if (d2 < dc) dc = d2, cover = cell;
  92. }
  93. }
  94. if (cover) {
  95. var v00 = [x0, y0], v01 = [x0, y1], v11 = [x1, y1], v10 = [x1, y0];
  96. cover.halfedges.push(
  97. edges.push(createBorderEdge(site = cover.site, v00, v01)) - 1,
  98. edges.push(createBorderEdge(site, v01, v11)) - 1,
  99. edges.push(createBorderEdge(site, v11, v10)) - 1,
  100. edges.push(createBorderEdge(site, v10, v00)) - 1
  101. );
  102. }
  103. }
  104. // Lastly delete any cells with no edges; these were entirely clipped.
  105. for (iCell = 0; iCell < nCells; ++iCell) {
  106. if (cell = cells[iCell]) {
  107. if (!cell.halfedges.length) {
  108. delete cells[iCell];
  109. }
  110. }
  111. }
  112. }