已知WGS84海拔高，求EGM96的数值

要将WGS84海拔高转换为EGM96的数值，需要进行大地水准面的转换。大地水准面（Geoid）是指在地球表面处处垂直于重力方向的平面，它与椭球面相比有一定的偏差。EGM96则是一种用于描述大地水准面形状的数学模型，通过计算可以得到某个位置上大地水准面高于WGS84椭球面的高度差（也称为大地高）。

在C语言中，可以通过使用EGM96的高斯正算公式来计算某个位置的大地高。以下是示例代码：

#include <math.h>

#define WGS84_A 6378137.0 // WGS84椭球体长半轴
#define WGS84_F (1.0 / 298.257223563) // WGS84椭球体扁率

double egm96_geoid_height(double lat, double lon, double height) {
    double sin_lat = sin(lat);
    double cos_lat = cos(lat);
    double sin_lon = sin(lon);
    double cos_lon = cos(lon);
    double N = WGS84_A / sqrt(1.0 - WGS84_F * (2.0 - WGS84_F) * sin_lat * sin_lat);
    double x = (N + height) * cos_lat * cos_lon;
    double y = (N + height) * cos_lat * sin_lon;
    double z = (N * (1.0 - WGS84_F * WGS84_F) + height) * sin_lat;
    double p = sqrt(x * x + y * y);
    double theta = atan2(z * WGS84_A, p * WGS84_B);
    double sin_theta = sin(theta);
    double cos_theta = cos(theta);
    double N_prime = WGS84_A / sqrt(1.0 - WGS84_F * (2.0 - WGS84_F) * sin_theta * sin_theta);
    double h_geoid = p * cos_theta + (z + N_prime * WGS84_F * WGS84_F * sin_theta * sin_theta - N * WGS84_F * WGS84_F) * sin_theta * sin_theta / cos_theta;
    return height - h_geoid;
}

这段代码中，`lat`为当前位置的纬度，`lon`为当前位置的经度，`height`为当前位置的WGS84椭球面高。函数返回当前位置的EGM96数值。

需要注意的是，这里使用的是EGM96模型，也可以使用其他的数学模型来描述大地水准面的形状，不同的模型之间可能会有一定的误差。