소방용수 접근성 분석, 공백 구역을 지도로 찾기

import geopandas as gpd

# 소방용수·건물을 GeoDataFrame으로 (원본은 보통 경위도 EPSG:4326)
water = gpd.read_file("fire_water.geojson").set_crs(4326, allow_override=True)
bldg  = gpd.read_file("buildings.geojson").set_crs(4326, allow_override=True)

# 거리 계산 전 반드시 미터 단위 투영좌표계로 — 한국 표준 UTM-K
water_m = water.to_crs(5179)
bldg_m  = bldg.to_crs(5179)

# 각 건물에 가장 가까운 소방용수를 붙이고, 그 거리(m)를 함께 계산
joined = gpd.sjoin_nearest(bldg_m, water_m, how="left", distance_col="용수거리m")

print(joined["용수거리m"].describe())          # 분포 — 중앙값·최댓값을 본다
공백 = joined[joined["용수거리m"] > 200]         # 도달 기준 초과(예: 200m)
print(f"공백 건물 {len(공백)}개 / 전체 {len(joined)}개")

covered = water_m.buffer(200).union_all()       # 모든 용수의 200m 도달 범위 합집합
# 관할 경계에서 도달 범위를 빼면 = 공백 구역
gap_area = boundary_m.geometry.union_all().difference(covered)

import folium

m = folium.Map(location=[37.55, 126.9], zoom_start=13, tiles="CartoDB positron")
# 공백 건물(빨강) · 소방용수(파랑)를 4326으로 되돌려 표시
for _, r in 공백.to_crs(4326).iterrows():
    folium.CircleMarker([r.geometry.y, r.geometry.x], radius=3, color="crimson").add_to(m)
m.save("소방용수_공백.html")

- [ ] 거리·버퍼 계산 전에 미터 단위 투영좌표계(EPSG:5179)로 변환했습니다.
- [ ] 소방용수 표준데이터의 좌표 정확성·최신성을 확인했습니다(폐지·이설 반영).
- [ ] 도달 기준 거리를 관할 특성에 맞게 정하고, 여러 단계로 나눠 봤습니다.
- [ ] 결과가 직선거리 기반이라는 한계를 보고서에 명시했습니다.
- [ ] 공백 구역은 현장 확인 우선순위로만 쓰고 단정하지 않았습니다.