题目：静止卫星遥感图像太阳及卫星天顶、方位角(SOZ/SOA/SAZ/SAA)计算方法
作者：白艺亭
日期：2020.02.25
摘要
遥感图像角度数据（太阳天顶角、太阳方位角、卫星天顶角、卫星方位角）是遥感图像基础数据的重要组成部分，在遥感图像预处理和遥感应用如气溶胶光学厚度反演中是必不可少的参数之一。国产遥感数据经常缺失角度数据支持，大多数产品都不给出成像时对应的逐像元角度数据，需要用户手动计算，而关于角度数据的定义及相应计算方法，网络上很多纷繁复杂的公式，且很多都没有参考来源，本文从遥感专业的角度，在查阅且推导了很多国际英文文献中角度计算公式后，给出四个角度的正确计算公式，为遥感图像的角度数据计算提供支持。
关键词 ： SOZ、SOA、SAZ、SAA计算; 遥感图像
1引言
一般遥感数据分为静止卫星数据和极轨卫星数据。国际上的静止卫星数据大多都带有角度数据（如H8 H9 GK2A ），位于波段数据之后，而例如我国的风云四号静止卫星则没有在0.5-2km空间分辨率数据后设置角度数据，因此如果要对其进行预处理（如大气校正）和气溶胶光学厚度反演时则需要逐像元的角度数据。对于小幅宽（小于250km）的极轨卫星数据，一般在进行预处理和气溶胶光学厚度反演应用时，我们一般选择成像时中心像元的太阳天顶角、太阳方位角、卫星天顶角、卫星方位角代表所有像元的该角度值，而在大幅宽的极轨卫星数据时，由于太阳及卫星观测角度大有不同，且受像元上空的大气空间异质性特征影响，无法利用中心像元代表所有像元，此时则需要计算极轨卫星逐像元角度数据。
2计算方法
网络上有关四个角度数据计算的方法成出不群，但大都没有参考文献来源，且对于方位角的定义千奇百怪，导致很多公式计算下来的结果并不可靠。本文从遥感专业角度，就以上问题于数月前开展相关研究，在查阅大量英文文献后，经推导筛选验证，最终选择确定了可靠的计算方法，现给出其定义及计算公式。如果小伙伴对推导感兴趣，我们可以就文章中的内容进行探讨。
2.1太阳天顶角
太阳天顶角(Solar zenith angle, SOZ)：像元位置处该时入射太阳光线与该像元地平面法线之间的夹角。取值范围：0°- 90°，有的地方写为0-180°是错误的，一定注意。
计算公式[1] [2]为：
【[1]中位于第一章15页，[2]中有详细推导过程，其计算本质上是利用纳皮尔准则Napier Rules在球形三角学中的应用推导出来的），关于纳皮尔准则球形三角学的知识参考[3] [4]仅供阅读了解 .】
式中，
为太阳天顶角，
为像元纬度（Latitude），
为像元太阳赤纬角，
为太阳时角。
太阳赤纬角为该时间太阳直射点与地心的连线 和 赤道平面所成的夹角，每年在-23.45°和23.45°之间变化，3.21 6.22 9.23 12.22分别为北半球一年中的春分，夏至，秋分和冬至日，直射点分别位于赤道，北回归线，赤道，南回归线。太阳赤纬角计算公式[2]（参考文献的第二章56页）为：
AST为像元当地时，一般如果知道成像时间北京时或者其他地方时的话可以根据经度差进行逐像元AST计算，例如为以北京时中心，成像为BT10：30，则计算AST
2.2太阳方位角
太阳方位角(Solar azimuth angle, SOA)：由正北方向起，顺时针转至太阳光线与像元连线在地表的投影所在直线，这个角度称为太阳方位角，取值范围为0° - 360°.
有的地方会定义为正南正北连线以东的角度为正，以西的角度为负，取值范围为-180° - 180°，此种定义也是不合适的，且常见于很多非遥感专业文献中，注意不要采用这种定义，例如此种表达也出现在参考文献【1】中，本文不予采用。
计算公式：
为太阳方位角，
为太阳高度角，取值范围为0-90，太阳入射光线与像元地平面的夹角。其他变量的含义同2.1.
2.3卫星天顶角
卫星天顶角（satellite zenith angle，SAZ）：卫星传感器与像元的连线 和 像元地平面法线所成的夹角，取值范围为0°到90°.
计算公式[6]为：
（图片来自该文章内容截图）
即为卫星天顶角SAZ.
式中，
为纬度，等于Lat.
为卫星星下点所在经度，
为像元经度，等于Lon（Longitude）.
R为地球半径，r为卫星到地心的距离，d为卫星到像元的距离.
为卫星高度角.
2.4卫星方位角
卫星方位角（satellite azimuth angle，SAA）：由正北方向起，顺时针转至卫星与像元连线在地表的投影所在直线，这个角度称为卫星方位角，取值范围为0° - 360°.
有的地方会定义为正南正北连线以东的角度为正，以西的角度为负，取值范围为-180° - 180°，此种定义也是不合适的.
根据纳皮尔准则，计算公式[6]为:
为球面角
Lon为像元经度，等于2.3中的
where if 
otherwise，
原文文章内容截图：
NOTE：
1.以上所有计算均假设地球是个球体，若是把地球当作椭球体，计算公式比较繁琐，但其计算精度会更高；
2.像元所代表地平面视为水平面（遥感当中一般都这么假设），若是视为垂直面或者有一定倾角的面，则上述天顶角和方位角的计算公式不普遍成立；
3.每一个角度求反函数前，若反函数为acos、asin等时，注意被反数组取值范围为[-1,1]，否则就会有错误产生。
3结果
以国产风云四号静止卫星为例，利用上述公式计算结果与官方给出4km的角度数据结果一致。
4总结
本文是作者2019年夏天完成的工作，一直没有总结整理，加上对于很多计算公式的正确与否一直处于判断状态，正直疫情在家时期，将其整理成篇，希望帮助到有需要的人。值得指出的是，在最终求解出来的角度值如有Nan值，则需要检查是否在上一步反函数的被反数组中出现了三角函数最大最小值之外的值，将其修正即可。
创作不易，耗费了整整几天的功夫整理、撰写，如蒙喜欢，麻烦点赞关注一键二联，不准白嫖~
《晴雯歌》红楼梦乐曲
南城有你
参考文献：
[1]Duffie J A, Beckman W A, Blair N. Solar Engineering of Thermal Processes, Photovoltaics and Wind[M]. John Wiley & Sons, 2020.
[2]William B. Stine, Michael Geyer. Power From The Sun[M].NET,2001.
[3]百度百科.球面三角学.
[4]百度文库.球面三角学.
[5]Kalogirou S A. Solar energy engineering: processes and systems[M]. Academic Press, 2013.
[6]Soler T, Eisemann D W. Determination of look angles to geostationary communication satellites[J]. Journal of surveying engineering, 1994, 120(3): 115-127.