经纬度转化为xy坐标系C++代码实现

经纬度坐标系和平面直角坐标系（xy坐标系）之间的转换需要用到地图投影算法。常见的地图投影算法有高斯-克吕格投影、墨卡托投影、等角圆柱投影等。下面给出一个简单的例子，使用墨卡托投影将经纬度坐标系转化为xy坐标系的C代码实现：

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

#define PI 3.14159265358979323846
#define EARTH_RADIUS 6378137.0

// 经纬度转墨卡托投影下的平面直角坐标系
void mercator(double lon, double lat, double *x, double *y)
{
    *x = lon * EARTH_RADIUS * PI / 180.0;
    *y = log(tan((90.0 + lat) * PI / 360.0)) * EARTH_RADIUS;
}

int main()
{
    double lon = 116.404; // 经度
    double lat = 39.915; // 纬度
    double x, y;

    mercator(lon, lat, &x, &y);

    printf("x = %lf, y = %lf\n", x, y);

    return 0;
}

在这个例子中，我们使用了墨卡托投影算法，将经度转化为x坐标，将纬度转化为y坐标。其中，EARTH_RADIUS是地球的半径，PI是圆周率。注意，这个例子中的结果是一个近似值，真实的投影坐标可能需要更加复杂的算法。