经纬度坐标与高斯坐标的转换代码
/功能说明:将绝对高斯坐标(y,x)转换成绝对的地理坐标(wd,jd)。/ // double y; 输入参数: 高斯坐标的横坐标,以米为单位
// double x; 输入参数: 高斯坐标的纵坐标,以米为单位
// short DH; 输入参数: 带号,表示上述高斯坐标是哪个带的
// double *L; 输出参数: 指向经度坐标的指针,其中经度坐标以秒为单位// double *B; 输出参数: 指向纬度坐标的指针,其中纬度坐标以秒为单位void GaussToGeo(double y, double x, short DH, double *L, double *B, double LP) {
double l0; // 经差
double tf; // tf = tg(Bf0),注意要将Bf转换成以弧度为单位
double nf ; // n = y * sqrt( 1 + etf ** 2) / c, 其中etf = e'**2 * cos(Bf0) ** 2 double t_l0; // l0,经差,以度为单位
double t_B0; // B0,纬度,以度为单位
double Bf0; // Bf0
double etf; // etf,其中etf = e'**2 * cos(Bf0) ** 2
double X_3 ;
double PI=3.14159265358979;
double b_e2=0.0067385254147;
double b_c=6399698.90178271;