开源的Proj4库

Proj4是一个开源的地图投影库，可用于将经纬度坐标转换为平面坐标，支持各种常见的地图投影方式。Proj4库提供了丰富的API，包括C语言、Python、Java、JavaScript等语言的接口，可以轻松实现各种投影转换需求。以下是使用Proj4库在C语言中实现将经纬度坐标转换为平面坐标的示例代码：

#include <proj_api.h>
#include <stdio.h>

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

    // 定义墨卡托投影和经纬度投影
    pj_merc = pj_init_plus("+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs");
    pj_latlong = pj_init_plus("+proj=latlong +ellps=WGS84 +datum=WGS84 +no_defs");

    // 设定经纬度坐标
    lon = 116.397458;   // 经度
    lat = 39.909715;    // 纬度

    // 将经纬度坐标转换为平面坐标
    pj_transform(pj_latlong, pj_merc, 1, 1, &lon, &lat, NULL);
    x = lon;
    y = lat;

    // 输出平面坐标
    printf("x=%f, y=%f\n", x, y);

    // 释放资源
    pj_free(pj_merc);
    pj_free(pj_latlong);

    return 0;
}

在上述示例代码中，我们使用Proj4库定义了墨卡托投影和经纬度投影，然后将经纬度坐标(116.397458, 39.909715)转换为平面坐标(x, y)输出。最后需要记得释放资源，避免内存泄漏。