mirror of
https://github.com/DOI-DO/j40-cejst-2.git
synced 2025-02-23 18:14:19 -08:00
543d147e61
* initial checkin * gitignore and docker-compose update * readme update and error on hud * encoding issue * one more small README change * data roadmap re-strcuture * pyproject sort * small update to score output folders * checkpoint * couple of last fixes
63 lines
2.2 KiB
Python
63 lines
2.2 KiB
Python
import os
|
|
from pathlib import Path
|
|
import shutil
|
|
|
|
from etl.sources.census.etl_utils import get_state_fips_codes
|
|
|
|
|
|
def generate_tiles(data_path: Path) -> None:
|
|
|
|
# remove existing mbtiles file
|
|
mb_tiles_path = data_path / "tiles" / "block2010.mbtiles"
|
|
if os.path.exists(mb_tiles_path):
|
|
os.remove(mb_tiles_path)
|
|
|
|
# remove existing mvt directory
|
|
mvt_tiles_path = data_path / "tiles" / "mvt"
|
|
if os.path.exists(mvt_tiles_path):
|
|
shutil.rmtree(mvt_tiles_path)
|
|
|
|
# remove existing score json files
|
|
score_geojson_dir = data_path / "score" / "geojson"
|
|
files_in_directory = os.listdir(score_geojson_dir)
|
|
filtered_files = [file for file in files_in_directory if file.endswith(".json")]
|
|
for file in filtered_files:
|
|
path_to_file = os.path.join(score_geojson_dir, file)
|
|
os.remove(path_to_file)
|
|
|
|
# join the state shape sqllite with the score csv
|
|
state_fips_codes = get_state_fips_codes()
|
|
for fips in state_fips_codes:
|
|
cmd = (
|
|
"ogr2ogr -f GeoJSON "
|
|
+ f"-sql \"SELECT * FROM tl_2010_{fips}_bg10 LEFT JOIN 'data/score/csv/data{fips}.csv'.data{fips} ON tl_2010_{fips}_bg10.GEOID10 = data{fips}.ID\" "
|
|
+ f"data/score/geojson/{fips}.json data/census/shp/{fips}/tl_2010_{fips}_bg10.dbf"
|
|
)
|
|
os.system(cmd)
|
|
|
|
# get a list of all json files to plug in the docker commands below
|
|
# (workaround since *.json doesn't seem to work)
|
|
geojson_list = ""
|
|
geojson_path = data_path / "score" / "geojson"
|
|
for file in os.listdir(geojson_path):
|
|
if file.endswith(".json"):
|
|
geojson_list += f"data/score/geojson/{file} "
|
|
|
|
if geojson_list == "":
|
|
logging.error(
|
|
"No GeoJson files found. Please run scripts/download_cbg.py first"
|
|
)
|
|
|
|
# generate mbtiles file
|
|
cmd = (
|
|
"tippecanoe --drop-densest-as-needed -zg -o /home/data/tiles/block2010.mbtiles --extend-zooms-if-still-dropping -l cbg2010 -s_srs EPSG:4269 -t_srs EPSG:4326 "
|
|
+ geojson_list
|
|
)
|
|
os.system(cmd)
|
|
|
|
# generate mvts
|
|
cmd = (
|
|
"tippecanoe --drop-densest-as-needed --no-tile-compression -zg -e /home/data/tiles/mvt "
|
|
+ geojson_list
|
|
)
|
|
os.system(cmd)
|