File name
Commit message
Commit date
import glob
from pyogrio import read_dataframe, write_dataframe
from tools.spatial_indexed_intersection import geom_overlay
cctv_list_by_region = glob.glob("DATA/refined/geopackage/CCTV_감시영역/*.gpkg")
grid = read_dataframe("DATA/refined/geopackage/100x100/100m격자총인구.gpkg")
for region in cctv_list_by_region:
region_gpkg = read_dataframe(region)
region_gpkg = region_gpkg.to_crs(grid.crs)
result = geom_overlay(region_gpkg, grid)
result = result.dissolve(by="GID").reset_index()
result["CCTV_감시영역"] = result['geometry'].area
write_dataframe(result, f"DATA/refined/geopackage/CCTV_감시영역_격자/{region.split('/')[-1].split(',')[0]}_격자.gpkg")