25:["$","title",null,{"children":"Air Travel Routing | Joe Weickert | Observable"}] 49:T430,airTravelEdges = { const edges = []; for (let r = 0; r < nestedVertices.length; r++) { for (let c = 0; c < nestedVertices[r].length; c++) { const rightNeighbor = nestedVertices[r][(c + 1) % nestedVertices[r].length]; edges.push([[nestedVertices[r][c].lng, nestedVertices[r][c].lat], [rightNeighbor.lng, rightNeighbor.lat]]); if (nestedVertices[(r + 1)]) { const bottomNeighbor = nestedVertices[(r + 1)][c]; edges.push([[nestedVertices[r][c].lng, nestedVertices[r][c].lat], [bottomNeighbor.lng, bottomNeighbor.lat]]); const bottomRightNeighbor = nestedVertices[(r + 1)][(c + 1) % nestedVertices[r].length]; edges.push([[nestedVertices[r][c].lng, nestedVertices[r][c].lat], [bottomRightNeighbor.lng, bottomRightNeighbor.lat]]); const bottomLeftNeighbor = nestedVertices[(r + 1)][(c - 1 < 0) ? nestedVertices[r].length - 1 : c - 1]; edges.push([[nestedVertices[r][c].lng, nestedVertices[r][c].lat], [bottomLeftNeighbor.lng, bottomLeftNeighbor.lat]]); } } } return edges; }4a:T78b,verticesWithNeighbors = { const vertices = []; for (let r = 0; r < nestedVertices.length; r++) { for (let c = 0; c < nestedVertices[r].length; c++) { const current = nestedVertices[r][c]; const rightNeighbor = nestedVertices[r][(c + 1) % nestedVertices[r].length]; const rDistance = Math.round(haversineDistance([current.lng, current.lat], [rightNeighbor.lng, rightNeighbor.lat], true) * 100) / 100; current.neighbors.push({id: rightNeighbor.id, distance: rDistance}); rightNeighbor.neighbors.push({id: current.id, distance: rDistance}); if (nestedVertices[(r + 1)]) { const bottomNeighbor = nestedVertices[(r + 1)][c]; const bDistance = Math.round(haversineDistance([current.lng, current.lat], [bottomNeighbor.lng, bottomNeighbor.lat], true) * 100) / 100; current.neighbors.push({id: bottomNeighbor.id, distance: bDistance}); bottomNeighbor.neighbors.push({id: current.id, distance: bDistance}); const bottomRightNeighbor = nestedVertices[(r + 1)][(c + 1) % nestedVertices[r].length]; const brDistance = Math.round(haversineDistance([current.lng, current.lat], [bottomRightNeighbor.lng, bottomRightNeighbor.lat], true) * 100) / 100; current.neighbors.push({id: bottomRightNeighbor.id, distance: brDistance}); bottomRightNeighbor.neighbors.push({id: current.id, distance: brDistance}); const bottomLeftNeighbor = nestedVertices[(r + 1)][(c - 1 < 0) ? nestedVertices[r].length - 1 : c - 1]; const blDistance = Math.round(haversineDistance([current.lng, current.lat], [bottomLeftNeighbor.lng, bottomLeftNeighbor.lat], true) * 100) / 100; current.neighbors.push({id: bottomLeftNeighbor.id, distance: blDistance}); bottomLeftNeighbor.neighbors.push({id: current.id, distance: blDistance}); } vertices.push(current); } } return vertices; }4b:[{"id":"2d2882bb161fcc99e2dad3a1a21830db287167e518fb87a5ded7a51624728538f03f4c16b24ef19b4ca56fe8ec4a81167c3aefdc4c1ca70fe9094e28614bcdc7","href":"https://static.observableusercontent.com/files/2d2882bb161fcc99e2dad3a1a21830db287167e518fb87a5ded7a51624728538f03f4c16b24ef19b4ca56fe8ec4a81167c3aefdc4c1ca70fe9094e28614bcdc7","name":"AirPorts@1.json","type":"application/json","lastModified":1701700782104.393,"size":236816},{"id":"ab156def1f2547ae49642efdfe7628de0164a62f991ddf5209ed57b3b7fde1eae55897207b044500274f4671b10e865f2d1e49e7f1e4a394ad11bfe1b0586815","href":"https://static.observableusercontent.com/files/ab156def1f2547ae49642efdfe7628de0164a62f991ddf5209ed57b3b7fde1eae55897207b044500274f4671b10e865f2d1e49e7f1e4a394ad11bfe1b0586815","name":"countries-50m.json","type":"application/json","lastModified":1690894307901.481,"size":756419}] 29:{"id":"0ad28c7c1b2c17e1","slug":null,"title":"Air Travel Routing","thumbnail":"c2d273a8edd07d210cf477031e730efee62b77ff8a9c97d9c3e6beac4bcac48c","access_level":"public","create_time":"$D2023-08-01T12:41:20.679Z","update_time":"$D2024-02-08T20:16:54.800Z","view_time":null,"first_public_time":"$D2023-08-01T12:41:20.679Z","has_unpublished_edits":null,"fork":{"id":"19cbc978f7478a7b","slug":"plot-world-map","title":"Plot World Map","owner":{"id":"f35c755083683fe5","login":"observablehq","name":"Observable"}},"stars":0,"starred":false,"owner":{"type":"team","id":"164dbf06be59d75e","login":"weickertjoe","name":"Joe Weickert","avatar_url":"https://avatars-next.observableusercontent.com/1169610d6fe2ae12a74ecea51bb14bd3a72502315e48c8bfdf27d06a3798b86e","tier":"starter_2024","bio":"","followed_by_you":false,"follows_you":false},"version":479,"publish_version":null,"role":null,"authored":false,"authors":[{"type":"individual","id":"f837ae4081a7fac2","login":"weickertjoe","name":"Joe Weickert","avatar_url":"https://avatars-next.observableusercontent.com/1169610d6fe2ae12a74ecea51bb14bd3a72502315e48c8bfdf27d06a3798b86e","tier":null,"bio":"","followed_by_you":false,"follows_you":false}],"members":[{"type":"individual","id":"f837ae4081a7fac2","login":"weickertjoe","name":"Joe Weickert","avatar_url":"https://avatars-next.observableusercontent.com/1169610d6fe2ae12a74ecea51bb14bd3a72502315e48c8bfdf27d06a3798b86e","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":0,"collection":null,"type":"notebook","body":{"titling":"auto","title":"Air Travel Routing","cells":[{"id":0,"value":"# Air Travel Routing","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":42,"value":"Plot.plot({\n projection: {type: projection, rotate: rotate},\n width,\n marks: [\n Plot.geo(land, { stroke: \"#000\"}),\n Plot.sphere(),\n Plot.dot(filteredVertices, {x: \"lng\", y: \"lat\", fill: 'blue'}),\n Plot.dot(ports, {x: \"lng\", y: \"lat\", title: \"name\", fill: 'red', tip: true})\n ]\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":152,"value":"viewof rotate = Inputs.form([\n Inputs.range([-180, 180], {step: 0.1, label: \"λ\"}),\n Inputs.range([-90, 90], {step: 0.1, label: \"φ\"}),\n])","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":214,"value":"countriesToBlock = []; //land.features.filter(f => [\"China\", \"Russia\", \"Brazil\"].includes(f.properties.name)).map(l => l.geometry)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":149,"value":"viewof projection = Inputs.select([\n \"orthographic\",\n \"equirectangular\"\n])","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":62,"value":"viewof pointDistance = Inputs.range([5, 50], {step: 1, label: \"Distance\", value: 5})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":10,"value":"land = topojson.feature(countries, countries.objects.countries)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":9,"value":"countries = FileAttachment(\"countries-50m.json\").json()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":95,"value":"$49","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":48,"value":"airTravelVertices = d3.range(-90, 95, pointDistance).flatMap(lat => d3.range(-180, 185, pointDistance).map(lng => {return {lat, lng, id: mutable nodeId++, neighbors: [], type: 'Waypoint'}}));","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":136,"value":"nestedVertices = {\n const nested = [];\n \n let currentLat = airTravelVertices[0].lat;\n let currentRow = [];\n for (let n = 0; n < airTravelVertices.length; n++) {\n if (airTravelVertices[n].lat !== currentLat) {\n nested.push(currentRow);\n currentRow = [airTravelVertices[n]];\n currentLat = airTravelVertices[n].lat\n } else {\n currentRow.push(airTravelVertices[n]);\n }\n }\n \n nested.push(currentRow);\n\n return nested;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":217,"value":"filteredVertices = airTravelVertices.filter(v => {\n return countriesToBlock.every(c => c.type === \"Polygon\" ? \n c.coordinates.every(polygon => isPointNotInPolygon(v, polygon)) \n : c.coordinates.every(c => c.every(polygon => isPointNotInPolygon(v, polygon))));\n});","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":383,"value":"$4a","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":406,"value":"ports = FileAttachment(\"AirPorts@1.json\").json()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":409,"value":"ports.forEach(p => {\n const nearestThree = airTravelVertices.sort((a, b) => haversineDistance([p.lng, p.lat], [a.lng, a.lat], true) - haversineDistance([p.lng, p.lat], [b.lng, b.lat], true)).slice(0, 3).map(v => { return {id: v.id, distance: Math.round((haversineDistance([p.lng, p.lat], [v.lng, v.lat], true) * 100)) / 100} });\n p.neighbors = nearestThree;\n});","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":325,"value":"function isPointNotInPolygon(v, polygon) {\n let intersectCount = 0;\n \n for (let i = 0, j = polygon.length - 1; i < polygon.length; j = i++)\n {\n if (((polygon[i][1] <= v.lat && v.lat < polygon[j][1]) || (polygon[j][1] <= v.lat && v.lat < polygon[i][1])) \n &&\n (v.lng < (polygon[j][0] - polygon[i][0]) * (v.lat - polygon[i][1]) / (polygon[j][1] - polygon[i][1]) + polygon[i][0])\n )\n {\n intersectCount++;\n }\n }\n\n return intersectCount % 2 === 0;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":373,"value":"function haversineDistance(coords1, coords2, isMiles) {\n function toRad(x) {\n return x * Math.PI / 180;\n }\n\n var lon1 = coords1[0];\n var lat1 = coords1[1];\n\n var lon2 = coords2[0];\n var lat2 = coords2[1];\n\n var R = 6371; // km\n\n var x1 = lat2 - lat1;\n var dLat = toRad(x1);\n var x2 = lon2 - lon1;\n var dLon = toRad(x2)\n var a = Math.sin(dLat / 2) * Math.sin(dLat / 2) +\n Math.cos(toRad(lat1)) * Math.cos(toRad(lat2)) *\n Math.sin(dLon / 2) * Math.sin(dLon / 2);\n var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a));\n var d = R * c;\n\n if(isMiles) d /= 1.60934;\n\n return d;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":81,"value":"mutable nodeId = 0","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":104,"value":"mutable edgeId = 0;","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":50,"value":"d3 = require(\"d3@5\");","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W4b","stdlib":"1"}} 24:null