crs-l3-tokyo-jgd-crossings

Rerun

import osmnx as ox

def get_tokyo_wards():
    # Query for Tokyo's 23 special wards
    # We want the boundaries and the OSM relation IDs.
    # In OSM, Tokyo's 23 wards are often part of the "Tokyo" administrative area
    # but we need the individual wards.
    # A good way is to search for administrative boundaries in Tokyo.
    
    # Let's try to get the wards via OSMnx/Overpass
    # 'Tokyo, Japan' might be too broad, so we look for administrative level 8 (wards in Japan)
    # within Tokyo.
    
    print("Fetching Tokyo ward boundaries from OSM...")
    try:
        # Using osmnx to get the boundary of Tokyo to use as a search area
        tokyo_boundary = ox.geocode_to_gdf("Tokyo, Japan")
        
        # Now we want the wards (admin_level=8) within this boundary
        # We'll use ox.features_from_polygon
        wards = ox.features_from_polygon(tokyo_boundary.geometry.iloc[0], tags={"boundary": "administrative", "admin_level": "8"})
        
        # Filter for wards that are actually part of the 23 special wards if possible, 
        # or just rely on the admin_level 8 and names.
        # The 23 wards are a specific set.
        
        return wards
    except Exception as e:
        print(f"Error fetching wards: {e}")
        return None

if __name__ == "__main__":
    wards = get_tokyo_wards()
    if wards is not None:
        print(f"Found {len(wards)} wards.")
        print(wards[['name', 'name:en', 'osm_id']].head())
        # Save to check
        wards.to_file("debug_wards.geojson", driver="GeoJSON")
    else:
        print("Failed to fetch wards.")

import osmnx as ox
import geopandas as gpd

def get_tokyo_wards():
    # Instead of fetching everything at once, let's try a more targeted approach.
    # We want admin_level 8 in Tokyo.
    # Let's use a more specific query if possible.
    print("Searching for Tokyo 23 wards...")
    
    # Tokyo's 23 wards are specifically designated.
    # Let's try getting them by name or via administrative level 8.
    # In OSM, Tokyo's special wards are admin_level 8.
    
    # We'll use ox.features_from_place to be more direct.
    try:
        # Try to get the wards directly
        # Note: "Tokyo, Japan" is quite large.
        # Let's refine the search.
        wards = ox.features_from_place("Tokyo, Japan", tags={"boundary": "administrative", "admin_level": "8"})
        
        # Many admin_level 8 features might be municipalities.
        # We need to ensure we only get the 23 special wards if possible.
        # For now, let's see what we get.
        return wards
    except Exception as e:
        print(f"Error: {e}")
        return None

if __name__ == "__main__":
    wards = get_tokyo_wards()
    if wards is not None:
        print(f"Found {len(wards)} features.")
        # print(wards[['name', 'admin_level']].head())
        wards.to_file("wards_raw.geojson", driver="GeoJSON")
    else:
        print("No wards found.")

import osmnx as ox
import pandas as pd

def test_ox():
    print("Testing OSMnx geocode...")
    try:
        gdf = ox.geocode_to_gdf("Tokyo, Japan")
        print(f"Geocode success: {len(gdf)} rows")
        print(gdf.geometry.iloc[0].bounds)
    except Exception as e:
        print(f"Geocode failed: {e}")

    print("\nTesting features_from_address...")
    try:
        gdf = ox.features_from_address("Tokyo, Japan", tags={"boundary": "administrative", "admin_level": "8"})
        print(f"Features success: {len(gdf)} rows")
    except Exception as e:
        print(f"Features failed: {e}")

if __name__ == "__main__":
    test_ox()

import osmnx as ox
import pandas as pd

def test_ox_v2():
    print("Testing ox.features_from_address with dist...")
    try:
        # ox.features_from_address(address, tags, dist=None)
        # Wait, looking at documentation, dist is usually optional but if it's used as a search radius...
        # Actually, let's use features_from_address with a small dist just to test.
        gdf = ox.features_from_address("Tokyo, Japan", tags={"boundary": "administrative", "admin_level": "8"}, dist=1000)
        print(f"Features success: {len(gdf)} rows")
    except Exception as e:
        print(f"Features failed: {e}")

if __name__ == "__main__":
    test_ox_v2()

import osmnx as ox
import geopandas as gpd
import pandas as pd

def get_tokyo_wards():
    print("Fetching Tokyo wards...")
    # Search for admin_level 8 in Tokyo
    # We'll use the bounding box of Tokyo to limit the search if needed, 
    # but features_from_address should work.
    # We'll fetch everything and then filter.
    wards = ox.features_from_address("Tokyo, Japan", tags={"boundary": "administrative", "admin_level": "8"}, dist=50000)
    
    # Filter for Polygons/MultiPolygons
    wards = wards[wards.geometry.type.isin(['Polygon', 'MultiPolygon'])]
    
    # We want to make sure we are in the 23 wards. 
    # The 23 wards are often named like "Shinjuku", "Shibuya", etc. 
    # They are in Tokyo.
    # Let's look at the names to see if we can identify them.
    
    # In Tokyo, the special wards are part of the 'Tokyo' prefecture.
    # We can check if they are within the Tokyo boundary.
    tokyo_boundary = ox.geocode_to_gdf("Tokyo, Japan").iloc[0]
    wards = wards[wards.intersects(tokyo_boundary.geometry)]
    
    # Standardize names
    # osm_id is the relation id in ox.features_from_address
    # Let's rename columns as requested.
    
    # We want: ward_id (osm_id), ward_name_en (name:en), ward_name (name)
    # Note: 'name' might be the native name if 'name:en' is missing.
    
    # Let's clean up the dataframe
    cols_to_keep = ['geometry', 'name', 'name:en', 'osm_id']
    # Check which of these columns exist
    existing_cols = [c for c in cols_to_keep if c in wards.columns]
    wards = wards[existing_cols]
    
    # Rename
    if 'osm_id' in wards.columns:
        wards = wards.rename(columns={'osm_id': 'ward_id'})
    if 'name:en' in wards.columns:
        wards = wards.rename(columns={'name:en': 'ward_name_en'})
    if 'name' in wards.columns:
        wards = wards.rename(columns={'name': 'ward_name'})
        
    # Ensure ward_name_en exists, if not, use ward_name
    if 'ward_name_en' not in wards.columns:
        wards['ward_name_en'] = wards['ward_name']
    
    return wards

if __name__ == "__main__":
    wards = get_tokyo_wards()
    print(f"Found {len(wards)} wards.")
    print(wards[['ward_id', 'ward_name', 'ward_name_en']].head())
    wards.to_file("debug_wards_final.geojson", driver="GeoJSON")