"use strict"; // get continuous paths (isolines) for all cells at once based on getType(cellId) comparison function getIsolines(graph, getType, options = {polygons: false, fill: false, halo: false, waterGap: false}) { const {cells, vertices} = graph; const isolines = {}; const checkedCells = new Uint8Array(cells.i.length); const addToChecked = cellId => (checkedCells[cellId] = 1); const isChecked = cellId => checkedCells[cellId] === 1; for (const cellId of cells.i) { if (isChecked(cellId) || !getType(cellId)) continue; addToChecked(cellId); const type = getType(cellId); const ofSameType = cellId => getType(cellId) === type; const ofDifferentType = cellId => getType(cellId) !== type; const onborderCell = cells.c[cellId].find(ofDifferentType); if (onborderCell === undefined) continue; // check if inner lake. Note there is no shoreline for grid features const feature = graph.features[cells.f[onborderCell]]; if (feature.type === "lake" && feature.shoreline?.every(ofSameType)) continue; const startingVertex = cells.v[cellId].find(v => vertices.c[v].some(ofDifferentType)); if (startingVertex === undefined) throw new Error(`Starting vertex for cell ${cellId} is not found`); const vertexChain = connectVertices({vertices, startingVertex, ofSameType, addToChecked, closeRing: true}); if (vertexChain.length < 3) continue; addIsoline(type, vertices, vertexChain); } return isolines; function addIsoline(type, vertices, vertexChain) { if (!isolines[type]) isolines[type] = {}; if (options.polygons) { if (!isolines[type].polygons) isolines[type].polygons = []; isolines[type].polygons.push(vertexChain.map(vertexId => vertices.p[vertexId])); } if (options.fill) { if (!isolines[type].fill) isolines[type].fill = ""; isolines[type].fill += getFillPath(vertices, vertexChain); } if (options.waterGap) { if (!isolines[type].waterGap) isolines[type].waterGap = ""; const isLandVertex = vertexId => vertices.c[vertexId].every(i => cells.h[i] >= 20); isolines[type].waterGap += getBorderPath(vertices, vertexChain, isLandVertex); } if (options.halo) { if (!isolines[type].halo) isolines[type].halo = ""; const isBorderVertex = vertexId => vertices.c[vertexId].some(i => cells.b[i]); isolines[type].halo += getBorderPath(vertices, vertexChain, isBorderVertex); } } } function getFillPath(vertices, vertexChain) { const points = vertexChain.map(vertexId => vertices.p[vertexId]); const firstPoint = points.shift(); return `M${firstPoint} L${points.join(" ")} Z`; } function getBorderPath(vertices, vertexChain, discontinue) { let discontinued = true; let lastOperation = ""; const path = vertexChain.map(vertexId => { if (discontinue(vertexId)) { discontinued = true; return ""; } const operation = discontinued ? "M" : "L"; const command = operation === lastOperation ? "" : operation; discontinued = false; lastOperation = operation; return ` ${command}${vertices.p[vertexId]}`; }); return path.join("").trim(); } // get single path for an non-continuous array of cells function getVertexPath(cellsArray) { const {cells, vertices} = pack; const cellsObj = Object.fromEntries(cellsArray.map(cellId => [cellId, true])); const ofSameType = cellId => cellsObj[cellId]; const ofDifferentType = cellId => !cellsObj[cellId]; const checkedCells = new Uint8Array(cells.c.length); const addToChecked = cellId => (checkedCells[cellId] = 1); const isChecked = cellId => checkedCells[cellId] === 1; let path = ""; for (const cellId of cellsArray) { if (isChecked(cellId)) continue; const onborderCell = cells.c[cellId].find(ofDifferentType); if (onborderCell === undefined) continue; const feature = pack.features[cells.f[onborderCell]]; if (feature.type === "lake" && feature.shoreline) { if (feature.shoreline.every(ofSameType)) continue; // inner lake } const startingVertex = cells.v[cellId].find(v => vertices.c[v].some(ofDifferentType)); if (startingVertex === undefined) throw new Error(`Starting vertex for cell ${cellId} is not found`); const vertexChain = connectVertices({vertices, startingVertex, ofSameType, addToChecked, closeRing: true}); if (vertexChain.length < 3) continue; path += getFillPath(vertices, vertexChain); } return path; } function getPolesOfInaccessibility(graph, getType) { const isolines = getIsolines(graph, getType, {polygons: true}); const poles = Object.entries(isolines).map(([id, isoline]) => { const multiPolygon = isoline.polygons.sort((a, b) => b.length - a.length); const [x, y] = polylabel(multiPolygon, 20); return [id, [rn(x), rn(y)]]; }); return Object.fromEntries(poles); } function connectVertices({vertices, startingVertex, ofSameType, addToChecked, closeRing}) { const MAX_ITERATIONS = vertices.c.length; const chain = []; // vertices chain to form a path let next = startingVertex; for (let i = 0; i === 0 || next !== startingVertex; i++) { const previous = chain.at(-1); const current = next; chain.push(current); const neibCells = vertices.c[current]; if (addToChecked) neibCells.filter(ofSameType).forEach(addToChecked); const [c1, c2, c3] = neibCells.map(ofSameType); const [v1, v2, v3] = vertices.v[current]; if (v1 !== previous && c1 !== c2) next = v1; else if (v2 !== previous && c2 !== c3) next = v2; else if (v3 !== previous && c1 !== c3) next = v3; if (next >= vertices.c.length) { ERROR && console.error("ConnectVertices: next vertex is out of bounds"); break; } if (next === current) { ERROR && console.error("ConnectVertices: next vertex is not found"); break; } if (i === MAX_ITERATIONS) { ERROR && console.error("ConnectVertices: max iterations reached", MAX_ITERATIONS); break; } } if (closeRing) chain.push(startingVertex); return chain; }