正轴等角双标准纬线圆锥投影算法代码

以下是正轴等角双标准纬线圆锥投影算法的 C++ 代码：

#include <iostream>
#include <cmath>

using namespace std;

const double PI = acos(-1.0); // 圆周率

// 计算经度差
double deltaLon(double lon1, double lon2) {
    double delta = lon1 - lon2;
    while (delta > PI) {
        delta -= 2 * PI;
    }
    while (delta < -PI) {
        delta += 2 * PI;
    }
    return delta;
}

// 正轴等角双标准纬线圆锥投影算法
void coneProjection(double lon, double lat, double lon0, double lat1, double lat2, double R, double &x, double &y) {
    double n = (sin(lat1) + sin(lat2)) / 2;
    double C = cos(lat1) * cos(lat1) + 2 * n * sin(lat1);
    double p = sqrt(2 * C - 2 * n * sin(lat));
    double theta = n * deltaLon(lon, lon0);

    x = R * p * sin(theta);
    y = R * (p * cos(theta) - sqrt(C - n * sin(lat))) * (lat1 < 0 ? -1 : 1);
}

int main() {
    double lon = 121.47; // 经度
    double lat = 31.23;  // 纬度
    double lon0 = 120;   // 中央经线
    double lat1 = 25;    // 标准纬线1
    double lat2 = 47;    // 标准纬线2
    double R = 6371000;  // 地球半径

    double x, y;
    coneProjection(lon / 180 * PI, lat / 180 * PI, lon0 / 180 * PI, lat1 / 180 * PI, lat2 / 180 * PI, R, x, y);

    cout << "x = " << x << endl;
    cout << "y = " << y << endl;

    return 0;
}

其中，`deltaLon(lon1, lon2)` 函数用于计算经度差，`coneProjection(lon, lat, lon0, lat1, lat2, R, x, y)` 函数为正轴等角双标准纬线圆锥投影算法的实现，`lon` 和 `lat` 分别为要投影的点的经度和纬度，`lon0` 为中央经线，`lat1` 和 `lat2` 分别为两个标准纬线，`R` 为地球半径，`x` 和 `y` 分别为投影后的 x 坐标和 y 坐标。