#!/bin/sh
# run from epsg5070/PERMANENT
g.region rast=ned_1_conus_clean -ap
r.mask ned_1_conus_clean --o
#v.db.addcolumn tl_2022_us_state col="STATEFP_INT int"
#v.db.update tl_2022_us_state col=STATEFP_INT qcol=STATEFP
v.to.rast tl_2022_us_state out=conus_statefp use=attr acol=STATEFP_INT
r.stats -cn conus_statefp > conus_statefp.cells
awk '{printf "%02d\n", $1}' conus_statefp.cells > conus_statefp.txt