js实现卫星天线方位角和仰角

参数含义：
SATLONG：经度
SATLAT：维度

function CalculateAZ(SATLONG,SATLAT) {var DTOR, RTOD, A, B, F, RSAT, ESHEIGHT, ESLAT, ESLONG, SATLAT, SATLONG, XS, YS, ZS, EE, BETA, RHO;var TXS, TYS, TZS, DIST, AZ, XS2, YS2, ZS2, TXS2, TYS2, TZS2, DIST2, RANGE;DTOR=Math.PI/180; RTOD=180/Math.PI; A=6378.3880;    B=6356.9120;    F=1/297;RSAT=Math.pow((6028.82*((24*60)-4)), 2/3);/*Check(window.document.Form.ESH, "0");Check(window.document.Form.ESLA, "0");Check(window.document.Form.ESLO, "0");Check(window.document.Form.SATLA, "0");Check(window.document.Form.SATLO, "0");*///以北京为观测点ESHEIGHT=0;//高度ESLAT=40;//维度ESLONG=116;//经度/*SATLAT=1.7282723;SATLONG=58.78332;*/XS=RSAT*Math.cos(DTOR*SATLAT)*Math.cos(DTOR*SATLONG);XS2=RSAT*Math.cos(-SATLAT*DTOR)*Math.cos(SATLONG*DTOR);YS=RSAT*Math.cos(DTOR*SATLAT)*Math.sin(DTOR*SATLONG);YS2=RSAT*Math.cos(-SATLAT*DTOR)*Math.sin(SATLONG*DTOR);ZS=RSAT*Math.sin(DTOR*SATLAT);ZS2=-RSAT*Math.sin(DTOR*SATLAT);EE=2*F - Math.pow(F, 2);BETA=Math.atan((1-EE)*Math.tan(ESLAT*DTOR));RHO=A*(1-F)/Math.sqrt(1-(2-F)*F*Math.cos(ESLAT*DTOR)*Math.cos(ESLAT*DTOR));TXS=(XS*Math.cos(ESLONG*DTOR)+YS*Math.sin(ESLONG*DTOR)-RHO*Math.cos(BETA))*Math.sin(ESLAT*DTOR)-(ZS-RHO*Math.sin(BETA))*Math.cos(ESLAT*DTOR);TXS2=(XS2*Math.cos(ESLONG*DTOR)+YS2*Math.sin(ESLONG*DTOR)-RHO*Math.cos(BETA))*Math.sin(ESLAT*DTOR)-(ZS2-RHO*Math.sin(BETA))*Math.cos(ESLAT*DTOR);TYS=-XS*Math.sin(ESLONG*DTOR)+YS*Math.cos(ESLONG*DTOR);TYS2=-XS2*Math.sin(ESLONG*DTOR)+YS2*Math.cos(ESLONG*DTOR);TZS=(XS*Math.cos(ESLONG*DTOR)+YS*Math.sin(ESLONG*DTOR)-RHO*Math.cos(BETA))*Math.cos(ESLAT*DTOR)+(ZS-RHO*Math.sin(BETA))*Math.sin(ESLAT*DTOR)-(ESHEIGHT/1000);TZS2=(XS2*Math.cos(ESLONG*DTOR)+YS2*Math.sin(ESLONG*DTOR)-RHO*Math.cos(BETA))*Math.cos(ESLAT*DTOR)+(ZS2-RHO*Math.sin(BETA))*Math.sin(ESLAT*DTOR)-(ESHEIGHT/1000);DIST=Math.sqrt(TXS*TXS+TYS*TYS+TZS*TZS);DIST2=Math.sqrt(TXS2*TXS2+TYS2*TYS2+TZS2*TZS2)AZ=180-RTOD*ATAN2(TXS,TYS);EL=Math.asin(TZS/DIST)*RTOD;RANGE=Math.max(DIST, DIST2);var m=[];m[0]=round(AZ, 3);//方位角m[1]=round(EL, 3);//仰角//round(RANGE, 3);//最大距离，单位kmreturn m;
}function Check(obRef, def) {if ((obRef.value == null) || (obRef.value == "")) {obRef.value = def;}
}function ATAN2(x, y) {var tmp;if (y == 0) {if (x > 0) {tmp = 0;} else {if (x < 0) {tmp = Math.PI;} else {tmp = 1/0;}}} else {if (y > 0) {tmp = (Math.PI/2) - Math.atan(x/y);} else {tmp = -(Math.PI/2) - Math.atan(x/y);}}return tmp
}function round(val, places) {val = val * Math.pow(10, places)val = Math.round(val)val = val / Math.pow(10, places)return val
}

调用方式:

var am = [];
am=CalculateAZ("经度","维度");
alert("仰角"+am[1]);
alert("方位角"+am[0]);

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！