[发明专利]高光谱卫片全谱段逐像元成像角度参数处理方法及介质

权利要求书

1.一种高光谱卫片全谱段逐像元成像角度参数处理方法，其特征在于，该方法包括：

S100，读取卫星影像，获取卫星影像的每个像元经纬度信息及高程信息；

S200，根据开窗位置、各景数据基准波段成像起始行计数及成像时间计算各波段逐像元成像时间；

S300，根据轨道信息文件通过内插获得对应时间的卫星定位信息；

S400，根据像元经纬度信息、高程信息、成像时间计算全谱段逐像元太阳高度角及方位角；根据像元经纬度信息、高程信息及卫星定位信息计算全谱段逐像元卫星高度角及方位角；

S500，输出太阳及卫星的全谱段逐像元高度角及方位角，存储为固定格式文件。

2.根据权利要求1所述的高光谱卫片全谱段逐像元成像角度参数处理方法，其特征在于，所述S100包括：根据RPC模型参数及DEM文件，获得各波段卫星影像的逐像元的经纬度、高程信息，RPC模型参数由CMOS传感器在焦面上的具体拼接位置等数据信息计算得到。

3.根据权利要求1所述的高光谱卫片全谱段逐像元成像角度参数处理方法，所述S200包括：

将卫星影像的整轨数据拆分为多景数据，并且将每一景数据中的分景时信息记录于分景元数据文件中，分景时信息包括各个波段每景影像的起始行的行计数。

4.根据权利要求3所述的高光谱卫片全谱段逐像元成像角度参数处理方法，所述S200包括：

读取每景图像基准波段起始行的行计数，根据各波段间开窗位置对应关系，依次加或减相邻波段间隔行数，获得各波段起始行计数，具体包括：

通过各波段的时间文件中获得对应行的扫描时间，顺序获取各波段整景5056行的扫描时间time_img具体为

其中，b为波段，line为各波段起始行行计数，pos为开窗位置。

5.根据权利要求4所述的高光谱卫片全谱段逐像元成像角度参数处理方法，所述S300包括：

根据各波段卫星的位置信息X、Y、Z的值，通过匹配查找每景数据的起始行扫描时间和轨道文件记录时间。

6.根据权利要求5所述的高光谱卫片全谱段逐像元成像角度参数处理方法，该方法还包括：

根据起始扫描行的秒数time_sec和下一秒time_gps的GPS位置信息线性插值，确定出700行的轨道位置信息，建立线性模型，计算逐行的卫星位置，其计算公式为

其中其中a表示斜率、b表示截距，X1、Y1、Z1表示扫描时间秒数对应的卫星位置，X2、Y2、Z2表示下一秒种卫星所在位置，time_img表示波段整景5056行的扫描时间。

7.根据权利要求1所述的高光谱卫片全谱段逐像元成像角度参数处理方法，所述计算全谱段逐像元太阳高度角及方位角包括：

计算太阳高度角α，计算方式为其中为当地纬度，δ为赤纬角，

其中Q为日角，n为日期在年内的顺序号，ω为时角，

ω＝15(t_s-12)，

以及，当t_s为真太阳时，其中t为当地平太阳时，L_st为标准时区的经度，L_loc为当地纬度；

计算太阳方位角γ，

其中

日出、日落时刻的太阳方位角计算公式为

波段像元经纬度通过所述S100读取得到，成像时间中的月、日、时、分、秒通过所述S200计算得到。