PyProj：如何在位置创建自定义NED坐标系？

Question

我有一个机器人的轨迹，以纬度/经度坐标给出，我想使用python将其转换为x / y坐标。 同样，x / y坐标的0点也应由我提供的另一条经/纬度定义（通常只是轨迹的起点）。

所以我一直在尝试使用PyProj库，该库成功将lat / long转换为x / y，但是给出结果的坐标系是英国的坐标系，但是理想情况下我希望0点是我的第一个经纬度组合。因此，尽管也许可以定义一个源自特定坐标系的自定义坐标系，但在Google上找不到任何内容。 因此：有人知道如何在特定的经纬度组合下定义NED坐标系吗？

到目前为止，这是我的代码：

import pyproj as proj

crs_wgs = proj.Proj(init='epsg:4326')  # assuming you're using WGS84 geographic
    crs_bng = proj.Proj(init='epsg:27700')  # use a locally appropriate projected CRS
    # then cast your geographic coordinate pair to the projected system
    x, y = proj.transform(crs_wgs, crs_bng, s['gps_lng'], s['gps_lat'])