publicdata/j40-cejst-2
1
0
Fork 0
mirror of https://github.com/DOI-DO/j40-cejst-2.git synced 2025-02-23 10:04:18 -08:00
Code Issues Projects Releases Packages Wiki Activity Actions
j40-cejst-2/data/data-pipeline/data_pipeline/ipython/tribal_and_tracts_overlap.ipynb

1423 lines
54 KiB
Text
Raw Normal View History

Backend release branch to main (#1822) * Create deploy_be_staging.yml (#1575) * Imputing income using geographic neighbors (#1559) Imputes income field with a light refactor. Needs more refactor and more tests (I spotchecked). Next ticket will check and address but a lot of "narwhal" architecture is here. * Adding HOLC indicator (#1579) Added HOLC indicator (Historic Redlining Score) from NCRC work; included 3.25 cutoff and low income as part of the housing burden category. * Update backend for Puerto Rico (#1686) * Update PR threshold count to 10 We now show 10 indicators for PR. See the discussion on the github issue for more info: https://github.com/usds/justice40-tool/issues/1621 * Do not use linguistic iso for Puerto Rico Closes 1350. Co-authored-by: Shelby Switzer <shelbyswitzer@gmail.com> * updating * Do not drop Guam and USVI from ETL (#1681) * Remove code that drops Guam and USVI from ETL * Add back code for dropping rows by FIPS code We may want this functionality, so let's keep it and just make the constant currently be an empty array. Co-authored-by: Shelby Switzer <shelbyswitzer@gmail.com> * Emma nechamkin/holc patch (#1742) Removing HOLC calculation from score narwhal. * updating ejscreen data, try two (#1747) * Rescaling linguistic isolation (#1750) Rescales linguistic isolation to drop puerto rico * adds UST indicator (#1786) adds leaky underground storage tanks * Changing LHE in tiles to a boolean (#1767) also includes merging / clean up of the release * added indoor plumbing to chas * added indoor plumbing to score housing burden * added indoor plumbing to score housing burden * first run through * Refactor DOE Energy Burden and COI to use YAML (#1796) * added tribalId for Supplemental dataset (#1804) * Setting zoom levels for tribal map (#1810) * NRI dataset and initial score YAML configuration (#1534) * update be staging gha * NRI dataset and initial score YAML configuration * checkpoint * adding data checks for release branch * passing tests * adding INPUT_EXTRACTED_FILE_NAME to base class * lint * columns to keep and tests * update be staging gha * checkpoint * update be staging gha * NRI dataset and initial score YAML configuration * checkpoint * adding data checks for release branch * passing tests * adding INPUT_EXTRACTED_FILE_NAME to base class * lint * columns to keep and tests * checkpoint * PR Review * renoving source url * tests * stop execution of ETL if there's a YAML schema issue * update be staging gha * adding source url as class var again * clean up * force cache bust * gha cache bust * dynamically set score vars from YAML * docsctrings * removing last updated year - optional reverse percentile * passing tests * sort order * column ordening * PR review * class level vars * Updating DatasetsConfig * fix pylint errors * moving metadata hint back to code Co-authored-by: lucasmbrown-usds <lucas.m.brown@omb.eop.gov> * Correct copy typo (#1809) * Add basic test suite for COI (#1518) * Update COI to use new yaml (#1518) * Add tests for DOE energy budren (1518 * Add dataset config for energy budren (1518) * Refactor ETL to use datasets.yml (#1518) * Add fake GEOIDs to COI tests (#1518) * Refactor _setup_etl_instance_and_run_extract to base (#1518) For the three classes we've done so far, a generic _setup_etl_instance_and_run_extract will work fine, for the moment we can reuse the same setup method until we decide future classes need more flexibility --- but they can also always subclass so... * Add output-path tests (#1518) * Update YAML to match constant (#1518) * Don't blindly set float format (#1518) * Add defaults for extract (#1518) * Run YAML load on all subclasses (#1518) * Update description fields (#1518) * Update YAML per final format (#1518) * Update fixture tract IDs (#1518) * Update base class refactor (#1518) Now that NRI is final I needed to make a small number of updates to my refactored code. * Remove old comment (#1518) * Fix type signature and return (#1518) * Update per code review (#1518) Co-authored-by: Jorge Escobar <83969469+esfoobar-usds@users.noreply.github.com> Co-authored-by: lucasmbrown-usds <lucas.m.brown@omb.eop.gov> Co-authored-by: Vim <86254807+vim-usds@users.noreply.github.com> * Update etl_score_geo.py Yikes! Fixing merge messup! * Create deploy_be_staging.yml (#1575) * Imputing income using geographic neighbors (#1559) Imputes income field with a light refactor. Needs more refactor and more tests (I spotchecked). Next ticket will check and address but a lot of "narwhal" architecture is here. * Adding HOLC indicator (#1579) Added HOLC indicator (Historic Redlining Score) from NCRC work; included 3.25 cutoff and low income as part of the housing burden category. * Update backend for Puerto Rico (#1686) * Update PR threshold count to 10 We now show 10 indicators for PR. See the discussion on the github issue for more info: https://github.com/usds/justice40-tool/issues/1621 * Do not use linguistic iso for Puerto Rico Closes 1350. Co-authored-by: Shelby Switzer <shelbyswitzer@gmail.com> * updating * Do not drop Guam and USVI from ETL (#1681) * Remove code that drops Guam and USVI from ETL * Add back code for dropping rows by FIPS code We may want this functionality, so let's keep it and just make the constant currently be an empty array. Co-authored-by: Shelby Switzer <shelbyswitzer@gmail.com> * Emma nechamkin/holc patch (#1742) Removing HOLC calculation from score narwhal. * updating ejscreen data, try two (#1747) * Rescaling linguistic isolation (#1750) Rescales linguistic isolation to drop puerto rico * adds UST indicator (#1786) adds leaky underground storage tanks * Changing LHE in tiles to a boolean (#1767) also includes merging / clean up of the release * added indoor plumbing to chas * added indoor plumbing to score housing burden * added indoor plumbing to score housing burden * first run through * Refactor DOE Energy Burden and COI to use YAML (#1796) * added tribalId for Supplemental dataset (#1804) * Setting zoom levels for tribal map (#1810) * NRI dataset and initial score YAML configuration (#1534) * update be staging gha * NRI dataset and initial score YAML configuration * checkpoint * adding data checks for release branch * passing tests * adding INPUT_EXTRACTED_FILE_NAME to base class * lint * columns to keep and tests * update be staging gha * checkpoint * update be staging gha * NRI dataset and initial score YAML configuration * checkpoint * adding data checks for release branch * passing tests * adding INPUT_EXTRACTED_FILE_NAME to base class * lint * columns to keep and tests * checkpoint * PR Review * renoving source url * tests * stop execution of ETL if there's a YAML schema issue * update be staging gha * adding source url as class var again * clean up * force cache bust * gha cache bust * dynamically set score vars from YAML * docsctrings * removing last updated year - optional reverse percentile * passing tests * sort order * column ordening * PR review * class level vars * Updating DatasetsConfig * fix pylint errors * moving metadata hint back to code Co-authored-by: lucasmbrown-usds <lucas.m.brown@omb.eop.gov> * Correct copy typo (#1809) * Add basic test suite for COI (#1518) * Update COI to use new yaml (#1518) * Add tests for DOE energy budren (1518 * Add dataset config for energy budren (1518) * Refactor ETL to use datasets.yml (#1518) * Add fake GEOIDs to COI tests (#1518) * Refactor _setup_etl_instance_and_run_extract to base (#1518) For the three classes we've done so far, a generic _setup_etl_instance_and_run_extract will work fine, for the moment we can reuse the same setup method until we decide future classes need more flexibility --- but they can also always subclass so... * Add output-path tests (#1518) * Update YAML to match constant (#1518) * Don't blindly set float format (#1518) * Add defaults for extract (#1518) * Run YAML load on all subclasses (#1518) * Update description fields (#1518) * Update YAML per final format (#1518) * Update fixture tract IDs (#1518) * Update base class refactor (#1518) Now that NRI is final I needed to make a small number of updates to my refactored code. * Remove old comment (#1518) * Fix type signature and return (#1518) * Update per code review (#1518) Co-authored-by: Jorge Escobar <83969469+esfoobar-usds@users.noreply.github.com> Co-authored-by: lucasmbrown-usds <lucas.m.brown@omb.eop.gov> Co-authored-by: Vim <86254807+vim-usds@users.noreply.github.com> * Update etl_score_geo.py Yikes! Fixing merge messup! * updated to fix linting errors (#1818) Cleans and updates base branch * Adding back MapComparison video * Add FUDS ETL (#1817) * Add spatial join method (#1871) Since we'll need to figure out the tracts for a large number of points in future tickets, add a utility to handle grabbing the tract geometries and adding tract data to a point dataset. * Add FUDS, also jupyter lab (#1871) * Add YAML configs for FUDS (#1871) * Allow input geoid to be optional (#1871) * Add FUDS ETL, tests, test-datae noteobook (#1871) This adds the ETL class for Formerly Used Defense Sites (FUDS). This is different from most other ETLs since these FUDS are not provided by tract, but instead by geographic point, so we need to assign FUDS to tracts and then do calculations from there. * Floats -> Ints, as I intended (#1871) * Floats -> Ints, as I intended (#1871) * Formatting fixes (#1871) * Add test false positive GEOIDs (#1871) * Add gdal binaries (#1871) * Refactor pandas code to be more idiomatic (#1871) Per Emma, the more pandas-y way of doing my counts is using np.where to add the values i need, then groupby and size. It is definitely more compact, and also I think more correct! * Update configs per Emma suggestions (#1871) * Type fixed! (#1871) * Remove spurious import from vscode (#1871) * Snapshot update after changing col name (#1871) * Move up GDAL (#1871) * Adjust geojson strategy (#1871) * Try running census separately first (#1871) * Fix import order (#1871) * Cleanup cache strategy (#1871) * Download census data from S3 instead of re-calculating (#1871) * Clarify pandas code per Emma (#1871) * Disable markdown check for link * Adding DOT composite to travel score (#1820) This adds the DOT dataset to the ETL and to the score. Note that currently we take a percentile of an average of percentiles. * Adding first street foundation data (#1823) Adding FSF flood and wildfire risk datasets to the score. * first run -- adding NCLD data to the ETL, but not yet to the score * Add abandoned mine lands data (#1824) * Add notebook to generate test data (#1780) * Add Abandoned Mine Land data (#1780) Using a similar structure but simpler apporach compared to FUDs, add an indicator for whether a tract has an abandonded mine. * Adding some detail to dataset readmes Just a thought! * Apply feedback from revieiw (#1780) * Fixup bad string that broke test (#1780) * Update a string that I should have renamed (#1780) * Reduce number of threads to reduce memory pressure (#1780) * Try not running geo data (#1780) * Run the high-memory sets separately (#1780) * Actually deduplicate (#1780) * Add flag for memory intensive ETLs (#1780) * Document new flag for datasets (#1780) * Add flag for new datasets fro rebase (#1780) Co-authored-by: Emma Nechamkin <97977170+emma-nechamkin@users.noreply.github.com> * Adding NLCD data (#1826) Adding NLCD's natural space indicator end to end to the score. * Add donut hole calculation to score (#1828) Adds adjacency index to the pipeline. Requires thorough QA * Adding eamlis and fuds data to legacy pollution in score (#1832) Update to add EAMLIS and FUDS data to score * Update to use new FSF files (#1838) backend is partially done! * Quick fix to kitchen or plumbing indicator Yikes! I think I messed something up and dropped the pctile field suffix from when the KP score gets calculated. Fixing right quick. * Fast flag update (#1844) Added additional flags for the front end based on our conversation in stand up this morning. * Tiles fix (#1845) Fixes score-geo and adds flags * Update etl_score_geo.py * Issue 1827: Add demographics to tiles and download files (#1833) * Adding demographics for use in sidebar and download files * Updates backend constants to N (#1854) * updated to show T/F/null vs T/F for AML and FUDS (#1866) * fix markdown * just testing that the boolean is preserved on gha * checking drop tracts works * OOPS! Old changes persisted * adding a check to the agvalue calculation for nri * updated with error messages * updated error message * tuple type * Score tests (#1847) * update Python version on README; tuple typing fix * Alaska tribal points fix (#1821) * Bump mistune from 0.8.4 to 2.0.3 in /data/data-pipeline (#1777) Bumps [mistune](https://github.com/lepture/mistune) from 0.8.4 to 2.0.3. - [Release notes](https://github.com/lepture/mistune/releases) - [Changelog](https://github.com/lepture/mistune/blob/master/docs/changes.rst) - [Commits](https://github.com/lepture/mistune/compare/v0.8.4...v2.0.3) --- updated-dependencies: - dependency-name: mistune dependency-type: indirect ... Signed-off-by: dependabot[bot] <support@github.com> Signed-off-by: dependabot[bot] <support@github.com> Co-authored-by: dependabot[bot] <49699333+dependabot[bot]@users.noreply.github.com> * poetry update * initial pass of score tests * add threshold tests * added ses threshold (not donut, not island) * testing suite -- stopping for the day * added test for lead proxy indicator * Refactor score tests to make them less verbose and more direct (#1865) * Cleanup tests slightly before refactor (#1846) * Refactor score calculations tests * Feedback from review * Refactor output tests like calculatoin tests (#1846) (#1870) * Reorganize files (#1846) * Switch from lru_cache to fixture scorpes (#1846) * Add tests for all factors (#1846) * Mark smoketests and run as part of be deply (#1846) * Update renamed var (#1846) * Switch from named tuple to dataclass (#1846) This is annoying, but pylint in python3.8 was crashing parsing the named tuple. We weren't using any namedtuple-specific features, so I made the type a dataclass just to get pylint to behave. * Add default timout to requests (#1846) * Fix type (#1846) * Fix merge mistake on poetry.lock (#1846) Signed-off-by: dependabot[bot] <support@github.com> Co-authored-by: Jorge Escobar <jorge.e.escobar@omb.eop.gov> Co-authored-by: Jorge Escobar <83969469+esfoobar-usds@users.noreply.github.com> Co-authored-by: dependabot[bot] <49699333+dependabot[bot]@users.noreply.github.com> Co-authored-by: Matt Bowen <83967628+mattbowen-usds@users.noreply.github.com> Co-authored-by: matt bowen <matthew.r.bowen@omb.eop.gov> * just testing that the boolean is preserved on gha (#1867) * updated with hopefully a fix; coercing aml, fuds, hrs to booleans for the raw value to preserve null character. * Adding tests to ensure proper calculations (#1871) * just testing that the boolean is preserved on gha * checking drop tracts works * adding a check to the agvalue calculation for nri * updated with error messages * tribal tiles fix (#1874) * Alaska tribal points fix (#1821) * tribal tiles fix * disabling child opportunity * lint * removing COI * removing commented out code * Pipeline tile tests (#1864) * temp update * updating with fips check * adding check on pfs * updating with pfs test * Update test_tiles_smoketests.py * Fix lint errors (#1848) * Add column names test (#1848) * Mark tests as smoketests (#1848) * Move to other score-related tests (#1848) * Recast Total threshold criteria exceeded to int (#1848) In writing tests to verify the output of the tiles csv matches the final score CSV, I noticed TC/Total threshold criteria exceeded was getting cast from an int64 to a float64 in the process of PostScoreETL. I tracked it down to the line where we merge the score dataframe with constants.DATA_CENSUS_CSV_FILE_PATH --- there where > 100 tracts in the national census CSV that don't exist in the score, so those ended up with a Total threshhold count of np.nan, which is a float, and thereby cast those columns to float. For the moment I just cast it back. * No need for low memeory (#1848) * Add additional tests of tiles.csv (#1848) * Drop pre-2010 rows before computing score (#1848) Note this is probably NOT the optimal place for this change; it might make more sense for each source to filter its own tracts down to the acceptable tract list. However, that would be a pretty invasive change, where this is central and plenty of other things are happening in score transform that could be moved to sources, so for today, here's where the change will live. * Fix typo (#1848) * Switch from filter to inner join (#1848) * Remove no-op lines from tiles (#1848) * Apply feedback from review, linter (#1848) * Check the values oeverything in the frame (#1848) * Refactor checker class (#1848) * Add test for state names (#1848) * cleanup from reviewing my own code (#1848) * Fix lint error (#1858) * Apply Emma's feedback from review (#1848) * Remove refs to national_df (#1848) * Account for new, fake nullable bools in tiles (#1848) To handle a geojson limitation, Emma converted some nullable boolean colunms to float64 in the tiles export with the values {0.0, 1.0, nan}, giving us the same expressiveness. Sadly, this broke my assumption that all columns between the score and tiles csvs would have the same dtypes, so I need to account for these new, fake bools in my test. * Use equals instead of my worse version (#1848) * Missed a spot where we called _create_score_data (#1848) * Update per safety (#1848) Co-authored-by: matt bowen <matthew.r.bowen@omb.eop.gov> * Add tests to make sure each source makes it to the score correctly (#1878) * Remove unused persistent poverty from score (#1835) * Test a few datasets for overlap in the final score (#1835) * Add remaining data sources (#1853) * Apply code-review feedback (#1835) * Rearrange a little for readabililty (#1835) * Add tract test (#1835) * Add test for score values (#1835) * Check for unmatched source tracts (#1835) * Cleanup numeric code to plaintext (#1835) * Make import more obvious (#1835) * Updating traffic barriers to include low pop threshold (#1889) Changing the traffic barriers to only be included for places with recorded population * Remove no land tracts from map (#1894) remove from map * Issue 1831: missing life expectancy data from Maine and Wisconsin (#1887) * Fixing missing states and adding tests for states to all classes * Removing low pop tracts from FEMA population loss (#1898) dropping 0 population from FEMA * 1831 Follow up (#1902) This code causes no functional change to the code. It does two things: 1. Uses difference instead of - to improve code style for working with sets. 2. Removes the line EXPECTED_MISSING_STATES = ["02", "15"], which is now redundant because of the line I added (in a previous pull request) of ALASKA_AND_HAWAII_EXPECTED_IN_DATA = False. * Add tests for all non-census sources (#1899) * Refactor CDC life-expectancy (1554) * Update to new tract list (#1554) * Adjust for tests (#1848) * Add tests for cdc_places (#1848) * Add EJScreen tests (#1848) * Add tests for HUD housing (#1848) * Add tests for GeoCorr (#1848) * Add persistent poverty tests (#1848) * Update for sources without zips, for new validation (#1848) * Update tests for new multi-CSV but (#1848) Lucas updated the CDC life expectancy data to handle a bug where two states are missing from the US Overall download. Since virtually none of our other ETL classes download multiple CSVs directly like this, it required a pretty invasive new mocking strategy. * Add basic tests for nature deprived (#1848) * Add wildfire tests (#1848) * Add flood risk tests (#1848) * Add DOT travel tests (#1848) * Add historic redlining tests (#1848) * Add tests for ME and WI (#1848) * Update now that validation exists (#1848) * Adjust for validation (#1848) * Add health insurance back to cdc places (#1848) Ooops * Update tests with new field (#1848) * Test for blank tract removal (#1848) * Add tracts for clipping behavior * Test clipping and zfill behavior (#1848) * Fix bad test assumption (#1848) * Simplify class, add test for tract padding (#1848) * Fix percentage inversion, update tests (#1848) Looking through the transformations, I noticed that we were subtracting a percentage that is usually between 0-100 from 1 instead of 100, and so were endind up with some surprising results. Confirmed with lucasmbrown-usds * Add note about first street data (#1848) * Issue 1900: Tribal overlap with Census tracts (#1903) * working notebook * updating notebook * wip * fixing broken tests * adding tribal overlap files * WIP * WIP * WIP, calculated count and names * working * partial cleanup * partial cleanup * updating field names * fixing bug * removing pyogrio * removing unused imports * updating test fixtures to be more realistic * cleaning up notebook * fixing black * fixing flake8 errors * adding tox instructions * updating etl_score * suppressing warning * Use projected CRSes, ignore geom types (#1900) I looked into this a bit, and in general the geometry type mismatch changes very little about the calculation; we have a mix of multipolygons and polygons. The fastest thing to do is just not keep geom type; I did some runs with it set to both True and False, and they're the same within 9 digits of precision. Logically we just want to overlaps, regardless of how the actual geometries are encoded between the frames, so we can in this case ignore the geom types and feel OKAY. I also moved to projected CRSes, since we are actually trying to do area calculations and so like, we should. Again, the change is small in magnitude but logically more sound. * Readd CDC dataset config (#1900) * adding comments to fips code * delete unnecessary loggers Co-authored-by: matt bowen <matthew.r.bowen@omb.eop.gov> * Improve score test documentation based on Lucas's feedback (#1835) (#1914) * Better document base on Lucas's feedback (#1835) * Fix typo (#1835) * Add test to verify GEOJSON matches tiles (#1835) * Remove NOOP line (#1835) * Move GEOJSON generation up for new smoketest (#1835) * Fixup code format (#1835) * Update readme for new somketest (#1835) * Cleanup source tests (#1912) * Move test to base for broader coverage (#1848) * Remove duplicate line (#1848) * FUDS needed an extra mock (#1848) * Add tribal count notebook (#1917) (#1919) * Add tribal count notebook (#1917) * test without caching * added comment Co-authored-by: lucasmbrown-usds <lucas.m.brown@omb.eop.gov> * Add tribal overlap to downloads (#1907) * Add tribal data to downloads (#1904) * Update test pickle with current cols (#1904) * Remove text of tribe names from GeoJSON (#1904) * Update test data (#1904) * Add tribal overlap to smoketests (#1904) * Issue 1910: Do not impute income for 0 population tracts (#1918) * should be working, has unnecessary loggers * removing loggers and cleaning up * updating ejscreen tests * adding tests and responding to PR feedback * fixing broken smoke test * delete smoketest docs * updating click * updating click * Bump just jupyterlab (#1930) * Fixing link checker (#1929) * Update deps safety says are vulnerable (#1937) (#1938) Co-authored-by: matt bowen <matt@mattbowen.net> * Add demos for island areas (#1932) * Backfill population in island areas (#1882) * Update smoketest to account for backfills (#1882) As I wrote in the commend: We backfill island areas with data from the 2010 census, so if THOSE tracts have data beyond the data source, that's to be expected and is fine to pass. If some other state or territory does though, this should fail This ends up being a nice way of documenting that behavior i guess! * Fixup lint issues (#1882) * Add in race demos to 2010 census pull (#1851) * Add backfill data to score (#1851) * Change column name (#1851) * Fill demos after the score (#1851) * Add income back, adjust test (#1882) * Apply code-review feedback (#1851) * Add test for island area backfill (#1851) * Fix bad rename (#1851) * Reorder download fields, add plumbing back (#1942) * Add back lack of plumbing fields (#1920) * Reorder fields for excel (#1921) * Reorder excel fields (#1921) * Fix formating, lint errors, pickes (#1921) * Add missing plumbing col, fix order again (#1921) * Update that pickle (#1921) * refactoring tribal (#1960) * updated with scoring comparison * updated for narhwal -- leaving commented code in for now * pydantic upgrade * produce a string for the front end to ingest (#1963) * wip * i believe this works -- let's see the pipeline * updated fixtures * Adding ADJLI_ET (#1976) * updated tile data * ensuring adjli_et in * Add back income percentile (#1977) * Add missing field to download (#1964) * Remove pydantic since it's unused (#1964) * Add percentile to CSV (#1964) * Update downloadable pickle (#1964) * Issue 105: Configure and run `black` and other pre-commit hooks (clean branch) (#1962) * Configure and run `black` and other pre-commit hooks Co-authored-by: matt bowen <matthew.r.bowen@omb.eop.gov> * Removing fixed python version for black (#1985) * Fixup TA_COUNT and TA_PERC (#1991) * Change TA_PERC, change TA_COUNT (#1988, #1989) - Make TA_PERC_STR back into a nullable float following the rules requestsed in #1989 - Move TA_COUNT to be TA_COUNT_AK, also add a null TA_COUNT_C for CONUS that we can fill in later. * Fix typo comment (#1988) * Issue 1992: Do not impute income for null population tracts (#1993) * Hotfix for DOT data source DNS issue (#1999) * Make tribal overlap set score N (#2004) * Add "Is a Tribal DAC" field (#1998) * Add tribal DACs to score N final (#1998) * Add new fields to downloads (#1998) * Make a int a float (#1998) * Update field names, apply feedback (#1998) * Add assertions around codebook (#2014) * Add assertion around codebook (#1505) * Assert csv and excel have same cols (#1505) * Remove suffixes from tribal lands (#1974) (#2008) * Data source location (#2015) * data source location * toml * cdc_places * cdc_svi_index * url updates * child oppy and dot travel * up to hud_recap * completed ticket * cache bust * hud_recap * us_army_fuds * Remove vars the frontend doesn't use (#2020) (#2022) I did a pretty rough and simple analysis of the variables we put in the tiles and grepped the frontend code to see if (1) they're ever accessed and (2) if they're used, even if they're read once. I removed everything I noticed was not accessed. * Disable file size limits on tiles (#2031) * Disable file size limits on tiles * Remove print debugs I know. * Update file name pattern (#2037) (#2038) * Update file name pattern (#2037) * Remove ETL from generation (2037) I looked more carefully, and this ETL step isn't used in the score, so there's no need to run it every time. Per previous steps, I removed it from constants so the code is there it won't run by default. * Round ALL the float fields for the tiles (#2040) * Round ALL the float fields for the tiles (#2033) * Floor in a simpler way (#2033) Emma pointed out that all teh stuff we're doing in floor_series is probably unnecessary for this case, so just use the built-in floor. * Update pickle I missed (#2033) * Clean commit of just aggregate burden notebook (#1819) added a burden notebook * Update the dockerfile (#2045) * Update so the image builds (#2026) * Fix bad dict (2026) * Rename census tract field in downloads (#2068) * Change tract ID field name (2060) * Update lockfile (#2061) * Bump safety, jupyter, wheel (#2061) * DOn't depend directly on wheel (2061) * Bring narwhal reqs in line with main * Update tribal area counts (#2071) * Rename tribal area field (2062) * Add missing file (#2062) * Add checks to create version (#2047) (#2052) * Fix failing safety (#2114) * Ignore vuln that doesn't affect us 2113 https://nvd.nist.gov/vuln/detail/CVE-2022-42969 landed recently and there's no fix in py (which is maintenance mode). From my analysis, that CVE cannot hurt us (famous last words), so we'll ignore the vuln for now. * 2113 Update our gdal ppa * that didn't work (2113) * Don't add the PPA, the package exists (#2113) * Fix type (#2113) * Force an update of wheel 2113 * Also remove PPA line from create-score-versions * Drop 3.8 because of wheel 2113 * Put back 3.8, use newer actions * Try another way of upgrading wheel 2113 * Upgrade wheel in tox too 2113 * Typo fix 2113 Signed-off-by: dependabot[bot] <support@github.com> Co-authored-by: Emma Nechamkin <97977170+emma-nechamkin@users.noreply.github.com> Co-authored-by: Shelby Switzer <shelby.c.switzer@omb.eop.gov> Co-authored-by: Shelby Switzer <shelbyswitzer@gmail.com> Co-authored-by: Emma Nechamkin <Emma.J.Nechamkin@omb.eop.gov> Co-authored-by: Matt Bowen <83967628+mattbowen-usds@users.noreply.github.com> Co-authored-by: Jorge Escobar <83969469+esfoobar-usds@users.noreply.github.com> Co-authored-by: lucasmbrown-usds <lucas.m.brown@omb.eop.gov> Co-authored-by: Jorge Escobar <jorge.e.escobar@omb.eop.gov> Co-authored-by: dependabot[bot] <49699333+dependabot[bot]@users.noreply.github.com> Co-authored-by: matt bowen <matthew.r.bowen@omb.eop.gov> Co-authored-by: matt bowen <matt@mattbowen.net>
2022-12-01 18:50:54 -08:00
{
"cells": [
{
"cell_type": "code",
"execution_count": 2,
"id": "f0b6f7e2",
"metadata": {},
"outputs": [],
"source": [
"import geopandas as gpd\n",
"import pyogrio\n",
"from data_pipeline.etl.sources.census.etl import CensusETL\n",
"from data_pipeline.etl.sources.tribal.etl import TribalETL\n",
"\n",
"import time\n",
"\n",
"begin = time.time()"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "1e3e65af",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>tribalId</th>\n",
" <th>landAreaName</th>\n",
" <th>Classification</th>\n",
" <th>geometry</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>LAR0001</td>\n",
" <td>Cheyenne River LAR</td>\n",
" <td>1</td>\n",
" <td>MULTIPOLYGON (((-100.49935 45.47125, -100.4993...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>LAR0002</td>\n",
" <td>Crow Creek LAR</td>\n",
" <td>1</td>\n",
" <td>POLYGON ((-99.42137 44.27733, -99.42138 44.273...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td>LAR0003</td>\n",
" <td>Flandreau LAR</td>\n",
" <td>1</td>\n",
" <td>MULTIPOLYGON (((-96.56655 44.08786, -96.57165 ...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td>LAR0004</td>\n",
" <td>Fort Berthold LAR</td>\n",
" <td>1</td>\n",
" <td>POLYGON ((-102.78362 47.99900, -102.78192 47.9...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>LAR0005</td>\n",
" <td>Lake Traverse (Sisseton) LAR</td>\n",
" <td>1</td>\n",
" <td>MULTIPOLYGON (((-97.28946 45.76084, -97.28955 ...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>...</th>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>592</th>\n",
" <td>{0886416F-643E-497E-89D3-E9CC0240158D}</td>\n",
" <td>Chilkat</td>\n",
" <td>None</td>\n",
" <td>POINT (-135.88440 59.40390)</td>\n",
" </tr>\n",
" <tr>\n",
" <th>593</th>\n",
" <td>{2029C35B-86D7-4751-A946-EA0772C81A80}</td>\n",
" <td>Chilkoot</td>\n",
" <td>None</td>\n",
" <td>POINT (-135.44500 59.23580)</td>\n",
" </tr>\n",
" <tr>\n",
" <th>594</th>\n",
" <td>{24DF6536-95CB-4964-94DF-16E440ABCA92}</td>\n",
" <td>Craig</td>\n",
" <td>None</td>\n",
" <td>POINT (-133.14830 55.47640)</td>\n",
" </tr>\n",
" <tr>\n",
" <th>595</th>\n",
" <td>{ACDE097A-9BDA-4FCA-9DB7-297DA6B73F88}</td>\n",
" <td>Douglas</td>\n",
" <td>None</td>\n",
" <td>POINT (-134.41970 58.30190)</td>\n",
" </tr>\n",
" <tr>\n",
" <th>596</th>\n",
" <td>{5E1D1895-FF41-4B11-9EDB-0C1254A360C4}</td>\n",
" <td>Agdaagux</td>\n",
" <td>None</td>\n",
" <td>POINT (-162.31030 55.06170)</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"<p>597 rows × 4 columns</p>\n",
"</div>"
],
"text/plain": [
" tribalId landAreaName \\\n",
"0 LAR0001 Cheyenne River LAR \n",
"1 LAR0002 Crow Creek LAR \n",
"2 LAR0003 Flandreau LAR \n",
"3 LAR0004 Fort Berthold LAR \n",
"4 LAR0005 Lake Traverse (Sisseton) LAR \n",
".. ... ... \n",
"592 {0886416F-643E-497E-89D3-E9CC0240158D} Chilkat \n",
"593 {2029C35B-86D7-4751-A946-EA0772C81A80} Chilkoot \n",
"594 {24DF6536-95CB-4964-94DF-16E440ABCA92} Craig \n",
"595 {ACDE097A-9BDA-4FCA-9DB7-297DA6B73F88} Douglas \n",
"596 {5E1D1895-FF41-4B11-9EDB-0C1254A360C4} Agdaagux \n",
"\n",
" Classification geometry \n",
"0 1 MULTIPOLYGON (((-100.49935 45.47125, -100.4993... \n",
"1 1 POLYGON ((-99.42137 44.27733, -99.42138 44.273... \n",
"2 1 MULTIPOLYGON (((-96.56655 44.08786, -96.57165 ... \n",
"3 1 POLYGON ((-102.78362 47.99900, -102.78192 47.9... \n",
"4 1 MULTIPOLYGON (((-97.28946 45.76084, -97.28955 ... \n",
".. ... ... \n",
"592 None POINT (-135.88440 59.40390) \n",
"593 None POINT (-135.44500 59.23580) \n",
"594 None POINT (-133.14830 55.47640) \n",
"595 None POINT (-134.41970 58.30190) \n",
"596 None POINT (-162.31030 55.06170) \n",
"\n",
"[597 rows x 4 columns]"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Load Tribal geojson\n",
"tribal_gdf = gpd.read_file(\n",
" TribalETL().NATIONAL_TRIBAL_GEOJSON_PATH,\n",
" # Use `pyogrio` because it's vectorized and faster.\n",
" engine=\"pyogrio\",\n",
")\n",
"\n",
"tribal_gdf"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "89fedd44",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>tribalId</th>\n",
" <th>landAreaName</th>\n",
" <th>Classification</th>\n",
" <th>geometry</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>LAR0001</td>\n",
" <td>Cheyenne River LAR</td>\n",
" <td>1</td>\n",
" <td>MULTIPOLYGON (((-100.49935 45.47125, -100.4993...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>LAR0002</td>\n",
" <td>Crow Creek LAR</td>\n",
" <td>1</td>\n",
" <td>POLYGON ((-99.42137 44.27733, -99.42138 44.273...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td>LAR0003</td>\n",
" <td>Flandreau LAR</td>\n",
" <td>1</td>\n",
" <td>MULTIPOLYGON (((-96.56655 44.08786, -96.57165 ...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td>LAR0004</td>\n",
" <td>Fort Berthold LAR</td>\n",
" <td>1</td>\n",
" <td>POLYGON ((-102.78362 47.99900, -102.78192 47.9...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>LAR0005</td>\n",
" <td>Lake Traverse (Sisseton) LAR</td>\n",
" <td>1</td>\n",
" <td>MULTIPOLYGON (((-97.28946 45.76084, -97.28955 ...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>...</th>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>365</th>\n",
" <td>TSA0354</td>\n",
" <td>Seminole TSA</td>\n",
" <td>None</td>\n",
" <td>POLYGON ((-96.49048 34.90423, -96.49146 34.903...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>366</th>\n",
" <td>TSA0355</td>\n",
" <td>Seneca Cayuga TSA</td>\n",
" <td>None</td>\n",
" <td>POLYGON ((-94.61803 36.62531, -94.62083 36.625...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>367</th>\n",
" <td>TSA0356</td>\n",
" <td>Tonkawa TSA</td>\n",
" <td>None</td>\n",
" <td>POLYGON ((-97.24698 36.68082, -97.24697 36.677...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>368</th>\n",
" <td>TSA0357</td>\n",
" <td>Wichita Caddo and Delaware TSA</td>\n",
" <td>None</td>\n",
" <td>POLYGON ((-97.99931 35.36425, -97.99948 35.360...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>369</th>\n",
" <td>TSA0358</td>\n",
" <td>Wyandotte TSA</td>\n",
" <td>None</td>\n",
" <td>POLYGON ((-94.61820 36.82030, -94.61821 36.816...</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"<p>370 rows × 4 columns</p>\n",
"</div>"
],
"text/plain": [
" tribalId landAreaName Classification \\\n",
"0 LAR0001 Cheyenne River LAR 1 \n",
"1 LAR0002 Crow Creek LAR 1 \n",
"2 LAR0003 Flandreau LAR 1 \n",
"3 LAR0004 Fort Berthold LAR 1 \n",
"4 LAR0005 Lake Traverse (Sisseton) LAR 1 \n",
".. ... ... ... \n",
"365 TSA0354 Seminole TSA None \n",
"366 TSA0355 Seneca Cayuga TSA None \n",
"367 TSA0356 Tonkawa TSA None \n",
"368 TSA0357 Wichita Caddo and Delaware TSA None \n",
"369 TSA0358 Wyandotte TSA None \n",
"\n",
" geometry \n",
"0 MULTIPOLYGON (((-100.49935 45.47125, -100.4993... \n",
"1 POLYGON ((-99.42137 44.27733, -99.42138 44.273... \n",
"2 MULTIPOLYGON (((-96.56655 44.08786, -96.57165 ... \n",
"3 POLYGON ((-102.78362 47.99900, -102.78192 47.9... \n",
"4 MULTIPOLYGON (((-97.28946 45.76084, -97.28955 ... \n",
".. ... \n",
"365 POLYGON ((-96.49048 34.90423, -96.49146 34.903... \n",
"366 POLYGON ((-94.61803 36.62531, -94.62083 36.625... \n",
"367 POLYGON ((-97.24698 36.68082, -97.24697 36.677... \n",
"368 POLYGON ((-97.99931 35.36425, -97.99948 35.360... \n",
"369 POLYGON ((-94.61820 36.82030, -94.61821 36.816... \n",
"\n",
"[370 rows x 4 columns]"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Drop the points from the Tribal data (because these cannot be joined to a (Multi)Polygon tract data frame)\n",
"tribal_gdf = tribal_gdf[tribal_gdf.geom_type != \"Point\"]\n",
"tribal_gdf"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "5940556f",
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>STATEFP10</th>\n",
" <th>COUNTYFP10</th>\n",
" <th>TRACTCE10</th>\n",
" <th>GEOID10</th>\n",
" <th>NAME10</th>\n",
" <th>NAMELSAD10</th>\n",
" <th>MTFCC10</th>\n",
" <th>FUNCSTAT10</th>\n",
" <th>ALAND10</th>\n",
" <th>AWATER10</th>\n",
" <th>INTPTLAT10</th>\n",
" <th>INTPTLON10</th>\n",
" <th>geometry</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>20</td>\n",
" <td>071</td>\n",
" <td>958100</td>\n",
" <td>20071958100</td>\n",
" <td>9581</td>\n",
" <td>Census Tract 9581</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>2016176814</td>\n",
" <td>0</td>\n",
" <td>+38.4804076</td>\n",
" <td>-101.8059837</td>\n",
" <td>POLYGON ((-101.79971 38.69806, -101.79097 38.6...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>20</td>\n",
" <td>175</td>\n",
" <td>965600</td>\n",
" <td>20175965600</td>\n",
" <td>9656</td>\n",
" <td>Census Tract 9656</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>1603575701</td>\n",
" <td>2204351</td>\n",
" <td>+37.1805849</td>\n",
" <td>-100.8547406</td>\n",
" <td>POLYGON ((-101.06766 37.20440, -101.06768 37.2...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td>20</td>\n",
" <td>175</td>\n",
" <td>965700</td>\n",
" <td>20175965700</td>\n",
" <td>9657</td>\n",
" <td>Census Tract 9657</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>9466451</td>\n",
" <td>358282</td>\n",
" <td>+37.0625361</td>\n",
" <td>-100.9131437</td>\n",
" <td>POLYGON ((-100.94250 37.06497, -100.94251 37.0...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td>20</td>\n",
" <td>043</td>\n",
" <td>020300</td>\n",
" <td>20043020300</td>\n",
" <td>203</td>\n",
" <td>Census Tract 203</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>211593206</td>\n",
" <td>7045771</td>\n",
" <td>+39.7881238</td>\n",
" <td>-094.9734666</td>\n",
" <td>POLYGON ((-94.95518 39.90129, -94.95475 39.901...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>20</td>\n",
" <td>043</td>\n",
" <td>020200</td>\n",
" <td>20043020200</td>\n",
" <td>202</td>\n",
" <td>Census Tract 202</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>352687026</td>\n",
" <td>2968059</td>\n",
" <td>+39.7540484</td>\n",
" <td>-095.1060098</td>\n",
" <td>POLYGON ((-95.02575 39.88295, -95.02585 39.883...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>...</th>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>74129</th>\n",
" <td>35</td>\n",
" <td>049</td>\n",
" <td>000600</td>\n",
" <td>35049000600</td>\n",
" <td>6</td>\n",
" <td>Census Tract 6</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>1629471</td>\n",
" <td>0</td>\n",
" <td>+35.6758519</td>\n",
" <td>-105.9446097</td>\n",
" <td>POLYGON ((-105.95207 35.67367, -105.95215 35.6...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>74130</th>\n",
" <td>35</td>\n",
" <td>049</td>\n",
" <td>000700</td>\n",
" <td>35049000700</td>\n",
" <td>7</td>\n",
" <td>Census Tract 7</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>1285597</td>\n",
" <td>0</td>\n",
" <td>+35.6802004</td>\n",
" <td>-105.9558818</td>\n",
" <td>POLYGON ((-105.96221 35.67223, -105.96245 35.6...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>74131</th>\n",
" <td>35</td>\n",
" <td>049</td>\n",
" <td>000800</td>\n",
" <td>35049000800</td>\n",
" <td>8</td>\n",
" <td>Census Tract 8</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>1916797</td>\n",
" <td>0</td>\n",
" <td>+35.6805095</td>\n",
" <td>-105.9703558</td>\n",
" <td>POLYGON ((-105.98159 35.67739, -105.98143 35.6...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>74132</th>\n",
" <td>35</td>\n",
" <td>049</td>\n",
" <td>000900</td>\n",
" <td>35049000900</td>\n",
" <td>9</td>\n",
" <td>Census Tract 9</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>2545563</td>\n",
" <td>0</td>\n",
" <td>+35.6692966</td>\n",
" <td>-105.9755351</td>\n",
" <td>POLYGON ((-105.96362 35.67616, -105.96365 35.6...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>74133</th>\n",
" <td>35</td>\n",
" <td>049</td>\n",
" <td>001001</td>\n",
" <td>35049001001</td>\n",
" <td>10.01</td>\n",
" <td>Census Tract 10.01</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>2617281</td>\n",
" <td>0</td>\n",
" <td>+35.6647341</td>\n",
" <td>-105.9468629</td>\n",
" <td>POLYGON ((-105.94510 35.65705, -105.94563 35.6...</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"<p>74134 rows × 13 columns</p>\n",
"</div>"
],
"text/plain": [
" STATEFP10 COUNTYFP10 TRACTCE10 GEOID10 NAME10 NAMELSAD10 \\\n",
"0 20 071 958100 20071958100 9581 Census Tract 9581 \n",
"1 20 175 965600 20175965600 9656 Census Tract 9656 \n",
"2 20 175 965700 20175965700 9657 Census Tract 9657 \n",
"3 20 043 020300 20043020300 203 Census Tract 203 \n",
"4 20 043 020200 20043020200 202 Census Tract 202 \n",
"... ... ... ... ... ... ... \n",
"74129 35 049 000600 35049000600 6 Census Tract 6 \n",
"74130 35 049 000700 35049000700 7 Census Tract 7 \n",
"74131 35 049 000800 35049000800 8 Census Tract 8 \n",
"74132 35 049 000900 35049000900 9 Census Tract 9 \n",
"74133 35 049 001001 35049001001 10.01 Census Tract 10.01 \n",
"\n",
" MTFCC10 FUNCSTAT10 ALAND10 AWATER10 INTPTLAT10 INTPTLON10 \\\n",
"0 G5020 S 2016176814 0 +38.4804076 -101.8059837 \n",
"1 G5020 S 1603575701 2204351 +37.1805849 -100.8547406 \n",
"2 G5020 S 9466451 358282 +37.0625361 -100.9131437 \n",
"3 G5020 S 211593206 7045771 +39.7881238 -094.9734666 \n",
"4 G5020 S 352687026 2968059 +39.7540484 -095.1060098 \n",
"... ... ... ... ... ... ... \n",
"74129 G5020 S 1629471 0 +35.6758519 -105.9446097 \n",
"74130 G5020 S 1285597 0 +35.6802004 -105.9558818 \n",
"74131 G5020 S 1916797 0 +35.6805095 -105.9703558 \n",
"74132 G5020 S 2545563 0 +35.6692966 -105.9755351 \n",
"74133 G5020 S 2617281 0 +35.6647341 -105.9468629 \n",
"\n",
" geometry \n",
"0 POLYGON ((-101.79971 38.69806, -101.79097 38.6... \n",
"1 POLYGON ((-101.06766 37.20440, -101.06768 37.2... \n",
"2 POLYGON ((-100.94250 37.06497, -100.94251 37.0... \n",
"3 POLYGON ((-94.95518 39.90129, -94.95475 39.901... \n",
"4 POLYGON ((-95.02575 39.88295, -95.02585 39.883... \n",
"... ... \n",
"74129 POLYGON ((-105.95207 35.67367, -105.95215 35.6... \n",
"74130 POLYGON ((-105.96221 35.67223, -105.96245 35.6... \n",
"74131 POLYGON ((-105.98159 35.67739, -105.98143 35.6... \n",
"74132 POLYGON ((-105.96362 35.67616, -105.96365 35.6... \n",
"74133 POLYGON ((-105.94510 35.65705, -105.94563 35.6... \n",
"\n",
"[74134 rows x 13 columns]"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Load Census tracts geojson\n",
"census_tract_gdf = gpd.read_file(\n",
" CensusETL.NATIONAL_TRACT_JSON_PATH,\n",
" # Use `pyogrio` because it's vectorized and faster.\n",
" engine=\"pyogrio\",\n",
")\n",
"\n",
"census_tract_gdf"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "595b2a2a",
"metadata": {
"scrolled": true
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/var/folders/lx/xmq8p65j71v9xq2bhsd2j5w40000gp/T/ipykernel_768/2956500515.py:2: UserWarning: Geometry is in a geographic CRS. Results from 'area' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.\n",
"\n",
" census_tract_gdf[\"area_tract\"] = census_tract_gdf.area\n"
]
},
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>STATEFP10</th>\n",
" <th>COUNTYFP10</th>\n",
" <th>TRACTCE10</th>\n",
" <th>GEOID10</th>\n",
" <th>NAME10</th>\n",
" <th>NAMELSAD10</th>\n",
" <th>MTFCC10</th>\n",
" <th>FUNCSTAT10</th>\n",
" <th>ALAND10</th>\n",
" <th>AWATER10</th>\n",
" <th>INTPTLAT10</th>\n",
" <th>INTPTLON10</th>\n",
" <th>geometry</th>\n",
" <th>area_tract</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>20</td>\n",
" <td>071</td>\n",
" <td>958100</td>\n",
" <td>20071958100</td>\n",
" <td>9581</td>\n",
" <td>Census Tract 9581</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>2016176814</td>\n",
" <td>0</td>\n",
" <td>+38.4804076</td>\n",
" <td>-101.8059837</td>\n",
" <td>POLYGON ((-101.79971 38.69806, -101.79097 38.6...</td>\n",
" <td>0.208156</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>20</td>\n",
" <td>175</td>\n",
" <td>965600</td>\n",
" <td>20175965600</td>\n",
" <td>9656</td>\n",
" <td>Census Tract 9656</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>1603575701</td>\n",
" <td>2204351</td>\n",
" <td>+37.1805849</td>\n",
" <td>-100.8547406</td>\n",
" <td>POLYGON ((-101.06766 37.20440, -101.06768 37.2...</td>\n",
" <td>0.162976</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td>20</td>\n",
" <td>175</td>\n",
" <td>965700</td>\n",
" <td>20175965700</td>\n",
" <td>9657</td>\n",
" <td>Census Tract 9657</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>9466451</td>\n",
" <td>358282</td>\n",
" <td>+37.0625361</td>\n",
" <td>-100.9131437</td>\n",
" <td>POLYGON ((-100.94250 37.06497, -100.94251 37.0...</td>\n",
" <td>0.000995</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td>20</td>\n",
" <td>043</td>\n",
" <td>020300</td>\n",
" <td>20043020300</td>\n",
" <td>203</td>\n",
" <td>Census Tract 203</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>211593206</td>\n",
" <td>7045771</td>\n",
" <td>+39.7881238</td>\n",
" <td>-094.9734666</td>\n",
" <td>POLYGON ((-94.95518 39.90129, -94.95475 39.901...</td>\n",
" <td>0.022990</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>20</td>\n",
" <td>043</td>\n",
" <td>020200</td>\n",
" <td>20043020200</td>\n",
" <td>202</td>\n",
" <td>Census Tract 202</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>352687026</td>\n",
" <td>2968059</td>\n",
" <td>+39.7540484</td>\n",
" <td>-095.1060098</td>\n",
" <td>POLYGON ((-95.02575 39.88295, -95.02585 39.883...</td>\n",
" <td>0.037373</td>\n",
" </tr>\n",
" <tr>\n",
" <th>...</th>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>74129</th>\n",
" <td>35</td>\n",
" <td>049</td>\n",
" <td>000600</td>\n",
" <td>35049000600</td>\n",
" <td>6</td>\n",
" <td>Census Tract 6</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>1629471</td>\n",
" <td>0</td>\n",
" <td>+35.6758519</td>\n",
" <td>-105.9446097</td>\n",
" <td>POLYGON ((-105.95207 35.67367, -105.95215 35.6...</td>\n",
" <td>0.000162</td>\n",
" </tr>\n",
" <tr>\n",
" <th>74130</th>\n",
" <td>35</td>\n",
" <td>049</td>\n",
" <td>000700</td>\n",
" <td>35049000700</td>\n",
" <td>7</td>\n",
" <td>Census Tract 7</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>1285597</td>\n",
" <td>0</td>\n",
" <td>+35.6802004</td>\n",
" <td>-105.9558818</td>\n",
" <td>POLYGON ((-105.96221 35.67223, -105.96245 35.6...</td>\n",
" <td>0.000128</td>\n",
" </tr>\n",
" <tr>\n",
" <th>74131</th>\n",
" <td>35</td>\n",
" <td>049</td>\n",
" <td>000800</td>\n",
" <td>35049000800</td>\n",
" <td>8</td>\n",
" <td>Census Tract 8</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>1916797</td>\n",
" <td>0</td>\n",
" <td>+35.6805095</td>\n",
" <td>-105.9703558</td>\n",
" <td>POLYGON ((-105.98159 35.67739, -105.98143 35.6...</td>\n",
" <td>0.000191</td>\n",
" </tr>\n",
" <tr>\n",
" <th>74132</th>\n",
" <td>35</td>\n",
" <td>049</td>\n",
" <td>000900</td>\n",
" <td>35049000900</td>\n",
" <td>9</td>\n",
" <td>Census Tract 9</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>2545563</td>\n",
" <td>0</td>\n",
" <td>+35.6692966</td>\n",
" <td>-105.9755351</td>\n",
" <td>POLYGON ((-105.96362 35.67616, -105.96365 35.6...</td>\n",
" <td>0.000253</td>\n",
" </tr>\n",
" <tr>\n",
" <th>74133</th>\n",
" <td>35</td>\n",
" <td>049</td>\n",
" <td>001001</td>\n",
" <td>35049001001</td>\n",
" <td>10.01</td>\n",
" <td>Census Tract 10.01</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>2617281</td>\n",
" <td>0</td>\n",
" <td>+35.6647341</td>\n",
" <td>-105.9468629</td>\n",
" <td>POLYGON ((-105.94510 35.65705, -105.94563 35.6...</td>\n",
" <td>0.000261</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"<p>74134 rows × 14 columns</p>\n",
"</div>"
],
"text/plain": [
" STATEFP10 COUNTYFP10 TRACTCE10 GEOID10 NAME10 NAMELSAD10 \\\n",
"0 20 071 958100 20071958100 9581 Census Tract 9581 \n",
"1 20 175 965600 20175965600 9656 Census Tract 9656 \n",
"2 20 175 965700 20175965700 9657 Census Tract 9657 \n",
"3 20 043 020300 20043020300 203 Census Tract 203 \n",
"4 20 043 020200 20043020200 202 Census Tract 202 \n",
"... ... ... ... ... ... ... \n",
"74129 35 049 000600 35049000600 6 Census Tract 6 \n",
"74130 35 049 000700 35049000700 7 Census Tract 7 \n",
"74131 35 049 000800 35049000800 8 Census Tract 8 \n",
"74132 35 049 000900 35049000900 9 Census Tract 9 \n",
"74133 35 049 001001 35049001001 10.01 Census Tract 10.01 \n",
"\n",
" MTFCC10 FUNCSTAT10 ALAND10 AWATER10 INTPTLAT10 INTPTLON10 \\\n",
"0 G5020 S 2016176814 0 +38.4804076 -101.8059837 \n",
"1 G5020 S 1603575701 2204351 +37.1805849 -100.8547406 \n",
"2 G5020 S 9466451 358282 +37.0625361 -100.9131437 \n",
"3 G5020 S 211593206 7045771 +39.7881238 -094.9734666 \n",
"4 G5020 S 352687026 2968059 +39.7540484 -095.1060098 \n",
"... ... ... ... ... ... ... \n",
"74129 G5020 S 1629471 0 +35.6758519 -105.9446097 \n",
"74130 G5020 S 1285597 0 +35.6802004 -105.9558818 \n",
"74131 G5020 S 1916797 0 +35.6805095 -105.9703558 \n",
"74132 G5020 S 2545563 0 +35.6692966 -105.9755351 \n",
"74133 G5020 S 2617281 0 +35.6647341 -105.9468629 \n",
"\n",
" geometry area_tract \n",
"0 POLYGON ((-101.79971 38.69806, -101.79097 38.6... 0.208156 \n",
"1 POLYGON ((-101.06766 37.20440, -101.06768 37.2... 0.162976 \n",
"2 POLYGON ((-100.94250 37.06497, -100.94251 37.0... 0.000995 \n",
"3 POLYGON ((-94.95518 39.90129, -94.95475 39.901... 0.022990 \n",
"4 POLYGON ((-95.02575 39.88295, -95.02585 39.883... 0.037373 \n",
"... ... ... \n",
"74129 POLYGON ((-105.95207 35.67367, -105.95215 35.6... 0.000162 \n",
"74130 POLYGON ((-105.96221 35.67223, -105.96245 35.6... 0.000128 \n",
"74131 POLYGON ((-105.98159 35.67739, -105.98143 35.6... 0.000191 \n",
"74132 POLYGON ((-105.96362 35.67616, -105.96365 35.6... 0.000253 \n",
"74133 POLYGON ((-105.94510 35.65705, -105.94563 35.6... 0.000261 \n",
"\n",
"[74134 rows x 14 columns]"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Create a measure for the entire census tract area\n",
"census_tract_gdf[\"area_tract\"] = census_tract_gdf.area\n",
"census_tract_gdf"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "0ea396ed",
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/var/folders/lx/xmq8p65j71v9xq2bhsd2j5w40000gp/T/ipykernel_768/1353983773.py:2: UserWarning: `keep_geom_type=True` in overlay resulted in 1123 dropped geometries of different geometry types than df1 has. Set `keep_geom_type=False` to retain all geometries\n",
" gdf_joined = gpd.overlay(census_tract_gdf, tribal_gdf, how=\"union\")\n"
]
}
],
"source": [
"# Performing overlay funcion\n",
"gdf_joined = gpd.overlay(census_tract_gdf, tribal_gdf, how=\"union\")"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "7fb3ef69",
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/var/folders/lx/xmq8p65j71v9xq2bhsd2j5w40000gp/T/ipykernel_768/2727120487.py:3: UserWarning: Geometry is in a geographic CRS. Results from 'area' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.\n",
"\n",
" gdf_joined['area_joined'] = gdf_joined.area\n"
]
},
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>STATEFP10</th>\n",
" <th>COUNTYFP10</th>\n",
" <th>TRACTCE10</th>\n",
" <th>GEOID10</th>\n",
" <th>NAME10</th>\n",
" <th>NAMELSAD10</th>\n",
" <th>MTFCC10</th>\n",
" <th>FUNCSTAT10</th>\n",
" <th>ALAND10</th>\n",
" <th>AWATER10</th>\n",
" <th>INTPTLAT10</th>\n",
" <th>INTPTLON10</th>\n",
" <th>area_tract</th>\n",
" <th>tribalId</th>\n",
" <th>landAreaName</th>\n",
" <th>Classification</th>\n",
" <th>geometry</th>\n",
" <th>area_joined</th>\n",
" <th>tribal_area_as_a_share_of_tract_area</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>20</td>\n",
" <td>043</td>\n",
" <td>020100</td>\n",
" <td>20043020100</td>\n",
" <td>201</td>\n",
" <td>Census Tract 201</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>454634616.0</td>\n",
" <td>2601186.0</td>\n",
" <td>+39.8206800</td>\n",
" <td>-095.2567279</td>\n",
" <td>0.048098</td>\n",
" <td>LAR0210</td>\n",
" <td>Iowa LAR</td>\n",
" <td>1</td>\n",
" <td>POLYGON ((-95.33994 39.97506, -95.33994 39.975...</td>\n",
" <td>4.998139e-04</td>\n",
" <td>0.010391</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>20</td>\n",
" <td>013</td>\n",
" <td>480600</td>\n",
" <td>20013480600</td>\n",
" <td>4806</td>\n",
" <td>Census Tract 4806</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>882293538.0</td>\n",
" <td>1376818.0</td>\n",
" <td>+39.8596443</td>\n",
" <td>-095.6255187</td>\n",
" <td>0.093019</td>\n",
" <td>LAR0210</td>\n",
" <td>Iowa LAR</td>\n",
" <td>1</td>\n",
" <td>POLYGON ((-95.45656 40.00025, -95.45528 40.000...</td>\n",
" <td>3.209294e-03</td>\n",
" <td>0.034502</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td>31</td>\n",
" <td>147</td>\n",
" <td>964500</td>\n",
" <td>31147964500</td>\n",
" <td>9645</td>\n",
" <td>Census Tract 9645</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>677848509.0</td>\n",
" <td>6076731.0</td>\n",
" <td>+40.1522236</td>\n",
" <td>-095.5858870</td>\n",
" <td>0.072289</td>\n",
" <td>LAR0210</td>\n",
" <td>Iowa LAR</td>\n",
" <td>1</td>\n",
" <td>MULTIPOLYGON (((-95.38162 40.02744, -95.38119 ...</td>\n",
" <td>1.476624e-03</td>\n",
" <td>0.020427</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td>29</td>\n",
" <td>087</td>\n",
" <td>960300</td>\n",
" <td>29087960300</td>\n",
" <td>9603</td>\n",
" <td>Census Tract 9603</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>412869716.0</td>\n",
" <td>6745159.0</td>\n",
" <td>+39.9730230</td>\n",
" <td>-095.1479701</td>\n",
" <td>0.044239</td>\n",
" <td>LAR0210</td>\n",
" <td>Iowa LAR</td>\n",
" <td>1</td>\n",
" <td>POLYGON ((-95.38119 40.02755, -95.38162 40.027...</td>\n",
" <td>1.965514e-07</td>\n",
" <td>0.000004</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>20</td>\n",
" <td>085</td>\n",
" <td>082600</td>\n",
" <td>20085082600</td>\n",
" <td>826</td>\n",
" <td>Census Tract 826</td>\n",
" <td>G5020</td>\n",
" <td>S</td>\n",
" <td>690868809.0</td>\n",
" <td>947758.0</td>\n",
" <td>+39.4553966</td>\n",
" <td>-095.6731404</td>\n",
" <td>0.072404</td>\n",
" <td>LAR0211</td>\n",
" <td>Kickapoo (Kansas) LAR</td>\n",
" <td>1</td>\n",
" <td>POLYGON ((-95.71031 39.65308, -95.69902 39.653...</td>\n",
" <td>5.285627e-06</td>\n",
" <td>0.000073</td>\n",
" </tr>\n",
" <tr>\n",
" <th>...</th>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>76317</th>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>TSA0353</td>\n",
" <td>Sac and Fox TSA</td>\n",
" <td>None</td>\n",
" <td>MULTIPOLYGON (((-96.62002 35.75143, -96.62001 ...</td>\n",
" <td>6.560647e-17</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>76318</th>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>TSA0354</td>\n",
" <td>Seminole TSA</td>\n",
" <td>None</td>\n",
" <td>MULTIPOLYGON (((-96.77536 35.03300, -96.77536 ...</td>\n",
" <td>7.207055e-18</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>76319</th>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>TSA0355</td>\n",
" <td>Seneca Cayuga TSA</td>\n",
" <td>None</td>\n",
" <td>POLYGON ((-94.61836 36.74340, -94.61836 36.743...</td>\n",
" <td>7.016721e-18</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>76320</th>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>TSA0356</td>\n",
" <td>Tonkawa TSA</td>\n",
" <td>None</td>\n",
" <td>MULTIPOLYGON (((-97.24698 36.69942, -97.24692 ...</td>\n",
" <td>2.612218e-17</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>76321</th>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>TSA0358</td>\n",
" <td>Wyandotte TSA</td>\n",
" <td>None</td>\n",
" <td>POLYGON ((-94.61828 36.78970, -94.61834 36.795...</td>\n",
" <td>1.555259e-18</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"<p>76322 rows × 19 columns</p>\n",
"</div>"
],
"text/plain": [
" STATEFP10 COUNTYFP10 TRACTCE10 GEOID10 NAME10 NAMELSAD10 \\\n",
"0 20 043 020100 20043020100 201 Census Tract 201 \n",
"1 20 013 480600 20013480600 4806 Census Tract 4806 \n",
"2 31 147 964500 31147964500 9645 Census Tract 9645 \n",
"3 29 087 960300 29087960300 9603 Census Tract 9603 \n",
"4 20 085 082600 20085082600 826 Census Tract 826 \n",
"... ... ... ... ... ... ... \n",
"76317 NaN NaN NaN NaN NaN NaN \n",
"76318 NaN NaN NaN NaN NaN NaN \n",
"76319 NaN NaN NaN NaN NaN NaN \n",
"76320 NaN NaN NaN NaN NaN NaN \n",
"76321 NaN NaN NaN NaN NaN NaN \n",
"\n",
" MTFCC10 FUNCSTAT10 ALAND10 AWATER10 INTPTLAT10 INTPTLON10 \\\n",
"0 G5020 S 454634616.0 2601186.0 +39.8206800 -095.2567279 \n",
"1 G5020 S 882293538.0 1376818.0 +39.8596443 -095.6255187 \n",
"2 G5020 S 677848509.0 6076731.0 +40.1522236 -095.5858870 \n",
"3 G5020 S 412869716.0 6745159.0 +39.9730230 -095.1479701 \n",
"4 G5020 S 690868809.0 947758.0 +39.4553966 -095.6731404 \n",
"... ... ... ... ... ... ... \n",
"76317 NaN NaN NaN NaN NaN NaN \n",
"76318 NaN NaN NaN NaN NaN NaN \n",
"76319 NaN NaN NaN NaN NaN NaN \n",
"76320 NaN NaN NaN NaN NaN NaN \n",
"76321 NaN NaN NaN NaN NaN NaN \n",
"\n",
" area_tract tribalId landAreaName Classification \\\n",
"0 0.048098 LAR0210 Iowa LAR 1 \n",
"1 0.093019 LAR0210 Iowa LAR 1 \n",
"2 0.072289 LAR0210 Iowa LAR 1 \n",
"3 0.044239 LAR0210 Iowa LAR 1 \n",
"4 0.072404 LAR0211 Kickapoo (Kansas) LAR 1 \n",
"... ... ... ... ... \n",
"76317 NaN TSA0353 Sac and Fox TSA None \n",
"76318 NaN TSA0354 Seminole TSA None \n",
"76319 NaN TSA0355 Seneca Cayuga TSA None \n",
"76320 NaN TSA0356 Tonkawa TSA None \n",
"76321 NaN TSA0358 Wyandotte TSA None \n",
"\n",
" geometry area_joined \\\n",
"0 POLYGON ((-95.33994 39.97506, -95.33994 39.975... 4.998139e-04 \n",
"1 POLYGON ((-95.45656 40.00025, -95.45528 40.000... 3.209294e-03 \n",
"2 MULTIPOLYGON (((-95.38162 40.02744, -95.38119 ... 1.476624e-03 \n",
"3 POLYGON ((-95.38119 40.02755, -95.38162 40.027... 1.965514e-07 \n",
"4 POLYGON ((-95.71031 39.65308, -95.69902 39.653... 5.285627e-06 \n",
"... ... ... \n",
"76317 MULTIPOLYGON (((-96.62002 35.75143, -96.62001 ... 6.560647e-17 \n",
"76318 MULTIPOLYGON (((-96.77536 35.03300, -96.77536 ... 7.207055e-18 \n",
"76319 POLYGON ((-94.61836 36.74340, -94.61836 36.743... 7.016721e-18 \n",
"76320 MULTIPOLYGON (((-97.24698 36.69942, -97.24692 ... 2.612218e-17 \n",
"76321 POLYGON ((-94.61828 36.78970, -94.61834 36.795... 1.555259e-18 \n",
"\n",
" tribal_area_as_a_share_of_tract_area \n",
"0 0.010391 \n",
"1 0.034502 \n",
"2 0.020427 \n",
"3 0.000004 \n",
"4 0.000073 \n",
"... ... \n",
"76317 NaN \n",
"76318 NaN \n",
"76319 NaN \n",
"76320 NaN \n",
"76321 NaN \n",
"\n",
"[76322 rows x 19 columns]"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Calculate overlap\n",
"# Calculating the areas of the newly-created geometries\n",
"gdf_joined[\"area_joined\"] = gdf_joined.area\n",
"\n",
"# Calculating the areas of the newly-created geometries in relation\n",
"# to the original grid cells\n",
"gdf_joined[\"tribal_area_as_a_share_of_tract_area\"] = (\n",
" gdf_joined[\"area_joined\"] / gdf_joined[\"area_tract\"]\n",
")\n",
"gdf_joined"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "042da05e",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" GEOID10 landAreaName tribal_area_as_a_share_of_tract_area\n",
"0 01051030800 Poarch Creek LAR 0.002467\n",
"1 01053970400 Poarch Creek LAR 0.002367\n",
"2 01053970500 Poarch Creek LAR 0.000682\n",
"3 01101005408 Poarch Creek LAR 0.001391\n",
"4 02130000100 Annette Island LAR 0.000038\n",
"... ... ... ...\n",
"2585 56013940300 Wind River LAR 0.204039\n",
"2586 56013940400 Wind River LAR 0.053289\n",
"2587 56017967900 Wind River LAR 0.191189\n",
"2588 56033000600 Crow LAR 0.000565\n",
"2589 56035000102 Wind River LAR 0.000140\n",
"\n",
"[2590 rows x 3 columns]\n"
]
}
],
"source": [
"# Aggregating the results\n",
"results = gdf_joined.groupby([\"GEOID10\", \"landAreaName\"]).agg(\n",
" {\"tribal_area_as_a_share_of_tract_area\": \"sum\"}\n",
")\n",
"\n",
"results = results.reset_index()\n",
"\n",
"results.to_csv(\n",
" \"~/Downloads/tribal_area_as_a_share_of_tract_area.csv\", index=False\n",
")\n",
"\n",
"# Printing results\n",
"print(results)"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "34524a94",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Time taken to execute the function is 140.10310292243958\n"
]
}
],
"source": [
"end = time.time()\n",
"\n",
"print(\"Time taken to execute the ETL is\", end - begin)"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.9.6"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Reference in a new issue Copy permalink