U OÙi„ ã@s‚dZddlZddlZddlmZdejd<ddlZedœdd„Z e d kr~e ej ƒd krpe d ƒe d ƒe d ¡e ej d ƒdS)a+ Loads extracted buildings GeoPackage and reprojects to meter-based CRS for analysis. This script complements load_gpkg.py by handling buildings data specifically. It reprojects buildings from WGS84 to a meter-based coordinate system (EPSG:5070) suitable for spatial analysis and grid aggregation. éN)ÚPathÚ0Z USE_PYGEOS)Ústatec Csô| ¡}t ¡}|d|}||›d}||›d}| ¡s`td|›ƒtdƒt d¡td|j›dƒzt  t |ƒ¡}td t |ƒd ›d ƒtd |j ›ƒ|j dkrÂtd ƒd|_ |j   ¡dkrîtdƒ|jdd}tdƒntdƒ|j|dddtd|j›ƒtdƒ|jD]8}|dkr$|| ¡ ¡}td|›d|d ›dƒq$tdƒt| ¡ƒWn<tk r°}ztd|›ƒt d¡W5d}~XYnXtdƒtd |›ƒtd!|j›ƒtd"|j›ƒtd#ƒdS)$NZOSM_State_Filesz_buildings.gpkgz_buildings_projected.gpkgu⌠Buildings file not found: z= Run extract_state_buildings.py first to extract buildings.éu!📦 Loading buildings data from z...u ✅ Loaded ú,z buildingsu📋 Current CRS: u/âš ï¸ No CRS found, assuming WGS84 (EPSG:4326)z EPSG:4326iÎu=🔄 Reprojecting to EPSG:5070 (Albers Equal Area for USA)...)Zepsgu✅ Reprojection completeu✅ Already in EPSG:5070Zbuildings_projectedZGPKG)ZlayerZdriveru#💾 Saved projected buildings to: u$ 📋 Building attributes available:Zgeometryz z: z non-null valuesu 🔠Sample buildings data:u⌠Processing failed: u# ✅ Buildings processing complete!u 📠Folder: u 📥 Input: u 💾 Output: uL🔄 Next: Run visualize_all_attribute_grid.py to generate building heatmaps)ÚlowerrÚcwdÚexistsÚprintÚsysÚexitÚnameÚgpdZ read_fileÚstrÚlenZcrsZto_epsgZto_crsZto_fileÚcolumnsZnotnaÚsumÚheadÚ Exception) rZbase_dirZ state_dirZ input_pathZ output_pathZ buildingsZcolZnon_null_countÚe©rú(/data2/heatmap/scripts/load_buildings.pyÚload_buildings_gpkgsN       rÚ__main__éz4Usage: python -m scripts.load_buildings z.Example: python -m scripts.load_buildings ohior)Ú__doc__r ÚosÚpathlibrÚenvironZ geopandasrrrÚ__name__rÚargvr r rrrrÚs  :