2a:null 2b:["$","title",null,{"children":"Deadliest roads in each borough (draft 2) | Zhi Keng He | Observable"}] 31:T906,map = { let container = DOM.element("div", { style: `width:${width}px;height:${width / 3}px` }); yield container; let map = L.map(container, { scrollWheelZoom: false }).setView( [40.743, -73.93], 14 ); let osmLayer = L.tileLayer( "https://{s}.basemaps.cartocdn.com/light_all/{z}/{x}/{y}x2.png", { attribution: "© OpenStreetMap contributors", maxZoom: 20, minZoom: 11 } ).addTo(map); const segments = wktParser(selected_lion.object(0)["wkt"]); const layer = L.geoJSON(segments, { style: { color: "yellow", weight: 5, opacity: 0.3 } }).addTo(map); const ids = selected_collisions.map((i) => i["COLLISION_ID"]); const points = collisions_table .select([ "NUMBER OF PERSONS KILLED", "ON STREET NAME", "OFF STREET NAME", "CROSS STREET NAME", "COLLISION_ID", "datetime", "location" ]) .objects() .map((crash) => ({ type: "Feature", properties: { deaths: crash["NUMBER OF PERSONS KILLED"], on_street: ids.includes(crash["COLLISION_ID"]), datetime: crash["datetime"], address: `Location` }, geometry: { type: "Point", coordinates: crash.location } })); function onEachFeature(feature, layer) { // does this feature have a property named popupContent? if (feature.properties) { layer.bindPopup( `

${feature.properties.datetime.toLocaleString("en-US", { dateStyle: "full" })}

${feature.properties.address}

` ); } } const pointsLayer = L.geoJSON(points, { onEachFeature, pointToLayer: function (feature, latlng) { return L.circleMarker(latlng, { radius: feature.properties.on_street ? feature.properties.deaths * 4 : 2, fillColor: feature.properties.on_street ? "red" : "#e5d3c3", color: "#000", weight: 0, opacity: 0.5, fillOpacity: feature.properties.on_street ? 0.8 : 0.6 }); } }).addTo(map); map.fitBounds(layer.getBounds()); }32:T535,{ const dates_count = selected_collisions.reduce((dates, collision) => { const date = new Date(collision.datetime.setHours(0, 0, 0, 0)).setDate(1); //check if date exists let find = dates.find((i) => i.date === date); if (!find) { find = { date, count: 0, killed: 0, injured: 0 }; dates.push(find); } find.count += 1; find.killed += collision["NUMBER OF PERSONS KILLED"]; //find.injured += collision["NUMBER OF PERSONS INJURED"]; return dates; }, []); //breakdown array const dates_count_breakdown = ["killed"].reduce((rows, key) => { return [ ...rows, ...dates_count.map((i) => ({ date: new Date(i.date), total: i[key], type: key })) ]; }, []); const height = 180; const ksi_plot = Plot.plot({ width, height, y: { grid: true, label: "# of Deaths" }, x: { type: "band", tickFormat: (d, i) => { const date = new Date(d); return date.toLocaleString("en", { year: "numeric" }); }, label: null, interval: d3.utcYear }, marks: [ Plot.barY(dates_count_breakdown, { x: "date", y: "total", fill: "type" }), Plot.ruleY([0]) ] }); //${ksi_plot.legend( "color" )} return html`

${ksi_plot}`; }33:[{"id":"b3514b1be7730891969e08e1209563d76502db9cc5b67574eed39f59dbd454fb78dbfa8e32e66542169956d1331c3061f17c8a051b98de5ea5022dc59b43cd7a","href":"https://static.observableusercontent.com/files/b3514b1be7730891969e08e1209563d76502db9cc5b67574eed39f59dbd454fb78dbfa8e32e66542169956d1331c3061f17c8a051b98de5ea5022dc59b43cd7a","name":"crashes_killed.parquet","type":"application/octet-stream","lastModified":1661969646156.937,"size":56387},{"id":"4b7258d8bce8c3364e9758be0ebcb99c698eb2c6bfafe31dfb6e87c80f333ea005a484617b6f6421c05c2368a20b2513de3f50b12746b1a51e82886f8fea1d9f","href":"https://static.observableusercontent.com/files/4b7258d8bce8c3364e9758be0ebcb99c698eb2c6bfafe31dfb6e87c80f333ea005a484617b6f6421c05c2368a20b2513de3f50b12746b1a51e82886f8fea1d9f","name":"lion_wcrashes.parquet","type":"application/octet-stream","lastModified":1661963605141.127,"size":7992827}] 2e:{"id":"9a83d3eb39cce81a","slug":"deadliest-roads-in-each-borough-draft-2","title":"Deadliest roads in each borough (draft 2)","thumbnail":"53e80b37fdfc978c71276657dc49b6b17f64a97a23c97527dc41d23f7b13db61","access_level":"public","create_time":"$D2022-08-31T14:31:24.937Z","update_time":"$D2022-10-12T18:11:51.189Z","view_time":null,"first_public_time":"$D2022-08-31T18:22:36.995Z","fork":{"id":"22f909d7967a1c2b","slug":"deadliest-roads-in-each-borough","title":"Deadliest roads in each borough","owner":{"id":"16506c8d98fe0d4f","login":"zhik","name":"Zhi Keng He"}},"stars":0,"starred":false,"owner":{"type":"team","id":"16506c8d98fe0d4f","login":"zhik","name":"Zhi Keng He","avatar_url":"https://avatars-next.observableusercontent.com/1921c3cb8ffa82cde953aa4417dc8b601dbd05384629c1f8ecf0544be068a204","tier":"starter_2024"},"version":1002,"publish_version":null,"role":null,"authored":false,"authors":[{"type":"individual","id":"215e59083f0e53f8","login":"zhik","name":"Zhi Keng He","avatar_url":"https://avatars-next.observableusercontent.com/1921c3cb8ffa82cde953aa4417dc8b601dbd05384629c1f8ecf0544be068a204","tier":null}],"members":[{"type":"individual","id":"215e59083f0e53f8","login":"zhik","name":"Zhi Keng He","avatar_url":"https://avatars-next.observableusercontent.com/1921c3cb8ffa82cde953aa4417dc8b601dbd05384629c1f8ecf0544be068a204","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"subscribed":null,"forks":0,"collection":null,"type":"notebook","body":{"title":"Deadliest roads in each borough (draft 2)","cells":[{"id":0,"value":"# Deadliest roads in each borough (draft 2)","pinned":false,"hidden":false,"mode":"md"},{"id":527,"value":"viewof search = {\n const o = roads\n .params({ cutoff })\n .select(\n \"Street\",\n \"ModFaceCode\",\n \"crash_count\",\n \"killed_total\",\n \"length\",\n \"d_per_mile\",\n \"weighted_d_per_mile\"\n );\n\n const searchObjects = cutoff ? o.filter((d, $) => d.length >= $.cutoff) : o;\n return Inputs.search(searchObjects.objects());\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":795,"value":"viewof cutoff = Inputs.range([0, 5], { step: 0.5, value: 0 })","pinned":false,"hidden":false,"mode":"ojs"},{"id":155,"value":"viewof select = Inputs.table(search, {\n format: {\n crash_count: sparkbar(d3.max(search, (d) => d.crash_count)),\n killed_total: sparkbar(d3.max(search, (d) => d.killed_total))\n },\n header: {\n length: \"Length (miles)\",\n crash_count: \"Crashes since 2018\",\n killed_total: \"# of Deaths\",\n d_per_mile: \"Deaths per mile\",\n weighted_d_per_mile: \"Deaths per mile (adj)\"\n },\n multiple: false,\n rows: 10,\n columns: [\n \"Street\",\n \"crash_count\",\n \"killed_total\",\n \"length\",\n \"d_per_mile\",\n \"weighted_d_per_mile\"\n ],\n value: search[0]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":909,"value":"total_deaths_since_2018_geocoded = lion_table\n .rollup({\n sum: (d) => op.sum(d.killed_total)\n })\n .object(0)[\"sum\"]","pinned":false,"hidden":false,"mode":"ojs"},{"id":763,"value":"{\n const { crash_count, Street, killed_total } = selected_lion.object(0);\n function titleCase(str) {\n return str.toLowerCase().replace(/\\b\\w/g, (s) => s.toUpperCase());\n }\n const format = d3.format(\",\");\n\n return html`

Since 2018, ${format(\n crash_count\n )} crashes has occurred on ${titleCase(\n Street\n )}. Of those, there were ${\n selected_collisions.length\n } crashes resulted in ${killed_total} deaths.

`;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":282,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":443,"value":"$32","pinned":false,"hidden":false,"mode":"ojs"},{"id":776,"value":"","pinned":false,"hidden":false,"mode":"html"},{"id":531,"value":"roads = lion_table.orderby(aq.desc(\"killed_total\"))\n","pinned":false,"hidden":false,"mode":"ojs"},{"id":193,"value":"selected_lion = lion_table\n .params({\n boro: select ? select.LBoro : 0,\n ModFaceCode: select ? select.ModFaceCode : 0\n })\n .filter((d, $) => d[\"ModFaceCode\"] === $.ModFaceCode)","pinned":false,"hidden":false,"mode":"ojs"},{"id":327,"value":"selected_lion.objects()","pinned":false,"hidden":false,"mode":"ojs"},{"id":173,"value":"selected_collisions = {\n const prep = selected_lion.object(0).collision_ids ?? [];\n\n const ids = eval(prep);\n\n const crashes = collisions_table\n .params({\n ids\n })\n .filter((d, $) => op.includes($.ids, d[\"COLLISION_ID\"]))\n .objects();\n\n return crashes;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":157,"value":"### Appendix / Imports","pinned":false,"hidden":false,"mode":"md"},{"id":585,"value":"boro_lookup = [\n [\"Manhattan\", 1],\n [\"Bronx\", 2],\n [\"Brooklyn\", 3],\n [\"Queens\", 4],\n [\"Staten Island\", 5]\n]","pinned":false,"hidden":false,"mode":"ojs"},{"id":5,"value":"parquetWasm = import(\"https://cdn.jsdelivr.net/npm/parquet-wasm@0.1.1/web.js\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":7,"value":"readParquet = parquetWasm.readParquet","pinned":false,"hidden":false,"mode":"ojs"},{"id":9,"value":"async function getDataTable(file) {\n await parquetWasm.default();\n const arrayBuffer = await file.arrayBuffer();\n const arr = new Uint8Array(arrayBuffer);\n const arrowIPC = readParquet(arr);\n const dt = aq.fromArrow(arrowIPC);\n return dt;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":23,"value":"import { aq, op } from \"@uwdata/arquero\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":4,"value":"lion = FileAttachment(\"lion_wcrashes.parquet\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":52,"value":"lion_table = {\n const dt = await getDataTable(lion);\n\n aq.addFunction(\"Number\", Number);\n aq.addFunction(\"wktParser\", wktParser);\n aq.addFunction(\"turfLength\", turf.length);\n const dt2 = dt\n .filter((d) => d.killed_total)\n .filter((d) => d.Street !== \"UNNAMED STREET\")\n .derive({\n killed_total: (d) => Number(d.killed_total),\n crash_count: (d) => Number(d.crash_count),\n collision_ids: (d) => d.collision_ids,\n wkt: (d) => d.wkt,\n length: (d, $) => $.turfLength(wktParser(d.wkt), { units: \"miles\" })\n })\n .derive({\n d_per_mile: (d) => d.killed_total / d.length\n });\n\n const { min, max, avg } = dt2\n .rollup({\n min: (d) => op.min(d.crash_count),\n max: (d) => op.max(d.crash_count),\n avg: (d) => op.average(d.crash_count)\n })\n .object(0);\n\n return dt2.params({ min, max, avg }).derive({\n weighted_d_per_mile: (d, $) =>\n ((d.crash_count - $.min) / ($.max - $.min)) * d.d_per_mile\n });\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":108,"value":"collisions = FileAttachment(\"crashes_killed.parquet\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":116,"value":"collisions_table = {\n const dt = await getDataTable(collisions);\n aq.addFunction(\"Number\", Number);\n return dt\n .derive({\n location: (d) => {\n const [lat, lng] = op.split(\n op.replace(op.replace(d[\"LOCATION\"], \"(\", \"\"), \")\", \"\"),\n \",\"\n );\n return [+lng, +lat];\n },\n datetime: (d) => {\n const [month, day, year] = op.split(d[\"CRASH DATE\"], \"/\");\n const [hour, min] = op.split(d[\"CRASH TIME\"], \":\");\n\n return op.datetime(+year, +month - 1, +day, +hour, +min);\n },\n COLLISION_ID: (d) => Number(d.COLLISION_ID)\n })\n .select(aq.not([\"LOCATION\", \"CRASH DATE\", \"CRASH TIME\"]));\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":422,"value":"collisions_table.view(10)","pinned":false,"hidden":false,"mode":"ojs"},{"id":268,"value":"wktParser = require(\"https://bundle.run/wellknown@0.5.0\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":304,"value":"L = require(\"leaflet\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":544,"value":"function sparkbar(max) {\n return (x) => htl.html`

${x.toLocaleString(\"en\")}`;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":335,"value":"html``","pinned":false,"hidden":false,"mode":"ojs"},{"id":673,"value":"turf = require(\"@turf/turf\")","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W33","stdlib":"1","titling":"auto"}} 2f:"$2e:body"