Verification method of spectral calibration accuracy for hyperspectral remote sensors
-
摘要: 鉴于精确的光谱辐射定标精度验证对于高光谱遥感器的现场定标和数据应用非常必要,以SVC光谱辐亮度计为例,利用多波段辐亮度标准传递探测器和新型的光谱可调积分球参考光源,设计了一种高光谱遥感器光谱定标精度的验证方法。该方法利用新型光谱可调积分球参考光源在待测波段内分别输出光谱形状单调上升和单调下降的光谱辐亮度状态,通过光谱匹配技术,即平移改变SVC光谱辐亮度计的波长,分析比对MRSTD和SVC光谱辐亮度计测量辐亮度的相对偏差。比对结果为光谱定标验证精度优于±0.2 nm,辐射定标验证精度小于5%。Abstract: Validation for radiometric and spectral calibration accuracy of hyperspectral remote sensor is necessary. Taking SVC spectral radiometer as an example, new type spectrum-tunable reference light source (STIS) and muti-band radiance standard transfer detector(MRSTD) are used to design validation method for radiometric and spectral calibration accuracy. In this method, spectral level of spectral shape is monotonically increasing and monotonically decreasing in measured band by new spectral adjustable integrating sphere reference light source. Wavelength of SVC spectral radiance is changed by spectral matching technique. Relative deviation of radiance measured by MRSTD and SVC spectral radiometer are calculated and contrasted. Contrast shows that spectral calibration accuracy is better than ±0.2 nm, and absolute radiometric calibration accuracy is less than 5%.
-
Keywords:
- standard transfer detector /
- calibration /
- reference light source /
- radiance /
- validation
-
引言
高光谱遥感器是指利用分光技术获取连续窄带光谱的遥感器,其通道光谱带宽Δλ与中心波长λ的关系为Δλ/λ≈0.01[1]。在高光谱遥感器使用前,需要进行光谱定标和辐射定标。光谱定标的目的是建立高光谱遥感器光谱通道与不同入瞳波长的对应关系,辐射定标时建立高光谱遥感器的输出值与相应已知的、用国际单位(SI)表示的辐射量基准之间定量关系[2-3]。光谱定标和辐射定标是保证高光谱遥感器获取数据的科学可信度和最大限度的定量化应用基础。
高光谱遥感器的光谱定标和辐射定标通常在实验室内进行,其中光谱定标通常利用基于元素特征谱的光谱定标参考光源[4](低压汞灯和高压钠灯等特征谱线灯安装在积分球里面)或经过波长校准的单色仪或可调谐激光器[5]等精细扫描的方法实现,辐射定标通常采用光谱平坦性较好的标准灯板系统或传统卤素灯积分球光源完成。由于高光谱遥感器工作环境复杂,通常需要在野外进行观测,尤其是长期观测,受本身老化和外界干扰等因素的影响,其光谱中心波长漂移会发生漂移,其辐射响应特性会不断衰减[6]。前期高光谱遥感器应用研究表明,中心波长漂移与通道光谱带宽比值在超出1%情况下,在光谱形状变化剧烈的谱段,辐射测量的偏差有可能超过5%,甚至达到20%[7]。因此有必要对高光谱遥感器光谱辐射定标精度在现场环境下进行精度验证,实现光谱定标精度优于±0.5 nm,确保辐射定标不确定小于5%。
目前在现场环境下,可利用大气廓线中的太阳夫琅禾费线以及氧气、二氧化碳和水汽等气体吸收光谱特征对高光谱遥感器进行光谱定标精度验证,例如Gao等人采用光谱匹配技术实现了机载高光谱成像光谱仪PHILLS的光谱定标[8-9]。通过选择大气干洁、地表开阔的场地,利用太阳光经过不同反射率漫反射参考板或靶标的反射信号对高光谱遥感器进行辐射定标精度验证[10-11]。这种定标方法较精确,但是在实际操作中,因为很难能找到稳定、晴好、无云的天气条件,进而导致定标数据质量下降。
现阶段尚未有高精度的光谱定标和辐射定标实现现场环境下光谱辐射定标验证方法。本文设计了一种高光谱遥感器现场光谱辐射定标精度的验证方法,以SVC光谱辐亮度计为例,利用新型具有高稳定性的光谱可调积分球参考光源STIS(a spectrum-tunable reference light source) [12],在待验证波段内输出光谱辐亮度形状单调上升或单调下降状态,通过比对多波段辐亮度标准传递探测器MRSTD(a muti-band radiance standard transfer detector)[11]和SVC光谱辐亮度计辐亮度测量值的相对偏差,实现SVC光谱辐亮度计光谱辐射定标的精度验证和校正,其中光谱定标验证精度小于±0.2 nm, 辐射定标验证精度优于5%。
1 精度验证的方法
本文使用STIS和MRSTD验证高光谱光学遥感器实验室光谱定标精度。STIS采用多种波长的LED密集阵列发光单元和卤钨灯发光单元,在400 nm~800 nm具有分立波段输出功能,可输出半高带宽约20 nm~40 nm近似高斯形状的光谱辐亮度状态。STIS采用循环水冷温控和高精度可编程直流电源驱动技术,保障了光源的稳定性。STIS角度均匀性在±60°以内优于99.4%,面均匀性优于99.5%,2 h内稳定性优于99.8%。STIS积分球内径250 mm,出光口直径为60 mm,如图 1所示。
MRSTD是基于低温辐射计传递的多波段辐亮度标准传递探测器[13],用于绝对辐射定标以及监测参考目标的辐射特性变化。MRSTD的工作波段有412 nm、443 nm、490 nm、520 nm、565 nm、670 nm、750 nm和865 nm,其通道带宽约为(10±2)nm。MRSTD采用激光导入积分球的单色面光源进行光谱辐射定标,其中心波长定标准确度优于±0.1 nm,其绝对不确定度 < 1.26%。
设计的高光谱遥感器光谱定标的精度验证方法如下:在MRSTD各通道的带内区域,利用STIS输出光谱形状单调上升和单调下降的光谱辐亮度状态。使用MRSTD和SVC光谱辐亮度计交替测量STIS输出的光谱辐亮度,如图 2所示。通过比较MRSTD和SVC光谱辐亮度计辐亮度测量结果的相对偏差,分析SVC光谱辐亮度计的光谱辐射定标精度。
利用标准传递探测器第i通道响应DNi和其带内响应度为Ri,测量光谱可调积分球参考光源输出辐亮度Lc(λi)作为参考值,如(1)式和(2)式所示:
$$ {{L}_{c}}\left( {{\lambda }_{i}} \right)={D{{N}_{i}}}/{{{R}_{i}}}\; $$ (1) $$ {{R}_{i}}=\int_{{{\lambda }_{l}}}^{{{\lambda }_{u}}}{R\left( {{\lambda }_{i}} \right)}\text{d}\lambda $$ (2) SVC光谱辐亮度计在350 nm~1 000 nm范围内的光谱分辨率,即半高带宽优于3.5 nm,小于MRSTD通道带宽。根据前期研究[14],SVC光谱辐亮度计通道光谱响应度一般接近于sinc函数,与MRSTD通道的类似高斯光谱响应度曲线形状不同,因此两种仪器测量的辐亮度无法直接进行比对。
解决方法:首先对SVC光谱辐亮度计测量的光谱辐亮度曲线进行线性插值,插值波长间隔为0.2 nm;其次卷积SVC光谱辐亮度计测量的光谱辐亮度Lsvc(λi)和MRSTD第i通道的相对光谱响应度R(λi),得到光谱辐亮度计对应MRSTD第i通道辐亮度测量值LBSR(λi):
$$ {{L}_{BSR}}\left( {{\lambda }_{i}} \right)=\int_{{{\lambda }_{l}}}^{{{\lambda }_{u}}}{{{L}_{SVC}}\left( {{\lambda }_{i}} \right)R\left( {{\lambda }_{i}} \right)}{\text{d}\lambda }/{{{R}_{i}}}\; $$ (3) 计算MRSTD的参考辐亮度Lc(λi)和SVC光谱辐亮度计辐亮度LBSR(λi)测量值之间的相对偏差Δi:
$$ {{\Delta }_{i}}=\frac{{{L}_{BSR}}\left( {{\lambda }_{i}} \right)-{{L}_{c}}\left( {{\lambda }_{i}} \right)}{{{L}_{c}}\left( {{\lambda }_{i}} \right)}\times 100% $$ (4) 当SVC光谱辐亮度计存在波长偏移时,参考辐亮度Lc(λi)和辐亮度LBSR (λi)测量值之间的相对偏差Δi差异明显。在MRSTD通道对应(10±2)nm光谱范围,SVC光谱辐亮度计的光谱带宽可默认为不发生变化。此时参考前期研究方法[8-9],利用光谱匹配技术,即对SVC光谱辐亮度计测量的光谱辐亮度LSVC(λi)进行波长平移,判断Δi的变化情况。
根据测量不确定度合格评定通用比对标准[15],当SVC光谱辐亮度计不存在波长偏移时,参考辐亮度Lc(λi)和辐亮度LBSR (λi)测量值的相对偏差Δi应满足下式:
$$ {\Delta _i} \le \sqrt {{u^2}\left( {{\rm{SVC}}} \right) + {u^2}\left( {{\rm{MRSTD}}} \right)}<5\% $$ (5) 式中:u(SVC)为SVC光谱辐亮度计的辐射定标不确定度;u(MRSTD)为MRSTD辐射定标不确定度。设计最大相对偏差Δi<5%为判据,认为此时平移后的波长为SVC光谱辐亮度计真实的波长,进而验证SVC光谱辐亮度计的光谱辐射定标精度。
2 实验及数据处理
利用STIS在MRSTD第2通道的带内区域输出4种光谱辐亮度度状态,对应的峰值中心波长分别为474.5 nm、475.9 nm、496.8 nm和499.7 nm。MRSTD测量的辐射亮度值分别为4.409μW·cm-2·sr-1、6.0377μW·cm-2·sr-1、13.950 7 μW·cm-2·sr-1和13.967 7 μW·cm-2·sr-1。图 3(a)为SVC光谱辐亮度计测量的光谱辐亮度曲线。图 3(b)为MRSTD第2通道490 nm波段的相对光谱响应度曲线,其辐亮度测量值为参考值,如表 1的第一行所示。
表 1 490 nm波段各波长相对偏移量(%) Table 1. Relative deviation of each wavelength in 490 nm band波长平移量/nm 相对偏差Δi 474.5 nm 475.9 nm 496.8 nm 499.7 nm 1.0 2.68 3.13 -4.65 -5.80 0.8 1.10 1.62 -3.66 -4.43 0.6 -0.47 0.14 -2.69 -3.08 0.4 -2.01 -1.32 -1.76 -1.74 0.2 -3.53 -2.76 -0.86 -0.43 0 -5.04 -4.18 0 0.85 -0.2 -6.53 -5.58 0.83 2.12 SVC光谱辐亮度计的波长不作平移时,计算MRSTD的参考辐亮度Lc(λi)和SVC光谱辐亮度计辐亮度LBSR (λi)测量值之间的相对偏差Δi最大为5.04%。将SVC波长向长波方向平移0.4 nm后,针对4种光谱辐亮度状态,计算的相对偏差Δi均小于2.01%。SVC光谱辐亮度计波长平移(长波方向为正)与相对偏差Δi的变化情况如图 4所示,相对偏差Δi变化的具体数值见表 1。
从表 1可以看出,针对STIS输出的474.5 nm和475.9 nm两种光谱辐亮度状态,SVC光谱辐亮度计波长平移与相对偏差Δi的变化规律趋向单调上升,对于STIS输出的496.8 nm和499.7 nm2种光谱辐亮度状态,SVC光谱辐亮度计波长平移与相对偏差Δi的变化规律呈单调下降趋势。通过两种STIS分别输出光谱形状单调上升和单调下降的两种光谱辐亮度状态,结合SVC光谱辐亮度计波长平移与MRSTD的参考辐亮度和SVC辐亮度测量值相对偏差Δi的变化规律和相对偏差Δi优于5%的判据,得到SVC光谱辐亮度计的波长偏移了0.4 nm。
SVC光谱辐亮度计760 nm波段的光谱和辐射定标精度验证情况如图 5所示,SVC的波长偏移了1.8 nm。SVC光谱辐亮度计波长平移与MRSTD的参考辐亮度和SVC辐亮度测量值相对偏差Δi的变化规律明显。760 nm波段测量与490 nm波段测量结果相比,SVC光谱辐亮度计波长平移与相对偏差Δi的变化规律呈单调上升和下降趋势较为平坦,主要是由于与峰值中心波长474.5 nm、475.9 nm、496.8 nm和499.7 nm的光谱辐亮度状态相比,STIS输出峰值中心波长748.3 nm和759.9 nm状态的光谱形状单调上升和单调下降的变化较为平缓。
同样我们验证了SVC光谱辐亮度计在350 nm~800 nm光谱范围内MRSTD其余波段的光谱辐射定标精度。受制于STIS内部LED的种类,在800 nm以后无法输出半高带宽约20 nm~40 nm近似高斯形状的光谱辐亮度状态。
3 结论
本文以SVC野外光谱辐亮度计模拟高光谱遥感器,通过STIS和MRSTD设计了一种高光谱光学遥感器光谱辐射定标的精度验证方法。借助2种STIS在待验证的光谱通道内分别输出光谱形状单调上升和单调下降的两种光谱辐亮度状态,利用SVC光谱辐亮度计波长平移与MRSTD的参考辐亮度和SVC辐亮度测量值相对偏差Δi的变化规律,可有效验证SVC光谱辐亮度计光谱和辐射定标精度,其中光谱定标验证精度优于0.2 nm,辐射定标精度优于5%。这种光谱辐射定标精度验证方法可用于高光谱遥感器350 nm~800 nm光谱范围内实验室定标和现场定标,同时也为光谱连续型光学遥感器与窄带滤光片型光学遥感器之间的辐亮度测量提供一种参考方法。
-
表 1 490 nm波段各波长相对偏移量
(%) Table 1 Relative deviation of each wavelength in 490 nm band
波长平移量/nm 相对偏差Δi 474.5 nm 475.9 nm 496.8 nm 499.7 nm 1.0 2.68 3.13 -4.65 -5.80 0.8 1.10 1.62 -3.66 -4.43 0.6 -0.47 0.14 -2.69 -3.08 0.4 -2.01 -1.32 -1.76 -1.74 0.2 -3.53 -2.76 -0.86 -0.43 0 -5.04 -4.18 0 0.85 -0.2 -6.53 -5.58 0.83 2.12 -
[1] 翟文超, 徐骏, 郑小兵, 等.应用于超光谱定标的光谱可调光源光学设计[J].红外与激光工程, 2014, 43(3):951-955. http://d.old.wanfangdata.com.cn/Periodical/hwyjggc201403050 Zhai Wenchao, Xu Jun, Zheng Xiaobing, et al. Optical design of a spectrally tunable source for hyperspectral calibration[J]. Infrared and Laser Engineering, 2014, 43(3):951-955. http://d.old.wanfangdata.com.cn/Periodical/hwyjggc201403050
[2] 乔延利, 郑小兵, 王先华, 等.卫星光学传感器全过程辐射定标[J].遥感学报, 2006, 10(5):616-623. http://d.old.wanfangdata.com.cn/Periodical/ygxb200605002 Qiao Yanli, Zheng Xiaobing, Wang Xianhua, et al. Whole-process radiometric calibration of optical remote sensors[J]. Journal of Remote Sensing, 2006, 10(5):616-623. http://d.old.wanfangdata.com.cn/Periodical/ygxb200605002
[3] 郑小兵.高精度卫星光学遥感器辐射定标技术[J].航天返回与遥感, 2011, 32 (05): 36-43. doi: 10.3969/j.issn.1009-8518.2011.05.009 Zheng Xiaobing. High-accuracy radiometric calibration of satellite optical remote sensors[J]. Spacecraft Recovery & Remote Sensing, 2011, 32 (05): 36-43. doi: 10.3969/j.issn.1009-8518.2011.05.009
[4] 李新, 张国伟, 寻丽娜, 等.短波红外平场光谱仪的波长定标[J].光学学报, 2008, 28 (5): 920-926 http://d.old.wanfangdata.com.cn/Periodical/gxxb200805016 Li Xin, Zhang Guowei, Xun Lina, et al. Wavelength calibration of shortwave infrared flat spectroradiometer[J]. Acta Optica Sinica, 2008, 28 (5): 920-926. http://d.old.wanfangdata.com.cn/Periodical/gxxb200805016
[5] 徐秋云, 郑小兵, 李健军, 等.外部导入激光的积分球辐射源的研制[J].光学精密工程, 2009, 17(4): 738-744. doi: 10.3321/j.issn:1004-924X.2009.04.008 Xu Qiuyun, Zheng Xiaobing, Li Jianjun, et al.. Development of laser illuminated integrating sphere source[J]. Optics and Precision Engineering, 2009, 17(4): 738-744. doi: 10.3321/j.issn:1004-924X.2009.04.008
[6] 陈洪耀, 张黎明, 李鑫, 等.高光谱遥感器飞行中基于大气特征谱线的光谱定标技术[J].光学学报, 2013, 33 (5): 287-293. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gxxb201305042 Chen Hongyao, Zhang Liming, Li Xin, et al. Hyperspectral sensor in flight spectral calibration based on characteristic spectra of atmosphere[J]. Actc Optica Sinica, 2013, 33 (5): 287-293. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gxxb201305042
[7] Green R O. Spectral calibration requirement for Earth-looking imaging spectrometers in the solar-reflected spectrum[J]. Applied Optics, 1998, 37(4): 683- 690. doi: 10.1364/AO.37.000683
[8] Gao Bocai, Montes M J, Davis C O.Refinement of wavelength calibrations of hyperspectral imaging datausing a spectrum-matching technique[J]. Remote Sensing of Environment, 2004, 90(4):424-433. doi: 10.1016/j.rse.2003.09.002
[9] Green R O. Determination of the in-flight spectral calibration of AVIRIS using atmospheric absorption features[M].US:JPL Publication, 1995.
[10] 陈伟, 晏磊, 勾志阳, 等.无人机多光谱传感器场地绝对辐射定标研究[J].光谱学与光谱分析, 2012, 32(12):3169-3174. doi: 10.3964/j.issn.1000-0593(2012)12-3169-06 Chen Wei, Yan Lei, Gou Zhiyang, et al. In-flight absolute radiometric calibration of UAV multispectral sensor[J]. Spectroscopy and Spectral Analysis, 2012, 32(12):3169-3174. doi: 10.3964/j.issn.1000-0593(2012)12-3169-06
[11] 勾志阳, 晏磊, 陈伟, 等.无人机高光谱成像仪场地绝对辐射定标及验证分析[J].光谱学与光谱分析, 2012, 32(5):430-434. http://d.old.wanfangdata.com.cn/Periodical/gpxygpfx201202032 Gou Zhiyang, Yan Lei, Chen Wei, et al. In-flight absolute radiometric calibration of UAV hyperspectral camera and its validation analysis[J]. Spectroscopy and Spectral Analysis, 2012, 32(5):430-434. http://d.old.wanfangdata.com.cn/Periodical/gpxygpfx201202032
[12] 袁银麟, 徐骏, 翟文超, 等.大孔径可调光谱积分球参考光源研制和检测[J].光学学报, 2013, 33(7):0712004-1-8. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gxxb201307034 Yuan Yinlin, Xu Jun, Zhai Wenchao, et al. Design and test of a spectrally tunable integrating sphere reference light source with large exit aperture[J]. Acta Optica Sinca, 2013, 33(7):0712004-1-8. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gxxb201307034
[13] 张艳娜, 郑小兵, 李健军, 等.多波段滤光片辐亮度计的系统级定标方法研究[J].光学学报, 2013, 33(8):0812004-1-7. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gxxb201308020 Zhang Yanna, Zheng Xiaobing, Li Jianjun, et al.Research on system level calibration method of multi-band filter radiometer[J]. Acta Optica Sinca, 2013, 33(8):0812004-1-7. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gxxb201308020
[14] Brown S W, Johnson B C, Feinholz M E, et al. Stray-light correction algorithm for spectrographs[J]. Metrologia, 2003m 40(1)81-84. http://cn.bing.com/academic/profile?id=dc2ebfca864649dbb0d25038e3e06d65&encoded=0&v=paper_preview&mkt=zh-cn
[15] 倪育才.实用测量不确定度评定[M].第三版.北京:中国计量出版社, 2009. Ni Yucai. Evaluation of practical measurement uncertainty[M].3rd ed.Beijing: Chinese Metrology Press, 2009.
-
期刊类型引用(8)
1. 张嘉显,王军,曹召良. 基于粒子群算法的高均匀性接触角测量光源研究. 激光与光电子学进展. 2025(03): 186-192 . 百度学术
2. 吴铭扬,张巧芬,王桂棠,吴黎明,邓耀华. 基于LED反光杯系统的优化研究. 应用光学. 2024(02): 329-336 . 本站查看
3. 何慧萌,杨婷,施会丽,王平,邴喆,王星,白勃. 考虑HPSA的IRS辅助室内VLC系统光源优化. 西安电子科技大学学报. 2024(02): 46-55 . 百度学术
4. 徐承军,李嘉群,张鹏. 港口高压开关柜局部放电信号自适应小波去噪方法. 起重运输机械. 2024(11): 67-75 . 百度学术
5. 李谱林,吴玲玲,陈靖,占春连,余欢. 大面积高均匀照明的LED阵列设计. 激光与光电子学进展. 2024(13): 368-377 . 百度学术
6. 汤雄,王鹏,吕玲玲. LED植物照明灯具技术发展及应用现状. 应用技术学报. 2023(02): 120-124 . 百度学术
7. 曹震,巢渊,徐魏,杜帅帅,张敏. 基于改进梯度加权的零件图像高精度聚焦方法. 电子测量与仪器学报. 2023(11): 132-142 . 百度学术
8. 樊华羽,李典,郝海兵,梁益华. 基于Alpha-stable的粒子群算法变异策略研究及气动优化应用. 西北工业大学学报. 2022(06): 1385-1393 . 百度学术
其他类型引用(5)