[发明专利]卫星影像像元观测天顶角和方位角的计算方法

说明书

技术领域

本发明涉及遥感卫星影像观测天顶角和方位角的计算，能够应用在农业、林业、气象、生态环境以及国防军事等遥感部门。

背景技术

观测角度(天顶角和方位角)信息对于遥感数据的处理和应用具有重要的作用。首先，卫星观测角度的不同会导致辐射传输中大气路径的不同，因此观测角度通常是遥感影像大气校正处理中很重要的输入参数。其次，在给定的入射条件下，传感器所接收到的反射辐射强度还受到观测角度的影响，因此通常利用双向反射分布函数(Bidirectional Reflectance Distribution Function，BRDF)来描述光线如何在物体表面进行反射。此外，当具有地形起伏时，观测角度还会对影像的几何特性产生较大的影响。

对于视场角较小的卫星影像，一景影像范围内观测角度的变化较小，数据提供商通常对一景影像提供一个观测天顶角和方位角(如Landsat-5)。对于视场角较大的影像，则需要针对不同的区域提供更多的观测角度(如MODIS，HJ-1A/B)。事实上，如果需要精确地利用卫星的观测角度信息，即使对于Landsat-5之类视场角较小的卫星影像，一景给定一组观测角度也是不够的，尤其是在高纬地区。另外，很多视场角较大的遥感影像(如GF-1宽幅影像)也只对一景影像给定一组观测角度，这样就更不精确了。

鉴于观测角度(天顶角和方位角)信息对于遥感数据的处理和应用的重要的作用，而现有的卫星影像产品所提供的观测角度信息通常不完整，因此研究卫星影像各个像素观测角度的计算方法具有重要的实用价值。

发明内容

本发明的目的在于解决现有技术的不足，提出了一种根据卫星影像几何成像模型计算像素观测天顶角和方位角的方法。只需要知道卫星影像的成像模型，就可以计算出影像中每一个像素的观测天顶角和方位角，在对影像进行几何校正、重投影等几何变换时，可将各像素的观测角度一起进行变换，因此不会因为几何变换而导致影像的观测角度信息丢失。

本发明的技术方案如下：

一种基于卫星影像几何成像模型的像素观测天顶角和方位角计算方法，所述方法中几何成像模型可以是用星历参数描述的严格成像模型，也可以是通用的三维成像模型(如有理函数模型)；包括以下步骤，

步骤1，恢复观测视线：对于影像上任意一点p(像素坐标为(x，y))，取两个不同的高程值h₁和h₂(h₁＜h₂)，利用卫星影像的成像模型得到p点对应的两个空间点P₁(φ₁，λ₁，h₁)和P₂(φ₂，λ₂，h₂)，恢复像平面点p的观测视线

步骤2，计算观测天顶角：连接地心O与P₁点，计算OP₁与P₂P₁的夹角，即为像平面点p的观测天顶角θ_zenith；

步骤3，计算观测方位角初值：计算P₂P₁与过P₁点的子午面的夹角，即为像平面点p观测方位角的初值θ_azimuth；

步骤4，对观测方位角进行南北向修正：如果φ₂＜φ₁，则取θ_azimuth＝180°-θ_azimuth；否则，转到步骤5；

步骤5，对观测方位角进行东西向修正：如果λ₂＜λ₁，则取θ_azimuth＝-θ_azimuth；否则，结束计算。

附图说明

图1是本发明实施例一流程图；

图2是卫星观测天顶角示意图；

图3是OP′₁P₂所在大圆平面示意图；

图4是卫星观测方位角示意图；

图5是平面O′IP₂示意图；

图6是子午面OIP′₁示意图。

具体实施方式

本发明提出了一种根据卫星影像几何成像模型计算像素观测天顶角和方位角的方法。以下结合一个实施例和附图1～6详细说明本发明技术方案。