下载后可任意编辑测量平差程序设计1
角度(度分秒)到弧度 AngleToRadian#define PI 3
14159265double AngleToRadian(double angle){int D,M;double S,radian,degree, angle,MS;D=int(angle+0
3);MS=angle-D;M=int((MS)*100+0
3);S=(MS*100-M)*100;degree=D+M/60
0+S/3600
0;radian=degree*PI/180
0;return radian;}注意:防止数据溢出,要加个微小量,例如 0
弧度换角度(度分秒) RadianToAngle#define PI 3
14159265double RadianToAngle(double radian){int D,M;double S,radian,degree,MS,angle;degree=radian*180/PI;D=int(degree);MS=degree-D;M=int(MS*60);S=(MS*60-M)*60;angle=D+M/100
0+S/10000
0;return angle;}3
已知两点求坐标方位角 Azimuth#include double Azimuth(double xi,double yi,double xj,double yj){double Dx,Dy,S,T;Dx=xj-xi;Dy=yj-yi;S=sqrt(Dx*Dx+Dy*Dy);if(S