Developing applications and services to make transit better for everyone - one train, bus, ferry, tram, metro ride at a time.
Public
Edited
Aug 12, 2024
1 fork
Step1_01_BufferedRoutes = {
const buffered_features = [];
const buffer_distance = 0.025;
for (const feature of Input_SydneyTransitLinesGeoJSON.features) {
if (_Imports.Geospatial.turf.length(feature, km) < 1000) {
const buffered = _Imports.Geospatial.turf.buffer(
feature,
buffer_distance,
km
);
buffered_features.push(buffered);
}
}
console.log("Finished Step1_01_BufferedRoutes");
return _Utils.Generate.GeoJSON("FeatureCollection", buffered_features, {});
}
JSTS_IO_HANDLERS = {
return {
gjreader: new _Imports.Geospatial.jsts.io.GeoJSONReader(),
gjwriter: new _Imports.Geospatial.jsts.io.GeoJSONWriter()
};
}
SydneyTrainRoutes_PreUnionsied_Geometries = {
const combined_routes = _Imports.Geospatial.turf.truncate(
Step1_01_BufferedRoutes,
{ precision: 5, coordinates: 2 }
);
return combined_routes.features.map((feature) =>
JSTS_IO_HANDLERS.gjreader.read(feature.geometry)
);
}
SydneyTrainRoutes_Unionised_Geometry = {
const geom = SydneyTrainRoutes_PreUnionsied_Geometries;
let ConstantGeometry = geom.shift();
for (const CurrentGeometry of geom) {
try {
ConstantGeometry = ConstantGeometry.union(CurrentGeometry);
console.log(
"Unionised of",
SydneyTrainRoutes_PreUnionsied_Geometries.length,
"geometries"
);
} catch (e) {}
}
const union_geometry = JSTS_IO_HANDLERS.gjwriter.write(ConstantGeometry);
const union_feature = _Utils.Generate.GeoJSON("Feature", union_geometry, {});
return union_feature;
}
SydneyTrainRoutes_SimplifiedUnionised_Geometry = _Imports.Geospatial.turf.simplify(
SydneyTrainRoutes_Unionised_Geometry,
{
tolerance: 0.0001,
highQuality: false
}
)
Step2_01_BboxUpperBounds = {
const UpperBounds = _Imports.Geospatial.turf.bboxPolygon(
_Imports.Geospatial.turf.bbox(SydneyTrainRoutes_SimplifiedUnionised_Geometry)
).geometry.coordinates[0][3];
console.log("Finished Step2_01_BboxUpperBounds");
return UpperBounds;
}
Step2_02_RoutesAsSVG = {
const RASVG = _Imports.Geospatial.geojson2svg
.default()
.styles(() => ({ fill: "black" }))
.projection(function (coord) {
return [
coord[0] + Step2_01_BboxUpperBounds[0] * -1,
coord[1] + Step2_01_BboxUpperBounds[1] * -1
];
})
.data(
_Imports.Geospatial.turf.rewind(
SydneyTrainRoutes_SimplifiedUnionised_Geometry
)
)
.render();
console.log("Finished Step2_02_RoutesAsSVG");
return RASVG;
}
Step2_03_RoutesPathCommands = {
const parser = new DOMParser();
const doc = parser.parseFromString(Step2_02_RoutesAsSVG, "application/xml");
console.log("Finished Step2_03_RoutesPathCommands");
return doc.querySelector(".polygon").getAttribute("d");
}
Step2_04_RoutesSAT = {
const paths = _Imports.Vector.mat.getPathsFromStr(
Step2_03_RoutesPathCommands
);
const mats = _Imports.Vector.mat.findMats(paths, 3);
const sats = mats.map((mat) => _Imports.Vector.mat.toScaleAxis(mat, 1.5));
console.log("Finished Step2_04_RoutesSAT");
return sats;
}
Step2_05_SATMultiLineStrings = {
const LsCoordinates = [];
const ReceivePathStr = (pathString) => {
if (pathString.length < 2) return;
const pathCoords = [pathString[0], pathString[pathString.length - 1]];
LsCoordinates.push([
[
pathCoords[0][0] - Step2_01_BboxUpperBounds[0] * -1,
pathCoords[0][1] - Step2_01_BboxUpperBounds[1] * -1
],
[
pathCoords[1][0] - Step2_01_BboxUpperBounds[0] * -1,
pathCoords[1][1] - Step2_01_BboxUpperBounds[1] * -1
]
]);
};
for (const MAT of Step2_04_RoutesSAT) {
if (!MAT.cpNode) return;
_Imports.Vector.mat.traverseEdges(MAT.cpNode, (cpNode) => {
if (_Imports.Vector.mat.isTerminating(cpNode)) return;
const bezier = _Imports.Vector.mat.getCurveToNext(cpNode);
if (bezier) ReceivePathStr(bezier);
});
}
console.log("Finished Step2_05_SATMultiLineStrings ");
return _Imports.Geospatial.turf.truncate(
_Imports.Geospatial.turf.multiLineString(LsCoordinates),
{ precision: 5, coordinates: 2 }
);
}
{
const container = yield htl.html`<div style="height: 500px;">`;
const map = L.map(container);
const layer = L.geoJSON(Step2_05_SATMultiLineStrings).addTo(map);
map.fitBounds(layer.getBounds(), { maxZoom: 9 });
L.tileLayer("https://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png", {
attribution:
"© <a href=https://www.openstreetmap.org/copyright>OpenStreetMap</a> contributors"
}).addTo(map);
}
Step3_01_CoordMapping = {
const coordMapping = new Map();
const coordUnmapping = new Map();
for (const coord of _Imports.Geospatial.turf.coordAll(
Step2_05_SATMultiLineStrings
)) {
if (!coordMapping.has(coord.toString())) {
const cid = coordMapping.size + 1;
coordMapping.set(coord.toString(), cid);
coordUnmapping.set(cid, coord);
}
}
console.log("Finished Step3_01_CoordMapping ");
return { coordMapping, coordUnmapping };
}
Step3_02_GraphNetwork = {
const graph = _Imports.Graph.ngraph.graph();
const coordMapping = Step3_01_CoordMapping.coordMapping;
for (const coord of [...coordMapping.values()]) {
graph.addNode(coord);
}
const dedupSegs = new Set();
const uniqueSegs = [];
for (const lstr of Step2_05_SATMultiLineStrings.geometry.coordinates) {
const segid1 = [
coordMapping.get(lstr[0].toString()),
coordMapping.get(lstr[1].toString())
].toString();
const segid2 = [
coordMapping.get(lstr[1].toString()),
coordMapping.get(lstr[0].toString())
].toString();
if (
!dedupSegs.has(segid1) &&
!dedupSegs.has(segid2) &&
coordMapping.get(lstr[0].toString()) !==
coordMapping.get(lstr[1].toString())
) {
dedupSegs.add(segid1);
dedupSegs.add(segid2);
uniqueSegs.push(lstr);
}
}
for (const lstr of uniqueSegs) {
graph.addLink(
coordMapping.get(lstr[0].toString()),
coordMapping.get(lstr[1].toString())
);
}
console.log("Finished Step3_02_GraphNetwork");
return _Imports.Graph.ngraph.todot.default(graph);
}
Step3_02_Hubs = {
const graph = _Imports.Graph.ngraph.fromdot.default(Step3_02_GraphNetwork);
const hubs = new Map();
const hubsToBeRemoved = [];
graph.forEachNode((node) => {
if (node.links.size > 2) {
hubs.set(node.id, new Set(node.links.entries()));
hubsToBeRemoved.push(node.id);
}
});
hubsToBeRemoved.map((h) => {
graph.removeNode(h);
});
console.log("Finished Step3_02_Hubs");
return { hubs, graph: _Imports.Graph.ngraph.todot.default(graph) };
}
Step3_025_Hub2Hub = {
const graph = _Imports.Graph.ngraph.fromdot.default(Step3_02_Hubs.graph);
const HubsRTG = new Map();
const connections = new Map();
const ConnectionsReplacement = new Map();
for (const hub of [...Step3_02_Hubs.hubs.entries()]) {
for (const link of [...hub[1].values()].map((link) =>
[link[0].fromId, link[0].toId].filter((v) => v !== hub[0])
)) {
if (HubsRTG.has(hub[0])) {
HubsRTG.set(hub[0], [...HubsRTG.get(hub[0]), link[0]]);
} else {
HubsRTG.set(hub[0], [link[0]]);
}
}
}
for (const hub of [...HubsRTG.entries()]) {
const connection = hub[1]
.map((v) => HubsRTG.has(v))
.map((v, i) => (v ? hub[1][i] : null))
.filter((v) => v !== null);
if (connection.length > 0) {
connections.set(hub[0], connection[0]);
graph.removeNode(connection[0]);
}
}
for (const connection of [...connections.entries()]) {
ConnectionsReplacement.set(
connection[0],
_Imports.Geospatial.turf.midpoint(
Step3_01_CoordMapping.coordUnmapping.get(connection[0]),
Step3_01_CoordMapping.coordUnmapping.get(connection[1])
).geometry.coordinates
);
}
console.log("Finished Step3_025_Hub2Hub");
return {
ConnectionsReplacement,
graph: _Imports.Graph.ngraph.todot.default(graph)
};
}
Step3_03_GraphWalked = {
const graph = _Imports.Graph.ngraph.fromdot.default(Step3_025_Hub2Hub.graph);
for (const hubId of [...Step3_02_Hubs.hubs.keys()]) {
graph.removeNode(hubId);
}
const paths = [];
const WalkToFindPath = (nodeIds) => {
const nextLink = [
...graph.getNode(nodeIds[nodeIds.length - 1]).links.values()
]
.map((link) => [link.fromId, link.toId])
.flat()
.filter((v) => !nodeIds.includes(v));
return nextLink[0];
};
graph.forEachNode((node) => {
if (node.links !== null && node.links.size == 1) {
if (
paths.filter((path) => path[path.length - 1] == node.id).length == 0
) {
const nodeIds = [node.id];
while (nodeIds[nodeIds.length - 1] !== undefined) {
nodeIds.push(WalkToFindPath(nodeIds));
}
paths.push(nodeIds.slice(0, nodeIds.length - 1));
}
}
});
console.log("Finished Step3_03_GraphWalked");
return paths;
}
Step3_04_AddHubs = {
const HubsRTG = new Map();
const paths = [];
for (const hub of [...Step3_02_Hubs.hubs.entries()]) {
for (const link of [...hub[1].values()]
.map((l) => l[0])
.map((l) => [l.fromId, l.toId])
.flat()
.filter((v) => v !== hub[0])) {
HubsRTG.set(link, hub[0]);
}
}
for (const incompletePath of Step3_03_GraphWalked) {
paths.push(
[
HubsRTG.get(incompletePath[0]),
...incompletePath,
HubsRTG.get(incompletePath[incompletePath.length - 1])
].filter((v) => v !== undefined)
);
}
console.log("Finished Step3_04_AddHubs");
return paths;
}
Step3_05_GetSegments = {
const LineStrings = [];
for (const path of Step3_04_AddHubs) {
LineStrings.push(
_Imports.Geospatial.turf.lineString(
path.map((id) => {
if (Step3_025_Hub2Hub.ConnectionsReplacement.has(id)) {
return Step3_025_Hub2Hub.ConnectionsReplacement.get(id);
} else {
return Step3_01_CoordMapping.coordUnmapping.get(id);
}
}),
{
id: LineStrings.length + 1
}
)
);
}
console.log("Finished Step3_05_GetSegments");
return _Imports.Geospatial.turf.truncate(
_Utils.Generate.GeoJSON("FeatureCollection", LineStrings, {}),
{ precision: 5, coordinates: 2 }
);
}
{
const container = yield htl.html`<div style="height: 500px;">`;
const map = L.map(container);
const layer = L.geoJSON(Step3_05_GetSegments, {
style: function () {
return {
color:
"#" +
Math.floor(Math.random() * 16777215)
.toString(16)
.padEnd(6, "0")
};
}
}).addTo(map);
map.fitBounds(layer.getBounds(), { maxZoom: 9 });
L.tileLayer("shttps://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png", {
attribution:
"© <a href=https://www.openstreetmap.org/copyright>OpenStreetMap</a> contributors"
}).addTo(map);
}
Purpose-built for displays of data
Observable is your go-to platform for exploring data and creating expressive data visualizations. Use reactive JavaScript notebooks for prototyping and a collaborative canvas for visual data exploration and dashboard creation.
