（20181109）角度转弧度——dms2rad/rad2dms

因为需要，补充C/C++版本的角度转弧度【参考  测绘程序设计】：

//#include "stdafx.h"
#include <iostream>
using namespace std;

//
//  将度分秒连写的(double型)角度化为弧度值
double dms2rad(double a)
{
	//提取角度值的符号
	double sign = (a<0.0) ? -1.0 : 1.0;
	a = fabs(a);

	//提取角度值的整度
	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 rad2dms(double a)
{
	a = a*206264.806247096363;

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

	int d = (int)(a / 3600.0 + 0.0000001);
	a = a - d*3600.0;

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

	int m = (int)(a / 60.0 + 0.0001);
	a = a - m*60.0;
	if (a<0.0){ m = m - 1; a = a + 60.0; }

	a = d*10000.0 + m*100.0 + a;

	return a*sign;
}

int main()
{
	double A = 330545.35;   //33°5'45.35''
	double B = 455145.25;   //45°51'45.25''
	printf("\n A=%lf", A);
	printf("\n B=%lf", B);

	double AB = dms2rad(A) + dms2rad(B);
	printf("\n A+B=%lf", rad2dms(AB));

	printf("\n");


	system("pause");
	return 0;
}

 

 

=========================================

角度转弧度matlab：

function hd=dms2rad(JD);%角度化弧度函数
%JD=99.3527;
  w=abs(JD);
  d=fix(w);
  df=(w-d)*100;
  f=fix(df);
  m=(df-f)*100;
  ddd=d+f/60+m/3600;
  hd=ddd*pi/180;
  if JD<0;
      hd=-hd;
  end
  return

弧度转角度：

function JD=rad2dms(hd)
%弧度化角度函数
   w=abs(hd);
   df=w*180/pi;
   d=fix(df);
   mf=(df-d)*60;
   m=fix(mf);
   sf=(mf-m)*60;
   JD=d+m/100+sf/10000;
   if hd<0;
       JD=-JD;
   end
return