目标光源差异对阵列式光谱辐射计测量的影响及紫外杂散辐射修正研究
2020-03-20吴志峰王彦飞代彩红
李 玲, 吴志峰, 王彦飞, 代彩红
中国计量科学研究院光学与激光计量研究所, 北京 100029
引 言
对地观测技术作为研究空间环境生态系统的重要手段[1],通过测量地球系统的光谱辐射度、反射率等物理量,反演计算地表碳水循环和云层覆盖率等信息[2-4]。无论星上在轨遥感器还是地面同步观测设备,全面、高质量的光谱辐射度定标是获得有效遥感观测数据的必要前提[5]。阵列式光谱辐射计因机械结构紧凑,以及全波段数据的快速采集,广泛应用于空间对地观测和遥感技术等领域[6-7]。然而,阵列式光谱辐射计内部元件在设计加工中,会出现由光栅刻槽精度误差、光学元件反射等原因造成的杂散辐射,严重影响光谱辐射度测量的准确性。
另一方面,外场光谱辐射计的目标光源是地表太阳辐射,其光谱分布与实验室的定标光源[8](光谱辐射照度标准灯)存在较大差异。实验室测的紫外杂散辐射比例不能直接应用于外场杂散辐射修正。因此,光谱辐射计的高精度杂散辐射修正,是地基验证场光谱辐射测量仪器获得真实可靠遥感数据的重要前提。
图1 实验室定标光源与地基验证场的光谱辐射亮度测量结果
Fig.1Theradiancemeasurementresultsofcalibrationlampinlaboratoryandsolarradianceinground-basedsite
1 光谱辐射计紫外杂散辐射测量
1.1 实验
为研究光谱辐射计紫外杂散辐射性能,本文选取多种典型的光谱辐射计进行实验,包括:AvaSpec-ULS2048(Avantes),Maya2000 Pro(Ocean Optics),HAAS-2000(Everfine Corp.),BLACK-Comet C-50(StellarNet),CR280(Colorimetry Research),PSR+3500(Spectral Evolution),和HR-1024i(SVC)。每台光谱辐射计定标后测量光谱辐射照度(Ec)或光谱辐射亮度(Lc),再在光路中增加截止滤光片(450 nm前截止),并保持其他实验条件和参数设置(积分时间、平均次数等)不变,测量光谱辐照度(Ex)或光谱辐射亮度(Lx),最后计算不同光谱辐射计在紫外波段的杂散辐射比例σ/%,其中σ=Ex/Ec×100或σ=Lx/Lc×100)。
1.2 结果分析
如图2所示,与扫描式光谱辐射计OL756(Gooch & Housego)测量结果对比,实验中的阵列式光谱辐射计在紫外波段都存在杂散辐射。尤其在250~260 nm波段杂散辐射比例高达80%左右。地基验证场同步观测设备中常用的光谱辐射度计CR280,在380 nm的杂散辐射比例为70%左右。
图2 (a) 光谱辐射计杂散辐射比例;(b) 地物光谱辐射计杂散辐射比例
Fig.2(a)Thestraylightradiosofarrayspectroradiometers;(b)Thestraylightradiosoffieldspectroradiometers
2 光谱辐射计紫外杂散辐射修正方法研究
杂散辐射是限制光谱辐射计高精度测量的重要原因,由于实验室定标光源与太阳辐射光谱分布不同,实验室测量的紫外杂散辐射比例不能直接应用于外场杂散辐射修正。本文针对光源的光谱分布特点,分别基于带通滤光片和可调谐激光器,研究光谱辐射计紫外杂散辐射修正方法。
2.1 基于带通滤光片
利用多组具有特定光谱透过率的带通滤光片,即420,550,650和800 nm带通(长波通短波截止)滤光片,分别测量不同滤光片下光谱辐射计CR280的杂散辐射信号比例。如图3所示,使用420 nm长通滤光片后,仪器测得380, 382, 384, …, 420 nm的信号值Ystray1,实际为420 nm至红外波段引起的杂散辐射信号。使用550 nm长通滤光片后,仪器测得的380, 382, 384, …, 420 nm的信号值Ystray2,实际为550 nm至红外波段引起的杂散辐射信号。将Ystray1-Ystray2=Ystray,即可获得420~550 nm之间光谱辐射信号总和对380, 382, 384, …, 420 nm造成的杂散辐射影响。
图3 不同带通滤光片下CR280的杂散辐射信号比例
用Ystray除以420~550 nm内光谱辐射总信号,得到杂散辐射比例
Rstray,380, 420~550 nm=(Ystray1,380-Ystray2,380)/Ymeas, 420~550 nm
以此类推,使用550,650和800 nm带通滤光片可以分别计算出550~650以及650~800 nm对380, 382, …, 420 nm测量值的杂散辐射比例。
Ytrue, 380=Ymeas, 380-Ystray, 380=Ymeas, 380-[Ystray,380, 420~550 nm+Ystray,380,550~650 nm+Ystray,380, 650~800 nm+
Ystray,380, 800~1 080 nm]=Ymeas, 380-[Rstray,380, 420~550 nm·
Y420~550 nm+Rstray,380, 550~650 nm·Y550~650 nm+
Rstray,380, 650~800 nm·Y650~800 nm+Rstray,380, 800~1 080 nm·
Y800~1 080 nm]
对于连续分布的宽谱段光源,相邻波长的杂散辐射影响接近,可利用上述数学模型近似计算出紫外各波长下的杂散辐射修正结果。如图4所示,对地基验证场地表太阳辐射的光谱辐射亮度测量结果进行修正后,可见光对紫外380~400 nm的杂散辐射影响明显减少。
图4 地基验证场的光谱辐射亮度测量结果和杂散辐射修正结果
由于短波(紫外)通长波截止的滤光片工艺非常困难,实验中无法使用。上述修正模型中只考虑了长波对短波的杂散辐射影响。另一方面,修正模型中对相邻波长的杂散辐射作了近似,因此,基于带通滤光片的修正方法受光源光谱分布的影响,很难实现杂散辐射的高精度修正。
2.2 基于可调谐激光器
任意波长的单色光入射到光谱辐射计时都会引起杂散辐射,光谱辐射计在某像素点(波长)下光谱辐射测量结果等于真实信号与杂散辐射信号之和。本文利用可调谐激光器(EKSPLA公司NL242-SH/SFG),改变输出波长,获得光谱辐射计各个像素点(波长)对不同波长入射光的杂散辐射线扩展函数(LSF),实验中光谱辐射计(AvaSpec-ULS2048, UV/Vis)的辐射信号由激光器的输出功率做归一化处理,如图5所示,光谱辐射计的杂散信号在5×10-3~5×10-4量级。
图5 光谱辐射计Avantes(UV/Vis)的LSF分布
入射波长为400 nm的杂散辐射线扩展函数LSF如图6所示,将带通内(IB)的辐射信号设置为零,再利用式(1)[9-10],计算出该波长LSF对应的杂散辐射比例di, J,得到光谱杂散辐射信号分布函数(SDF)。
(1)
由此,各个像素点(波长)对不同入射波长的杂散辐射比例di, J可写为矩阵D
(2)
仪器光谱辐射度测量值及真实值用矩阵形式表示为
Ym=[I+D]YIB=AYIB
(3)
YIB=A-1YIB=CYIB
(4)
根据杂散辐射修正式(1)—式(4),用MATLAB软件反演计算出光谱辐射计的真实值YIB。
图6 400 nm的线扩展函数LSF(黑色实线)及杂散辐射分布函数(SDF)(红色实线)
图7 Avantes光谱辐射计(UV/Vis)测量405 nm窄带光源的杂散辐射修正结果
Fig.7Themeasuredsignalandstraylight-correctedresultfor405nmnarrow-bandlightsourceofAvantes(UV/Vis)
为了验证杂散辐射修正方法,用Avantes光谱辐射计(UV/Vis)测量了405 nm窄带光源的光谱辐射信号,如图7所示,利用式(1)—式(4)进行杂散辐射修正后,杂散辐射信号由1×10-4左右降低到1×10-5左右,并且谱线两边的杂散宽峰也被消除。
图8 光谱辐射计CR280测量405 nm窄带光源的杂散辐射修正结果
Fig.8Themeasuredsignalandstraylight-correctedresultfor405nmnarrow-bandlightsourceofCR208
此外,利用可调谐激光器修正方法,对地基验证场光谱辐射计CR280进行了紫外杂散辐射修正,如图8所示,修正后测量结果在窄带光源两侧的杂散信号降低了1个数量级。
3 结 论
针对不同光谱分布的光源研究了紫外杂散辐射修正方法,结论如下:
(1)利用带通滤光片获得可见及红外光谱辐射引起的紫外杂散辐射比例,建立了紫外杂散辐射数学修正模型。对于连续分布的宽谱段光源,带通滤光片修正法具有实验简便易行、测试过程高效等优点;
(2)基于可调谐激光器建立了高精度杂散辐射测量系统,解决了各个像素点杂散辐射线扩展函数的测量难题,实现了光谱辐射计单像素点的高精度修正。适用于阵列式光谱辐射计对不同光源紫外测量结果的修正,修正后杂散辐射信号均降低了一个数量级,谱线两边的杂散宽峰也被显著消除。