使用逐飞库

输入参数：GPS获取的经度维度（直接获得的原始数据）

/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
//@函数：GPS坐标系WGS84转UTM坐标系
//输入参数：GPS获取的经度维度
/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
double angle_rad_umt = 3.1415926535897932 / 180.0;
double utm_longitude,utm_latitude;//函数输出结果，经过转化后的utm坐标点
void LonLat2UTM(double longitude, double latitude)
{
    double lon = longitude;
    double lat = latitude;
    double a = 6378.137;
    double e = 0.0818192;
    double k0 = 0.9996;
    double E0 = 500;
    double N0 = 0;

  
    double zoneNumber = floor(lon/6) + 31;
   
    double lambda0 = (zoneNumber - 1) * 6 - 180 + 3;
    lambda0 = lambda0 * angle_rad_umt;
    
    double phi = lat * angle_rad_umt;
    double lambda = lon * angl