太阳天顶角和方位角计算

---

计算太阳天顶角方位角，不知道在线计算器算得准不准，今天自己来验证一下。

• 在线计算器：太阳高度角在线计算器
• 参考博文：太阳高度角/方位角计算公式

---

参数

位置：东经115.78388，北纬40.34924
时间：2019/10/08 -- 10:20AM
海拔：481米
时间制式：24h
时区：H(GMT+8:00)
基准时间：太阳时
零方位角：北

真太阳时

先根据北京时间计算得到东经115.78388的平太阳时

10:20-(120-115.78388)*4=10:03:08

再根据10月08号的真太阳时时差为+12分44秒计算得到对应的真太阳时

10:03:08+00:12:44=10:15:52=10.264h

太阳时角

t=(10.264-12)*15=-26.04°

太阳赤纬

这里补充说明一下。

太阳赤纬，是地球赤道平面与太阳和地球中心的连线之间的夹角。赤纬角以年为周期，在 ±23°26′ 的范围内移动，成为季节的标志。
赤纬值日变化很小，一年内任何一天的赤纬角δ可用下式计算：

sinδ=0.39795*cos[0.98563*(N-173)/180*PI]

式中N为日数，自每年1月1日开始计算,δ单位为弧度。

更准确的太阳赤纬计算公式则为：

δ=0.006918-0.399912*cos(b)+0.070257*sin(b)-0.006758*cos(2*b)+0.000907*sin(2*b)-0.002697*cos(3*b)+0.00148*sin(3*b)

其中delta的单位为度(deg)，b单位为弧度。

下面采用简化公式进行计算：

2019年10月08号对应的N=281
sinδ=0.39795*cos[0.98563*(108)/180*PI]≈-0.112678
δ=-6.4697°（这里为了计算方便保留度的单位）

太阳高度角

sinHs=sin(40.34924)*sin(-6.4697)+cos(40.34924)*cos(-6.4697)*cos(-26.04)≈0.6067
Hs≈37.367°

太阳天顶角与高度角互余，那么太阳天顶角就等于：

SZA=90°-Hs≈52.633°

太阳方位角

cosAs=(sinHs*sin(40.3924)-sin(-6.4697))/(cosHs*cos(40.3924))≈0.83574
As≈33.307°

此处的方位角是以南为零方位角的，转换为以北为零方位角：

SAA=180°-As=146.693°

结果对比与总结

相同参数下在线计算器计算结果：

整理至一张表中：

-	SZA	SAA
手算	52.633°	146.693°
在线算	51.49°	147.4°

如果前面的计算没有什么问题，结合在线计算器中给出的偏差，同时考虑到太阳赤纬计算公式的不同，以及计算过程中多次的四舍五入，可以认为两者的偏差是可接受的。
如果后续研究中并不需要很高的精度，可以直接采用在线计算器计算得到，上述计算过程可以为结果提供理论依据。