#!/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_county col="COUNTYNS_INT int"
#v.db.update tl_2022_us_county col=COUNTYNS_INT qcol=COUNTYNS
v.to.rast tl_2022_us_county out=conus_countyns use=attr acol=COUNTYNS_INT
r.stats -cn conus_countyns > conus_countyns.cells
awk '{printf "%08d\n", $1}' conus_countyns.cells > conus_countyns.txt