31:[] 2f:{"title":"PM2.5 to AQI Conversion","cells":[{"id":0,"value":"# PM2.5 to AQI Conversion\n\nMy home weather station measures PM2.5 concentration in micrograms per cubic meter (μg/m³). But if I’m concerned about air quality, I want to know the Air Quality Index (AQI). AQI is a policy tool and as such has a complicated definition. So, here’s a tool to convert between PM2.5 and AQI (assuming that PM2.5 is the only pollutant you care about).","pinned":false,"hidden":false,"mode":"md"},{"id":299,"value":"Drag either slider below to choose the selected PM2.5 or AQI.","pinned":false,"hidden":false,"mode":"md"},{"id":204,"value":"viewof pm25 = Inputs.range([0, 500], {label: \"PM2.5 (μg/m³)\", value: 50, step: 0.1})","pinned":false,"hidden":false,"mode":"ojs"},{"id":206,"value":"bind(Inputs.range([0, 500], {label: \"AQI\", step: 1}), viewof pm25, {transform: pm25_aqi, invert: aqi_pm25})","pinned":false,"hidden":false,"mode":"ojs"},{"id":62,"value":"Plot.plot({\n grid: true,\n x: {\n label: \"PM2.5 →\"\n },\n y: {\n label: \"↑ AQI\"\n },\n color: {\n type: \"identity\"\n },\n marks: [\n Plot.rect(categories, {\n x1: 0,\n x2: d => aqi_pm25(d.max),\n y1: 0,\n y2: \"max\",\n stroke: \"color\",\n strokeWidth: 3,\n fill: \"color\",\n reverse: true\n }),\n Plot.line(d3.range(501), {\n x: d => d, \n y: pm25_aqi\n }),\n Plot.dot([[pm25, pm25_aqi(pm25)]]),\n Plot.text(categories, {\n x: 0,\n y: \"max\",\n text: \"name\",\n textAnchor: \"start\",\n stroke: \"white\",\n strokeWidth: 3,\n strokeLinejoin: \"round\",\n dx: 4,\n dy: 12\n }),\n Plot.text(categories, {\n x: 0,\n y: \"max\",\n text: \"name\",\n textAnchor: \"start\",\n dx: 4,\n dy: 12\n })\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":268,"value":"The cell below defines the six official [Air Quality Index (AQI) categories](https://www.airnow.gov/aqi/aqi-basics/): “Each category corresponds to a different level of health concern. Each category also has a specific color. The color makes it easy for people to quickly determine whether air quality is reaching unhealthy levels in their communities.”","pinned":false,"hidden":false,"mode":"md"},{"id":78,"value":"categories = [\n {max: 50, color: \"green\", name: \"Good\"},\n {max: 100, color: \"yellow\", name: \"Moderate\"},\n {max: 150, color: \"orange\", name: \"Unhealthy for sensitive groups\"},\n {max: 200, color: \"red\", name: \"Unhealthy\"},\n {max: 300, color: \"purple\", name: \"Very unhealthy\"},\n {max: 500, color: \"maroon\", name: \"Hazardous\"}\n]","pinned":true,"hidden":false,"mode":"ojs"},{"id":262,"value":"The `pm25_aqi` function converts from a PM2.5 concentration in micrograms per cubic meter to the corresponding AQI value (assuming that PM2.5 is the only contributor to AQI). It is the inverse of `aqi_pm25`.","pinned":false,"hidden":false,"mode":"md"},{"id":6,"value":"function pm25_aqi(pm25) {\n const c = Math.floor(10 * pm25) / 10;\n const a = c < 0 ? 0 // values below 0 are considered beyond AQI\n : c < 12.1 ? lerp( 0, 50, 0.0, 12.0, c)\n : c < 35.5 ? lerp( 51, 100, 12.1, 35.4, c)\n : c < 55.5 ? lerp(101, 150, 35.5, 55.4, c)\n : c < 150.5 ? lerp(151, 200, 55.5, 150.4, c)\n : c < 250.5 ? lerp(201, 300, 150.5, 250.4, c)\n : c < 350.5 ? lerp(301, 400, 250.5, 350.4, c)\n : c < 500.5 ? lerp(401, 500, 350.5, 500.4, c)\n : 500; // values above 500 are considered beyond AQI\n return Math.round(a);\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":258,"value":"The `aqi_pm25` function converts from an AQI value to the corresponding PM2.5 concentration in micrograms per cubic meter (assuming that PM2.5 is the only contributor to AQI). It is the inverse of `pm25_aqi`.","pinned":false,"hidden":false,"mode":"md"},{"id":91,"value":"function aqi_pm25(aqi) {\n const a = Math.round(aqi);\n const c = a < 0 ? 0 // values below 0 are considered beyond AQI\n : a <= 50 ? lerp( 0.0, 12.0, 0, 50, a)\n : a <= 100 ? lerp( 12.1, 35.4, 51, 100, a)\n : a <= 150 ? lerp( 35.5, 55.4, 101, 150, a)\n : a <= 200 ? lerp( 55.5, 150.4, 151, 200, a)\n : a <= 300 ? lerp(150.5, 250.4, 201, 300, a)\n : a <= 400 ? lerp(250.5, 350.4, 301, 400, a)\n : a <= 500 ? lerp(350.5, 500.4, 401, 500, a)\n : 500; // values above 500 are considered beyond AQI\n return Math.floor(10 * c) / 10;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":253,"value":"The `lerp` function is like d3.scaleLinear, redux.","pinned":false,"hidden":false,"mode":"md"},{"id":8,"value":"function lerp(ylo, yhi, xlo, xhi, x) {\n return ((x - xlo) / (xhi - xlo)) * (yhi - ylo) + ylo;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":247,"value":"The `bind` function is like [Inputs.bind](https://github.com/observablehq/inputs/blob/main/README.md#bind), except it applies a transform to convert between units, and only propagates on trusted events. This way when the user interacts with the target, it’ll propagate to the source, but the source won’t propagate back to the target. This transform needs to be invertible so that you can drag either range input to affect the other.","pinned":false,"hidden":false,"mode":"md"},{"id":211,"value":"function bind(target, source, {\n invalidation = Inputs.disposal(target),\n transform = d => d,\n invert = d => d\n} = {}) {\n const onsource = (event) => {\n if (!event.isTrusted) return;\n target.value = transform(source.value);\n };\n const ontarget = (event) => {\n if (!event.isTrusted) return;\n source.value = invert(target.value);\n source.dispatchEvent(new Event(\"input\", {bubbles: true}));\n };\n target.value = transform(source.value);\n target.addEventListener(\"input\", ontarget);\n source.addEventListener(\"input\", onsource);\n invalidation.then(() => source.removeEventListener(\"input\", onsource));\n return target;\n}","pinned":true,"hidden":false,"mode":"ojs"}],"files":"$W31","stdlib":"1","titling":"auto"}