31:T632d,map3 = { lmc; lspin; let cluster_sizes = []; let container = DOM.element("div", { class: "mainmap", style: `width:${width}px;height:${width * 0.7}px;` }); yield container; var visited = {}; let map = L.map(container, { tap: false }).setView([40, -75.12], 11); esriLeaflet .tiledMapLayer({ url: "https://tiles.arcgis.com/tiles/fLeGjb7u4uXqeF9q/arcgis/rest/services/CityBasemap/MapServer?f=jsapi&cacheKey=8c0da44bee92edf8", maxZoom: 18 }) .addTo(map); esriLeaflet .tiledMapLayer({ url: "https://tiles.arcgis.com/tiles/fLeGjb7u4uXqeF9q/arcgis/rest/services/CityBasemap_Labels/MapServer?f=jsapi&cacheKey=b5e156f30610659b", maxZoom: 18 }) .addTo(map); let target_scale = "# Visits weighted by % Unvaccinated"; const blocksTooltipContent = l => ``; let blocks = L.geoJSON(aq_tracts_ranked.objects()); let removedblocks = L.layerGroup([]); let removeNoPopBlocks = false; const vsites = L.layerGroup( vaccine_sites_aq.objects().map(v => { return L.marker([v.lat, v.long], { title: v["Site_Name"], vsite: v }); }) ); //L.geoJSON(enrolled_Vaccination_Sites_PublicView); let highlightblocks; //= L.geoJSON(in_filterCT); let selectedblocks_info = []; let selectedblocks = L.geoJson([]); let tractsInfo_places; let isFocused = false; let cluster_num_mean = 0; let selectedplaces = L.layerGroup([]); let hoverclusters; let selectedclusters; let highlightedMarker; //let keepHighlightedMarker = false; let mean_places_num = aq_places .rollup({ mean: `d => op.mean(d['total_visits']) ` }) .get("mean", 0); const placesIcon = place => { let color = clusterColorFunction( logisticSize( 1, 0, place["geoid_visits"], mean_places_num, 2.0 / mean_places_num ) ); return L.divIcon({ color: color, className: "default", iconAnchor: [8, 24], tooltipAnchor: [0, -22], popupAnchor: [0, -24], html: `${svg_marker}

` }); }; const placesMarker = p => { let location = p["Location"] ? "Location" : "location_name"; return L.marker([p.latitude, p.longitude], { title: p[location], place: p, icon: placesIcon(p) }).bindTooltip( p[location] + "
" + Math.round(p.geoid_visits * 100) / 100 + " average visits", { direction: "top" } ); }; const rerfreshBlocksColor = () => { blocks.eachLayer(l => { l.setTooltipContent(blocksTooltipContent(l)); l.setStyle( styleFunction( l.feature[target_scale], aq_tracts_collect.get(target_scale + "_mean", 0), 3.0 / aq_tracts_collect.get(target_scale + "_mean", 0), false ) ); }); }; const refreshClusterColor = () => { let num_arr = []; map.eachLayer(layer => { if (layer.getChildCount) { let info = layer .getAllChildMarkers() .map(a => { let a_num = a.options.place["geoid_visits"]; a_num = a_num || 0; return a_num; }) .map(Number); let num = info.reduce((a, b) => a + b, 0); num_arr.push(num); } }); if (num_arr.length != 0) cluster_num_mean = num_arr.reduce((a, b) => a + b) / num_arr.length; if (hoverclusters) hoverclusters.refreshClusters(); try { selectedclusters.refreshClusters(); } catch (TypeError) {} }; const clusterIconCreateFunction = size => { return cluster => { let info = cluster .getAllChildMarkers() .map(a => { let a_num = a.options.place["geoid_visits"]; a_num = a_num || 0; return a_num; }) .map(Number); let raw_num = info.reduce((a, b) => a + b, 0); let short_num = raw_num > 1000 ? Math.round(raw_num / 1000) + "k" : Math.round(raw_num); let clusterColor = getClusterColor( raw_num, cluster_num_mean, 1.0 / cluster_num_mean ); return L.divIcon({ html: `

${short_num}

`, className: "cluster", iconSize: L.point(size, size) }); }; }; map.createPane("paneA").style.zindex = 1100; map.createPane("paneB").style.zIndex = 200; let clusterOptions = { showCoverageOnHover: false, zoomToBoundsOnClick: true, maxClusterRadius: 40, animate: true, disableClusteringAtZoom: 18, spiderfyOnMaxZoom: true, chunkedLoading: true, singleMarkerMode: false, iconCreateFunction: clusterIconCreateFunction(25), clusterPane: "paneB" }; let legendBar; const createLegendBar = data => { L.Control.LegendBar = L.Control.extend({ onAdd: map => { let div = L.DomUtil.create("div"); let divbody = L.DomUtil.create("div"); let divtitle = L.DomUtil.create("h1"); div.className = "legend-bar"; divtitle.className = "legend-title"; divbody.className = "legend-body"; let cscale = L.DomUtil.create("div"); cscale.className = "legend-scale"; let topscale = L.DomUtil.create("h1"); let bottomscale = L.DomUtil.create("h1"); topscale.innerHTML = formatNum( aq_tracts_collect.get(target_scale + "_max", 0), 0 ); bottomscale.innerHTML = formatNum( aq_tracts_collect.get(target_scale + "_min", 0), 0 ); let rampcscale = L.DomUtil.create("div"); rampcscale.appendChild( vramp(colorFunction, undefined, undefined, true) ); cscale.appendChild(topscale); cscale.appendChild(rampcscale); cscale.appendChild(bottomscale); divbody.appendChild(cscale); div.appendChild(divtitle); div.appendChild(divbody); return div; }, onRemove: map => {} }); L.control.legendBar = opts => { return new L.Control.LegendBar(opts); }; legendBar = L.control.legendBar({ position: "topleft" }); legendBar.addTo(map); }; let tractsInfo; let searchBar; const createSearchBar = data => { L.Control.SearchBar = L.Control.extend({ onAdd: map => { let div = L.DomUtil.create("div"); div.className = "search-box"; div.id = "searchbox"; div.innerHTML = ``; L.DomEvent.on(div, "keyup", () => { if (highlightblocks) highlightblocks.remove(); let in_searchbox_text = div.firstChild.value; let highlightblocks_table = aq_tracts_ranked .params({ in_searchbox_text }) .filter(d => op.indexof(d["CT Number"], in_searchbox_text) !== -1) .objects(); highlightblocks = L.geoJSON(highlightblocks_table); in_searchbox_text.length > 1 ? highlightblocks.addTo(map) : highlightblocks.remove(); }); L.DomEvent.disableScrollPropagation(div); L.DomEvent.disableClickPropagation(div); return div; }, onRemove: map => {} }); L.control.searchBar = opts => { return new L.Control.SearchBar(opts); }; searchBar = L.control.searchBar({ position: "topright" }); searchBar.addTo(map); }; let render_table; const updateRenderTable = () => { if (tractsInfo_places) { render_table = tractsInfo_places.groupby("safegraph_place_id").rollup({ Location: r => op.any(r["Location"]), Category: r => op.any(r["Category"]), Visits: r => op.sum(r["Visits"]), "Total Visits": r => op.sum(r["Total Visits"]), "Median Dwell": r => op.mean(r["Median Dwell"]), origin_geoid: r => op.any(r["origin_geoid"]), latitude: r => op.any(r.latitude), longitude: r => op.any(r.longitude), geoid_visits: r => op.sum(r["Visits"]) }); } }; const createTractsInfo = data => { let show_table = render_table .derive({ "Lookup ID": d => d.safegraph_place_id }) .orderby(aq.desc("Visits")) .select( "Location", "Category", "Visits", "Median Dwell", "Total Visits", "origin_geoid", "Lookup ID" ); L.Control.TractsInfo = L.Control.extend({ onAdd: map => { let div = L.DomUtil.create("div"); let titleBox = L.DomUtil.create("div"); let detailsBox = L.DomUtil.create("div"); let tableBox = L.DomUtil.create("div"); tableBox.className = "table-box"; tableBox.id = "tableBox"; div.className = "tracts-info"; div.id = "infobox"; titleBox.innerHTML = `

Census Tract ${tractsInfo_places .select("origin_geoid") .dedupe() .columnArray("origin_geoid")}

`; detailsBox.innerHTML = `

Estimated Outgoing: ${Math.round( tractsInfo_places .rollup({ sum: d => op.sum(d["Visits"]) }) .get("sum", 0) )} people

`; tableBox.className = "places"; tableBox.innerHTML = show_table .toHTML({ limit: 100, maxdigits: 2 }) .replaceAll("", ``) .replaceAll("og", ` og`); let trs = Array.prototype.slice.call( tableBox.getElementsByTagName("table")[0].rows ); trs.shift(); trs.forEach(a => { let split = a.lastChild.innerHTML.split("_||_"); const getDisplayMarkers = () => { let target_marker, hover_marker; selectedclusters.eachLayer(l => { if (l.options.place.safegraph_place_id == split[0]) { target_marker = l; hover_marker = L.marker(l._latlng, { zIndexOffset: 1000, icon: L.divIcon({ iconAnchor: [12, 12], className: "default", html: `

` }) }); } }); return [target_marker, hover_marker]; }; L.DomEvent.on(a, "mouseover", e => { highlightedMarker = getDisplayMarkers()[1]; highlightedMarker.addTo(map); }); L.DomEvent.on(a, "mouseout", e => { if (highlightedMarker) highlightedMarker.remove(); }); L.DomEvent.on(a, "click", e => { map.flyTo(getDisplayMarkers()[0]._latlng, 18); }); }); div.appendChild(titleBox); div.appendChild(detailsBox); div.appendChild(tableBox); L.DomEvent.disableScrollPropagation(div); L.DomEvent.disableClickPropagation(div); return div; }, onRemove: map => {} }); L.control.tractsInfo = opts => { return new L.Control.TractsInfo(opts); }; if (tractsInfo) tractsInfo.remove(); tractsInfo = L.control.tractsInfo({ position: "bottomright" }); tractsInfo.addTo(map); }; let selectedClusterOptions = Object.assign({}, clusterOptions); selectedClusterOptions.clusterPane = "paneA"; selectedclusters = new L.markerClusterGroup(selectedClusterOptions); const removeSelectedMenu = () => { if (selectedclusters) selectedclusters.remove(); if (selectedblocks) selectedblocks.remove(); if (tractsInfo) tractsInfo.remove(); selectedclusters.clearLayers(); selectedblocks_info = []; selectedblocks = L.geoJson([]); selectedplaces = L.layerGroup([]); tractsInfo_places = undefined; }; const addSelectedBlocks = arr => { map.spin(true); setTimeout(() => { if (selectedblocks) selectedblocks.remove(); selectedblocks_info = selectedblocks_info.concat(arr); selectedblocks = L.geoJSON(selectedblocks_info, { weight: 5, opacity: 1, fillColor: "black", pane: "paneA", onEachFeature: (f, layer) => { layer.on("click", e => { map.spin(true); setTimeout(() => { const index = selectedblocks_info.indexOf(layer.feature); if (index > -1) { selectedblocks_info.splice(index, 1); } layer.remove(); tractsInfo_places = tractsInfo_places .params({ layer }) .filter(d => d.origin_geoid != layer.feature.origin_geoid); selectedplaces.eachLayer(selectedlayer => { if ( selectedlayer.options.place.origin_geoid == layer.feature.origin_geoid ) { selectedplaces.removeLayer(selectedlayer); } }); updateRenderTable(); selectedclusters.clearLayers(); selectedplaces = L.layerGroup([]); render_table.objects().forEach(p => { selectedplaces.addLayer(placesMarker(p)); }); tractsInfo.remove(); if (tractsInfo_places.numRows() > 0) createTractsInfo(); selectedclusters.addLayers(selectedplaces); selectedclusters.refreshClusters(); if (selectedclusters.getLayers().length < 1) isFocused = false; map.spin(false); }, 1); }); } }); selectedclusters.clearLayers(); selectedplaces = L.layerGroup([]); selectedclusters.addLayers(selectedplaces); //Repeat, but important for speed. selectedclusters.addTo(map); let arr_places = []; arr.forEach(a => { arr_places = arr_places.concat(a.safegraph_place_info); }); let raw_table = aq.from(arr_places).derive({ Location: d => d.location_name, Category: d => d.category, Visits: d => Math.round(d.geoid_visits), "Median Dwell": d => Math.round(d.median_dwell), "Total Visits": d => Math.round(d.total_visits) }); if (tractsInfo_places) { tractsInfo_places = tractsInfo_places.concat(raw_table); } else { tractsInfo_places = raw_table; } updateRenderTable(); render_table.objects().forEach(p => { selectedplaces.addLayer(placesMarker(p)); }); selectedclusters.addLayers(selectedplaces); refreshClusterColor(); selectedclusters.on("animationend", e => { refreshClusterColor(); }); if (hoverclusters) hoverclusters.remove(); selectedblocks.addTo(map); createTractsInfo(); map.spin(false); }, 1); }; let highlightStyle = { weight: 1, opacity: 0.5, fillOpacity: 0.7, color: "purple", fillColor: "purple" }; let optionsInfo; const createOptionsInfo = data => { L.Control.OptionsInfo = L.Control.extend({ onAdd: map => { let div = L.DomUtil.create("div"); let shownSection = L.DomUtil.create("div"); let hiddenSection = L.DomUtil.create("div"); let hiddenSectionWrapper = L.DomUtil.create("div"); // shownSection Elements, shownSection.className = "upperSection-area"; // Delesect all currently selected census tracts button let deselectCTbtn = L.DomUtil.create("div"); //deselectCTbtn.className = "toggle-btn"; deselectCTbtn.innerHTML = `Deselect All Census Tracts`; L.DomEvent.on(deselectCTbtn, "pointerdown", () => { removeSelectedMenu(); }); shownSection.appendChild(deselectCTbtn); // hiddenSection Elements hiddenSectionWrapper.className = "lowerSection-area"; hiddenSection.className = "lowerSection-area-inside"; hiddenSectionWrapper.style.maxHeight = "0px"; let hideBtnBar = L.DomUtil.create('div'); hideBtnBar.className = 'hide-btn-bar'; // Hide/Show Census Tract Overlays let hideCTbtn = L.DomUtil.create("div"); hideCTbtn.className = "lowersec-btn"; hideCTbtn.innerHTML = "Hide Census Tracts"; L.DomEvent.on(hideCTbtn, "pointerdown", () => { if (map.hasLayer(blocks)) { hideCTbtn.innerHTML = "Show Census Tracts"; blocks.remove(); } else { hideCTbtn.innerHTML = "Hide Census Tracts"; blocks.addTo(map); } }); hideBtnBar.appendChild(hideCTbtn); // Hide/Show Vaccine Site Markers let hideVSites = L.DomUtil.create("div"); hideVSites.className = "lowersec-btn"; hideVSites.innerHTML = "Hide Vaxx Sites"; L.DomEvent.on(hideVSites, "pointerdown", () => { if (map.hasLayer(vsites)) { hideVSites.innerHTML = "Show Vaxx Sites"; vsites.remove(); } else { hideVSites.innerHTML = "Hide Vaxx Sites"; vsites.addTo(map); } }); hideBtnBar.appendChild(hideVSites); hiddenSection.appendChild(hideBtnBar); let selectColorBar = L.DomUtil.create('div'); selectColorBar.className = 'select-color-bar'; //bookmark // Select CT Color Scale let selectCTColor = L.DomUtil.create("div"); selectCTColor.className = "scale-select"; let selectCTColorLbl = L.DomUtil.create('div'); selectCTColorLbl.innerHTML = `CT Color: `; let selectCTColorOp = L.DomUtil.create('div'); selectCTColorOp.innerHTML = ` `; L.DomEvent.on(selectCTColorOp, "change", () => { target_scale = selectCTColor.lastChild.value; rerfreshBlocksColor(); legendBar.remove(); createLegendBar(); //refreshClusterColor(); //refreshMenus(); }); selectCTColor.appendChild(selectCTColorLbl); selectCTColor.appendChild(selectCTColorOp); selectColorBar.appendChild(selectCTColor); // Select POI Color Scale let selectPOIColor = L.DomUtil.create('div'); selectPOIColor.className = 'scale-select'; let selectPOIColorLbl = L.DomUtil.create('div'); selectPOIColorLbl.innerHTML = `POI Color: `; let selectPOIColorOp = L.DomUtil.create('div'); selectPOIColorOp.innerHTML = ` `; L.DomEvent.on(selectPOIColorOp, "change", () => { // target_scale = selectCTColor.lastChild.value; // rerfreshBlocksColor(); // legendBar.remove(); // createLegendBar(); //refreshClusterColor(); //refreshMenus(); }); selectPOIColor.appendChild(selectPOIColorLbl); selectPOIColor.appendChild(selectPOIColorOp); selectColorBar.appendChild(selectPOIColor); hiddenSection.appendChild(selectColorBar); // Selection Range for bottom-%-multi-census tract let rangeBox = L.DomUtil.create("div"); let confirmBox = L.DomUtil.create("div"); rangeBox.className = "range-slider"; rangeBox.innerHTML = ``; let rangeVal = rangeBox.firstChild.value; confirmBox.className = "toggle-btn"; confirmBox.innerHTML = `Select bottom ${rangeVal}% Tracts`; L.DomEvent.on(rangeBox, "input", () => { removeSelectedMenu(); rangeVal = rangeBox.firstChild.value; confirmBox.innerHTML = `Select bottom ${rangeVal}% Tracts`; }); L.DomEvent.on(confirmBox, "pointerdown", () => { removeSelectedMenu(); let totalIndex = aq_tracts_ranked.numRows(); let criticalIndex = Math.round((totalIndex * rangeVal) / 100); let blocksLayersToAdd = []; blocks.eachLayer(l => { if (l.feature[target_scale + "_rank"] <= criticalIndex) { blocksLayersToAdd.push(l.feature); } }); addSelectedBlocks(blocksLayersToAdd); refreshClusterColor(); }); hiddenSection.appendChild(rangeBox); hiddenSection.appendChild(confirmBox); // Expand button for hidden section let expandHiddenSectionBtn = L.DomUtil.create("div"); let isHiddenSectionExpanded = false; expandHiddenSectionBtn.className = "toggle-btn upperSection-area expand-section-btn"; expandHiddenSectionBtn.innerHTML = "--Expand--"; L.DomEvent.on(expandHiddenSectionBtn, "pointerdown", () => { if (isHiddenSectionExpanded) { hiddenSectionWrapper.style.maxHeight = "0px"; isHiddenSectionExpanded = false; } else { hiddenSectionWrapper.style.maxHeight = "1000px"; isHiddenSectionExpanded = true; } }); hiddenSectionWrapper.appendChild(hiddenSection); div.appendChild(shownSection); div.appendChild(hiddenSectionWrapper); div.appendChild(expandHiddenSectionBtn); L.DomEvent.disableScrollPropagation(div); L.DomEvent.disableClickPropagation(div); return div; }, onRemove: map => {} }); L.control.optionsInfo = opts => { return new L.Control.OptionsInfo(opts); }; optionsInfo = L.control.optionsInfo({ position: "topright" }); optionsInfo.addTo(map); }; //let focusBlock; hoverclusters = new L.markerClusterGroup(clusterOptions); blocks.eachLayer(l => { let target_table = l.feature.safegraph_place_info.map(p => { p.origin_geoid = l.feature.origin_geoid; return placesMarker(p); }); // console.log(target_table); let places = L.layerGroup(target_table, { title: l.feature.origin_geoid }); l.bindTooltip(blocksTooltipContent(l), { sticky: true, direction: "top" }); l.feature.places_info = places; l.feature.places_info.from = l.feature.properties.NAME10; l.setStyle( styleFunction( l.feature[target_scale], aq_tracts_collect.get(target_scale + "_mean", 0), 3.0 / aq_tracts_collect.get(target_scale + "_mean", 0), false ) ); l.on("mouseover", e => { if (!isFocused) { blocks.eachLayer(l => { let hoverstyle = l.options; hoverstyle["fillOpacity"] = 0.2; l.setStyle(hoverstyle); }); e.target.setStyle(highlightStyle); cluster_sizes = []; hoverclusters.clearLayers(); hoverclusters.addLayers(places); hoverclusters.addTo(map); refreshClusterColor(); hoverclusters.on("animationend", e => { refreshClusterColor(); }); hoverclusters.bringToBack(); blocks.bringToFront(); } }); l.on("mouseout", e => { if (!isFocused) { rerfreshBlocksColor(); l.setStyle(l.feature["Style"]); hoverclusters.remove(); blocks.bringToFront(); } }); l.on("click", e => { isFocused = true; //e.target.setStyle(l.feature["Style"]); rerfreshBlocksColor(); addSelectedBlocks([l.feature]); }); }); vsites.eachLayer(l => { l.options.zIndexOffset = -500; l.bindTooltip(html``, { direction: "top" }); l.options.icon = L.divIcon({ iconAnchor: [10, 10], iconSize: 20, html: `

`, className: "default" }); }); createLegendBar(); createSearchBar(map); createOptionsInfo(); vsites.addTo(map); blocks.addTo(map); }32:Te84,html``33:T537,aq_tracts = aq_tracts_unroll .derive({ safegraph_place_info: d => { return { safegraph_place_id: d.safegraph_place_id, location_name: d.location_name, latitude: d.latitude, longitude: d.longitude, longitude: d.longitude, area: d.area_square_feet, category: d.top_category, geoid_visits: d.geoid_visits, median_dwell: d.median_dwell, total_visits: d.total_visits, total_sq_ft: d.place_total_sq_ft, p_vaccinated: d.percent_vaccinated_visits, p_nvaccinated: d.percent_not_vaccinated_visits, vaccinated: d.total_vaccinated_visits }; } }) .groupby("origin_geoid") .rollup({ type: op.any("type"), properties: op.any("properties"), geometry: op.any("geometry"), safegraph_place_info: op.array_agg("safegraph_place_info"), total_visits: op.sum("geoid_visits"), total_weighted_visits: op.sum("weighted_visits") }) .join(demographic_data_alt, ["origin_geoid", "GEOID"]) .join( clean_vaccination_data .rename({ census_tract_geoid: "origin_geoid" }) .select("origin_geoid", "vaccine_percent") ) .join(processed_under_18_table, "origin_geoid") .derive({ visits_vaccinated: d => d.total_visits * d.vaccine_percent }) .select(aq.not("GEOID"))34:Tbd7,demographic_data_alt = demog_data .derive({ female: (d) => d[" female"], male: (d) => d[" male"], age_info: (d) => { return { "Under 20": op.parse_float(d[" <20"]), "20 to 34": op.parse_float(d[" 20-34"]), "35 to 54": op.parse_float(d[" 35-54"]), "55 to 75": op.parse_float(d[" 55-74"]), "Over 75": op.parse_float(d[" 75+"]) }; }, ethnicity_info: (d) => { return { "American Indian": op.parse_float( d[" American Indian and Alaska Native alone"] ), Asian: op.parse_float(d[" Asian alone"]), Black: op.parse_float(d[" Black or African American alone"]), White: op.parse_float(d[" White alone"]), "Pacific Islander": op.parse_float( d[" Native Hawaiian and Other Pacific Islander alone"] ), Other: op.parse_float(d[" Some other race alone"]), "Two or More": op.parse_float(d[" Two or more races:"]) }; }, income_info: (d) => { return { "Under $10k": op.parse_float(d["Less than $10,000"]), "$10k to $15k": op.parse_float(d["$10,000 to $14,999"]), "$15k to $20k": op.parse_float(d["$15,000 to $19,999"]), "$20k to $25k": op.parse_float(d["$20,000 to $24,999"]), "$25k to $30k": op.parse_float(d["$25,000 to $29,999"]), "$30k to $35k": op.parse_float(d["$30,000 to $34,999"]), "$35k to $40k": op.parse_float(d["$35,000 to $39,999"]), "$40k to $45k": op.parse_float(d["$40,000 to $44,999"]), "$45k to $50k": op.parse_float(d["$45,000 to $49,999"]), "$50k to $60k": op.parse_float(d["$50,000 to $59,999"]), "$60k to $75k": op.parse_float(d["$60,000 to $74,999"]), "$75k to $100k": op.parse_float(d["$75,000 to $99,999"]), "$100k to $130k": op.parse_float(d["$100,000 to $124,999"]), "$130k to $150k": op.parse_float(d["$125,000 to $149,999"]), "$150k to $200k": op.parse_float(d["$150,000 to $199,999"]), "$200k or more": op.parse_float(d["$200,000 or more"]) }; }, // poverty_info: (d) => { // return { // Above: op.parse_float(d[" abovePovertyLevel"]), // Below: op.parse_float(d[" belewPovertyLevel"]) // }; // }, family_info: (d) => { return { "Have Child": op.parse_float(d[" own child under18"]), "No Child": op.parse_float(d[" no child under18"]) }; } // transportation_info: (d) => { // return { // "Car/Truck": op.parse_float(d[" Car truck or van"]), // "Public Transportation": op.parse_float(d[" Public transportation"]), // Walking: op.parse_float(d[" walked"]), // "Taxi/Bike/Other": op.parse_float( // d[" Taxicab motorcycle bicycle or other means"] // ) // }; // } }) .select( "GEOID", "female", "male", "age_info", "ethnicity_info", "income_info", // "poverty_info", "family_info" // "transportation_info" )35:[{"id":"4b7abc566bf2e3804be16a86912c3861b47cd81941ecccef7f711bd82e271c3e71246f0d12874432e29bae0c9ecd877541a1392710648190ff3365668f8ed036","href":"https://static.observableusercontent.com/files/4b7abc566bf2e3804be16a86912c3861b47cd81941ecccef7f711bd82e271c3e71246f0d12874432e29bae0c9ecd877541a1392710648190ff3365668f8ed036","name":"Shot_0.png","type":"image/png","lastModified":1626487931326.223,"size":17372},{"id":"c24f236f8d41b7ba6ec6c5397480eae55a086bd0235782ad8a771e3528ef5efa9b90380362d85e2da66ec8ae0de4d5fc75c86f888af500229956b7135bbf4c87","href":"https://static.observableusercontent.com/files/c24f236f8d41b7ba6ec6c5397480eae55a086bd0235782ad8a771e3528ef5efa9b90380362d85e2da66ec8ae0de4d5fc75c86f888af500229956b7135bbf4c87","name":"cluster-yellow.png","type":"image/png","lastModified":1624030021132.619,"size":7407},{"id":"1e73259a9ae7c756f70eae018e5de940b237eee54b0fbc2b2f3524349402144832061f7f2d41f7f8f6af9b72514fc11c4f8192eb72de07b16f161c8933d02222","href":"https://static.observableusercontent.com/files/1e73259a9ae7c756f70eae018e5de940b237eee54b0fbc2b2f3524349402144832061f7f2d41f7f8f6af9b72514fc11c4f8192eb72de07b16f161c8933d02222","name":"cluster-green.png","type":"image/png","lastModified":1624030021132.619,"size":7387},{"id":"f48a4b4b32e7de96136e48e3a2e49c993571edcbed283800aad2a3dd701bced7bfc6e299acee90417352fd7a68cda5ab7ed802952bea05f851351b384c414d22","href":"https://static.observableusercontent.com/files/f48a4b4b32e7de96136e48e3a2e49c993571edcbed283800aad2a3dd701bced7bfc6e299acee90417352fd7a68cda5ab7ed802952bea05f851351b384c414d22","name":"cluster-red.png","type":"image/png","lastModified":1624030021132.619,"size":7258},{"id":"c2c399b14abb2ecafcb73c7d73dd785122fff86e6ddcbcf4994cd09c42c04eca32cc8fdc300abcc7f1384345fae618fc0efb1bf64a6449238c4793c06cb36751","href":"https://static.observableusercontent.com/files/c2c399b14abb2ecafcb73c7d73dd785122fff86e6ddcbcf4994cd09c42c04eca32cc8fdc300abcc7f1384345fae618fc0efb1bf64a6449238c4793c06cb36751","name":"cluster-orange@2.png","type":"image/png","lastModified":1624030021132.619,"size":7981},{"id":"56aaee97d3327d111d242489d73959ec785007707f1a383bd3d595806405061364640044f5b9f2730b8f8f3b4da425010217e0d82b950982f986dea599494051","href":"https://static.observableusercontent.com/files/56aaee97d3327d111d242489d73959ec785007707f1a383bd3d595806405061364640044f5b9f2730b8f8f3b4da425010217e0d82b950982f986dea599494051","name":"cluster-teal@2.png","type":"image/png","lastModified":1624030021132.619,"size":7782},{"id":"81e997672a1505fcbf98e5144904cb9191d0fe2775b4bf6bea560d3c02d12325d3726249d90b42c45869a042dcef089285e03e52148c3cdacfdfb364eb3fd412","href":"https://static.observableusercontent.com/files/81e997672a1505fcbf98e5144904cb9191d0fe2775b4bf6bea560d3c02d12325d3726249d90b42c45869a042dcef089285e03e52148c3cdacfdfb364eb3fd412","name":"new-marker.png","type":"image/png","lastModified":1624030021132.619,"size":3719},{"id":"015c6d5badf2e0bb9fbc96715e4369326ba7a3f8b65669560601bcd96afe56a5a3072cb03e580200dc35d0d72278ff7c1463ce31954f9034ed133cf47ec37849","href":"https://static.observableusercontent.com/files/015c6d5badf2e0bb9fbc96715e4369326ba7a3f8b65669560601bcd96afe56a5a3072cb03e580200dc35d0d72278ff7c1463ce31954f9034ed133cf47ec37849","name":"new-marker-highlighted.png","type":"image/png","lastModified":1624030021132.619,"size":3822},{"id":"34145651f31d1414a80b4238c0fcfb44fc32e06bccc2e6d10e204c042f3466ad0a225a638310b59508af60f92b763d473ad5edfc91bf08ff0661c64bd062d69a","href":"https://static.observableusercontent.com/files/34145651f31d1414a80b4238c0fcfb44fc32e06bccc2e6d10e204c042f3466ad0a225a638310b59508af60f92b763d473ad5edfc91bf08ff0661c64bd062d69a","name":"color-marker-orange.png","type":"image/png","lastModified":1624030021132.619,"size":3091},{"id":"33ebf2c62cd6e278f44fd811b20815eb54d85087ed4a4c3ae5e1b2e77efc36f459ae61b791a46facdf597461c11187401e399b748bcc6959b2bd5b3f32d63c6c","href":"https://static.observableusercontent.com/files/33ebf2c62cd6e278f44fd811b20815eb54d85087ed4a4c3ae5e1b2e77efc36f459ae61b791a46facdf597461c11187401e399b748bcc6959b2bd5b3f32d63c6c","name":"color-marker-teal.png","type":"image/png","lastModified":1624030021132.619,"size":3157},{"id":"afb21452bf07bc05f3d352d60d6e8ed73a6cdb4ebb7fc578383ce9930a3057adc3914522894c880aa931e49c3d0ba86d829df5b7730444c22f0854c59cc63f93","href":"https://static.observableusercontent.com/files/afb21452bf07bc05f3d352d60d6e8ed73a6cdb4ebb7fc578383ce9930a3057adc3914522894c880aa931e49c3d0ba86d829df5b7730444c22f0854c59cc63f93","name":"color-marker-green.png","type":"image/png","lastModified":1624030021132.619,"size":3261},{"id":"b4d384c3660e03be2605fdac7a6f5d83e6242cc93501b44faa967d690af001394355ad6a7e5f874c8b940dc9a9ccd27a812801a752d0493443e92f22383aa14f","href":"https://static.observableusercontent.com/files/b4d384c3660e03be2605fdac7a6f5d83e6242cc93501b44faa967d690af001394355ad6a7e5f874c8b940dc9a9ccd27a812801a752d0493443e92f22383aa14f","name":"color-marker-red.png","type":"image/png","lastModified":1624030021132.619,"size":3323},{"id":"49d561492937e3cdb92998ac0dfefa393375dfb5e3122dd8429f0ad8cc776cb71024fcdd5e26e39a1503fb48abe01dfab530d5c8181c6ccddb8799a60c5adde8","href":"https://static.observableusercontent.com/files/49d561492937e3cdb92998ac0dfefa393375dfb5e3122dd8429f0ad8cc776cb71024fcdd5e26e39a1503fb48abe01dfab530d5c8181c6ccddb8799a60c5adde8","name":"color-marker-yellow.png","type":"image/png","lastModified":1624030021132.619,"size":3585},{"id":"db7e70eb2941ea9f69fd9fea947eadb638b9f622a7cbc0c399811f897f36c2088c8b1d45979c8402d181101a0c7b4b73f266d302dbdd424005662dc2a74611dc","href":"https://static.observableusercontent.com/files/db7e70eb2941ea9f69fd9fea947eadb638b9f622a7cbc0c399811f897f36c2088c8b1d45979c8402d181101a0c7b4b73f266d302dbdd424005662dc2a74611dc","name":"Census_Tracts_2010.geojson","type":"application/geo+json","lastModified":1624030021132.619,"size":1298632},{"id":"fb638891b1fe1acde7a23ef391f770835a67610a4ab95b7495c09a8c6f1f71717dd6d66f294325b3f3fd63ef6b40ec7e3e3b28cc6633e057b84bd803d88a6b61","href":"https://static.observableusercontent.com/files/fb638891b1fe1acde7a23ef391f770835a67610a4ab95b7495c09a8c6f1f71717dd6d66f294325b3f3fd63ef6b40ec7e3e3b28cc6633e057b84bd803d88a6b61","name":"Census_Block_Groups_2010.geojson","type":"application/geo+json","lastModified":1624030021132.619,"size":2482102},{"id":"449e620e7479e41a75d2150a5179d02ecc61aeb6e7cf7512770871c1493d67d9da92e5ad8175de33dba82799dbb58672b41e613a99b62a4f86b07d5a2d54ca9f","href":"https://static.observableusercontent.com/files/449e620e7479e41a75d2150a5179d02ecc61aeb6e7cf7512770871c1493d67d9da92e5ad8175de33dba82799dbb58672b41e613a99b62a4f86b07d5a2d54ca9f","name":"ACSDT5Y2019.B09001_data_with_overlays_2021-06-28T160411.csv","type":"text/csv","lastModified":1624911049671.281,"size":72655},{"id":"621d3de3e73f985a02641ff1b195fddbfe648bd7cb9704deada5a0b98aa202b6e8a6ebc704bd8aee62006a454394ed764dd9b9c1f12df263c6519bc0bd78d4dc","href":"https://static.observableusercontent.com/files/621d3de3e73f985a02641ff1b195fddbfe648bd7cb9704deada5a0b98aa202b6e8a6ebc704bd8aee62006a454394ed764dd9b9c1f12df263c6519bc0bd78d4dc","name":"ACSDP5Y2019.DP05_data_with_overlays_2021-06-28T173533.csv","type":"text/csv","lastModified":1624916228874.559,"size":866081}] 2f:{"title":"Philly mapping — COVID-19 Vaccinations Tracts View","cells":[{"id":0,"value":"md`# Philly mapping — COVID-19 Vaccinations Tracts View`","pinned":false,"hidden":false,"mode":"ojs"},{"id":9391,"value":"Select a CT or a percentile (based on the variable in the dropdown in the top right of the map below), and see the places that most people from that selection are visiting with an ordered table that can be used to quickly navigate.","pinned":false,"hidden":false,"mode":"md"},{"id":1735,"value":"$31","pinned":true,"hidden":false,"mode":"ojs"},{"id":9992,"value":"md`*Note: Visit categories, such as unvaccinated visitors, are estimates. Certain businesses and areas may have different policies that could affect the number of visits per category. For example, a university could require vaccination, which would place real unvaccinated visits much lower than our estimates.*`","pinned":false,"hidden":false,"mode":"ojs"},{"id":9602,"value":"import { aq_places_vsites } from \"@wattslab/philly-covid19-optimal-vsites\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":2498,"value":"md`## Settings`","pinned":false,"hidden":false,"mode":"ojs"},{"id":9606,"value":"vaccine_sites_aq_old = aq.from(\n enrolled_Vaccination_Sites_PublicView.features.map(a => {\n const b = a.properties;\n b.lat = a.geometry.coordinates[1];\n b.long = a.geometry.coordinates[0];\n return b;\n })\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":10158,"value":"vaccine_sites_aq = aq\n .from(\n await d3.json(\n `https://data.cdc.gov/resource/5jp2-pgaw.json?loc_admin_city=Philadelphia`\n )\n )\n .rename({\n loc_admin_street1: \"Street_Address\",\n loc_name: \"Site_Name\",\n loc_admin_city: \"City\",\n loc_admin_state: \"State\",\n loc_admin_zip: \"ZIP_Code\",\n latitude: \"lat\",\n longitude: \"long\"\n })\n .derive({\n Type: d => \"Pharmacy\",\n ObjectId: d => op.row_number()\n })","pinned":true,"hidden":false,"mode":"ojs"},{"id":9149,"value":"enrolled_Vaccination_Sites_PublicView = d3.json(\n \"https://services.arcgis.com/fLeGjb7u4uXqeF9q/ArcGIS/rest/services/Enrolled_Vaccination_Sites_PublicView/FeatureServer/0/query?where=State%3D%27PA%27&objectIds=&time=&geometry=&geometryType=esriGeometryEnvelope&inSR=&spatialRel=esriSpatialRelIntersects&resultType=none&distance=0.0&units=esriSRUnit_Meter&returnGeodetic=false&outFields=*&returnGeometry=true&featureEncoding=esriDefault&multipatchOption=xyFootprint&maxAllowableOffset=&geometryPrecision=&outSR=&datumTransformation=&applyVCSProjection=false&returnIdsOnly=false&returnUniqueIdsOnly=false&returnCountOnly=false&returnExtentOnly=false&returnQueryGeometry=false&returnDistinctValues=false&cacheHint=false&orderByFields=&groupByFieldsForStatistics=&outStatistics=&having=&resultOffset=&resultRecordCount=&returnZ=false&returnM=false&returnExceededLimitFeatures=true&quantizationParameters=&sqlFormat=none&f=pgeojson&token=\"\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":3932,"value":"date = \"2021-06-14\"//\"2021-03-29\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":1768,"value":"$32","pinned":true,"hidden":false,"mode":"ojs"},{"id":9613,"value":"vsite_color = ({\n pharmacy: \"rgb(237, 100, 100)\",\n hospital: \"rgb(100, 180, 237)\",\n pdph: \"rgb(60, 199, 190)\",\n clinic: \"rgb(237, 207, 100)\",\n community: \"rgb(100, 237, 132)\",\n other: \"rgb(219, 219, 219)\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":4,"value":"md`## Resources`","pinned":false,"hidden":false,"mode":"ojs"},{"id":10072,"value":"formatNum = (num, decimalplaces = 0) => {\n let decimalplacesfactor = Math.pow(10, decimalplaces);\n return num < 1\n ? Math.round((num * 100) / decimalplacesfactor) + '%'\n : num >= 1000\n ? Math.round((num / 1000) * decimalplacesfactor) / decimalplacesfactor + 'k'\n : Math.round(num * decimalplacesfactor) / decimalplacesfactor;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":10120,"value":"formatNum(0.01, 0)","pinned":true,"hidden":false,"mode":"ojs"},{"id":9728,"value":"luminanceBW = color => {\n let rgb = color.startsWith('#')\n ? hexToRGB(color)\n .split(\"(\")[1]\n .split(\")\")[0]\n .split(',')\n : color\n .split(\"(\")[1]\n .split(\")\")[0]\n .split(',');\n let lum = 0.299 * rgb[0] + 0.587 * rgb[1] + 0.114 * rgb[2];\n return lum > 130 ? 'rgb(0, 0, 0)' : 'rgb(255, 255, 255)';\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":9758,"value":"hexToRGB = (h) => {\n let r = 0, g = 0, b = 0;\n if (h.length == 4) {\n r = \"0x\" + h[1] + h[1];\n g = \"0x\" + h[2] + h[2];\n b = \"0x\" + h[3] + h[3];\n } else if (h.length == 7) {\n r = \"0x\" + h[1] + h[2];\n g = \"0x\" + h[3] + h[4];\n b = \"0x\" + h[5] + h[6];\n }\n return \"rgb(\"+ +r + \",\" + +g + \",\" + +b + \")\";\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":7855,"value":"colorFunction = d3.interpolatePlasma","pinned":false,"hidden":false,"mode":"ojs"},{"id":9428,"value":"styleFunction = (x, c, k, reverse) => ({\n fillColor: colorFunction(\n Math.abs(1 / (1 + Math.pow(Math.E, -k * (x - c))) - (reverse ? 1 : 0))\n ),\n weight: 0.3,\n color: \"black\",\n fillOpacity: 0.3\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":9178,"value":"clusterColorFunction = colorFunction //d3.interpolateHsl(\"rgb(52, 235, 159)\",\"rgb(235, 64, 52)\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":9201,"value":"getClusterColor = (x, c, k) =>\n clusterColorFunction(1 / (1 + Math.pow(Math.E, -k * (x - c))))","pinned":false,"hidden":false,"mode":"ojs"},{"id":9181,"value":"svg_marker = ``","pinned":false,"hidden":false,"mode":"ojs"},{"id":6626,"value":"places_scale_focus = 'total_visits'","pinned":false,"hidden":false,"mode":"ojs"},{"id":3021,"value":"tile_mode = 'light'","pinned":false,"hidden":false,"mode":"ojs"},{"id":5819,"value":"clusterSize = (x, arr) => logisticSize(30, 15, x, 50, 1 / 50)","pinned":false,"hidden":false,"mode":"ojs"},{"id":6013,"value":"tractsOpacity = (x, med, max, min) => logisticSize(max, min, x, med, 4.0 / med)","pinned":false,"hidden":false,"mode":"ojs"},{"id":5999,"value":"logisticSize = (max, min, x, c, k) =>\n (max - min) / (1 + Math.pow(Math.E, -k * (x - c))) + min","pinned":false,"hidden":false,"mode":"ojs"},{"id":3320,"value":"md`## Tracts Data`","pinned":false,"hidden":false,"mode":"ojs"},{"id":8085,"value":"CT_options = aq_tracts_col\n .columnNames()\n .filter(name => !['origin_geoid'].includes(name))\n// .filter((name) => name.includes(\"_\"))\n// .map((name) => name.replace(\"_\", \" \"))","pinned":false,"hidden":false,"mode":"ojs"},{"id":10047,"value":"aq_tracts_collect.view()","pinned":true,"hidden":false,"mode":"ojs"},{"id":9409,"value":"aq_tracts_collect = aq_tracts_alt\n .params({ CT_means_collect })\n .rollup(CT_means_collect)","pinned":true,"hidden":false,"mode":"ojs"},{"id":9411,"value":"CT_means_collect = {\n let obj = {};\n CT_options.forEach(option => {\n obj[option + '_mean'] = `d => op.mean(d['${option}'])`;\n obj[option + '_max'] = `d => op.max(d['${option}'])`;\n obj[option + '_min'] = `d => op.min(d['${option}'])`;\n });\n return obj;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":9440,"value":"CT_rank_collect = {\n let obj = {};\n CT_options.forEach(option => {\n obj[option + `_rank`] = `d => op.rank()`;\n });\n return obj;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":9458,"value":"aq_tracts_ranked = {\n let table = aq_tracts_alt;\n CT_options.forEach(option => {\n let key = option + '_rank';\n let newCollect = {};\n newCollect[key] = `d => op.rank()`;\n table = table.orderby(`${option}`).derive(newCollect);\n });\n return table;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":6653,"value":"aq_tracts_unroll = ct_data\n .join(poi_data, \"safegraph_place_id\")\n .join(aq_tracts_raw, [\"origin_geoid\", \"geoid\"])\n .join(vaccination_ct_data, \"safegraph_place_id\")\n .select(aq.not(\"geoid\")) //, 'type', 'geometry'\n .rename({\n estimated_visits_1: \"geoid_visits\",\n estimated_visits_2: \"place_total_visits\",\n visits_sq_ft_1: \"geoid_sq_ft\",\n visits_sq_ft_2: \"place_total_sq_ft\"\n })\n .derive({\n geoid_latitude: (d) => op.parse_float(d.properties.INTPTLAT10),\n geoid_longitude: (d) => op.parse_float(d.properties.INTPTLON10)\n })\n .derive({\n geoid_place_distance: (d) =>\n op.sqrt(\n op.pow(d.latitude - d.geoid_latitude, 2) +\n op.pow(d.longitude - d.geoid_longitude, 2)\n )\n })\n .derive({\n weighted_visits: (d) => d.geoid_visits * d.geoid_place_distance\n })","pinned":false,"hidden":false,"mode":"ojs"},{"id":6668,"value":"$33","pinned":false,"hidden":false,"mode":"ojs"},{"id":9230,"value":"aq_tracts.get('safegraph_place_info', 0)","pinned":false,"hidden":false,"mode":"ojs"},{"id":7307,"value":"aq_tracts_alt = aq_tracts\n .derive({\n \"CT Number\": (d) => \"CT\" + d.properties.NAME10\n })\n .select(\n aq.not(\n \"age_info\",\n \"ethnicity_info\",\n \"income_info\",\n //\"poverty_info\",\n \"family_info\"\n //\"transportation_info\"\n )\n )\n .join(aq_tracts_col, \"origin_geoid\")\n .relocate(\"CT Number\", { before: \"type\" })\n .filter((r) => r.total_visits > 200)","pinned":false,"hidden":false,"mode":"ojs"},{"id":8228,"value":"aq_tracts_col = aq_tracts\n .derive({\n \"% Under 12\": (d) => d[\"% under 12\"],\n \"% Under 20\": (d) => d.age_info[\"Under 20\"],\n \"% Over 65\": (d) => d.age_info[\"Over 75\"],\n \"% Have Child\": (d) => d.family_info[\"Have Child\"],\n // \"% In Poverty\": d => d.poverty_info[\"Below\"],\n \"% Vaccinated\": (d) => d.vaccine_percent,\n \"# Visits weighted by % Under 12\": (d) => d[\"Total under 12\"],\n \"# Visits weighted by % Unvaccinated\": (d) =>\n d.total_visits - d.visits_vaccinated\n })\n .select(aq.not(aq_tracts.columnNames().filter((e) => e !== \"origin_geoid\")))","pinned":false,"hidden":false,"mode":"ojs"},{"id":4297,"value":"aq_tracts_raw = aq\n .from((await FileAttachment(\"Census_Tracts_2010.geojson\").json()).features)\n .derive({\n geoid: x => x.properties.GEOID10\n })","pinned":false,"hidden":false,"mode":"ojs"},{"id":8251,"value":"aq_tracts_raw.sample(5).get('properties', 0)","pinned":false,"hidden":false,"mode":"ojs"},{"id":51,"value":"phillyCensusBlockGroups = await FileAttachment(\"Census_Block_Groups_2010.geojson\").json()","pinned":false,"hidden":false,"mode":"ojs"},{"id":3927,"value":"FIPS = 42101","pinned":false,"hidden":false,"mode":"ojs"},{"id":8850,"value":"md`## Places Data`","pinned":false,"hidden":false,"mode":"ojs"},{"id":7686,"value":"POI_options = ({\n \"Total Visits\": 'estimated_visits',\n \"% Not Vaccinated\": 'total_unvaccinated_visits'\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":8855,"value":"ct_data_weighted_collect = ct_data_weighted\n .derive({\n geoid_info: d => ({\n origin_geoid: d.origin_geoid,\n estimate_visits: d.estimated_visits,\n visits_sq_ft: d.visits_sq_ft,\n under_20_visits: d.under_20_visits,\n unvaccinated_vists: d.unvaccinated_visits,\n under_12_visits: d.under_12_visits\n })\n })","pinned":false,"hidden":false,"mode":"ojs"},{"id":8911,"value":"ct_data_weighted = ct_data\n .join(demographic_data_alt.rename({ GEOID: 'origin_geoid' }), 'origin_geoid')\n .join(vaccination_ct_data, 'safegraph_place_id')\n .join(processed_under_18_table, 'origin_geoid')\n .derive({\n under_20_visits: d => d.age_info['Under 20'] * d.estimated_visits,\n unvaccinated_visits: d =>\n d.percent_not_vaccinated_visits * d.estimated_visits,\n under_12_visits: d => d['% under 12'] * d.estimated_visits\n })\n .select(\n aq.not(\n demographic_data_alt\n .columnNames()\n .concat(\n vaccination_ct_data.columnNames(),\n processed_under_18_table.columnNames()\n )\n .filter(e => ['origin_geoid', 'safegraph_place_id'].indexOf(e) < 0)\n )\n )","pinned":false,"hidden":false,"mode":"ojs"},{"id":8865,"value":"aq_places = ct_data_weighted_collect\n //.filter(`p => p.total_visits > ${places_filter}`)\n .groupby(\"safegraph_place_id\")\n .rollup({\n origin_geoids_info: op.array_agg('geoid_info'),\n total_visits: op.sum('estimated_visits'),\n total_unvaccinated_visits: op.sum('unvaccinated_visits')\n })\n .join_left(poi_data, 'safegraph_place_id')","pinned":false,"hidden":false,"mode":"ojs"},{"id":4233,"value":"ct_data = aq.from(\n await d3.csv(\n `${base_URL}networks_${FIPS}/bipartite_CT_${date}.csv`\n )\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":4401,"value":"poi_data = aq.from(\n await d3.csv(`${base_URL}networks_${FIPS}/poi_summary_${date}.csv`)\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":8777,"value":"## Demographic Data","pinned":false,"hidden":false,"mode":"md"},{"id":6830,"value":"import { demographic_data } from \"@wattslab/covid-philly-census-tracts\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":10245,"value":"import { demog_data } from \"@wattslab/philadelphia-census-tract-data\"","pinned":true,"hidden":false,"mode":"ojs"},{"id":8915,"value":"demog_data.view()","pinned":true,"hidden":false,"mode":"ojs"},{"id":6960,"value":"$34","pinned":false,"hidden":false,"mode":"ojs"},{"id":8657,"value":"processed_under_18_table = raw_under_18_table\n .rename(\n raw_under_18_table\n .filter((r) => r.geoid == \"id\")\n .select(aq.matches(\"B0900\"))\n .fold(aq.all())\n .filter((r) => op.includes(r.value, \"Estimate\"))\n .derive({\n value: (r) =>\n op.replace(r.value, \"Estimate!!Total:!!In households:!!\", \"\")\n })\n .filter((r) => !op.includes(r.value, \"Estimate\"))\n .pivot(\"key\", \"value\")\n .object()\n )\n .filter((r) => r.geoid != \"id\")\n .select(aq.not(\"NAME\", aq.matches(\"B09001\")))\n .join(raw_population_data)\n .derive({\n \"Total under 12\": (r) =>\n op.parse_int(r[\"Under 3 years\"]) +\n op.parse_int(r[\"3 and 4 years\"]) +\n op.parse_int(r[\"5 years\"]) +\n op.parse_int(r[\"6 to 8 years\"]) +\n op.parse_int(r[\"9 to 11 years\"])\n })\n .derive({ \"% under 12\": (r) => r[\"Total under 12\"] / r.total_population })\n .select(\"geoid\", \"Total under 12\", \"% under 12\")\n .rename({ geoid: \"origin_geoid\" })","pinned":false,"hidden":false,"mode":"ojs"},{"id":8625,"value":"raw_under_18_table = aq\n .from(\n await FileAttachment(\n \"ACSDT5Y2019.B09001_data_with_overlays_2021-06-28T160411.csv\"\n ).csv()\n )\n .rename({ GEO_ID: \"geoid\" })\n .derive({ geoid: (r) => op.replace(r.geoid, \"1400000US\", \"\") })","pinned":false,"hidden":false,"mode":"ojs"},{"id":8810,"value":"raw_population_data = aq\n .from(\n await FileAttachment(\n \"ACSDP5Y2019.DP05_data_with_overlays_2021-06-28T173533.csv\"\n ).csv()\n )\n .rename({ GEO_ID: \"geoid\", DP05_0001E: \"total_population\" })\n .filter((r) => r.geoid != \"id\")\n .derive({ geoid: (r) => op.replace(r.geoid, \"1400000US\", \"\") })\n .select(\"geoid\", \"total_population\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":7478,"value":"md`## Vaccine Data`","pinned":false,"hidden":false,"mode":"ojs"},{"id":7483,"value":"raw_vaccination_data = aq.from(\n (\n await d3.json(\n \"https://phl.carto.com/api/v2/sql?q=SELECT * FROM covid_vaccines_by_census_tract\"\n )\n ).rows\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":7479,"value":"clean_vaccination_data = raw_vaccination_data\n .select(\"census_tract_geoid\", \"dose1_percent_range\", \"etl_timestamp\")\n .derive({\n vaccine_ordinal: (r) =>\n op.indexof(\n [\n null,\n \"<20%\",\n \"20-29%\",\n \"30-39%\",\n \"40-49%\",\n \"50-59%\",\n \"60-69%\",\n \"70%+\"\n ],\n r.dose1_percent_range\n ),\n vaccine_percent: (r) =>\n op.recode(r.dose1_percent_range, {\n null: 0,\n \"<20%\": 0.1,\n \"20-29%\": 0.25,\n \"30-39%\": 0.35,\n \"40-49%\": 0.45,\n \"50-59%\": 0.55,\n \"60-69%\": 0.65,\n \"70%+\": 0.85\n }),\n etl_timestamp: (r) => op.parse_date(r.etl_timestamp)\n })","pinned":false,"hidden":false,"mode":"ojs"},{"id":7503,"value":"vaccination_ct_data = ct_data\n .join(\n clean_vaccination_data\n .rename({ census_tract_geoid: \"origin_geoid\" })\n .select(\"origin_geoid\", \"vaccine_percent\")\n )\n .join(\n demographic_data_alt\n .rename({ GEOID: \"origin_geoid\" })\n .derive({ under_20_percent: r => r.age_info[\"Under 20\"] })\n .select(\"origin_geoid\", \"under_20_percent\")\n )\n .derive({\n vaccine_weighted_visits: d => d.estimated_visits * d.vaccine_percent\n })\n .groupby(\"safegraph_place_id\")\n .rollup({\n percent_vaccinated_visits: r =>\n op.sum(r.vaccine_weighted_visits) / op.sum(r.estimated_visits),\n total_vaccinated_visits: r => op.sum(r.vaccine_weighted_visits),\n total_visits: r => op.sum(r.estimated_visits)\n })\n .derive({\n percent_not_vaccinated_visits: d => 1 - d.percent_vaccinated_visits\n })\n .ungroup()","pinned":false,"hidden":false,"mode":"ojs"},{"id":3323,"value":"md`## Libraries and utilities`","pinned":false,"hidden":false,"mode":"ojs"},{"id":10238,"value":"base_URL = \"https://upenn-seas-wattscovid19lab.s3.us-east-2.amazonaws.com/paco/stats/data/\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":3806,"value":"sortTimes = (a, b) => a.getTime() - b.getTime()","pinned":false,"hidden":false,"mode":"ojs"},{"id":15,"value":"d3 = require('d3@5', 'd3-array@2', 'd3-regression@1')","pinned":false,"hidden":false,"mode":"ojs"},{"id":10032,"value":"vramp = (color, n = 512, maxHeight = 150, reversed = false) => {\n const canvas = DOM.canvas(1, n);\n const context = canvas.getContext(\"2d\");\n canvas.style.maxHeight = maxHeight + \"px\";\n canvas.style.width = \"20px\";\n if (reversed) canvas.style.transform = \"scaleY(-1)\";\n for (let i = 0; i < n; i++) {\n context.fillStyle = color(i / (n - 1));\n context.fillRect(0, i, 1, 1);\n }\n return canvas;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":890,"value":"import { vl } from '@vega/vega-lite-api'","pinned":false,"hidden":false,"mode":"ojs"},{"id":4039,"value":"import { aq, op } from '@uwdata/arquero'","pinned":false,"hidden":false,"mode":"ojs"},{"id":33,"value":"L = require('leaflet')","pinned":false,"hidden":false,"mode":"ojs"},{"id":41,"value":"html``","pinned":false,"hidden":false,"mode":"ojs"},{"id":5513,"value":"lmc = {\n L;\n const leafletmarkercluster = await require('leaflet.markercluster@1');\n if (!leafletmarkercluster._style) {\n const href = await require.resolve(\n \"leaflet.markercluster@1/dist/MarkerCluster.css\"\n );\n const href2 = await require.resolve(\n \"leaflet.markercluster@1/dist/MarkerCluster.Default.css\"\n );\n document.head.appendChild(\n (leafletmarkercluster._style = html``)\n );\n }\n return leafletmarkercluster;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":9833,"value":"lspin = require(\"leaflet-spin@1.1.2/leaflet.spin.js\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":3725,"value":"scope = (value, min, max) =>\n value > min ? (max ? (value < max ? value : max) : value) : min","pinned":false,"hidden":false,"mode":"ojs"},{"id":3139,"value":"stats = require('simple-statistics')","pinned":false,"hidden":false,"mode":"ojs"},{"id":2682,"value":"import { summarize, sigFigs, set } from \"@markwhiting/stats\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":188,"value":"import { radio, select, slider } from \"@jashkenas/inputs\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":5308,"value":"import {\n Button,\n Checkbox,\n Toggle,\n Radio,\n Range,\n Select,\n Text,\n Textarea,\n Search,\n Table\n} from \"@observablehq/inputs\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":861,"value":"import { Scrubber } from \"@mbostock/scrubber\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":3880,"value":"import { SiteList } from '@markwhiting/covid-philly'","pinned":false,"hidden":false,"mode":"ojs"},{"id":9863,"value":"import { ora } from \"@victordidenko/ora\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":3899,"value":"title = \"Mapping dashboard\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":10140,"value":"esriLeaflet = require(\"https://unpkg.com/esri-leaflet@2.5.3/dist/esri-leaflet.js\")","pinned":true,"hidden":false,"mode":"ojs"}],"files":"$W35","stdlib":"1","titling":"auto"}