将GPS坐标转换为地图平面坐标

import pyproj

class 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, y

if __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)