python将经纬高转换为enu坐标，并且以第一帧数据为原点

可以使用pyproj库来进行经纬高到enu坐标的转换，以下是示例代码：

import pyproj

# 定义WGS84坐标系
wgs84 = pyproj.CRS('EPSG:4326')

# 定义ENU坐标系
enu = pyproj.CRS(proj='enu', ellps='WGS84', datum='WGS84')

# 定义原点经纬高
origin_lon = 120.0
origin_lat = 30.0
origin_alt = 0.0

# 定义转换器
transformer = pyproj.Transformer.from_crs(wgs84, enu, always_xy=True)
x, y, z = transformer.transform(origin_lon, origin_lat, origin_alt)

# 输出结果
print(f"ENU坐标系下的原点坐标为：({x}, {y}, {z})")

注意，以上代码只是将一个点从经纬高坐标系转换到ENU坐标系，如果需要将一组数据转换，需要对每个点进行单独的转换。

相关问题

python将经纬高转换为enu坐标

可以使用pyproj库中的Transformer类来实现经纬高到ENU坐标的转换，具体代码如下：

import pyproj

# 定义WGS84坐标系
wgs84 = pyproj.CRS('EPSG:4326')

# 定义ENU坐标系
enu = pyproj.CRS('ENU')

# 定义转换器
transformer = pyproj.Transformer.from_crs(wgs84, enu)

# 经纬高坐标
lon,

经纬高坐标系-ECEF坐标系-ENU坐标系 转换 C

以下是经纬高（WGS84）坐标系到ECEF坐标系的转换代码：

#include <stdio.h>
#include <math.h>

#define PI 3.14159265358979323846
#define a 6378137.0        // 地球长半轴
#define b 6356752.314245   // 地球短半轴
#define e ((a*a-b*b)/(a*a)) // 地球椭球体第一偏心率

void WGS84_to_ECEF(double lat, double lon, double alt, double *x, double *y, double *z)
{
    double N = a / sqrt(1 - e*sin(lat)*sin(lat));
    *x = (N + alt) * cos(lat) * cos(lon);
    *y = (N + alt) * cos(lat) * sin(lon);
    *z = (N*(1-e) + alt) * sin(lat);
}

int main()
{
    double lat = 39.909187;
    double lon = 116.397451;
    double alt = 50.0;
    double x, y, z;
    WGS84_to_ECEF(lat*PI/180.0, lon*PI/180.0, alt, &x, &y, &z);
    printf("ECEF坐标系下的坐标为：(%lf, %lf, %lf)\n", x, y, z);
    return 0;
}

以下是ECEF坐标系到ENU坐标系的转换代码：

#include <stdio.h>
#include <math.h>

#define PI 3.14159265358979323846

void ECEF_to_ENU(double x, double y, double z, double lat0, double lon0, double alt0, double *e, double *n, double *u)
{
    double slat0 = sin(lat0);
    double clat0 = cos(lat0);
    double slon0 = sin(lon0);
    double clon0 = cos(lon0);

    double dx = x - (clon0*x + slon0*y);
    double dy = y - (-slon0*x + clon0*y);
    double dz = z - alt0;

    *e = -slon0*dx + clon0*dy;
    *n = -slat0*clon0*dx - slat0*slon0*dy + clat0*dz;
    *u = clat0*clon0*dx + clat0*slon0*dy + slat0*dz;
}

int main()
{
    double x = -2339108.3;
    double y = 5004615.5;
    double z = 3224800.9;
    double lat0 = 39.909187*PI/180.0;
    double lon0 = 116.397451*PI/180.0;
    double alt0 = 50.0;
    double e, n, u;
    ECEF_to_ENU(x, y, z, lat0, lon0, alt0, &e, &n, &u);
    printf("ENU坐标系下的坐标为：(%lf, %lf, %lf)\n", e, n, u);
    return 0;
}