//提取角度值的整度

int d=(int)((a+0.00001)/10000.0);

a=a-d*10000.0;

if(a<0.0){ d=d-1; a=a+10000;}

//提取角度值的整分及秒值

int m=(int)((a+0.00001)/100.0);

a=a-m*100;

if(a<0.0){ m=m-1; a=a+100.0;}

a=sign*(d*3600.0+m*60.0+a)/206264.806247096363;

return a;

}

// 将角度的弧度值化为度分秒连写的角度（double 型）

double arc_dms(double a)

{

a=a*206264.806247096363;

double sign=(a<0.0) ? -1.0 : 1.0;