,6it SrSSKrSSKrSSKrSSKJr SSKrSrSr Sr Sr \ S:XGa\ "5r \ RR5r\ "\ R 5r\S -S 3r\"S 5S- \- r\\S3- R/5(a\rO \"S 5\- r\\S3- r\\S\S3- r\\S\S3- r\\4H7r\R/5(aM\"S\35 \R("S 5 M9 \"S\R;5S\S - SS35 \"S\R<35 \"S\R<35 \R>"\SS9r \R>"\S\3S9r!\ RD(d\!RD(a\"S5 \R("S 5 \ "\ \!5r#\#RI\S\3SS 9 \"S!\RK535 gg!\a(r\"S \35 \R("S 5 SrCGNSrCff=f)"u Calculates total road length per grid cell for a given state. ✅ Features: - Reads {state}_roads_projected.gpkg and {state}_grid_.gpkg - Computes total road length per tile (in meters) - Writes {state}_road_density_.gpkg N)PathGRID_CELL_SIZE_METERScU=(d [R"[5nUc[S[S35e[ U5$![an[SUS[S35UeSnAff=f)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 +/data2/heatmap/scripts/calculate_density.pyresolve_cell_sizers}>")),=">IFGXFYYZ [  9~ !),IJ[I\\] ^  s A A,A''A,cN[S5 UR5RSS0S9n[R"US/USS/SS9n[S [U5S S 35 URRUS '[S5 URS5S R5nUSRU5RS5US 'URS/S9RUR 5n[S5 U$![ an[ SU35UeS nAff=f)z@Calculate total road length per grid cell using in-cell lengths.u5📊 Computing road segments within each grid cell...indexcell_id)columnsgeometry intersection)howz+Failed to compute road/grid intersections: Nu ✅ Created ,z road pieces clipped to cells.length_mu-🧮 Summing clipped lengths per grid cell...ru,✅ Density calculated with in-cell lengths.)print reset_indexrenamegpdoverlay Exception RuntimeErrorlenrlengthgroupbysummapfillnadrop set_indexr)roadsgrid grid_local intersectionsrgrouped grid_results rcalculate_densityr/#s2 AB!!#**GY3G*HJY :,   :. /  L]+A..L MN - 6 6 = =M* 9:##I.z:>>@G' 266w?FFqIJz//9+/6@@LK 89  YHNOUXXYsD D$DD$c[R"SS9nURSSS9 URSSS[S 3S 9 UR 5$) Nz6Calculate road length per grid cell for a given state.) descriptionstatezState name (e.g., alabama).)help cell_size?z.Grid cell size in meters. Defaults to env var z if unset.)nargsr3)argparseArgumentParser add_argumentr parse_args)parsers rr:r:@sh  $ $LF &CD %%6$7zC    __main__u❌ ikmoutputgermanyz_roads_projected.gpkg_grid_z.gpkg_road_density_u❌ Missing input file: u📂 Working on z (z.1fz km grid)u 📥 Roads: u 📥 Grid: roads_projected)layergrid_u1⚠️ One of the input files is empty. Aborting. road_density_GPKG)rEdriveru💾 Saved density results to: )&__doc__r7rsyspathlibr geopandasrr rr/r:__name__argsr2lowerr4r errrexitgrid_size_label germany_direxistsbase_dir roads_path grid_path output_pathpath capitalizename read_filer)r*emptyresultto_fileresolver<rrcsI + :  z E)ug%:;;JeWF?*;5AAIwn_4EUKKKY'{{}} ,TF3 4 HHQK(  U--/09T>#2Fi PQ L) *+ L( )* MM*,= >E ==E/1B*C DD {{djj AB  ud +F NN; o5F&GPVNW +K,?,?,A+B CDU  SEl  s GH %HH