Published
Edited
May 14, 2020
Insert cell
md`# Choropleth Mapping [GEOG 3540]
OWI's Per County in Iowa`
Insert cell
d3 = require("d3@5")
Insert cell
import {legend} from "@d3/color-legend"
Insert cell
simple = require("simple-statistics@7.0.7/dist/simple-statistics.min.js")
Insert cell
format = d => `${d}%`
Insert cell
topojson = require("topojson-client@3")
Insert cell
Iowa = FileAttachment("iowa_counties_topo.json").json()
Insert cell
counties = topojson.feature(Iowa, Iowa.objects.iowa_counties)
Insert cell
csv_data = d3.csvParse(await FileAttachment("data4.csv").text(),({FIPS, DRUNKDRIVING,POP2010 }) => [FIPS, [+DRUNKDRIVING/+POP2010]])
Insert cell
data = Object.assign(new Map(csv_data), {title: "Percentage of OWI's Per Iowa County"})
Insert cell
DrunkDrive = Array.from(csv_data.values(), function(d){
if(d[1][0])
return d[1][0];
else
return 0;
})
Insert cell
YlGnBu = [d3.color("#ffffcc"), d3.color("#a1dab4"), d3.color("#41b6c4"), d3.color("#2c7fb8"),d3.color("#253494")]
Insert cell
naturalbreaks = simple.ckmeans(DrunkDrive, YlGnBu.length).map(v => v.pop())
Insert cell
//more information on sequential scales: https://observablehq.com/@d3/sequential-scales
// color = d3.scaleSequentialQuantile([...data.values()], d3.interpolateBlues)

// color = d3.scaleQuantile()
// .domain(med_age)
// .range()

color = d3.scaleThreshold()
.domain(naturalbreaks)
.range(YlGnBu)
Insert cell
width = 975
Insert cell
height = 610
Insert cell
margin = 100
Insert cell
//Rotate the map sets the longitude of origin for our UTM Zone 15N projection.
projection = d3.geoTransverseMercator().rotate([94,0]).fitExtent([[80, 80], [width, height]], counties);
//d3 reference for projections: https://github.com/d3/d3-geo/blob/master/README.md

//use the following url for specific projection settings: https://github.com/veltman/d3-stateplane
//Use this code to set up the map projection (if different than geographic projection)

//projection = d3.geoAlbers().fitExtent([[margin, margin], [width - margin, height - margin]], counties)

//projection = d3.geoMercator().fitExtent([[margin, margin], [width - margin, height - margin]], counties)
Insert cell
//Using a path generator to project geometry onto the map
path = d3.geoPath().projection(projection);
Insert cell
choropleth = {
const svg = d3.create("svg")
.attr("viewBox", [0, 0, width, height]);

svg.append("g")
.attr("transform", "translate(360,20)")
.append(() =>
legend({
color: color,
title: data.title,
width: 260,
tickFormat: ".3f"
})
);

svg.append("g")
.selectAll("path")
.data(counties.features)
.join("path")
.attr("stroke", "black")
.attr("stroke-linejoin", "round")
.attr("stroke-width", 0.5)
// .attr("fill", function(d){
// console.log(color(data.get(d.properties.FIPS)[0]))
// return color(data.get(d.properties.FIPS)[0]);
// })
.attr("fill", function(d){
if(data.get(d.properties.FIPS))
return color(data.get(d.properties.FIPS)[0]);
else
return "black";
})
.attr("d", path)
.append("title")
.text(d => " Percent of OWI's: " + data.get(d.properties.FIPS));

return svg.node();
}
Insert cell

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.
Learn more