为何转换投影？如何将等时线多边形转回经纬度（lon/lat）？

我正在运行OSMnx等时线示例代码，注意到数据最初以经纬度（如-122.271, 37.871）格式提取，但随后被转换为其他投影格式。想请教两个问题：

1. 该投影转换的必要性是什么？
2. 如何将生成的isochrone_polys转回经纬度（lon/lat）格式？

一、投影转换的必要性 #

• 保证距离计算精准：经纬度是球面坐标，直接用它计算距离会受地球曲率影响产生误差，区域越大误差越明显。转成UTM投影后，坐标单位是米，平面上的距离、长度计算完全符合真实地理尺度，这对等时线依赖的行走时间计算（时间=距离/速度）来说是核心前提。
• 让缓冲区操作有实际意义：代码里对道路和节点做了buffer（比如edge_buff=25），这个操作必须在平面坐标系下进行才有地理意义。如果用经纬度做缓冲区，单位是度，既不符合真实空间逻辑，生成的形状也会因为球面扭曲完全失真。
• 优化可视化效果：UTM投影能让局部区域的地图保持正确的形状和比例，避免经纬度投影带来的拉伸变形，让等时线和道路网络的可视化更贴合真实地理样貌。

二、将isochrone_polys转回经纬度格式的方法 #

借助GeoPandas的坐标转换功能，结合原始图的坐标系就能实现，步骤如下：

1. 确认坐标系：原始未投影的图G用的是WGS84坐标系（EPSG:4326，即经纬度标准），可以通过G.graph['crs']查看具体参数；投影后的图Gb用的是UTM坐标系，参数存在Gb.graph['crs']里。
2. 执行坐标转换：

# 把isochrone_polys封装成带UTM坐标系的GeoDataFrame
gdf_utm = gpd.GeoDataFrame(geometry=isochrone_polys, crs=Gb.graph['crs'])
# 转换为WGS84经纬度坐标系
gdf_lonlat = gdf_utm.to_crs(G.graph['crs'])
# 提取转回经纬度的多边形列表
isochrone_polys_lonlat = gdf_lonlat.geometry.tolist()

3. 验证结果：打印转换后多边形的坐标，确认已经变成（经度，纬度）的格式：

print(isochrone_polys_lonlat[0].exterior.coords[:5])

内容的提问来源于stack exchange，提问作者Ben