+6i HSrSSKrSSKrSSKJr SSKrSSKJr SS\S\ S\R4Sjjr \ S :XGaG\R"S S 9r\RS S S9 \R!5r\R$R'5r\"S5S- \- r\\S3- R+5(a\rO \"S5\- r\\S3- r\\S3- r\R+5(d\"S\35 \R4"S5 \ "\S5r\"\5r\R>\:wa,\RA\5r\"S\RC5S35 O \"S\35 \RE\SSS9 \"S\RG535 gg!\ar\"S5 \ "\5rSrCNSrCff=f)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_crspathlayerreturncU(a'[SUSUS35 [R"XS9nO%[SUS35 [R"U5n[S[U5SS35 [S UR35 U$) z,Load a GeoPackage layer into a GeoDataFrame.u📦 Loading layer 'z' from z...)ru 📦 Loading u ✅ Loaded ,z features.zCRS: )printgpd read_filelencrs)rrgdfs #/data2/heatmap/scripts/load_gpkg.py load_gpkgrsy  $UG74&<=mmD. dV3'(mmD! KC|: ./ E#''  J__main__z;Load a location's roads GeoPackage and reproject if needed.) descriptionlocationz/Location name (e.g., alabama, belgium, france).)helpoutputgermanyz _roads.gpkgz_roads_projected.gpkgu❌ File not found: roadsu<⚠️ Could not load 'roads' layer, trying default layer...u🔁 Reprojected to local UTM (z).u✅ Already in target CRS: roads_projectedGPKG)rdriveru💾 Saved projected layer to: )N)$__doc__argparsesyspathlibr geopandasr scripts.crs_utilsrstr GeoDataFramer__name__ArgumentParserparser add_argument parse_argsargsrlower germany_direxistsbase_dir input_path output_pathr exitr Exceptione target_crsrto_crs to_stringto_fileresolverrr;s  / D  0@0@  z  $ $QF  )Z[    D}}""$Hx.9,x7K +..6688>H,xj 44Jz)>??K      $ZL12  $ G,#3'J ww*jj$ / 0D0D0F/GrJK +J<89KK #4VKD +K,?,?,A+B CDU6 $ LN #$s= FF!FF!