Skip to content

Instantly share code, notes, and snippets.

@mbostock
Last active June 11, 2024 21:02
Show Gist options
  • Save mbostock/2522624ada2c1f9e0fafb75cca09442b to your computer and use it in GitHub Desktop.
Save mbostock/2522624ada2c1f9e0fafb75cca09442b to your computer and use it in GitHub Desktop.
U.S. Population Density
license:gpl-3.0
height:680
border:no
cb_*
topo.json
<!DOCTYPE html>
<img src="topo.svg" width="960" height="680">
Display the source blob
Display the rendered blob
Raw
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
{
"private": true,
"license": "gpl-3.0",
"author": {
"name": "Mike Bostock",
"url": "https://bost.ocks.org/mike"
},
"scripts": {
"prepublish": "bash prepublish"
},
"devDependencies": {
"d3-scale": "^1.0.4",
"d3-scale-chromatic": "^1.1.0",
"d3-geo-projection": "^1.2.1",
"ndjson-cli": "^0.3.0",
"shapefile": "^0.5.9",
"topojson-server": "^2.0.0",
"topojson-client": "^2.1.0",
"topojson-simplify": "^2.0.0"
}
}
#!/bin/bash
# U.S. Albers
PROJECTION='d3.geoAlbersUsa().scale(1280).translate([480, 300])'
# The state FIPS codes.
STATES="01 02 04 05 06 08 09 10 11 12 13 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 44 45 46 47 48 49 50 51 53 54 55 56"
# The ACS 5-Year Estimate vintage.
YEAR=2014
# The display size.
WIDTH=960
HEIGHT=680
# Download the census tract boundaries.
# Extract the shapefile (.shp) and dBASE (.dbf).
# Download the census tract population estimates.
for STATE in ${STATES}; do
if [ ! -f cb_${YEAR}_${STATE}_tract_500k.shp ]; then
curl -o cb_${YEAR}_${STATE}_tract_500k.zip \
"http://www2.census.gov/geo/tiger/GENZ${YEAR}/shp/cb_${YEAR}_${STATE}_tract_500k.zip"
unzip -o \
cb_${YEAR}_${STATE}_tract_500k.zip \
cb_${YEAR}_${STATE}_tract_500k.shp \
cb_${YEAR}_${STATE}_tract_500k.dbf
fi
if [ ! -f cb_${YEAR}_${STATE}_tract_B01003.json ]; then
curl -o cb_${YEAR}_${STATE}_tract_B01003.json \
"http://api.census.gov/data/${YEAR}/acs5?get=B01003_001E&for=tract:*&in=state:${STATE}&key=${CENSUS_KEY}"
fi
done
# Construct TopoJSON.
if [ ! -f topo.json ]; then
geo2topo -n \
tracts=<(for STATE in ${STATES}; do \
ndjson-join 'd.id' \
<(shp2json -n cb_${YEAR}_${STATE}_tract_500k.shp \
| geoproject -n "${PROJECTION}" \
| ndjson-map 'd.id = d.properties.GEOID, d') \
<(ndjson-cat cb_${YEAR}_${STATE}_tract_B01003.json \
| ndjson-split 'd.slice(1)' \
| ndjson-map '{id: d[1] + d[2] + d[3], B01003: +d[0]}') \
| ndjson-map -r d3=d3-array 'd[0].properties = {density: d3.bisect([1, 10, 50, 200, 500, 1000, 2000, 4000], (d[1].B01003 / d[0].properties.ALAND || 0) * 2589975.2356)}, d[0]'; \
done) \
| topomerge -k 'd.id.slice(0, 5)' counties=tracts \
| topomerge -k 'd.id.slice(0, 2)' states=counties \
| topomerge --mesh -f 'a !== b' counties=counties \
| topomerge --mesh -f 'a !== b' states=states \
| topomerge -k 'd.properties.density' tracts=tracts \
| toposimplify -p 1 -f \
| topoquantize 1e5 \
> topo.json
fi
# Convert to SVG (while dropping the last line).
cat \
<(topo2geo -n < topo.json tracts=- \
| ndjson-map -r d3=d3-scale-chromatic '(d.properties.title = d.id, d.properties.fill = d3.schemeOrRd[9][d.id], d)') \
<(topo2geo -n < topo.json counties=- \
| ndjson-map '(d.properties.stroke = "black", d.properties.strokeWidth = 0.5, d.properties.strokeOpacity = 0.3, d)') \
<(topo2geo -n < topo.json states=- \
| ndjson-map '(d.properties.stroke = "black", d.properties.strokeOpacity = 0.3, d)') \
| geo2svg --stroke=none -n -p 1 -w ${WIDTH} -h ${HEIGHT} \
| sed '$d' \
> topo.svg
# Insert the legend.
tail -n +4 \
< legend.svg \
>> topo.svg
Display the source blob
Display the rendered blob
Raw
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@zhx995
Copy link

zhx995 commented Jun 11, 2024

Hello, I want to ask if this has GeoJSON?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment