### 绘制基于UTM坐标系的运动轨迹或覆盖范围图 要在 UTM 坐标系下使用 Python 的 `matplotlib` 或 `geopandas` 绘制运动轨迹或覆盖范围图，可以通过以下方法实现： #### 数据准备 假设已知一组经纬度数据表示运动轨迹点的位置。这些数据需要转换到指定的 UTM 投影坐标系（如 EPSG:32649）。可以利用 `pyproj` 库完成投影变换。 ```python import geopandas as gpd from shapely.geometry import Point import pyproj import pandas as pd import matplotlib.pyplot as plt # 创建原始地理坐标 (经度, 纬度) data = { 'lon': [106.561203, 106.571203, 106.581203], 'lat': [29.558078, 29.568078, 29.578078] } df = pd.DataFrame(data) # 定义 WGS84 和 UTM 投影对象 wgs84 = pyproj.CRS('EPSG:4326') utm_zone_49n = pyproj.CRS('EPSG:32649') project = pyproj.Transformer.from_crs(wgs84, utm_zone_49n, always_xy=True).transform # 将经纬度转换为 UTM 坐标 df['easting'], df['northing'] = project(df['lon'].values, df['lat'].values) # 转换为 GeoDataFrame 并设置 CRS gdf = gpd.GeoDataFrame( df, geometry=gpd.points_from_xy(df.easting, df.northing), crs=utm_zone_49n ) ``` #### 可视化部分 通过 `matplotlib` 和 `geopandas` 结合绘制运动轨迹和覆盖区域。 ```python fig, ax = plt.subplots(figsize=(10, 8)) # 如果有背景地图，则加载并绘制 background_map = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres')) background_map.to_crs(utm_zone_49n).plot(ax=ax, color='lightgray', edgecolor='black', alpha=0.5) # 绘制运动轨迹线 trajectory_line = gpd.GeoSeries([gpd.LineString(list(gdf.geometry))], crs=utm_zone_49n) trajectory_line.plot(ax=ax, color='blue', linewidth=2, label="Trajectory") # 绘制单个点作为起点/终点标记 start_point = gpd.GeoSeries([Point(gdf.iloc[0].geometry.x, gdf.iloc[0].geometry.y)], crs=utm_zone_49n) end_point = gpd.GeoSeries([Point(gdf.iloc[-1].geometry.x, gdf.iloc[-1].geometry.y)], crs=utm_zone_49n) start_point.plot(ax=ax, color='green', markersize=50, marker='o', label="Start Point") end_point.plot(ax=ax, color='red', markersize=50, marker='x', label="End Point") plt.legend() plt.title("Motion Trajectory and Coverage Area in UTM Zone", fontsize=16) plt.xlabel("Easting (m)", fontsize=14) plt.ylabel("Northing (m)", fontsize=14) plt.grid(True) plt.show() ``` 上述代码实现了以下几个功能： - 使用 `pyproj` 进行坐标转换[^1]。 - 利用 `shapely` 构建几何对象，并存储于 `GeoPandas` 中以便进一步操作。 - 配合 `matplotlib` 实现高质量的地图可视化[^2]。 #### 注意事项 - 当前使用的 UTM 区域为 EPSG:32649（WGS84 / UTM zone 49N），实际应用中需根据目标位置调整至合适的 UTM 子区。 - 若涉及较大面积或多子区的情况，可能需要额外处理边界条件。 ---