c0i (dZddlZddlZddlZddlmZddlZddlZ ddl m Z dZ dZ d dZdZed k(r.eZej$j'Z e ej(Zed ez Zeed z Zeededzdz Zej9sedeej0d edeej:edZej>redeej0d edej@ededz ddeeeZ!e!jEededzddedejGyy#e$r(Zed eej0d YdZ[#dZ[wwxYw)!u Creates a regular grid (square tiles) covering a state's projected road data. ✅ Features: - Reads each state's {state}_roads_projected.gpkg - Generates grid tiles (default 2 km) - Saves as {state}_grid_2km.gpkg - Logs extent and CRS info N)Path)boxGRID_CELL_SIZE_METERSc|xstjt}|tdtd t |S#t$r}td|dtd|d}~wwxYw)z8Resolve cell size from CLI args or environment variable.Nz8Cell size must be passed as an argument or provided via .zInvalid cell size 'z'. Provide an integer or set )osgetenvCELL_SIZE_ENV_VAR ValueErrorint)explicit_value raw_valueexcs scripts/create_grid.pyresolve_cell_sizers}>")),=">IFGXFYYZ [  9~ !),IJ[I\\] ^  s ? A$AA$c t|j\}}}}tj|||z|}tj|||z|}|ddD cgc] }|ddD]} t|| ||z| |z"} }} t j d| i|j } tdt| dd| Scc} }w)zQGenerate a regular grid (square tiles) over the extent of the input GeoDataFrame.Ngeometry)crsu ✅ Created ,z grid cells.) total_boundsnparangergpd GeoDataFramerprintlen) gdf cell_sizexminyminxmaxymaxcolsrowsxypolygonsgrids r make_gridr*&s --D$d 99T4)+Y 7D 99T4)+Y 7Dcr7"AD"IqAq!i-Y7,58" 7   Z2 @D LT1 \ 23 K 7s%B4ctjd}|jdd|jdddtd  |j S) Nz>Create regular grid tiles for a given state's projected roads.) descriptionstatezState name (e.g., alabama).)helpr?z.Grid cell size in meters. Defaults to env var z if unset.)nargsr.)argparseArgumentParser add_argumentr parse_args)parsers rr4r44sh  $ $TF &CD %%6$7zC    __main__u❌ outputz_roads_projected.gpkg_grid_izkm.gpkgu$❌ Projected roads file not found: u📦 Loading roads from roads_projected)layeru⚠️ No features found in u 🗺️ CRS: u📏 Generating grid with z.1fz km cells...)rgrid_kmGPKG)r<driveru💾 Saved grid to: )i)$__doc__r1rsyspathlibr geopandasrnumpyrshapely.geometryrr rr*r4__name__argsr-lowerrr errrexitbase_dir input_path output_pathexists read_fileremptyrr)to_fileresolver6rrUs  +    z