import pyprojclass GPSToXYConverter:def __init__(self, lat0, lon0):"""地图原点经纬度(lat0, lon0)。"""self.projection = pyproj.Proj(f'+proj=tmerc +lat_0={lat0} +lon_0={lon0} +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs')def gps_to_xy(self, lat, lon):"""将GPS坐标转换为地图平面坐标系（米）。"""x, y = self.projection(lon, lat)return x, yif __name__ == '__main__':# 创建GPSToXYConverter对象，指定地图原点经纬度gps_converter = GPSToXYConverter(39.976702, 116.303945)  # 地图原点：北京市东城区王府井大街# 定义测试用GPS坐标lat, lon = 39.958601, 116.299656  # 北京市朝阳区三元桥地铁站# 转换GPS坐标为地图平面坐标x, y = gps_converter.gps_to_xy(lat, lon)# 打印地图平面坐标print(f'({x:.2f}m, {y:.2f}m)')

这个示例代码使用了Proj库的Proj方法，创建了一个GPSToXYConverter类，用于将GPS坐标转换为地图平面坐标系（米）。在初始化GPSToXYConverter对象时，需要传入地图原点经纬度。在gps_to_xy方法中，我们将GPS坐标作为参数传入projection方法中，得到对应的平面坐标值（单位：米）。

在if name == ‘main’:语句块中，我们创建了一个GPSToXYConverter对象，并指定了地图原点为北京市东城区王府井大街的GPS坐标。然后，我们定义了要测试的GPS坐标为北京市朝阳区三元桥地铁站的GPS坐标。接着，我们调用gps_to_xy方法，将这个GPS坐标转换为地图平面坐标，并打印输出结果

(-366.48m, -2009.82m)