点を面に割り当てる ── 行政区・商圏との突き合わせ

import geopandas as gpd
points = gpd.read_file("stores.geojson")areas = gpd.read_file("admin_areas.geojson")
# 1. まず確認する。None なら定義されていない（要注意）print(f"点の座標系: {points.crs}")print(f"面の座標系: {areas.crs}")
if points.crs is None or areas.crs is None:    raise ValueError("座標系が定義されていません。元データの仕様を確認してください")
# 2. 揃える（面の側に合わせる例）if points.crs != areas.crs:    print(f"揃えます: {points.crs} → {areas.crs}")    points = points.to_crs(areas.crs)
# 3. 揃ったことを確認assert points.crs == areas.crs, "座標系が揃っていません"

import geopandas as gpd

def assign_area(points, areas, area_key="area_name"):    """点を面に割り当て、取りこぼしと重複を確認する。"""    before = len(points)
    joined = gpd.sjoin(points, areas[[area_key, "geometry"]],                       how="left", predicate="within")
    # 1. どこにも割り当てられなかった点    unassigned = joined[joined[area_key].isna()]
    # 2. 複数に割り当てられた点（面が重なっている）    dup_ids = joined.index[joined.index.duplicated()].unique()
    print(f"  入力 {before:,} 件 → 結果 {len(joined):,} 件")    print(f"  割り当てなし: {len(unassigned):,} 件 "          f"({len(unassigned) / before:.1%})")    print(f"  複数に割り当て: {len(dup_ids):,} 件")
    if len(joined) != before:        print("  ※ 件数が変わっています。面の重なりを確認してください")
    return joined, unassigned, dup_ids

import geopandas as gpd

def assign_nearest(unassigned, areas, area_key="area_name",                   max_distance=100.0):    """最も近い面に割り当てる。ただし距離の上限を設ける。
    距離を扱うので、事前に適切な座標系へ変換しておくこと。    """    if unassigned.empty:        return unassigned
    result = gpd.sjoin_nearest(        unassigned, areas[[area_key, "geometry"]],        how="left", max_distance=max_distance,        distance_col="distance_m",    )
    still_none = result[result[area_key].isna()]    print(f"  最寄りに割り当て: {len(result) - len(still_none):,} 件")    print(f"  それでも割り当て不可: {len(still_none):,} 件")    if not result.empty and "distance_m" in result:        print(f"  割り当てた距離の最大: {result['distance_m'].max():.1f} m")
    return result