【不動産情報ライブラリAPI】駅名を入力するだけでXYZ座標を自動変換！地図APIを簡単に活用するPythonプログラム

import pandas as pd
from pyproj import Transformer
import math
import requests

# 駅名から対応するXYZ座標を計算する関数
def lonlat_to_xyz(lon, lat, z=15):
    """
    経度・緯度をXYZタイル座標に変換
    """
    transformer = Transformer.from_crs("EPSG:4326", "EPSG:3857")  # WGS84 -> Web Mercator
    x_mercator, y_mercator = transformer.transform(lat, lon)
    
    # タイルの解像度 (Web Mercator の全体サイズを 256px * 2^z タイルで分割)
    tile_size = 256
    initial_resolution = 40075016.68557849 / tile_size  # 地球周囲長 / 初期タイルサイズ
    resolution = initial_resolution / (2**z)
    
    # XYZタイル座標の計算
    x_tile = int((x_mercator + 20037508.342789244) / (tile_size * resolution))
    y_tile = int((20037508.342789244 - y_mercator) / (tile_size * resolution))
    
    return x_tile, y_tile, z

# 駅名を基にAPIを実行する関数
def execute_api(station_name, csv_file):
    """
    駅名を基にXYZ座標を計算し、APIを実行
    """
    # 駅データの読み込み
    df = pd.read_csv(csv_file)
    
    # 駅名で検索
    station_row = df[df["station_name"] == station_name]
    if station_row.empty:
        raise ValueError(f"駅名 '{station_name}' はデータベースに存在しません。")
    
    # 経度・緯度の取得
    lon = station_row["longitude"].values[0]
    lat = station_row["latitude"].values[0]
    
    # XYZ座標を計算
    x, y, z = lonlat_to_xyz(lon, lat)
    
    # APIリクエストの設定
    API_KEY = '****'  # APIキーを設定してください
    url = 'https://www.reinfolib.mlit.go.jp/ex-api/external/XKT013?'
    headers = {"Ocp-Apim-Subscription-Key": API_KEY}
    params = {
        'response_format': 'geojson',
        'z': z,
        'x': x,
        'y': y
    }
    
    # APIリクエスト
    response = requests.get(url, headers=headers, params=params)
    if response.status_code == 200:
        return response.json()
    else:
        raise RuntimeError(f"APIリクエスト失敗: {response.status_code}, {response.text}")

# プログラムの実行例
if __name__ == "__main__":
    station_name = "天王寺"  # 駅名を指定
    csv_file = "station_data.csv"  # 駅データのCSVファイル
    try:
        geojson_data = execute_api(station_name, csv_file)
        print("API実行成功！取得データ:")
        print(geojson_data)
    except Exception as e:
        print(f"エラー: {e}")

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