Delphi_utils geomapper: update to fips 2020 population estimates

Author: dshemetov

_delphi_utils_python/data_proc/geomap/consistency_checks.ipynb

@@ -24,7 +24,7 @@
 FIPS_MSA_URL = "https://www2.census.gov/programs-surveys/metro-micro/geographies/reference-files/2018/delineation-files/list1_Sep_2018.xls"
 JHU_FIPS_URL = "https://raw.githubusercontent.com/CSSEGISandData/COVID-19/master/csse_covid_19_data/UID_ISO_FIPS_LookUp_Table.csv"
 STATE_CODES_URL = "http://www2.census.gov/geo/docs/reference/state.txt?#"
-FIPS_POPULATION_URL = "https://www2.census.gov/programs-surveys/popest/datasets/2010-2019/counties/totals/co-est2019-alldata.csv"
+FIPS_POPULATION_URL = "https://www2.census.gov/programs-surveys/popest/datasets/2010-2020/counties/totals/co-est2020-alldata.csv"
 FIPS_PUERTO_RICO_POPULATION_URL = (
     "https://www2.census.gov/geo/docs/maps-data/data/rel/zcta_county_rel_10.txt?"
 )
@@ -360,32 +360,29 @@ def create_fips_population_table():
     )
     census_pop["pop"] = census_pop["POPESTIMATE2019"]
     census_pop = census_pop[["fips", "pop"]]
+
+    # Set population for Dukes and Nantucket combo county
+    dukes_pop = int(census_pop.loc[census_pop["fips"] == "25007", "pop"])
+    nantu_pop = int(census_pop.loc[census_pop["fips"] == "25019", "pop"])
     census_pop = pd.concat(
         [
             census_pop,
             pd.DataFrame(
                 {
-                    "fips": ["70002", "70003"],
-                    "pop": [0, 0],
+                    "fips": [
+                        "70002",  # Dukes and Nantucket combo county
+                        "70003"   # Kansas City
+                    ],
+                    "pop": [
+                        dukes_pop + nantu_pop,
+                        491918    # via Google
+                    ],
                 }
             ),
         ]
     )
     census_pop = census_pop.reset_index(drop=True)
 
-    # Set population for Dukes and Nantucket
-    dn_fips = "70002"
-    dukes_fips = "25007"
-    nantu_fips = "25019"
-
-    census_pop.loc[census_pop["fips"] == dn_fips, "pop"] = (
-        census_pop.loc[census_pop["fips"] == dukes_fips, "pop"].values
-        + census_pop.loc[census_pop["fips"] == nantu_fips, "pop"].values
-    )
-
-    # Set population for Kansas City
-    census_pop.loc[census_pop["fips"] == "70003", "pop"] = 491918  # via Google
-
     # Get the file with Puerto Rico populations
     df_pr = pd.read_csv(FIPS_PUERTO_RICO_POPULATION_URL)
     df_pr["fips"] = df_pr["STATE"].astype(str).str.zfill(2) + df_pr["COUNTY"].astype(

_delphi_utils_python/data_proc/geomap/geo_data_proc.py

@@ -71,6 +71,8 @@ fips,hhs
 02020,10
 02050,10
 02060,10
+02063,10
+02066,10
 02068,10
 02070,10
 02090,10
@@ -91,7 +93,6 @@ fips,hhs
 02220,10
 02230,10
 02240,10
-02261,10
 02275,10
 02282,10
 02290,10
@@ -3234,60 +3235,60 @@ fips,hhs
 69100,9
 69110,9
 69120,9
-78000,2
-33000,1
-09000,1
-32000,9
-04000,9
-30000,8
-15000,9
-69000,9
-25000,1
-19000,7
 01000,4
-51000,3
-41000,10
 66000,9
-22000,6
-06000,9
-39000,5
-38000,8
-53000,10
-12000,4
-46000,8
-17000,5
-48000,6
-34000,2
 08000,8
+12000,4
+37000,4
+09000,1
 36000,2
-40000,6
-47000,4
+29000,7
+21000,4
+35000,6
 49000,8
-45000,4
-55000,5
-11000,3
-02000,10
+69000,9
+18000,5
+24000,3
+78000,2
+26000,5
+46000,8
+38000,8
+10000,3
 72000,2
-23000,1
-16000,10
-44000,1
+05000,6
+40000,6
+22000,6
+02000,10
+17000,5
+42000,3
+55000,5
+31000,7
+32000,9
+53000,10
 54000,3
-10000,3
-21000,4
+60000,9
+44000,1
+19000,7
+30000,8
+27000,5
 28000,4
-18000,5
-05000,6
+33000,1
+34000,2
 56000,8
-24000,3
+45000,4
+48000,6
+15000,9
+06000,9
+41000,10
 50000,1
-37000,4
-20000,7
-31000,7
-29000,7
-27000,5
+23000,1
 13000,4
-42000,3
-35000,6
-60000,9
-26000,5
 70000,9
+16000,10
+20000,7
+04000,9
+51000,3
+11000,3
+39000,5
+25000,1
+47000,4