用C语言写一段墨卡托投影的程序

墨卡托投影是一种地图投影方式，在海洋导航和GIS应用中非常常见，它保持了等角特性，即地球上任意两点间的距离在投影到平面上后仍然保持不变。在C语言中编写墨卡托投影的程序需要对数学和地理坐标有深入理解，并使用适当的数据结构和库。

由于墨卡托投影的具体算法比较复杂，这里提供一个简化的伪代码示例，实际编写时可能需要用到像`proj4`这样的库：

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

// 墨卡托转换函数
void mercator(double lat, double lon, double *x, double *y) {
    double phi = lat * M_PI / 180; // 地球半径和转换因子
    double k = 6378137; // 墨卡托比例系数

    *x = k * lon;
    *y = k * tan(phi);
}

int main() {
    double lon, lat, x, y;
    
    // 输入经度和纬度
    printf("请输入经度（度）: ");
    scanf("%lf", &lon);
    printf("请输入纬度（度）: ");
    scanf("%lf", &lat);

    // 转换到墨卡托坐标
    mercator(lat, lon, &x, &y);

    printf("墨卡托坐标 (x, y): %.2f, %.2f\n", x, y);

    return 0;
}