写出经纬度转平面直角坐标的墨卡托投影C代码

在C语言中，经纬度转换为墨卡托投影坐标的过程可以使用一些基础的数学公式来实现。下面是一个简单的例子，使用经纬度 (经度, 纬度) 转换为 x 和 y 的墨卡托坐标。假设纬度范围在 [-90, 90]，经度在 [-180, 180]：

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

// 经纬度转墨卡托坐标的转换函数
void mercator(double lon, double lat, double *x, double *y) {
    const double R = 6378137; // 地球半径，单位米
    double phi = lat * M_PI / 180; // 将纬度转换为弧度
    double lambda = lon * M_PI / 180; // 将经度转换为弧度

    *x = R * cos(phi) * tan(lambda); // 墨卡托投影x坐标
    *y = R * sin(phi);              // 墨卡托投影y坐标
}

int main() {
    double lat, lon;
    double x, y;

    // 示例输入经纬度
    lat = 51.5074;   // 北纬51.5074度
    lon = -0.1278;  // 西经0.1278度

    // 转换并打印结果
    mercator(lon, lat, &x, &y);
    printf("Mercator coordinates: (%lf, %lf)\n", x, y);

    return 0;
}

注意，这个示例只是一个基础版本，实际应用可能还需要考虑偏移和比例尺调整以适应具体的墨卡托投影需求。同时，这个代码并未包含误差处理和地理区域的边界检查。