要使用 Python 的 `matplotlib` 库结合贝塞尔曲线，根据每秒 1 个点且每个点间隔 0.5m 的经纬度信息生成 300m * 2km 四角边缘圆滑、无中间填充且形状接近 4km 常规跑道的轨迹二维图，可以按照以下步骤实现： ```python import math import matplotlib.pyplot as plt from scipy.special import comb import numpy as np # 初始经纬度 start_lat = 30.0 start_lon = 120.0 # 计算经纬度每 0.5m 的变化量 lat_per_m = 1 / (111 * 1000) lon_per_m = 1 / (111 * 1000 * math.cos(math.radians(start_lat))) # 方形区域的长和宽（单位：米） length = 300 width = 2000 # 计算步数 length_steps = int(length / 0.5) width_steps = int(width / 0.5) # 生成方形边缘的经纬度数据 coordinates = [] # 上边 for j in range(width_steps): new_lat = start_lat new_lon = start_lon + j * 0.5 * lon_per_m coordinates.append((new_lat, new_lon)) # 右边 for i in range(length_steps): new_lat = start_lat + i * 0.5 * lat_per_m new_lon = start_lon + (width_steps - 1) * 0.5 * lon_per_m coordinates.append((new_lat, new_lon)) # 下边 for j in range(width_steps - 1, -1, -1): new_lat = start_lat + (length_steps - 1) * 0.5 * lat_per_m new_lon = start_lon + j * 0.5 * lon_per_m coordinates.append((new_lat, new_lon)) # 左边 for i in range(length_steps - 1, 0, -1): new_lat = start_lat + i * 0.5 * lat_per_m new_lon = start_lon coordinates.append((new_lat, new_lon)) # 定义贝塞尔曲线函数 def bezier_curve(points, t): n = len(points) - 1 result = np.zeros(2) for i in range(n + 1): result += comb(n, i) * ((1 - t) ** (n - i)) * (t ** i) * np.array(points[i]) return result # 对四角进行圆滑处理 smoothed_coordinates = [] num_points = 20 # 每个角插入的点数 for i in range(len(coordinates)): if i % length_steps == 0 or i % width_steps == 0: prev_point = coordinates[i - 1] current_point = coordinates[i] next_point = coordinates[(i + 1) % len(coordinates)] for t in np.linspace(0, 1, num_points): smoothed_point = bezier_curve([prev_point, current_point, next_point], t) smoothed_coordinates.append(smoothed_point) else: smoothed_coordinates.append(coordinates[i]) # 分离纬度和经度 lats = [coord[0] for coord in smoothed_coordinates] lons = [coord[1] for coord in smoothed_coordinates] # 绘制二维图片 plt.figure(figsize=(10, 6)) plt.plot(lons, lats, linewidth=2) plt.title('300m x 2000m Square Trajectory with Rounded Corners') plt.xlabel('Longitude') plt.ylabel('Latitude') plt.grid(True) plt.axis('equal') # 显示图形 plt.show() ```