本申请公开了一种高光谱卫星的在轨数据快速预处理方法，涉及遥感卫星影像在轨数据预处理技术领域，包括以下步骤：卫星根据任务指令执行目标拍摄任务，获取原始数据和辅助数据；对辅助数据解析，获取时间、轨道、姿态有效辅助数据；结合有效辅助数据和任务点经纬度对所选谱段影像进行目标区域剪裁，获取粗略目标区域数据；读取辐射定标参数，对目标区域数据各谱段进行快速辐射纠正；读取几何定标参数，利用优化线性严密模型最小二乘求解法进行地理坐标快速解算；基于地理坐标和线性内插，生成高精度的附带地理编码影像产品，完成在轨预处理任务，能够满足有限资源的在轨星上处理的高时效需求，实现高光谱影像在轨数据的快速、精确的预处理。
1.一种高光谱卫星的在轨数据快速预处理方法，其特征在于，包括以下步骤：步骤1，卫星接收并解析任务指令，根据指令要求执行拍摄任务，获取原始数据，同时获取星上测量器件输出平台的相关辅助数据；
步骤2，根据辅助数据的编码格式，分别解析出时间、轨道、姿态有效辅助数据，其中：时间数据包括行时和UTC时间；轨道数据是按照UTC时间解析对应的WGS84坐标系下的卫星位置三维坐标；姿态数据是按照UTC时间解析对应的四元数据；
步骤3，根据步骤2中的有效辅助数据和任务点经纬度对所选谱段影像进行目标区域剪裁，获取目标区域数据；
步骤4，根据辐射定标参数，对目标区域数据各谱段进行快速辐射纠正；
步骤5，结合几何定标参数和严密成像几何模型，采用优化线性严密模型最小二乘求解法进行地理坐标快速解算；
步骤6，四角点、中心点地面坐标求解之后，可通过线性内插法计算该剪裁区域任一点坐标，进而生成高精度的附带地理编码影像产品，从而完成在轨预处理任务；
步骤3.1，获取时间数据中的中间行对应的UTC时间；
步骤3.2，根据步骤3.1中的UTC时间，在轨道数据和姿态数据中获取对应UTC时间的卫星位置GNSS坐标 和姿态四元数 ；
步骤3.3，根据式（1），利用姿态四元数获得J2000坐标系下的旋转矩阵R为：式（1）；
其中， 分别为岁差矩阵、章动矩阵、地球自转旋转矩阵、极移矩阵，进而可利用式（2）获得卫星的俯仰pitch、翻滚roll的转角为：式（2）；
其中，为相机载荷焦距， 为探测器每个探元的物理尺寸；
步骤3.5，以目标点为中心画正方形，剪裁有效目标区域，其中，边长可设为的1.5 2.5倍之间；
步骤3.6，同理，将其他选中的谱段按照步骤3.5方法裁剪。
2.如权利要求1所述的一种高光谱卫星的在轨数据快速预处理方法，其特征在于，所述步骤4具体包括以下步骤：步骤4.1，建立辐射定标参数的m*n*c的3维查找表，其中，m表示列，列号与影像的列号相对应；n表示行，行号与谱段号相对应；c表示表号，与积分级数组号相对应；
此外，辐射定标参数包括偏移参数和增益参数，每个存储单元采用4个Byte存储对应的
2个定标参数，前2个Byte存储偏移、后2个Byte存储增益，每个存储单元可用 表示；
步骤4.2，每次成像指定特定的积分级数，根据积分级数可确定查找表号确认定标参数，进而可以对剪裁区域执行多线程协同快速辐射纠正处理；以局部影像列为基本单元，设置二维线程组，线程组包括nn个线程组，每个线程组包括mm个线程，mm对应影像的列、nn对应谱段，一般可分别设置为列数、谱段数，即为m*n；每个线程按照下式进行计算，获取辐射纠正处理后的结果；第n组中第m个线程的纠正结果 的计算方法如式（4）：式（4）。
3.如权利要求1所述的一种高光谱卫星的在轨数据快速预处理方法，其特征在于，所述步骤5具体包括以下步骤：步骤5.1，考虑几何定标偏移矩阵描述了相机载荷、星敏、GNSS天线的安装及偏移误差，对严密成像几何模型进行简化，如下式（5）所示：式（5）；
其中， 偏移矩阵可从定标参数中获取；、为相机内畸变参数，可从定标参数中获取；x、y为图像像素坐标；为比例系数；
步骤5.2，选取四角点和中心点像素，带入式（5）求解地面坐标，式（5）改写为式（6）：式（6）；
其中，下标1、2、3、4、5分别代表四角点和中心点，式（6）为线性方程，在求解过程中，高程值分别选取最大、最小、平均，因此根据式（6）可建立9个方程，进而通过最小二乘求解最终地面坐标。
4.如权利要求1所述的一种高光谱卫星的在轨数据快速预处理方法，其特征在于，所述步骤6中线性内插按照下式（7）进行计算：式（7）；
其中，x为图像坐标，X为地理坐标，下标i表示待求解、下标e、b表示已知的2个点。
一种高光谱卫星的在轨数据快速预处理方法
技术领域
[0001] 本发明涉及遥感卫星影像在轨数据预处理技术领域，尤其涉及一种基于星上智能处理系统的高光谱数据快速预处理的技术方法。
背景技术
[0002] 卫星遥感在轨处理是在获取拍摄数据后在星上平台上执行计算分析任务，直接下传影像产品或信息情报供用户使用，显著减少传统地面处理模式的大量耗时，提升任务执行的时效性，使得卫星遥感在应急管理、国防军事等高时效要求的领域发挥了更大作用。目前在轨预处理技术方法受限于星上资源、系统设计与时效性约束，一般仅是针对普通光学影像的处理及试验，如单波段或3波段影像的局部处理，在特定任务场景取得了较好的效果。但现有技术方法在面向高光谱影像数据的处理及应用的适配性不佳，高光谱数据信息量更加丰富，在地物识别和目标判读方面具有更好的优势，而更丰富的信息也伴随更大的数据量，对星上处理的资源配置、处理时效和精度提出了巨大挑战。
发明内容
[0003] 本发明要解决的技术问题是针对以上不足，提供一种高光谱卫星的在轨数据快速预处理方法，能够满足有限资源的在轨星上处理的高时效需求，实现了高光谱影像在轨数据的快速、精确的预处理。
[0004] 为解决以上技术问题，本发明采用以下技术方案：
[0005] 一种高光谱卫星的在轨数据快速预处理方法，包括以下步骤：
[0006] 步骤1，卫星接收并解析任务指令，根据指令要求执行拍摄任务，获取原始数据，同时获取星上测量器件输出平台的相关辅助数据；
[0007] 步骤2，根据辅助数据的编码格式，分别解析出时间、轨道、姿态有效辅助数据，其中：时间数据包括行时和UTC时间；轨道数据是按照UTC时间解析对应的WGS84坐标系下的卫星位置三维坐标；姿态数据是按照UTC时间解析对应的四元数据；
[0008] 步骤3，根据步骤2中的有效辅助数据和任务点经纬度对所选谱段影像进行目标区域剪裁，获取目标区域数据；
[0009] 步骤4，根据辐射定标参数，对目标区域数据各谱段进行快速辐射纠正；
[0010] 步骤5，结合几何定标参数和严密成像几何模型，采用优化线性严密模型最小二乘求解法进行地理坐标快速解算；
[0011] 步骤6，四角点、中心点地面坐标求解之后，可通过线性内插法计算该剪裁区域任一点坐标，进而生成高精度的附带地理编码影像产品，从而完成在轨预处理任务。
[0012] 进一步的，所述步骤3具体包括以下步骤：
[0013] 步骤3.1，获取时间数据中的中间行对应的UTC时间；
[0014] 步骤3.2，根据步骤3.1中的UTC时间，在轨道数据和姿态数据中获取对应UTC时间的卫星位置GNSS坐标 和姿态四元数 。
[0015] 进一步的，所述步骤3具体还包括以下步骤：
[0016] 步骤3.3，根据式（1），利用姿态四元数获得J2000坐标系下的旋转矩阵R为：
[0017] 式（1）；
[0018] 将R转到本体坐标系下为：
[0019] ；
[0020] 其中， 、 、 、 分别为岁差矩阵、章动矩阵、地球自转旋转矩阵、极移矩阵，进而可利用式（2）获得卫星的俯仰pitch、翻滚roll的转角为：
[0021] 式（2）；
[0022] 其中，atan2为反正切函数。
[0023] 进一步的，所述步骤3具体还包括以下步骤：
[0024] 步骤3.4，获取水平、垂直方向的像素指向偏差；
[0025] 水平、垂直方向的像素指向偏差 为：
[0026] 式（3）；
[0027] 其中， 为相机载荷焦距， 为探测器每个探元的物理尺寸；
[0028] 步骤3.5，以目标点为中心画正方形，剪裁有效目标区域，其中，边长可设为的1.5 2.5倍之间；~
[0029] 步骤3.6，同理，将其他选中的谱段按照步骤3.5方法裁剪。
[0030] 进一步的，所述步骤4具体包括以下步骤：
[0031] 步骤4.1，建立辐射定标参数的m*n*c的3维查找表，其中，m表示列，列号与影像的列号相对应；n表示行，行号与谱段号相对应；c表示表号，与积分级数组号相对应；
[0032] 此外，辐射定标参数包括偏移参数和增益参数，每个存储单元采用4个Byte存储对应的2个定标参数，前2个Byte存储偏移、后2个Byte存储增益，每个存储单元可用 表示；
[0033] 步骤4.2，每次成像指定特定的积分级数，根据积分级数可确定查找表号确认定标参数，进而可以对剪裁区域执行多线程协同快速辐射纠正处理；以局部影像列为基本单元，设置二维线程组，线程组包括nn个线程组，每个线程组包括mm个线程，mm对应影像的列、nn对应谱段，一般可分别设置为列数、谱段数，即为m*n。每个线程按照下式进行计算，获取辐射纠正处理后的结果；第n组中第m个线程的纠正结果 的计算方法如式（4）：
[0034] 式（4）。
[0035] 进一步的，所述步骤5具体包括以下步骤：
[0036] 步骤5.1，考虑几何定标偏移矩阵描述了相机载荷、星敏、GNSS天线的安装及偏移误差，对严密成像几何模型进行简化，如下式（5）所示：
[0037] 式（5）；
[0038] 其中， 偏移矩阵可从定标参数中获取； 为相机内畸变参数，可从定标参数中获取；x、y为图像像素坐标； 为比例系数；
[0039] 步骤5.2，选取四角点和中心点像素，带入式（5）求解地面坐标，式（5）改写为式（6）： 式（6）
[0040] 其中，下标1、2、3、4、5分别代表四角点和中心点，式（6）为线性方程，在求解过程中，高程值分别选取最大、最小、平均，因此根据式（6）可建立9个方程，进而通过最小二乘求解最终地面坐标。
[0041] 进一步的，所述步骤6中线性内插按照下式（7）进行计算：
[0042] 式（7）；
[0043] 其中，x为图像坐标，X为地理坐标，下标i表示待求解、下标e、b表示已知的2个点。
[0044] 本发明采用以上技术方案，与现有技术相比，具有如下技术效果：
[0045] 本发明能够满足有限资源的在轨星上处理的高时效需求，实现了高光谱影像在轨数据的快速、精确的预处理。
附图说明
[0046] 为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍。在所有附图中，类似的元件或部分一般由类似的附图标记标识。附图中，各元件或部分并不一定按照实际的比例绘制。
[0047] 图1为本发明的三维查找表结构图；
[0048] 图2为本发明的中高光谱卫星的在轨数据快速预处理方法的运行流程图。
具体实施方式
[0049] 实施例1，如图2所示，一种高光谱卫星的在轨数据快速预处理方法，包括以下步骤：
[0050] 步骤1，卫星接收并解析任务指令，根据指令要求执行拍摄任务，获取原始数据，同时获取星上测量器件输出平台的相关辅助数据；
[0051] 步骤2，根据辅助数据的编码格式，分别解析出时间、轨道、姿态有效辅助数据。其中：
[0052] 时间数据包括行时和UTC时间；
[0053] 轨道数据是按照UTC时间解析对应的WGS84坐标系下的卫星位置三维坐标，间隔为1或2秒；
[0054] 姿态数据是按照UTC时间解析对应的四元数据，间隔为1或2秒。
[0055] 步骤3，根据步骤2中的有效辅助数据和任务点经纬度对所选谱段影像进行目标区域剪裁，获取目标区域数据。
[0056] 步骤3.1，获取时间数据中的中间行对应的UTC时间；
[0057] 步骤3.2，根据步骤3.1中的UTC时间，在轨道数据和姿态数据中获取对应UTC时间的卫星位置GNSS坐标 和姿态四元数 ；
[0058] 步骤3.3，根据式（1），利用姿态四元数获得J2000坐标系下的旋转矩阵R为：
[0059] 式（1）；
[0060] 将R转到本体坐标系下为：
[0061] ；
[0062] 其中， 、 、 、 分别为岁差矩阵、章动矩阵、地球自转旋转矩阵、极移矩阵，进而可利用式（2）获得卫星的俯仰pitch、翻滚roll的转角为：
[0063] 式（2）；
[0064] 其中，atan2为反正切函数。
[0065] 步骤3.4，获取水平、垂直方向的像素指向偏差。
[0066] 水平、垂直方向的像素指向偏差 为：
[0067] 式（3）；
[0068] 其中， 为相机载荷焦距， 为探测器每个探元的物理尺寸；
[0069] 步骤3.5，以目标点为中心画正方形，剪裁有效目标区域，其中，边长可设为的1.5 2.5倍之间。~
[0070] 步骤3.6，同理，将其他选中的谱段按照步骤3.5方法裁剪。
[0071] 步骤4，根据辐射定标参数，对目标区域数据各谱段进行快速辐射纠正。
[0072] 步骤4.1，建立辐射定标参数的m*n*c的3维查找表，见图1，其中，m表示列，列号与影像的列号相对应；n表示行，行号与谱段号相对应；c表示表号，与积分级数组号相对应。高光谱成像时，为了更好适配各种地物特征，包括多组成像积分级数。
[0073] 此外，辐射定标参数包括偏移参数和增益参数，每个存储单元采用4个Byte存储对应的2个定标参数，前2个Byte存储偏移、后2个Byte存储增益，每个存储单元可用 表示。
[0074] 步骤4.2，每次成像指定特定的积分级数，根据积分级数可确定查找表号确认定标参数，进而可以对剪裁区域执行多线程协同快速辐射纠正处理。以局部影像列为基本单元，设置二维线程组，线程组包括nn个线程组，每个线程组包括mm个线程，mm对应影像的列、nn对应谱段，一般可分别设置为列数、谱段数，即为m*n。每个线程按照下式进行计算，获取辐射纠正处理后的结果。第n组中第m个线程的纠正结果 的计算方法如式（4）：
[0075] 式（4）。
[0076] 步骤5，结合几何定标参数和严密成像几何模型，采用优化线性严密模型最小二乘求解法进行地理坐标快速解算。
[0077] 步骤5.1，考虑几何定标偏移矩阵描述了相机载荷、星敏、GNSS天线的安装及偏移误差，对严密成像几何模型进行简化，如下式（5）所示：
[0078] 式（5）；
[0079] 其中， 偏移矩阵可从定标参数中获取； 为相机内畸变参数，可从定标参数中获取；x、y为图像像素坐标； 为比例系数；可利用地球椭球性质，将严密模型改写成关于 的二次方程求解获取，为常规计算方法，这里不赘述。
[0080] 步骤5.2，选取四角点和中心点像素，带入式（5）求解地面坐标，式（5）改写为式（6）：
[0081] 式（6）；
[0082] 其中，下标1、2、3、4、5分别代表四角点和中心点。式（6）为线性方程，在求解过程中，高程值分别选取最大、最小、平均，因此根据式（6）可建立9个方程，进而通过最小二乘求解最终地面坐标。
[0083] 步骤6，四角点、中心点地面坐标求解之后，可通过线性内插法计算该剪裁区域任一点坐标，进而生成高精度的附带地理编码影像产品，从而完成在轨预处理任务。线性内插按照下式（7）进行计算：
[0084] 式（7）；
[0085] 其中，x为图像坐标，X为地理坐标，下标i表示待求解、下标e、b表示已知的2个点，例如四角点中的2个。
[0086] 步骤7，可选的步骤，根据实际需要，对其他谱段进行同样处理，并根据要求执行基于地理坐标的相对几何配准。
[0087] 本发明的描述是为了示例和描述起见而给出的，而并不是无遗漏的或者将本发明限于所公开的形式。很多修改和变化对于本领域的普通技术人员而言是显然的。选择和描述实施例是为了更好的说明本发明的原理和实际应用，并且使本领域的普通技术人员能够理解本发明从而设计适于特定用途的带有各种修改的各种实施例。
