### 使用 Python 根据经纬度和时间绘制带有底图并进行渲染的航迹图 在轨迹可视化中，除了绘制简单的经纬度点，通常还需要将轨迹叠加在地图底图上，并结合时间信息进行渲染，以增强可视化效果和信息表达能力。可以使用 `folium`、`plotly` 和 `cartopy` 等库实现这一目标。 #### 使用 `folium` 绘制带有底图的轨迹图 `folium` 是基于 Leaflet.js 的 Python 地图可视化库，支持多种地图底图（如 OpenStreetMap、Mapbox、Stamen 等），适合绘制交互式轨迹图： ```python import folium # 假设 data 是包含 latitude、longitude 和 timestamp 的 DataFrame m = folium.Map(location=[data.iloc[0]['latitude'], data.iloc[0]['longitude']], zoom_start=13) # 添加轨迹线 folium.PolyLine( locations=data[['latitude', 'longitude']].values.tolist(), color='blue', weight=2.5, opacity=1 ).add_to(m) # 添加轨迹点标记（可选） for _, row in data.iterrows(): folium.Marker( location=[row['latitude'], row['longitude']], popup=row['timestamp'].strftime('%Y-%m-%d %H:%M:%S') ).add_to(m) m.save("trajectory_with_map.html") ``` 该方法可以在地图上清晰展示轨迹路径和时间信息，适合导航地图中信息展示的需求[^2]。 #### 使用 `plotly` 实现时间维度的轨迹渲染 `plotly` 支持将时间作为动态维度，实现轨迹的动画播放效果： ```python import plotly.express as px fig = px.line_mapbox( data, lat="latitude", lon="longitude", animation_frame="timestamp", mapbox_style="open-street-map", zoom=10, height=600 ) fig.show() ``` 该方法将时间信息与轨迹路径结合，有助于动态展示用户活动区域[^1]。 #### 使用 `cartopy` 进行地理投影和底图叠加 如果需要更专业的地理投影支持，例如高斯-克吕格投影等，可以使用 `cartopy` 库进行地图绘制： ```python import cartopy.crs as ccrs import matplotlib.pyplot as plt ax = plt.axes(projection='mercator') ax.set_extent([data['longitude'].min() - 0.1, data['longitude'].max() + 0.1, data['latitude'].min() - 0.1, data['latitude'].max() + 0.1]) ax.stock_img() plt.plot(data['longitude'], data['latitude'], transform=ccrs.PlateCarree(), color='blue', linewidth=2) plt.scatter(data['longitude'], data['latitude'], transform=ccrs.PlateCarree(), color='red', s=10) plt.show() ``` 该方法支持多种地图投影方式，适合我国标准地形图中常用的高斯-克吕格投影需求[^3]。 #### 使用 `pyecharts` 绘制交互式轨迹图 `pyecharts` 支持动态轨迹展示，并可结合时间维度进行动画渲染： ```python from pyecharts import options as opts from pyecharts.charts import Map, Timeline tl = Timeline() for _, group in data.groupby(data['timestamp'].dt.date): map_chart = ( Map() .add("轨迹点", list(zip(group['longitude'], group['latitude'])), "china") .set_global_opts(title_opts=opts.TitleOpts(title=f"轨迹 {group.iloc[0]['timestamp'].date()}")) ) tl.add(map_chart, group.iloc[0]['timestamp'].date().strftime("%Y-%m-%d")) tl.render("trajectory_timeline.html") ``` 该方法可以为用户提供调节地图视图的功能，如缩放和倾斜，用户可以根据自身需求进行自定义设置。