CN112284537B - 一种推扫式高光谱成像仪光谱弯曲和梯形畸变检测方法 - Google Patents
一种推扫式高光谱成像仪光谱弯曲和梯形畸变检测方法 Download PDFInfo
- Publication number
- CN112284537B CN112284537B CN202010966364.0A CN202010966364A CN112284537B CN 112284537 B CN112284537 B CN 112284537B CN 202010966364 A CN202010966364 A CN 202010966364A CN 112284537 B CN112284537 B CN 112284537B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- hyperspectral imager
- trapezoidal distortion
- detection target
- imaging data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 44
- 238000001228 spectrum Methods 0.000 title claims abstract description 41
- 238000005452 bending Methods 0.000 title claims abstract description 32
- 238000003384 imaging method Methods 0.000 claims abstract description 35
- 230000003595 spectral effect Effects 0.000 claims abstract description 34
- 238000005259 measurement Methods 0.000 claims abstract description 11
- 238000012935 Averaging Methods 0.000 claims abstract description 10
- 238000000034 method Methods 0.000 claims abstract description 9
- 238000012545 processing Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 238000004519 manufacturing process Methods 0.000 claims abstract 2
- 238000010521 absorption reaction Methods 0.000 claims description 8
- 238000012937 correction Methods 0.000 claims description 6
- 230000001678 irradiating effect Effects 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000000701 chemical imaging Methods 0.000 description 2
- 230000032683 aging Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical group [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2823—Imaging spectrometer
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
- G01M11/0242—Testing optical properties by measuring geometrical properties or aberrations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J2003/2866—Markers; Calibrating of scan
- G01J2003/2873—Storing reference spectrum
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Geometry (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本发明公开了一种推扫式高光谱成像仪光谱弯曲和梯形畸变检测方法,该方法步骤如下:(1)制作适用于光谱弯曲检测和梯形畸变检测的参考靶标;(2)暗背景数据测量;(3)光谱弯曲检测靶标成像数据测量;(4)梯形畸变检测靶标成像数据测量;(5)成像数据处理;(6)时间维方向平均;(7)光谱弯曲计算;(8)梯形畸变计算。本发明的以推扫式高光谱成像仪的光谱弯曲和梯形畸变的产生机理为基础,以参考靶标测量的方式,通过一定数据处理获取了高光谱成像仪全幅的光谱弯曲和梯形畸变,在保证测量精度的同时,提高了工作效率。
Description
技术领域
本发明涉及深空探测技术领域,特别是涉及一种推扫式高光谱成像仪光谱弯曲和梯形畸变检测方法。
背景技术
高光谱成像技术是一种图谱合一的成像技术,既可以获得成像目标的二维图像信息,又可以获得目标的光谱信息。高光谱成像技术不仅广泛应用于对地观测领域、深空探测领域。在这些领域,国内外先后实施了多个星载高光谱观测计划,如Hyperion、高分五号、CRISM等。
光栅分光和傅里叶变换是高光谱成像仪常用的两种分光方式,其中光栅分光不涉及严密的运动机构,更加常用。线阵推扫高光谱成像仪是一种最常见的高光谱成像仪,是高光谱成像仪的主流。线阵推扫高光谱成像仪的焦平面行和列分别构成光谱维和一个空间维,通过传感器运动形成另一个空间维,所形成的数据称为高光谱数据立方体。
然而,受光学相差、装配错位、部件老化等因素的影响,线阵推扫式高光谱成像仪经常存在光谱和空间错位问题。光谱维上的错位形成了同一波段(焦平面上的行)不同像元的中心波长和半高全宽不同,称为光谱弯曲。空间维上的错位,形成了同一像元不同波段(焦平面上的列)所对应空间位置不同,称为梯形畸变。两种类型的误配准都会扭曲光谱特征,从而降低了高光谱分类识别、定量化反演等应用的精度。据相关研究表明,高精度的高光谱数据应用需要光谱维错位不确定度小于高光谱成像仪半高全宽的1%,空间错配不确定度小于5%像元。因此,需要在发射前,通过一定的测量手段全面掌握高光谱成像仪的光谱弯曲和梯形畸变特性,为在轨数据的处理提供参考和输入参数支持。通过单色仪和点光源分别对高光谱成像仪进行光谱定标和几何定标是确定光谱弯曲和梯形畸变的常用方式。虽然这种测量方式是金标准,但是要掌握所有像元的畸变情况,采用这种方式耗时耗力。本发明提出了一种采用参考靶标实现快速测量,并通过一定的数据处理方法获得高光谱成像仪的全幅光谱弯曲和梯形畸变的方法。
发明内容
针对现有的技术空白和缺点,本发明所要解决的技术问题是提供一种快速的适用于推扫式高光谱成像仪的全幅光谱弯曲和梯形畸变测量方法。
为了解决上述技术问题,本发明提供的一种推扫式高光谱成像仪光谱弯曲和梯形畸变检测方法,其特点是:
(1)制作适用于光谱弯曲检测和梯形畸变检测的参考靶标.光谱弯曲检测靶标的光谱吸收特征应丰富。吸收特征在光高光谱成像仪响应光谱范围内应分布均匀,保证每5个波段所对应的光谱范围内有一个吸收特征。靶标表面物质分布均匀,反射各项均匀。梯形畸变检测靶标应具有丰富的表面纹理,且纹理相似性低,每5个像元不少于1个特征点。组成梯形畸变检测靶标的物质在高光谱成像仪的不同波段的反射光谱形状应近似直线,每个位置的反射光谱的峰峰值小于0.05;
(2)暗背景数据测量。在模拟的高光谱成像仪在轨运行时的工作环境下,遮住仪器镜头,测量工作环境温度下的暗背景数据;
(3)光谱弯曲检测靶标成像数据测量。将光谱弯曲检测靶标放入高光谱成像仪视场内,并充满视场。用均匀的平行光照射靶标,开启高光谱成像仪,记录光谱弯曲检测靶标的成像数据;
(4)梯形畸变检测靶标成像数据测量。将梯形畸变检测靶标放入高光谱成像仪视场内,并充满视场。用均匀的平行光照射靶标,开启高光谱成像仪,记录梯形畸变检测靶标的成像数据;
(5)成像数据处理。根据暗背景数据、高光谱成像仪的非均匀性校正系数和高光谱成像仪的辐射定标系数对测量得到的光谱弯曲检测靶标成像数据和梯形畸变检测靶标成像数据做非均匀性校正和辐射定标,得到辐射亮度数据;
(6)时间维方向平均。对成像光谱数据在时间维方向上求平均,得到一个经过多帧平均后的空间维×光谱维的图像;
(7)光谱弯曲计算。以空间维×光谱维的图像的中间列为参考,采用同名点搜索算法搜索成像数据每一个像元在光谱维方向上的同名点,所采用的同名点搜索算法精度可以达到亚像元级。
(8)梯形畸变计算。以空间维×光谱维的图像的中间行为参考,采用同名点搜索算法搜索成像数据每一个像元在空间维方向上的同名点,所采用的同名点搜索算法精度可以达到亚像元级。
本发明的有益效果如下:通过参考靶标的方式,通过若干次对靶标的测量,通过一定数据处理可以计算高光谱成像仪全幅的光谱弯曲和梯形畸变,在保证测量精度的同时,提高了工作效率。
附图说明
图1推扫式高光谱成像仪光谱弯曲和梯形畸变检测流程。
具体实施方式
下面结合附图对本发明做进一步的详细描述,但本实施例并不限于本发明,凡是采用本发明的相似方法及其相似变化,均应列入本发明的保护范围。
(1)首先根据高光谱成像仪的光谱响应范围、瞬时视场角、总视场角大小等指标制作尺寸和特性适用的光谱弯曲检测和梯形畸变检测的参考靶标.光谱弯曲检测靶标的光谱吸收特征应丰富。吸收特征在光高光谱成像仪响应光谱范围内应分布均匀,保证每5个波段所对应的光谱范围内有一个吸收特征。靶标表面物质分布均匀,反射各项均匀。梯形畸变检测靶标应具有丰富的表面纹理,且纹理相似性低,每5个像元不少于1个特征点。组成梯形畸变检测靶标的物质在高光谱成像仪的不同波段的反射光谱形状应近似直线,每个位置的反射光谱的峰峰值小于0.05;
(2)暗背景数据测量。在模拟的高光谱成像仪在轨运行时的工作环境下,遮住仪器镜头,测量工作环境温度下的暗背景数据;
(3)光谱弯曲检测靶标成像数据测量。将光谱弯曲检测靶标放入高光谱成像仪视场内,并充满视场。用均匀的平行光照射靶标,开启高光谱成像仪,记录光谱弯曲检测靶标的成像数据;
(4)梯形畸变检测靶标成像数据测量。将梯形畸变检测靶标放入高光谱成像仪视场内,并充满视场。用均匀的平行光照射靶标,开启高光谱成像仪,记录梯形畸变检测靶标的成像数据;
(5)成像数据处理。根据暗背景数据、高光谱成像仪的非均匀性校正系数和高光谱成像仪的辐射定标系数对测量得到的光谱弯曲检测靶标成像数据和梯形畸变检测靶标成像数据做非均匀性校正和辐射定标,得到辐射亮度数据;
(6)时间维方向平均。对成像光谱数据在时间维方向上求平均,得到一个经过多帧平均后的空间维×光谱维的图像;
(7)光谱弯曲计算。以空间维×光谱维的图像的中间列为参考,采用同名点搜索算法搜索成像数据每一个像元在光谱维方向上的同名点,所采用的同名点搜索算法精度可以达到亚像元级。
(8)梯形畸变计算。以空间维×光谱维的图像的中间行为参考,采用同名点搜索算法搜索成像数据每一个像元在空间维方向上的同名点,所采用的同名点搜索算法精度可以达到亚像元级。
Claims (1)
1.一种推扫式高光谱成像仪光谱弯曲和梯形畸变检测方法,其特征在于,包括以下步骤:
(1)制作适用于光谱弯曲检测和梯形畸变检测的参考靶标,光谱弯曲检测靶标的光谱吸收特征应丰富,吸收特征在光高光谱成像仪响应光谱范围内应分布均匀,保证每5个波段所对应的光谱范围内有一个吸收特征,靶标表面物质分布均匀,反射各项均匀,梯形畸变检测靶标应具有丰富的表面纹理,且纹理相似性低,每5个像元不少于1个特征点,组成梯形畸变检测靶标的物质在高光谱成像仪的不同波段的反射光谱形状应近似直线,每个位置的反射光谱的峰峰值小于0.05;
(2)暗背景数据测量,在模拟的高光谱成像仪在轨运行时的工作环境下,遮住仪器镜头,测量工作环境温度下的暗背景数据;
(3)光谱弯曲检测靶标成像数据测量,将光谱弯曲检测靶标放入高光谱成像仪视场内,并充满视场,用均匀的平行光照射靶标,开启高光谱成像仪,记录光谱弯曲检测靶标的成像数据;
(4)梯形畸变检测靶标成像数据测量,将梯形畸变检测靶标放入高光谱成像仪视场内,并充满视场; 用均匀的平行光照射靶标,开启高光谱成像仪,记录梯形畸变检测靶标的成像数据;
(5)成像数据处理,根据暗背景数据、高光谱成像仪的非均匀性校正系数和高光谱成像仪的辐射定标系数对测量得到的光谱弯曲检测靶标成像数据和梯形畸变检测靶标成像数据做非均匀性校正和辐射定标,得到辐射亮度数据;
(6)时间维方向平均,对成像光谱数据在时间维方向上求平均,得到一个经过多帧平均后的空间维×光谱维的图像;
(7)光谱弯曲计算,以空间维×光谱维的图像的中间列为参考,采用同名点搜索算法搜索成像数据每一个像元在光谱维方向上的同名点,所采用的同名点搜索算法精度可以达到亚像元级;
(8)梯形畸变计算,以空间维×光谱维的图像的中间行为参考,采用同名点搜索算法搜索成像数据每一个像元在空间维方向上的同名点,所采用的同名点搜索算法精度达到亚像元级。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010966364.0A CN112284537B (zh) | 2020-09-15 | 2020-09-15 | 一种推扫式高光谱成像仪光谱弯曲和梯形畸变检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010966364.0A CN112284537B (zh) | 2020-09-15 | 2020-09-15 | 一种推扫式高光谱成像仪光谱弯曲和梯形畸变检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112284537A CN112284537A (zh) | 2021-01-29 |
CN112284537B true CN112284537B (zh) | 2022-07-08 |
Family
ID=74420960
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010966364.0A Active CN112284537B (zh) | 2020-09-15 | 2020-09-15 | 一种推扫式高光谱成像仪光谱弯曲和梯形畸变检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112284537B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1990002928A1 (en) * | 1988-09-07 | 1990-03-22 | Sira Limited | Imaging spectrometer |
US6127673A (en) * | 1998-07-09 | 2000-10-03 | Lj Laboratories, L.L.C. | Apparatus and method for detecting curvature |
CN101294867A (zh) * | 2008-06-03 | 2008-10-29 | 清华大学 | Ccd成像设备相对光谱响应特性的标定方法及系统 |
CN102944308A (zh) * | 2012-10-17 | 2013-02-27 | 北京航空航天大学 | 一种时空联合调制干涉成像光谱仪姿态误差校正方法 |
CN207443053U (zh) * | 2017-11-07 | 2018-06-01 | 丹阳丹耀光学有限公司 | 一种tv畸变检测装置 |
CN108709641A (zh) * | 2018-06-19 | 2018-10-26 | 中国科学院上海技术物理研究所 | 一种成像光谱仪畸变补偿结构 |
CN109238991A (zh) * | 2018-11-14 | 2019-01-18 | 中国科学院合肥物质科学研究院 | 一种高光谱大视场成像光谱仪光谱弯曲校正方法 |
CN110794387A (zh) * | 2019-11-28 | 2020-02-14 | 中国科学院合肥物质科学研究院 | 一种机载高光谱成像激光雷达系统的辐射标定方法 |
CN111076816A (zh) * | 2019-12-27 | 2020-04-28 | 中国科学院西安光学精密机械研究所 | 大孔径静态干涉成像光谱仪全视场光谱定标误差校正方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10825152B2 (en) * | 2017-09-14 | 2020-11-03 | Canon U.S.A., Inc. | Distortion measurement and correction for spectrally encoded endoscopy |
-
2020
- 2020-09-15 CN CN202010966364.0A patent/CN112284537B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1990002928A1 (en) * | 1988-09-07 | 1990-03-22 | Sira Limited | Imaging spectrometer |
US6127673A (en) * | 1998-07-09 | 2000-10-03 | Lj Laboratories, L.L.C. | Apparatus and method for detecting curvature |
CN101294867A (zh) * | 2008-06-03 | 2008-10-29 | 清华大学 | Ccd成像设备相对光谱响应特性的标定方法及系统 |
CN102944308A (zh) * | 2012-10-17 | 2013-02-27 | 北京航空航天大学 | 一种时空联合调制干涉成像光谱仪姿态误差校正方法 |
CN207443053U (zh) * | 2017-11-07 | 2018-06-01 | 丹阳丹耀光学有限公司 | 一种tv畸变检测装置 |
CN108709641A (zh) * | 2018-06-19 | 2018-10-26 | 中国科学院上海技术物理研究所 | 一种成像光谱仪畸变补偿结构 |
CN109238991A (zh) * | 2018-11-14 | 2019-01-18 | 中国科学院合肥物质科学研究院 | 一种高光谱大视场成像光谱仪光谱弯曲校正方法 |
CN110794387A (zh) * | 2019-11-28 | 2020-02-14 | 中国科学院合肥物质科学研究院 | 一种机载高光谱成像激光雷达系统的辐射标定方法 |
CN111076816A (zh) * | 2019-12-27 | 2020-04-28 | 中国科学院西安光学精密机械研究所 | 大孔径静态干涉成像光谱仪全视场光谱定标误差校正方法 |
Non-Patent Citations (2)
Title |
---|
孙红胜.紫外-真空紫外成像光谱仪校准技术研究.《应用光学》.2020, * |
杨振杰.推扫式高光谱影像的相对辐射定标研究.《青海大学学报》.2018, * |
Also Published As
Publication number | Publication date |
---|---|
CN112284537A (zh) | 2021-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Childress et al. | PyWiFeS: a rapid data reduction pipeline for the Wide Field Spectrograph (WiFeS) | |
Lenhard et al. | Independent laboratory characterization of NEO HySpex imaging spectrometers VNIR-1600 and SWIR-320m-e | |
Fomalont et al. | The radio/optical catalog of the SSA 13 field | |
CN106768391B (zh) | 成像仪不同焦面谱段配准精度测试方法 | |
CN103323114B (zh) | 一种傅里叶变换红外成像光谱仪快速光谱定标校正方法 | |
Lin et al. | NPP VIIRS geometric performance status | |
CN115265803A (zh) | 一种远距离点目标红外光谱辐射特性测量方法 | |
CN112798013A (zh) | 一种对光学载荷在轨绝对辐射定标结果进行验证的方法 | |
Barducci et al. | Developing a new hyperspectral imaging interferometer for earth observation | |
CN112284537B (zh) | 一种推扫式高光谱成像仪光谱弯曲和梯形畸变检测方法 | |
US10656017B2 (en) | On-board processing of hyperspectral data | |
US5675513A (en) | Method of calibrating an interferometer and reducing its systematic noise | |
CN112834051B (zh) | 一种基于多光谱相机的辐射测温法 | |
Nieke et al. | Uniformity of imaging spectrometry data products | |
CN110702228B (zh) | 一种航空高光谱影像的边缘辐射校正方法 | |
CN113405660B (zh) | 一种水下光谱成像设备下潜在轨光谱辐射定标装置及其定标方法 | |
Kouyama et al. | One-year lunar calibration result of Hodoyoshi-1, Moon as an ideal target for small satellite radiometric calibration | |
Cai et al. | Precise reduction of solar spectra observed by the one-meter new vacuum solar telescope | |
Wang et al. | On-orbit modulation transfer function characterization of terra MODIS using the moon | |
Bakker et al. | Determining smile and keystone of lab hyperspectral line cameras | |
Zhou et al. | Smile effect detection for dispersive hypersepctral imager based on the doped reflectance panel | |
Ungermann et al. | Quantification and mitigation of the instrument effects and uncertainties of the airborne limb imaging FTIR GLORIA | |
Pieper et al. | Sensitivity of temperature and emissivity separation to atmospheric errors in LWIR hyperspectral imagery | |
Mouroulis et al. | Spectral response evaluation and computation for pushbroom imaging spectrometers | |
Bali et al. | Comparisons of IASI-A and AATSR measurements of top-of-atmosphere radiance over an extended period |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |