I0idZddlZddlZddlmZddlmZddlZde fdZ e dk(rIejd Z e jdd e jZe ej"yy) z Extracts building features from an OpenStreetMap .pbf file. This script complements extract_roads.py by extracting building polygons with attributes like building type, height, and building:levels for heatmap analysis. N)Path)OSMlocationc|j}tj}|dz |z }||dz }|jsht |j d}|r|d}t d|jn.t d|t dtjd||d z }t d |jd  tt|}|jd d idddgd}| |jrt dt dyt dt|dd|jdt dt |j d |j vr0|d j#j%d}t d|d|j vr!|dj'} t d|  j+d}t d  j-|d"d#$t d%|t d't d(|t d)|jt d*|jt d+tdt d,t d-t d.t d/y#t($r.} t d| tjdYd} ~ d} ~ wwxYw#t($r} t d!| Yd} ~ d} ~ wwxYw#t($r.} t d&| tjdYd} ~ d} ~ wwxYw)0Noutputz.osm.pbfz*.pbfru📦 Using PBF file: u❌ No PBF file found in z7 Run extract_roads.py first to download the OSM data.z_buildings.gpkgu 🏢 Extracting buildings from z...buildingTkeepF)r heightzbuilding:levelszaddr:housenumberz addr:streetamenityshopoffice) custom_filter filter_type keep_nodeskeep_relationsextra_attributesu/⚠️ No building data found in this OSM file.z< This location may not have detailed building data in OSM.u✅ Extracted ,z buildings from z data.u📋 Columns available: u🏘️ Top building types: r u📏 Height stats: u ❌ Building extraction failed: i)epsgu&🌐 Reprojected to WGS84 (EPSG:4326).u⚠️ Reprojection failed: buildingsGPKG)layerdriveru#💾 Saved extracted buildings to: u❌ Saving failed: u" ✅ Building extraction complete!u 📁 Folder: u 📥 Input: u 💾 Output: u🏢 Buildings extracted: z<------------------------------------------------------------z Next steps:z91. Run load_buildings.py to reproject buildings to meterszD2. Run visualize_all_attribute_grid.py to generate building heatmaps)lowerrcwdexistslistglobprintnamesysexitrstrget_data_by_custom_criteriaemptylen capitalizecolumns value_countsheaddescribe Exceptionto_crsto_file) rbase_dir location_dirpbf_file pbf_files output_pathosmrbuilding_types height_statses scripts/extract_buildings.pyextract_buildingsr:s&~~HxxzHh&1L (33H ?? **734  |H )(--9 : -l^< = K M HHQKH:_!==K  -hmm_C @A'#h- 33%t,   4 $    C D P Q  s9~a00@ATATAV@WW]^_ (i.?.?)@(ABC ** *&z2??AFFrJN 1.1AB C y(( ($X.779L (7 82$$$$/  67+[H 3K=AB /0 M, () L () M+**+ ,- &s9~a&8 9: (O - EF PQE  045   2 ,QC0112  #A3'(  sPAJ)"B8J)K#9"L) K 2#KK # L,LL L?#L::L?__main__zrLs _R_RD z $X $ $RF  )Z[    Ddmm$rK