【不動産情報ライブラリAPI】国土交通省のデータから駅名と座標を抽出し、CSVファイルを作成する方法

pip install geopandas pandas

import geopandas as gpd
import pandas as pd

def extract_station_data(input_file, output_csv):
    """
    駅データから駅名と緯度経度を抽出してCSVに保存する関数
    """
    # GeoJSONファイルを読み込む
    gdf = gpd.read_file(input_file)
    
    # 利用可能な列名とジオメトリタイプを確認
    print("利用可能な列名:", gdf.columns)
    print("ジオメトリタイプ:", gdf.geometry.geom_type.value_counts())
    
    extracted_data = []
    for _, row in gdf.iterrows():
        if row.geometry.geom_type == "Point":
            # Point型の場合は緯度経度を直接取得
            station_name = row["N02_005"]
            lon, lat = row.geometry.x, row.geometry.y
            extracted_data.append({"station_name": station_name, "latitude": lat, "longitude": lon})
        elif row.geometry.geom_type == "MultiLineString":
            # MultiLineString型の場合は重心を使用
            centroid = row.geometry.centroid
            station_name = row["N02_005"]
            lon, lat = centroid.x, centroid.y
            extracted_data.append({"station_name": station_name, "latitude": lat, "longitude": lon})
        else:
            # サポートされていないジオメトリタイプを表示
            print(f"Unsupported geometry type: {row.geometry.geom_type}")
    
    # PandasデータフレームとしてCSVファイルに保存
    df = pd.DataFrame(extracted_data)
    df.to_csv(output_csv, index=False, encoding="utf-8-sig")
    print(f"CSVファイルにデータを保存しました: {output_csv}")

# 入力ファイル（GeoJSON）と出力先（CSV）の指定
input_file = "N02-20_Station.geojson"
output_csv = "station_data.csv"

# 関数を実行
extract_station_data(input_file, output_csv)

station_name,latitude,longitude
天王寺,34.64675,135.51562
新大阪,34.73028,135.49825
大阪,34.70249,135.50017