Corrective method for spectral offset error caused by radial distortion in the large aperture static imaging spectrometer
AN Ling-ping, WANG Shuang, ZHANG Geng, LI Juan, LIU Xue-bin
引用本文:
安玲坪,王爽,张耿,李娟,刘学斌. 大孔径静态干涉成像光谱仪径向畸变导致的谱线偏移误差的校正[J]. 中国光学, 2021,14(2): 382-389. doi: 10.37188/CO.2020-0084
AN Ling-ping, WANG Shuang, ZHANG Geng, LI Juan, LIU Xue-bin. Corrective method for spectral offset error caused byradial distortion in the large aperture static imaging spectrometer[J]. Chinese Optics, 2021, 14(2): 382-389. doi:10.37188/CO.2020-0084
在线阅读 View online: https://doi.org/10.37188/CO.2020-0084
您可能感兴趣的其他文章Articles you may be interested in大偏离度非球面检测畸变校正方法
Distortion correcting method when testing large-departure asphere中国光学. 2017, 10(3): 383 https://doi.org/10.3788/CO.20171003.0383基于中阶梯光栅的波长定标方法研究Spectral calibration based on echelle
中国光学. 2017, 10(3): 376 https://doi.org/10.3788/CO.20171003.0376大随机相位误差下条带模式合成孔径激光雷达成像实验
Stripmap mode synthetic aperture ladar imaging under large random phase errors condition中国光学. 2019, 12(1): 130 https://doi.org/10.3788/CO.20191201.0130低阶梯多级微反射镜高度误差分析及制作研究Error analysis and fabrication of low-stepped mirrors
中国光学. 2019, 12(4): 791 https://doi.org/10.3788/CO.20191204.0791Seya-Namioka单色仪中光栅曲率半径误差的影响及补偿
Effect and compensate of grating curvature radius error in Seya-Namioka monochromator中国光学. 2018, 11(4): 623 https://doi.org/10.3788/CO.20181104.0623光谱成像技术在海域目标探测中的应用
Application of spectral imaging technology in maritime target detection中国光学. 2017, 10(6): 708 https://doi.org/10.3788/CO.20171006.0708
第 14 卷 第 2 期2021年3月
中国光学Chinese Optics
Vol. 14 No. 2
Mar. 2021
文章编号 2095-1531（2021）02-0382-08
大孔径静态干涉成像光谱仪径向畸变导致的
谱线偏移误差的校正
安玲坪1,2，王 爽1 *，张 耿1，李 娟1，刘学斌1（1. 中国科学院 西安光学精密机械研究所 光谱成像技术重点实验室，陕西 西安 710119；2. 中国科学院大学 北京 100049）摘要：为提高大孔径静态干涉成像光谱仪在视场增大时的光谱定标精度，减小径向畸变对光谱精度的影响，本文提出一种基于光谱——畸变关联模型的光谱定标系数修正方法，给出了波数和波长修正公式。采用594.1 nm和632.8 nm气体激光器对成像光谱仪进行了光谱成像实验，并对数据进行了处理和分析。结果表明，当存在0.3%的桶形畸变时，边缘视场的反演光谱存在2 nm左右的偏移，利用本文方法校正后，谱线偏移减小到0.1 nm左右。该方法仅需根据镜头畸变参数即可完成修正，简化了实验室光谱定标流程，提高了工作效率，也可应用于星载干涉光谱数据的在轨参数校正。关 键 词：大孔径静态干涉成像光谱仪；光谱定标；误差校正；畸变模拟中图分类号：TP751.1 文献标志码：A doi：10.37188/CO.2020-0084
Corrective method for spectral offset error caused by radialdistortion in the large aperture static imaging spectrometer
AN Ling-ping1,2，WANG Shuang1 *，ZHANG Geng1，LI Juan1，LIU Xue-bin1
（1. Key Laboratory of Spectral Imaging Technology, Xi’an Institute of Optics and Precision Mechanics,
Chinese Academy of Sciences, Xi’an 710119, China；
2. University of Chinese Academy of Sciences, Beijing 100049, China）
* Corresponding author，E-mail: wangshuang@opt.ac.cn
Abstract: In order to improve the spectral calibration accuracy of the large aperture static imaging spectro-meter when its field of view is increased, and to reduce the influence of radial distortion on its spectral accur-acy, we propose a corrective method for spectral calibration coefficients based on a spectral distortion correl-ation model. To begin the process, the wave number and wavelength correction formulas are given. Using594.1 nm and 632.8 nm gas lasers, a spectroscopic imaging experiment was performed on the imaging spec-trometer, and the data was processed and analyzed. The results show that when there is a barrel distortion of0.3%, the inversion spectrum at the edge of the field of view shifts approximately 2 nm. After implementingthe corrective method of this paper, the line shift is reduced to approximately 0.1 nm. This method only needsto be corrected according to the lens distortion parameters, which simplifies the laboratory spectral calibra-tion process and improves work efficiency. It can also be applied to the orbit parameter correction of space-
收稿日期：2020-05-08；修订日期：2020-06-15
基金项目：国家重点研发计划项目（No. 2018YFB0504900，No. 2018YFB0504901）
Supported by National Key R&D Program of China (No. 2018YFB0504900, No. 2018YFB0504901)
第 2 期
安玲坪, 等: 大孔径静态干涉成像光谱仪径向畸变导致的谱线偏移误差的校正
383
borne interference spectral data.
Key words: Large Aperture Static Imaging Spectrometer (LASIS)；spectral calibration；error correction；dis-tortion simulation
1 引 言
成像光谱仪可同时获取目标的二维空间信息和一维光谱信息，具有“图谱合一”的特点，在目标特性检测、空间遥感以及气象探测等领域有重要应用[1-5]。20世纪80年代后期出现的无动镜干涉成像光谱技术[6-7]，为实现轻小型、高稳定性干涉成像光谱仪打下了基础。大孔径静态干涉成像光谱仪（Large Aperture Static Imaging Spectro-metry, LASIS）出现于90年代末[8]
，无狭缝和运动部件，高通量和高稳定性的优势并存，是实现高光谱成像的重要途经[9]
，有着广阔的应用前景。为获得高质量的光学信息，光学系统通常只在近轴区域，即小视场角成像。随着遥感仪器的进一步发展及数据实时性要求的不断提高，增大光谱仪器的视场角已成为一种趋势。
LASIS面阵探测器包含多个探测单元，当视场角增大时，单次成像包含更多目标信息，可有效减少扫描次数，提高工作效率的同时还可以减少目标随时间变化导致的误差。但是图像边缘的畸变效应也会更加明显，探测单元的离轴角会使光程差改变[10]，此时干涉图像存在径向畸变导致的误差，使得最终复原后光谱图两边相对于中间谱线位置发生偏移，有研究表明这种现象会影响辐射定标精度[11]
。谱线偏移后需要对光谱定标系数进行修正，且光谱定标对辐射定标影响很大。因此，校正径向畸变导致谱线偏移产生的误差非常有必要。
在实验室光谱定标过程中采用小视场多次定标的方法可实现精准定标，但是该方法大幅度降低了定标的工作效率，且随着成像光谱仪视场的增大，该方法很难满足实际工作需求。现有的畸变校正方法多为直接校正探测器采集到的信号[12-14]，这些方法并不适用于LASIS探测器接收到的干涉图，但是通过校正反演光谱可以实现间接校正。本文在分析LASIS的工作原理及光谱反演过程的基础上，提出了通过探测器畸变数据获得反演光谱偏移率，进而修正光谱定标系数的
方法。最后结合实验室定标数据证明了本方法能
够有效校正LASIS光谱定标误差，改善数据质量。
2 误差校正方法的理论分析
LASIS主要有前置光学系统，干涉系统，傅立叶成像系统和探测器4个部分[8,15]，结构如图1所示。LASIS光学系统是普通照相系统的前置镜，没有狭缝限制，能够实现高光通量，能量利用率高，具有大孔径优势。利用Sagnac横向剪切干涉仪产生光程差，不需要运动部件，系统简单，稳定性高。面阵探测器可获得一维空间的干涉信息强度，再通过主动或被动推扫即可获得另一维空间信息。
前置镜
Sagnac 横向剪切干涉仪准直镜
傅氏镜
探测器
图 1 LASIS结构示意图
Fig. 1 Schematic diagram of LASIS structure
LASIS干涉系统原理如图2所示。傅氏镜焦距为fF，某一像点S被Sagnac干涉仪剪切成在垂直于光轴方向相距d的2个虚像S1、S2，视场角为α的两束平行光通过傅氏镜后，在成像面处汇聚
到与光轴相距y的一点，此时存在光程差x，从而发生干涉[16]
。成像面
α
傅氏镜
S1d
yo
S2
x
fF
fF
图 2 LASIS干涉原理示意图Fig. 2 Schematic diagram of LASIS interference
384
中国光学
第 14 卷
x=d·sinα,（1）sinα≈tanα=yf,（2）Fx=dyf.（3）FLASIS测量目标光谱辐射的干涉信息，可通过数学计算方法利用变换的方式反演出目标光谱。光谱强度与干涉强度是傅立叶变换对，关系式如下：
I(x)=+∞B(υ)exp(2πiυx)dυ,（4）B(υ)=−∞+∞−∞I(x)exp(−2πiυx)dx,（5）式中，x为光程差，υ为波数，I为干涉强度，B为光谱强度。理想情况下光程差x可以取到无穷大来反演完整的目标光谱B，但是受光学系统和探测器尺寸限制，在干涉图采样过程中无法将光程差取到无穷大，令截止频率为L，得到的干涉图为：
I)=+L(x)L(x−LB(υ)exp(2πiυx)dυ=I(x)·rect2L,（6）式中，IL(x)相当于理想干涉图I(x)与矩形函rect(数
x)2L
相乘，此时得到的反演光谱相当于理想
反演光谱卷积了一个sinc函数，出现了展宽[17]。
B′(υ)=+∞−∞I(x)rect(x)2Lexp(−2πiυx)dx=B(υ)·[2L·sinc(2πυL)].（7）对于每个探测器像元，得到的信号强度是一定范围的积分，探测器像元(i,j)采集到的待测目标像元i的第j个光谱通道的辐射信号为：
I(i)=νj2νB′(υ)exp(2πiυx)dυ,（8）j1式中，νj1和νj2表示探测器像元(i,j)包含的最小和最大波数。
当视场角增大时，垂轴放大率会随之改变，成像面干涉图的径向畸变无法忽略，且畸变程度与像元位置到畸变中心的距离的平方成正比。对于图像来说，径向畸变不会改变清晰度，但是会改变图像的几何位置[18]，图2中两束相干光的汇聚点变为：
y′=y(1+kr2)=yK(r),（9）式中，k为径向畸变二次项系数，k=0时无畸变，
k>0时为正畸变，即桶形畸变，k<0时为负畸变，
即枕形畸变，畸变效果如图3所示。
(a) k<0(b) k=0(c) k>0
图 3 径向畸变效果图Fig. 3 Diagram of the radial distortion effects
此时，两束光线的光程差为：
′x′=dyf=dy(1+kr2)=x(1+kr2)=xK(r).FfF（10）光谱强度与干涉强度是傅立叶变换对，根据傅立叶变换的尺度变换性质，发生径向畸变后的干涉图对应的反演光谱谱线位置会发生偏移，且偏移程度与畸变程度和像元位置有关：
I(x)↔B(υ)I(xK(r))↔1(,υ)|K(r)|BK(r).（11）若在进行光谱定标时，将谱线偏移产生的影响包含在定标系数之中，便可以减弱畸变导致的误差，修正谱线位置，改善数据质量。谱线偏移带来的影响可用偏移率表示，将每一条反演光谱原本的定标系数除以该条谱线的偏移率便可校正径向畸变所导致的光谱定标误差。
3 仿真分析与校正模型
选取脉冲谱作为目标光谱，验证本文误差模型的正确性。脉冲信号δ(υ)与常数1是傅立叶变换对，根据傅立叶变换频移性质，中心波数为υ的脉冲谱δ(υ−υ0)的干涉曲线为e−jxυ0
的实数部分，即
脉冲信号的干涉图为余弦曲线。
δ(υ)↔1,δ(υ−υjxυ00)↔e−.（12）本文首先利用Matlab软件对探测器采集信号的过程进行仿真，然后编写程序对数据进行处理，即将干涉数据反演为光谱数据，具体流程如图
4所示。第 2 期
安玲坪, 等: 大孔径静态干涉成像光谱仪径向畸变导致的谱线偏移误差的校正
385
干涉图模拟光谱反演一列干涉曲线(m×1)光谱图复制选取区间待采样干涉图傅立叶逆变换(m×n)径向畸变补零畸变后的干涉图(m×n)切趾矩形采样采样后的干涉图(M×N)直流基线去除
图 4 仿真试验流程图
Fig. 4 Flow chart of the simulation experiment
余弦曲线由Matlab软件生成，面阵探测器可同时采集多个目标的干涉曲线，将余弦曲线进行复制，并代入式（9）模拟径向畸变过程，得到待采样干涉图。由于Matlab的数据是由离散的点组成的，因此m与n的取值要尽可能大才能更接近实际的连续信号。接着，将图像分割成整齐分布的矩形区域，并根据式（8）对矩形区域进行积分，获得待处理干涉数据，为与下文实际的探测器像元数目一致，M取256，N取2 048。根据干涉曲
线光谱反演的过程[19]编写数据处理程序，遍历所有干涉曲线得到最终的光谱图。
由于LASIS的有效谱段为可见—近红外范围，同时为了与实验室光谱定标过程相匹配，本文选用了中心波长分别为400、543、632、694和900 nm的5种脉冲信号模拟目标光源。一般成像光谱仪的畸变都在4%以内[10]
，本文分别模拟了无畸变以及畸变率K为1%，2%和3%的桶形畸变与枕形畸变，最终得到每一列谱线峰值波数相对于中心列谱线（无畸变）峰值波数的偏移率P，如图5（彩图见期刊电子版）所示，6种线型分别代表6种脉冲信号，不同的畸变率用不同颜色
表示。
可以看出，偏移率P随畸变率K和像元位置的改变而改变，同种颜色的曲线近似重合在一起，偏移率P与波长相关性较小，可认为与波长无关。理论上，偏移率为连续曲线，但是经过矩形采样后，采集到的信号为离散信号，加之Matlab程序计算偏移率的精度为0.001，因此出现了跳变。仿真实验结果显示，径向畸变对反演谱线的影响可以用谱线偏移率表征，偏移率是与波长无关，与像
元位置有关的二次曲线，5种脉冲信号在畸变率相同时，偏移率变化趋势基本一致，通过拟合可得：
P(i)=K×10−8×R(i)2−K×10−8×R(i)+1（13,）式中，K表示畸变程度，K%为畸变率，R(i)为探测器像元i与畸变中心O之间的距离，R(i)=|O−i|
。1.04400 nm543 nm1.02632 nm−3%694 nm−2%900 nm−1%P1.000K1%0.982%3%0.965001 0001 5002 000
Pixel图 5 Matlab模拟结果
Fig. 5 Matlab simulation results
根据所得畸变模型，任意目标点i的第j个谱段中心波数υ(i,j)以及中心波长λ(i,j)的径向畸变误差可通过如下方法进行校正：
υ′(i,j)=υ(i,j)÷P(i)λ′(i,j)=λ(i,j)×P(i).（14）4 验证实验与结果分析
为了验证本文所提误差校正方法的性能，本文进行了实验室光谱定标实验，采集了LASIS探测器的干涉信息，通过进一步的光谱反演处理，确定了不同像元的光谱定标参数（单位光程差和最大光程差），为谱线校正提供了初始光谱定标依据。
图6为实验室光谱定标实验原理图，定标光源为标准激光器，利用激光器获得特定中心波长的光束，单色光经过扩束整形后照亮平行光管靶面处的毛玻璃，通过平行光管转换成光谱定标所需的平行光。平行光进入LASIS光学系统后最终在探测器上形成干涉图，利用数据采集软件实现图像采集
毛激光器扩束器
玻平行光管
大孔径静态
璃
干涉成像光谱仪
图 6 光谱定标实验示意图Fig. 6 Schematic diagram of the spectral calibration exper-iment
386
中国光学
1.0Normalized DN0.80.60.40.20560580600620λ/nm640Left (594.1 nm)Middle (594.1 nm)Right (594.1 nm)Left (632.8 nm)Middle (632.8 nm)Right (632.8 nm)第 14 卷
和存储。最后使用数据处理软件根据图4中光谱反演部分的流程对干涉图进行数据处理并分析。
LASIS探测器面阵大小为256（光谱维）×2 048（空间维），标准激光器发出的光束无法覆盖探测器空间维整个视场，所以分3次分别在左、中、右3个位置进行成像实验。畸变数据利用平行光管对LASIS探测器进行一维内方位元素标定（空间维）获得。经过成像和拟合，得到探测器采集到的图像为正畸变，畸变中心在第1 070列，即O=1070，畸变率为0.3%，即K=0.3。代入式(13)得到该探测器各列像元反演光谱的偏移率为：
P(i)=3×10−9×(1070−i)2−3×10−9×(1070−i)+1,i=1,2,···,2048.（15）660 图 8 激光三视场归一化反演光谱Fig. 8 The normalized inversion spectra of the three fields
of view of the lasers
对每一条反演光谱通过高斯拟合得到峰值波长，并利用式（14）、式（15）进行谱线校正。图9（彩图见期刊电子版）绘制了误差校正前后有效成像区域反演光谱的中心波长，表1列出了部分数据结果。
实验分别记录了594.1 nm激光器和632.8 nm激光器成像实验的干涉图像，以及有效像元的反演谱线。图7展示了实验获得的干涉图像，图8（彩图见期刊电子版）绘出了不同视场归一化后的反演谱线。从图8可以看出当存在桶形畸变时，两侧谱线的峰值波长向波长增大的方向发生了偏移。
660650640630620610600590200600Before (594.1 nm)After (594.1 nm)Reference (594.1 nm)λ/nmBefore (632.8 nm)After (632.8 nm)Reference (632.8 nm)
(a) 594.1 nm(b) 632.8 nm 1 000Pixel1 4001 800图 7 激光干涉图Fig. 7 Interferogram of the lasers
图 9 光谱定标误差校正Fig. 9 Correction of spectral calibration
表 1 激光中心波长对比
Tab. 1 Comparison of the lasers’ center wavelength
中心波长/nm
像元 100 4001 0001 1001 7002 000 100 4001 0001 1001 7002 000
632.8参考值λ0594.1
校正前λ595.576 2594.760 3594.134 2594.134 2594.718 6595.450 6634.542 9633.569 3632.835 1632.811 4633.379 6634.115 1
误差|λ−λ0|1.476 20.660 30.034 20.034 20.618 61.350 61.742 90.769 30.035 10.011 40.579 61.315 1
校正后λ′594.009 7594.012 6594.126 0594.132 8594.059 4594.013 6632.822 1632.736 8632.822 1632.811 0632.710 7632.633 7
误差|λ′−λ0|0.090 30.087 40.026 00.032 80.040 60.086 40.022 10.063 20.022 10.011 00.089 30.166 3
第 2 期
安玲坪, 等: 大孔径静态干涉成像光谱仪径向畸变导致的谱线偏移误差的校正
较好的校正效果。
387
从图9和表1可以看出，随着与畸变中心距离的增大，中心波长向长波方向移动，但偏移量与距离不成线性关系。校正前中间误差较小，两侧误差较大，序号为100和2 000的像元中心波长的误差均大于1 nm。利用本文的校正方法进行处理后，整体的误差得到了显著降低，图9中红线与绿线基本重合，由表1可知，除632.8 nm右视场边缘外，误差均小于0.1 nm。632.8 nm在右视场的校正效果略差的原因可能是由于探测器在安装过程中出现了倾斜、旋转等，使得实际畸变分布不是理想的桶形畸变[20]，进而导致了误差。实验证明，本文提出的谱线位置误差校正方法具有
为了验证该误差校正方法在红外波段的有效性，本文利用LASIS进行了室外推扫实验，并选取一帧图像进行光谱反演，分析了不同像元位置处氧气在红外波段的吸收带。分别选取了左中右3个视场序号为100、400、1 000、1 100、1 700和2 000的复原光谱，将746~779 nm光谱曲线进行归一化处理后，得到的氧气吸收带如图10（彩图见期刊电子版）所示。可以看出两侧视场存在纳米级的偏移，经过校正后，吸收带曲线更加聚集，表明该误差校正方法在近红外波段依然有效。
1.00.8Normalized DN0.60.40.21004001 0001 1001 7002 000 1.00.8Normalized DN0.60.40.207507551004001 0001 1001 7002 000 760765770λ/nm(a) 校正前(a) Before correction7750750755760765770λ/nm(b) 校正后(b) After correction775图 10 氧气红外吸收带Fig. 10 Absorption band of oxygen in infrared wavelengths
5 结 论
本文提出了一种径向畸变误差校正模型以校正光谱反演时出现的谱线偏移，用于大孔径静态干涉成像光谱仪在使用大面阵探测器时的光谱定标。根据干涉图与光谱图之间的对应关系和径向畸变的理论和数学表达式推导出了径向畸变与谱线偏移的模型，用计算机仿真实验验证模型，并给出了误差校正模型，最后实验验证了模型的准确性和可靠性。实验结果表明，本文所提出的校正参考文献：
[1]
方法可以有效减少谱线误差，对于594.1 nm激光，边缘视场的谱线误差从1.5 nm减少到0.1 nm，对于632.8 nm激光，尽管右边缘误差校正效果相对较差，但是也减少到了0.2 nm，在近红外波段也有明显校正效果，表明该方法准确可靠，提高了数据质量。尽管该方法是基于干涉型成像光谱仪实验室光谱定标提出的，但当仪器投入使用后，畸变数据可以通过载荷的外方位元素进行修正，从而光谱定标误差也可以不断校正，实现再定标，该方法具有长久的参考价值。
刘红婕, 王凤蕊, 耿峰, 等. 荧光成像技术无损探测光学元件亚表面缺陷[J]. 光学 精密工程，2020，28（1）：50-59.
LIU H J, WANG F R, GENG F, et al.. Nondestructive detection of optics subsurface defects by fluorescence imagetechnique[J]. Optics and Precision Engineering, 2020, 28（1）: 50-59. （in Chinese）[2]
何志平, 李春来, 吕刚, 等. 月球表面原位光谱探测技术研究与应用(特约)[J]. 红外与激光工程, 2020, 49(5): 62-69.
HE ZH P, LI CH L, LÜ G, et al.. Research and applications of In-Situ lunar surface spectral detection
388
中国光学
technology(Invited)[J]. Infrared and Laser Engineering, 2020, 49(5): 62-69. (in Chinese)
第 14 卷
[3]
郑光辉, 焦彩霞, 上官晨曦, 等. 基于成像光谱技术的土壤剖面发生层划分初探[J]. 光谱学与光谱分析，2019，
39（3）：882-885.
ZHENG G H, JIAO C X, SHANGGUAN CH X, et al.. Horizon classification in soil profile using imagingspectroscopy[J]. Spectroscopy and Spectral Analysis, 2019, 39（3）: 882-885. （in Chinese）[4]
朱梦远, 杨红兵, 李志伟. 高光谱图像和叶绿素含量的水稻纹枯病早期检测识别[J]. 光谱学与光谱分析，2019，
39（6）：1898-1904.
ZHU M Y, YANG H B, LI ZH W. Early detection and identification of rice sheath blight disease based on hyperspectralimage and chlorophyll content[J]. Spectroscopy and Spectral Analysis, 2019, 39（6）: 1898-1904. （in Chinese）[5]
谭翠媚, 许廷发, 马旭, 等. 图-谱结合的压缩感知高光谱视频图像复原[J]. 中国光学，2018，11（6）：949-957.
TAN C M, XU T F, MA X, et al.. Graph-spectral hyperspectral video restoration based on compressive sensing[J].Chinese Optics, 2018, 11（6）: 949-957. （in Chinese）[6][7][8]
MAILHES C, VERMANDE P, CASTANIE F. Spectral image compression[J]. Journal of Optics, 1990, 21（3）: 121-LUCEY P G, HORTON K A, WILLIAMS T J, et al.. SMIFTS: a cryogenically cooled, spatially modulated imaging相里斌, 赵葆常, 薛鸣球. 空间调制干涉成像光谱技术[J]. 光学学报，1998，18（1）：18-22.
132.
infrared interferometer spectrometer[J]. Proceedings of SPIE, 1993, 1937: 130-141.
XIANG L B, ZHAO B CH, XUE M Q. Spatially modulated imaging interferometry[J]. Acta Optica Sinica, 1998,18（1）: 18-22. （in Chinese）[9]
高泽东, 郝群, 刘宇, 等. 高光谱成像与应用技术发展[J]. 计测技术，2019，39（4）：24-34.
GAO Z D, HAO Q, LIU Y, et al.. Hyperspectral imaging and application technology development[J]. Metrology &Measurement Technology, 2019, 39（4）: 24-34. （in Chinese）
[10]吕群波, 相里斌, 姚涛, 等. 光学畸变对大孔径静态干涉成像光谱仪影响的建模与仿真[J]. 光谱学与光谱分析，
2010，30（1）：142-145.
LÜ Q B, XIANG L B, YAO T, et al.. Modeling and simulation of effect of optical distortion on the large aperture staticimaging spectrometer[J]. Spectroscopy and Spectral Analysis, 2010, 30（1）: 142-145. （in Chinese）[11]郑玉权. 超光谱成像仪的精细光谱定标[J]. 光学 精密工程，2010，18（11）：2347-2354.
ZHENG Y Q. Precise spectral calibration for hyperspectral imager[J]. Optics and Precision Engineering, 2010, 18（11）:2347-2354. （in Chinese）
[12]郑子鹏, 邱波, 魏诗雅, 等. 曲线距离法的二维光纤光谱弯曲校正[J]. 光谱学与光谱分析，2019，39（10）：3051-3055.
ZHENG Z P, QIU B, WEI SH Y, et al.. Two-dimensional fiber spectral bending correction based on curve distancemethod[J]. Spectroscopy and Spectral Analysis, 2019, 39（10）: 3051-3055. （in Chinese）
[13]于丙文, 金伟, 金钦汉, 等. 基于能量重分配的波长偏移校正方法[J]. 高等学校化学学报，2019，40（8）：1600-1605.
YU B W, JIN W, JIN Q H, et al.. Wavelength drift correction method based on energy redistribution[J]. ChemicalJournal of Chinese Universities, 2019, 40（8）: 1600-1605. （in Chinese）
[14]朱丹彤, 沈宏海, 杨名宇, 等. 编码孔径成像光谱仪中编码元形变的分析校正[J]. 激光与光电子学进展，2018，
55（6）：061201.
ZHU D T, SHEN H H, YANG M Y, et al.. Analysis and correction of coded pixel distortion in coded aperture imagingspectrometer[J]. Laser & Optoelectronics Progress, 2018, 55（6）: 061201. （in Chinese）
[15]董瑛, 相里斌, 赵葆常. 大孔径静态干涉成像光谱仪的干涉系统分析[J]. 光学学报，2001，21（3）：330-334.
DONG Y, XIANG L B, ZHAO B CH. Analysis of interferometer systerm in a large aperture static imagingspectrometer[J]. Acta Optica Sinica, 2001, 21（3）: 330-334. （in Chinese）
[16]董瑛, 相里斌, 赵葆常. 大孔径静态干涉成象光谱仪中的横向剪切干涉仪[J]. 光子学报，1999，28（11）：991-995.
DONG Y, XIANG L B, ZHAO B CH. Lateral shearing interferometer in large aperture static imaging spectrometer[J].Acta Photonica Sinica, 1999, 28（11）: 991-995. （in Chinese）
[17]王爽. 大孔径静态干涉光谱成像仪信噪比研究[D]. 西安: 中国科学院研究生院(西安光学精密机械研究所), 2013.
WANG SH. Research on the signal-to-noise ratio of large aperture static imaging spectrometer[D]. Xi’an: Xi’anInstitute of Optics and Precision Mechanics, Chinese Academy of Sciences, 2013. (in Chinese).
[18]SUN Q CH, HOU Y Q, CHEN J. Lens distortion correction for improving measurement accuracy of digital image
第 2 期
安玲坪, 等: 大孔径静态干涉成像光谱仪径向畸变导致的谱线偏移误差的校正
389
correlation[J]. Optik, 2015, 126（21）: 3153-3157.
[19]邹曜璞. 星载傅里叶光谱仪星上数据处理研究[D]. 上海: 中国科学院研究生院(上海技术物理研究所), 2016.
ZOU Y P. Research of on-board data processing techniques for Fourier transform spectrometer[D]. Shanghai: ShanghaiInstitute of Technical Physics of the Chinese Academy of Sciences, 2016. (in Chinese).
[20]蔡萍, 李潇雁, 唐玉俊, 等. 改进的空间大口径追踪相机畸变校正[J]. 光学精密工程，2019，27（10）：2272-2279.
CAI P, LI X Y, TANG Y J, et al.. Improved distortion correction method for spacial large aperture tracking cameras[J].Optics and Precision Engineering, 2019, 27（10）: 2272-2279. （in Chinese）作者简介：
安玲坪(1996—)，女，山西忻州人，博士研究生，2018年于太原理工大学获得学士学位，主要从事光谱成像仪器定量化技术研究。E-mail：anlingping@opt.cn
王 爽(1980—)，男，黑龙江东宁人，博士，研究员，硕士生导师，2002年、2006年于北京理工大学分别获得学士、硕士学位，2014年于中国科学院西安光学精密机械研究所获得博士学位，主要从事光谱成像理论与仪器开发、光谱成像数据处理与定量化研究。E-mail：wangshuang@opt.ac.cn
因篇幅问题不能全部显示，请点此查看更多更全内容