OSMnx中interpolate_points函数详解及街道细分与图构建实践

import osmnx as ox
from shapely.geometry import LineString

# 示例LineString，从(0,0)到(10,0)
line = LineString([(0, 0), (10, 0)])

# interpolate_points返回一个生成器，每隔2个单位插值
points_generator = ox.utils_geo.interpolate_points(line, 2)

print("生成器对象:", points_generator)

# 将生成器转换为列表以查看所有点
points_list = list(points_generator)
print("插值点列表:", points_list)
# 预期输出: [(0.0, 0.0), (2.0, 0.0), (4.0, 0.0), (6.0, 0.0), (8.0, 0.0), (10.0, 0.0)]

# 安装所需库
!pip install osmnx geopandas momepy networkx matplotlib pandas

import osmnx as ox
import pandas as pd
import geopandas as gpd
import momepy
import networkx as nx
import matplotlib.pyplot as plt
from shapely.geometry import LineString # 用于创建LineString对象

# 定义感兴趣的地点
place_name = "Macon, Macon County, Illinois, United States"

# 生成步行网络图
graph = ox.graph_from_place(place_name, network_type="walk")

# 绘制原始图
fig, ax = ox.plot_graph(graph, bgcolor='w', node_size=0, edge_linewidth=0.5, edge_color='#999999')
plt.title("原始街道网络图")
plt.show()

# 将OSMnx图转换为GeoDataFrames
nodes, edges = ox.graph_to_gdfs(graph)

# 提取边的几何列，并转换为LineString对象列表
geometry_list = edges['geometry'].tolist()

print(f"原始街道线段数量: {len(geometry_list)}")
# print(geometry_list[:3]) # 打印前3个几何对象示例

# 定义插值距离（此处为经纬度单位，实际应用建议先进行投影转换）
# 例如，0.0005经纬度单位在大约40度纬度时，东西方向约为38米，南北方向约为55米。
interpolation_distance = 0.0005

new_subdivided_linestrings = [] # 初始化空列表，用于存储新的细分线段

for original_line in geometry_list:
    # 对每条原始LineString进行插值，获取点生成器
    interpolated_points_generator = ox.utils_geo.interpolate_points(original_line, interpolation_distance)

    # 将生成器转换为点列表
    interpolated_points = list(interpolated_points_generator)

    # 从插值点列表创建新的细分LineString
    # 每个新的LineString由相邻的两个插值点构成
    for i in range(len(interpolated_points)