Matlab的BLH转XYZ的坐标转换

BLH坐标系指的是大地坐标系，其中B表示纬度（latitude），L表示经度（longitude），H表示海拔高度（height）。将大地坐标系转换为空间直角坐标系（也就是XYZ坐标系），需要进行如下的转换：

1. 确定椭球体参数

在转换之前，需要先确定使用的参考椭球体的参数，包括长半轴a、扁率f等。例如，如果使用WGS84椭球体，则a=6378137m，f=1/298.257223563。

2. 计算子午线曲率半径N

根据椭球体参数和纬度B，可以计算出子午线曲率半径N，公式为：

N = a / sqrt(1 - e^2 * sin^2(B))

其中e为椭球体的第一偏心率，可计算为：

e^2 = f * (2 - f)

3. 计算空间直角坐标系中的位置向量

根据公式，可以将BLH坐标系中的位置向量转换为空间直角坐标系中的位置向量：

X = (N + H) * cos(B) * cos(L)
Y = (N + H) * cos(B) * sin(L)
Z = (N * (1 - e^2) + H) * sin(B)

其中，X、Y、Z就是空间直角坐标系中的位置向量，单位为米。