F1idZddlZddlZddlmZddlZddlmZdede dejfdZ e d k(r ejd Zejd d ej!Zej$j'Zeddz ez Zeedz Zeedz Zej/sedeej2de edZeeZej8ek7r,ej;eZedej=dn edeej?eddedejAyy)uP Loads a location's roads GeoPackage, checks CRS, and reprojects to an appropriate local UTM coordinate system for analysis. ✅ Features: - Works for any location (US states, countries, regions) - Auto-checks for missing files - Chooses a WGS84 UTM zone based on the location's centroid - Saves as {location}_roads_projected.gpkg N)Path)determine_utm_crspathlayerreturnctd|d|dtj||}tdt|ddtd|j|S) z,Load a GeoPackage layer into a GeoDataFrame.u📦 Loading layer 'z' from z...)ru ✅ Loaded ,z features.zCRS: )printgpd read_filelencrs)rrgdfs #/data2/heatmap/scripts/load_gpkg.py load_gpkgrsY wtfC 89 --E *C KC|: ./ E#''  J__main__z;Load a location's roads GeoPackage and reproject if needed.) descriptionlocationz/Location name (e.g., alabama, belgium, france).)helpz..outputz _roads.gpkgz_roads_projected.gpkgu❌ File not found: roadsu🔁 Reprojected to local UTM (z).u✅ Already in target CRS: roads_projectedGPKG)rdriveru💾 Saved projected layer to: )!__doc__argparsesyspathlibr geopandasr scripts.crs_utilsrstr GeoDataFramer__name__ArgumentParserparser add_argument parse_argsargsrlowerbase_dir input_path output_pathexistsr exitr target_crsrto_crs to_stringto_fileresolverrr7s  /D)9)9 z $X $ $QF  )Z[    D}}""$HDzH$x/Hxj 44Jz)>??K     $ZL12  J (C#3'J ww*jj$ / 0D0D0F/GrJK +J<89KK #4VKD +K,?,?,A+B CD=r