一种机载高光谱遥感数据光度校正方法
【专利说明】一种机载高光谱遥感数据光度校正方法 【技术领域】
[0001] 本发明涉及一种批量化的机载高光谱遥感数据光度校正方法，属于高光谱遥感光 度校正领域，它适用于对机载高光谱遥感数据进行光度校正，以消除由于地表二向性效应 引起的辐亮度值异常，实现机载高光谱影像数据的无缝拼接。 【【背景技术】】
[0002] 遥感技术发展经历了全色、彩色摄影、多光谱扫描阶段之后，高光谱遥感技术的出 现成为20世纪80年代遥感技术的最大成就之一。目前，机载高光谱设备已在国内矿产资源 勘查领域进行广泛的应用，并取得了丰硕的成果。然而由于地表目标的非朗伯体特性，反射 光在空间中的分布往往是不均一的，即光谱特性随太阳与传感器高度角与方位角变化，使 得沿垂直飞行方向上辐射亮度值往往差异较大，短时间内获取影像在航带内部，尤其航带 间呈现出明显的波浪状异常，直接影响影像数据的拼接和表达，同时造成一些含量较小的 矿物弱异常会淹没在地表二向性效应所致的条带噪声中。
[0003] 目前，光度校正在行星探测中应用较多，在对地观测中也有类似的应用，但主要以 卫星数据为主。然而，对于机载高光谱数据来说，相关研究和应用案例较少。为此，本发明开 发了一种光度校正方法，对机载高光谱数据辐亮度值进行校正，实现影像的无缝拼接，降低 矿物检出限。 【
【发明内容】
】
[0004] 1.目的：本发明目的是提供一种机载高光谱数据的光度校正方法，针对地表二向 性引起的多航带影像之间产生的波浪状辐亮度值异常进行校正，实现机载高光谱影像数据 的无缝拼接。
[0005] 2.技术方案:本发明一种机载高光谱遥感数据光度校正方法，包括以下七个步骤：
[0006] 步骤一:原始数据读入
[0007] 原始数据主要包括三种，第一种是影像数据，即原始获取的未经几何校正的辐亮 度数据;第二种是信息类数据，包括飞行高度和传感器GPS航迹;第三种是影像坐标数据。以 上三种数据均由传感器获取数据时一并获取，为已知数据，均用于后续计算。
[0008] 步骤二:影像角度参数计算
[0009] 某一单航带高光谱数据具有m行，η列，ai*SS方向上任意行数据的平均传感器观 测角，i = l……m。通过传感器航迹的GPS文件和影像坐标数据，根据三角函数定律计算出第 i行数据上任意象元与传感器间的距离h，传感器的飞行高度h已知，可根据公式（1)计算出 任意行的传感器观测角cu。
[0010]
⑴
[0011] 步骤三:辐亮度与传感器观测角拟合
[0012] 设定AFcuS单航带内平行SS方向上任意行数据的平均辐亮度值，则单航带内m行 数据的平均辐亮度数据与对应传感器观测角之间存在以下关系：
[0013] AFai = Aai2+Bai+C (2)
[0014] 方程中i = l……m，AFai已知，可通过反演求出任意观测角度下的A，B，C方程系 数。
[0015] 步骤四：校正系数计算
[0016] 设AF1(a=())为星下点处任意行数据的拟合平均辐亮度值，即通过公式(2)计算出当 传感器观测角为〇时，对应任意行数据的平均辐亮度值，通过步骤三中求得的m行数据对应 的方程系数A，B，C，可求出m行数据的AFi( a=Q)值。
[0017] 由于传感器在获取数据时是实时运动的过程，传感器观测角的不断变化造成获取 数据辐亮度值的不同。在此，通过将获取数据的传感器观测统一校正到星下点处<^ = 0,来 实现航空高光谱数据的光度校正。
[0018] Ki = AFai/AFi(a=〇) (3)
[0019] K为校正系数，Ki则表示在任意给定的传感器观测角的行数据的光度校正系数。
[0020] 步骤五:光度校正
[0021] 第i行数据中，任意像元点处的辐亮度值SF1(x,y)，该点的光度校正后的辐亮度值 为FCi( X,y)，则满足关系式：
[0022] FCi(x,y)=Fi(x,y)/Ki (4)
[0023]其中Fi(x,y)已知，带入步骤四中求出的Ki，可对单航带机载高光谱数据进行逐点光 度校正。
[0024]步骤六:多航带数据无缝拼接
[0025]重复步骤二到步骤五，可求出各个航带光度校正后的辐亮度数据，进行几何校正、 镶嵌拼接(无需做匀色处理)后，可获得多航带影像的无缝拼接数据。
[0026]本发明优点及功效:本发明一种机载高光谱遥感数据光度校正方法，基于航空高 光谱数据辐亮度值与传感器天顶角间的函数关系，对原始数据进行逐行光度校正，具有运 算量小，精度高等优点，可实现对于航空高光谱数据的无缝拼接。 【【附图说明】】
[0027]图1为本发明的流程图。
[0028]图2为单航带数据辐亮度值与观测角拟合曲线图。
[0029]图3为单航带数据光度校正前后效果对比图。
[0030]图4为多航带数据光度校正前拼接效果图。
[0031 ]图5为多航带数据光度校正后无缝拼接效果图。 【【具体实施方式】】
[0032]图1说明了本发明的基本方法流程，为了更好地说明该方法的实现，以新疆东天山 地区HyMap机载高光谱数据为例，进行光度校正实验。
[0033] (1)试验所用的设备为图形工作站，规格型号为Dell Precision t7600,操作系统 为Windows7(64位），CPU为2.66GHz，内容为32GB，硬盘为IT。
[0034] (2)具体步骤如下：
[0035]步骤一:原始数据读入
[0036] 实验数据采用中国国土资源航空物探遥感中心于2002年10月15日使用HyMap机载 高光谱成像光谱获取的东天山地区高光谱数据，实验数据共包含四个航带。除搜集影像数 据外，还准备了包括飞行高度、传感器GPS航迹影像坐标数据。
[0037]步骤二:影像角度参数计算
[0038]首先，通过传感器航迹的GPS文件和影像坐标信息，根据三角函数定律计算出第i 行数据上任意象元与传感器间的距离li，传感器的飞行高度h已知，可根据公式（1)计算出 任意行的传感器观测角CU。 _9]
⑴
[0040] 步骤三:辐亮度与传感器观测角拟合
[0041] 由于单航带内平行SS方向上任意行数据的平均辐亮度值八?€[1与€[1满足如下关系 式：
[0042] AFai = Aai2+Bai+C (2)
[0043] 通过拟合求出任意观测角度〇1下对应的A，B，C方程系数，附图2展示了实验数据中 第一行带的厶?€4与€ [1函数拟合情况。
[0044]步骤四:校正系数计算
[0045] 将cti = 0，带入步骤三中的公式(2 )，可计算出当传感器观测角为0时，对应任意行 数据的平均辐亮度值AFi(a=0)。
[0046] 通过将获取数据的传感器观测统一校正到星下点处(αι = 0)，来实现航空高光谱 数据的光度校正。
[0047] Ki = AFai/AFi(a=0) (3)
[0048] 将求得的AF1(a=())带入公式(3)中，可求出任意给定的传感器观测角的行数据的光 度校正系数K 1<3
[0049] 步骤五:光度校正
[0050] 第i行数据中，任意像元点处的辐亮度值SF1(x,y)，将步骤四种求得的1带入公式 (4)中，求得任意像元点处光度校正后的辐亮度值为Fc 1(x,y)。
[0051] FCi(x,y)=Fi(x,y)/Ki (4)
[0052] 根据公式(4)对单航带内高光谱数据进行逐点计算，可求得单航带光度校正后的 辐亮度数据。如附图3所示，经过光度校正的数据在垂直SS方向上的辐亮度值差异变小，在 视觉上显示更平滑。
[0053]步骤六:多航带数据无缝拼接
[0054]重复步骤二到步骤五，分别求出四个航带光度校正后的辐亮度数据，利用步骤一 中获取的git文件，对各航带辐亮度数据进行几何校正，最后通过镶嵌拼接(无需做匀色处 理)可获得多航带影像的无缝拼接数据通过镶嵌拼接形成最终的多航带无缝拼接数据。附 图4展示了未经光度校正处理的多航带数据拼接后的效果，从中可明显的看出由于辐亮度 值异常产生的条带间的波浪状条纹。附图5展示了经过光度校正处理后多航带数据无缝拼 接的效果，与附图4对比，其拼接效果有了极大的改善。
【主权项】
1. 一种机载高光谱遥感数据光度校正方法，包括以下七个步骤： 步骤一:原始数据读入 原始数据包括三种，第一种是影像数据，即原始获取的未经几何校正的辐亮度数据;第 二种是信息类数据，包括飞行高度和传感器GPS航迹;第三种是影像坐标数据；以上三种数 据均由传感器获取数据时一并获取，为已知数据，均用于后续计算； 步骤二:影像角度参数计算 某一单航带高光谱数据具有m行，η列，ai*SS方向上任意行数据的平均传感器观测角，i =1……m;通过传感器航迹的GPS文件和影像坐标数据，根据三角函数定律计算出第i行数 据上任意象元与传感器间的距离h，传感器的飞行高度h已知，根据公式（1)计算出任意行 的传感器观测角cu;步骤三:辐亮度与传感器观测角拟合 设定AFcu为单航带内平行SS方向上任意行数据的平均辐亮度值，则单航带内m行数据 的平均辐亮度数据与对应传感器观测角之间存在以下关系： AFai=Aai2+Bai+C (2) 公式⑵中i = l……m，AFai已知，通过反演求出任意观测角度下的A，B，C方程系数； 步骤四：校正系数计算 设AF1(a=())为星下点处任意行数据的拟合平均辐亮度值，通过公式(2)计算出当传感器 观测角为〇时，对应任意行数据的平均辐亮度值，通过步骤三中求得的m行数据对应的方程 系数A，B，C，求出m行数据的AFi( a=〇)值； 由于传感器在获取数据时是实时运动的过程，传感器观测角的不断变化造成获取数据 辐亮度值的不同；通过将获取数据的传感器观测统一校正到星下点处~=〇,来实现航空高 光谱数据的光度校正； Ki = AF〇i/AFi(a=〇) (3) K为校正系数，L则表示在任意给定的传感器观测角的行数据的光度校正系数； 步骤五:光度校正 第i行数据中，任意像元点处的辐亮度值为F i ( X , y )，该点的光度校正后的辐亮度值为 FCi(x,y)，则满足关系式： FCi(x,y)=Fi(x,y)/Ki (4) 其中Fi(X,y)已知，带入步骤四中求出的Ki，对单航带机载高光谱数据进行逐点光度校正； 步骤六:多航带数据无缝拼接 重复步骤二到步骤五，求出各个航带光度校正后的辐亮度数据，进行几何校正、镶嵌拼 接后，获得多航带影像的无缝拼接数据。
【专利摘要】本发明提供一种机载高光谱数据的光度校正方法，针对地表二向性引起的多航带影像之间产生的波浪状辐亮度值异常进行校正，实现机载高光谱影像数据的无缝拼接，具体包括以下七个步骤：步骤一：原始数据读入；步骤二：影像角度参数计算；步骤三：辐亮度与传感器观测角拟合；步骤四：校正系数计算；步骤五：光度校正；步骤六：多航带数据无缝拼接。所述光度校正方法是基于航空高光谱数据辐亮度值与传感器天顶角间的函数关系，对原始数据进行逐行光度校正，具有运算量小，精度高等优点，可实现对于航空高光谱数据的无缝拼接。
【IPC分类】G01J3/28, G01S17/89, G01S7/497
【公开号】CN105676208
【申请号】CN201610007130
【发明人】于峻川, 刘文良, 闫柏琨
【申请人】中国国土资源航空物探遥感中心
【公开日】2016年6月15日
【申请日】2016年1月6日