% 度分秒转弧度
function jrad=dtor(deg)
degd=fix(deg);
degf=(deg-degd)*100;
degm=fix(degf);
degs=(degf-degm)/36.0;
jrad=(degd+degm*60.00+degs)*pi/180;