CN117168619B - 一种星载高光谱成像仪光谱定标方法及*** - Google Patents
一种星载高光谱成像仪光谱定标方法及*** Download PDFInfo
- Publication number
- CN117168619B CN117168619B CN202311443081.8A CN202311443081A CN117168619B CN 117168619 B CN117168619 B CN 117168619B CN 202311443081 A CN202311443081 A CN 202311443081A CN 117168619 B CN117168619 B CN 117168619B
- Authority
- CN
- China
- Prior art keywords
- satellite
- hyperspectral imager
- wavelength
- calibrated
- calibration
- 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
- 238000001228 spectrum Methods 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 58
- 230000003595 spectral effect Effects 0.000 claims abstract description 121
- 238000005316 response function Methods 0.000 claims abstract description 42
- 238000003384 imaging method Methods 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 7
- 238000005286 illumination Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 6
- 230000003287 optical effect Effects 0.000 description 9
- 238000004891 communication Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000006467 substitution reaction Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 238000000701 chemical imaging Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
Landscapes
- Spectrometry And Color Measurement (AREA)
Abstract
本发明提供一种星载高光谱成像仪光谱定标方法及***，属于遥感光学技术领域，包括：通过单色仪单色平行光照射获取待标定星载高光谱成像仪的相对光谱辐照度分布，采集待标定星载高光谱成像仪的定标影像；由待标定星载高光谱成像仪的起始波长和终止波长，得到多个标定影像；利用多个定标影像，确定各中心波长与定标影像中的最大值所对应行号；基于相对光谱辐照度分布、多个定标影像和最大值所对应行号，计算待标定星载高光谱成像仪的相对光谱响应函数；根据相对光谱响应函数，确定待标定星载高光谱成像仪的各中心波长对应的半最大值全波。本发明获得遥感卫星定量化应用必需的光谱定标参数，填补了线性渐变滤光片式高光谱成像仪光谱定标方法的空白。
Description
技术领域
本发明涉及遥感光学技术领域，尤其涉及一种星载高光谱成像仪光谱定标方法及***。
背景技术
遥感星载高光谱成像仪通常具有几十至上百个光谱通道，其光谱分辨率可以达到中心波长的百分之一量级，具有信息量大、识别效果好等优点。就现阶段而言，常用的几种光谱分光方式主要包括色散型、干涉型以及滤光片型。相较于色散型或干涉型，滤光片型的优势在于光路紧凑、光学***中所包含的光学元器件数量少，重量轻，成本可控。其中线性渐变滤光片型具有以下两个显著的优点：光谱通道数多，通常可达几十至上百个；光谱分辨率高，通常可优于10nm，因此得到了广泛的应用。
现有的高光谱成像仪光谱定标方法缺乏对中心波长与成像器件位置对应关系的标定流程，且对于相对光谱响应函数的测算未涉及光谱通道数多且中心波长可变的情况。因此现有方法适用于色散型或干涉型高光谱成像仪，但无法满足光谱通道数多、通道中心与宽度可变的线性渐变滤光片型高光谱成像仪。
因此，针对线性渐变滤光片型高光谱成像仪，需要提出一种新的光谱定标方法。
发明内容
本发明提供一种星载高光谱成像仪光谱定标方法及***，用以解决现有技术中缺乏针对线性渐变滤光片型高光谱成像仪光谱定标方法的缺陷。
第一方面，本发明提供一种星载高光谱成像仪光谱定标方法，包括：
S1，在暗背景成像模式下采集待标定星载高光谱成像仪的预设数量的暗背景影像数据；
S2，通过单色仪的单色平行光照射获取所述待标定星载高光谱成像仪的相对光谱辐照度分布，采集所述待标定星载高光谱成像仪的定标影像；
S3，由所述待标定星载高光谱成像仪的起始波长和终止波长，得到多个标定影像；
S4，利用所述多个定标影像，确定各中心波长与定标影像中的最大值所对应行号；
S5，基于所述相对光谱辐照度分布、所述多个定标影像和所述最大值所对应行号，计算所述待标定星载高光谱成像仪的相对光谱响应函数；
S6，根据所述相对光谱响应函数，确定所述待标定星载高光谱成像仪的各中心波长对应的半最大值全波。
根据本发明提供的一种星载高光谱成像仪光谱定标方法，步骤S1之前还包括：
获取待标定星载高光谱成像仪的原始器件参数。
根据本发明提供的一种星载高光谱成像仪光谱定标方法，步骤S2包括：
所述单色仪接收宽谱照明光源之后输出单色光；
将所述单色光经过平行光管照射得到所述单色平行光；
由所述单色平行光照射所述待标定星载高光谱成像仪得到所述相对光谱辐照度分布；
以预设次数采集测量所述相对光谱辐照度分布，获得所述定标影像。
根据本发明提供的一种星载高光谱成像仪光谱定标方法，步骤S3包括：
从所述待标定星载高光谱成像仪的起始波长开始，以所述单色仪的预设中心波长变化步距重复步骤S2，直至达到所述待标定星载高光谱成像仪的终止波长，得到所述多个定标影像。
根据本发明提供的一种星载高光谱成像仪光谱定标方法，步骤S4包括：
S401，确定所述多个定标影像中任一中心波长定标影像，逐行计算所述任一中心波长定标影像的平均遥感影像像元亮度值，提取所述平均遥感影像像元亮度值中最大值所对应行号，基于所述最大值所对应行号确定任一中心波长输出行；
S402，逐个中心波长重复步骤S401，通过最小二乘拟合得到不同中心波长对应的中心波长输出行。
根据本发明提供的一种星载高光谱成像仪光谱定标方法，步骤S5包括：
S501，计算任一中心波长定标影像的平均遥感影像像元亮度值与对应单色光的相对光谱辐照度分布的比值，由所述比值除以所述比值在任一中心波长光谱通道波长下限与任一中心波长光谱通道波长上限之间的最大值，得到任一中心波长的相对光谱响应函数；
S502，采用最小二乘法对所述任一中心波长光谱通道的相对光谱响应函数进行拟合，得到所述待标定星载高光谱成像仪在任一中心波长对应光谱通道的相对光谱响应函数；
S503，逐个中心波长重复步骤S501和步骤S502，得到不同中心波长对应光谱通道的相对光谱响应函数。
根据本发明提供的一种星载高光谱成像仪光谱定标方法，步骤S6包括：
S601，确定任一中心波长对应光谱通道的相对光谱响应函数中的峰值响应二分之一且位于短波处的短波波长，以及位于长波处的长波波长，由所述长波波长减去所述短波波长，得到任一中心波长对应光谱通道的半最大值全波；
S602，逐个中心波长重复步骤S601，得到不同中心波长对应光谱通道的半最大值全波。
第二方面，本发明还提供一种星载高光谱成像仪光谱定标***，包括：
采集模块，用于在暗背景成像模式下采集待标定星载高光谱成像仪的预设数量的暗背景影像数据；
获取模块，用于通过单色仪的单色平行光照射获取所述待标定星载高光谱成像仪的相对光谱辐照度分布，采集所述待标定星载高光谱成像仪的定标影像；
标定模块，用于由所述待标定星载高光谱成像仪的起始波长和终止波长，得到多个标定影像；
确定模块，用于利用所述多个定标影像，确定各中心波长与定标影像中的最大值所对应行号；
计算模块，用于基于所述相对光谱辐照度分布、所述多个定标影像和所述最大值所对应行号，计算所述待标定星载高光谱成像仪的相对光谱响应函数；
综合模块，用于根据所述相对光谱响应函数，确定所述待标定星载高光谱成像仪的各中心波长对应的半最大值全波。
第三方面，本发明还提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上述任一种所述星载高光谱成像仪光谱定标方法。
第四方面，本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上述任一种所述星载高光谱成像仪光谱定标方法。
本发明提供的星载高光谱成像仪光谱定标方法及***，通过获得遥感卫星定量化应用必需的光谱定标参数，填补了线性渐变滤光片式高光谱成像仪光谱定标方法的空白，避免在定标影像数据获取过程中出现过曝等现象，确保定标数据的有效性；并且通过多次采集与计算避免偶然误差，确保定标数据的准确性，并采用中心波长变化步距，减少了重复操作的次数，提高了光谱定标的效率；定标过程中均使用常见光学分析仪器，如平行光管、宽谱照明光源等，降低了实验室光谱定标的难度。因此，本发明所获得的实验室光谱定标参数是后续开展发射前实验室绝对辐射定标与卫星在轨替代定标的重要基础，且对高光谱遥感卫星定量化应用具有重要意义。
附图说明
为了更清楚地说明本发明或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明提供的星载高光谱成像仪光谱定标方法的流程示意图之一；
图2是本发明提供的星载高光谱成像仪光谱定标方法的流程示意图之二；
图3是本发明提供的星载高光谱成像仪光谱定标***的结构示意图；
图4是本发明提供的电子设备的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
由于现有技术中普遍缺乏针对线性渐变滤光片式高光谱成像仪的光谱标定方法，本发明拟定了一个完整的实验室光谱定标方法。
图1是本发明实施例提供的星载高光谱成像仪光谱定标方法的流程示意图之一，如图1所示，包括：
S1，在暗背景成像模式下采集待标定星载高光谱成像仪的预设数量的暗背景影像数据；
S2，通过单色仪的单色平行光照射获取所述待标定星载高光谱成像仪的相对光谱辐照度分布，采集所述待标定星载高光谱成像仪的定标影像；
S3，由所述待标定星载高光谱成像仪的起始波长和终止波长，得到多个标定影像；
S4，利用所述多个定标影像，确定各中心波长与定标影像中的最大值所对应行号；
S5，基于所述相对光谱辐照度分布、所述多个定标影像和所述最大值所对应行号，计算所述待标定星载高光谱成像仪的相对光谱响应函数；
S6，根据所述相对光谱响应函数，确定所述待标定星载高光谱成像仪的各中心波长对应的半最大值全波。
具体地，如图2所示，本发明实施例首先使待定标高光谱成像仪工作在暗背景成像模式，采集若干次暗背景影像数据，例如50次，用于消除高光谱成像仪由于暗电流所产生的背景噪声；
由宽谱照明光源经单色仪输出单色光，经平行光管照射待定标高光谱成像仪，使用光谱辐射计测量单色平行光的相对光谱辐照度分布，并采集待定标高光谱成像仪输出的定标影像；
从待定标高光谱成像仪的起始波长开始，以1nm作为单色仪输出中心波长的变化步距并重复前述步骤，直至待定标高光谱成像仪的终止波长；
利用前述步骤中所采集的待定标高光谱成像仪的影像，确定各中心波长与其定标影像上最大值所对应的行号；
再利用前述步骤中所采集的待定标高光谱成像仪的影像与相对光谱辐照度分布，计算待定标高光谱成像仪的相对光谱响应函数；
最后根据得到的相对光谱响应函数，确定各中心波长所对应的半最大值全波（Full Width at Half Maximum，FWHM）。
本发明针对线性渐变滤光片式高光谱成像仪光谱通道数多、通道中心与宽度可变的特点，设计了一种实验室光谱定标方法，能够获得遥感卫星定量化应用必需的光谱定标参数，填补了线性渐变滤光片式高光谱成像仪光谱定标方法的空白。
在上述实施例的基础上，步骤S1之前还包括：
获取待标定星载高光谱成像仪的原始器件参数。
具体地，在正式标定之前，需要先获取线性渐变滤光片式高光谱成像仪的诸多原始器件参数，如带通透过率等光谱定标参数。
在上述实施例的基础上，步骤S2包括：
所述单色仪接收宽谱照明光源之后输出单色光；
将所述单色光经过平行光管照射得到所述单色平行光；
由所述单色平行光照射所述待标定星载高光谱成像仪得到所述相对光谱辐照度分布；
以预设次数采集测量所述相对光谱辐照度分布，获得所述定标影像。
具体地，由宽谱照明光源经单色仪输出单色光，然后经平行光管照射待定标高光谱成像仪。使用光谱辐射计测量单色平行光的相对光谱辐照度分布，每次应至少测量50次，并至少采集50次待定标高光谱成像仪输出的定标影像。
本发明提出的高光谱成像仪成像参数的设置方法，可避免在定标影像数据获取过程中出现过曝等现象，确保定标数据的有效性；并且通过多次采集与计算避免偶然误差，确保定标数据的准确性。
在上述实施例的基础上，步骤S3包括：
从所述待标定星载高光谱成像仪的起始波长开始，以所述单色仪的预设中心波长变化步距重复步骤S2，直至达到所述待标定星载高光谱成像仪的终止波长，得到所述多个定标影像。
具体地，本发明实施例为覆盖所有波长段的场景，设定从待定标高光谱成像仪的起始波长开始，按照预设中心波长变化步距，通常为1nm作为单色仪输出中心波长的变化步距重复步骤S2，直至达到待定标高光谱成像仪的终止波长。
本发明所设计的实验室光谱定标方法将1nm作为单色仪输出中心波长的变化步距，减少了重复操作的次数，提高了光谱定标的效率；定标过程中均使用常见光学分析仪器，如平行光管、宽谱照明光源等，降低了实验室光谱定标的难度。
在上述实施例的基础上，步骤S4包括：
S401，确定所述多个定标影像中任一中心波长定标影像，逐行计算所述任一中心波长定标影像的平均遥感影像像元亮度值，提取所述平均遥感影像像元亮度值中最大值所对应行号，基于所述最大值所对应行号确定任一中心波长输出行；
S402，逐个中心波长重复步骤S401，通过最小二乘拟合得到不同中心波长对应的中心波长输出行。
具体地，本发明实施例利用步骤S2、S3中所采集的待定标高光谱成像仪的影像，确定各中心波长与其定标影像上最大值所对应的行号，具体包括：
S401，取中心波长为的定标影像，逐行计算平均遥感影像像元亮度值值（Digital Number，DN），并比较平均DN值大小，取平均DN值最大的一行作为中心波长/>的输出行/>；
S402，逐中心波长重复步骤S401，得到不同中心波长对应的输出行；并通过最小二乘拟合得到二者的关系为：
= P（/>）
利用该关系即可计算得到任意中心波长所对应的成像器件位置。
在上述实施例的基础上，步骤S5包括：
S501，计算任一中心波长定标影像的平均遥感影像像元亮度值与对应单色光的相对光谱辐照度分布的比值，由所述比值除以所述比值在任一中心波长光谱通道波长下限与任一中心波长光谱通道波长上限之间的最大值，得到任一中心波长的相对光谱响应函数；
S502，采用最小二乘法对所述任一中心波长光谱通道的相对光谱响应函数进行拟合，得到所述待标定星载高光谱成像仪在任一中心波长对应光谱通道的相对光谱响应函数；
S503，逐个中心波长重复步骤S501和步骤S502，得到不同中心波长对应光谱通道的相对光谱响应函数。
具体地，本发明实施例利用步骤S2、S3中所采集的待定标高光谱成像仪的影像与相对光谱辐照度分布，计算待定标高光谱成像仪的相对光谱响应函数，具体包括：
S501，根据单色平行光的相对光谱辐照度分布，按照下式计算待定标高光谱成像仪当前光谱通道的相对光谱响应函数R（λ）：
R（λ） =
式中λ为中心波长，R（λ）为对应的归一化光谱响应，DN（λ）为对应的输出DN值，将每一行上所有像素的平均DN值作为输出DN值，RL（λ）为对应单色光的相对光谱辐照度分布，为待定标高光谱成像仪该光谱通道波长下限，/>为待定标高光谱成像仪该光谱通道波长上限；
S502，利用最小二乘法进行拟合，得到待定标高光谱成像仪该光谱通道的相对光谱响应函数R（λ）；
S503，逐中心波长重复步骤S501、S502，得到不同中心波长所对应的相对光谱响应函数R（λ）。
在上述实施例的基础上，步骤S6包括：
S601，确定任一中心波长对应光谱通道的相对光谱响应函数中的峰值响应二分之一且位于短波处的短波波长，以及位于长波处的长波波长，由所述长波波长减去所述短波波长，得到任一中心波长对应光谱通道的半最大值全波；
S602，逐个中心波长重复步骤S601，得到不同中心波长对应光谱通道的半最大值全波。
具体地，本发明实施例根据步骤S5中计算得到的相对光谱响应函数，确定各中心波长所对应的FWHM，具体包括：
S601，中心波长的相对光谱响应函数为R（/>），将响应值为峰值响应二分之一且位于短波处的波长记为/>、位于长波处的波长记为/>，则此时半波宽FWHM可以由下式计算：
FWHM =-/>
S602，逐中心波长重复步骤S601，即可得到不同中心波长所对应的FWHM。
本发明中所获得的实验室光谱定标参数是后续开展发射前实验室绝对辐射定标与卫星在轨替代定标的重要基础，且对高光谱遥感卫星定量化应用具有重要意义。
下面对本发明提供的星载高光谱成像仪光谱定标***进行描述，下文描述的星载高光谱成像仪光谱定标***与上文描述的星载高光谱成像仪光谱定标方法可相互对应参照。
图3是本发明实施例提供的星载高光谱成像仪光谱定标***的结构示意图，如图3所示，包括：采集模块31、获取模块32、标定模块33、确定模块34、计算模块35和综合模块36，其中：
采集模块31用于在暗背景成像模式下采集待标定星载高光谱成像仪的预设数量的暗背景影像数据；获取模块32用于通过单色仪的单色平行光照射获取所述待标定星载高光谱成像仪的相对光谱辐照度分布，采集所述待标定星载高光谱成像仪的定标影像；标定模块33用于由所述待标定星载高光谱成像仪的起始波长和终止波长，得到多个标定影像；确定模块34用于利用所述多个定标影像，确定各中心波长与定标影像中的最大值所对应行号；计算模块35用于基于所述相对光谱辐照度分布、所述多个定标影像和所述最大值所对应行号，计算所述待标定星载高光谱成像仪的相对光谱响应函数；综合模块36用于根据所述相对光谱响应函数，确定所述待标定星载高光谱成像仪的各中心波长对应的半最大值全波。
图4示例了一种电子设备的实体结构示意图，如图4所示，该电子设备可以包括：处理器(processor)410、通信接口(Communications Interface)420、存储器(memory)430和通信总线440，其中，处理器410，通信接口420，存储器430通过通信总线440完成相互间的通信。处理器410可以调用存储器430中的逻辑指令，以执行星载高光谱成像仪光谱定标方法，该方法包括：在暗背景成像模式下采集待标定星载高光谱成像仪的预设数量的暗背景影像数据；通过单色仪的单色平行光照射获取所述待标定星载高光谱成像仪的相对光谱辐照度分布，采集所述待标定星载高光谱成像仪的定标影像；由所述待标定星载高光谱成像仪的起始波长和终止波长，得到多个标定影像；利用所述多个定标影像，确定各中心波长与定标影像中的最大值所对应行号；基于所述相对光谱辐照度分布、所述多个定标影像和所述最大值所对应行号，计算所述待标定星载高光谱成像仪的相对光谱响应函数；根据所述相对光谱响应函数，确定所述待标定星载高光谱成像仪的各中心波长对应的半最大值全波。
此外，上述的存储器430中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器（ROM，Read-Only Memory）、随机存取存储器（RAM，Random Access Memory）、磁碟或者光盘等各种可以存储程序代码的介质。
另一方面，本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各方法提供的星载高光谱成像仪光谱定标方法，该方法包括：在暗背景成像模式下采集待标定星载高光谱成像仪的预设数量的暗背景影像数据；通过单色仪的单色平行光照射获取所述待标定星载高光谱成像仪的相对光谱辐照度分布，采集所述待标定星载高光谱成像仪的定标影像；由所述待标定星载高光谱成像仪的起始波长和终止波长，得到多个标定影像；利用所述多个定标影像，确定各中心波长与定标影像中的最大值所对应行号；基于所述相对光谱辐照度分布、所述多个定标影像和所述最大值所对应行号，计算所述待标定星载高光谱成像仪的相对光谱响应函数；根据所述相对光谱响应函数，确定所述待标定星载高光谱成像仪的各中心波长对应的半最大值全波。
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (9)
1.一种星载高光谱成像仪光谱定标方法，其特征在于，包括：
S1，在暗背景成像模式下采集待标定星载高光谱成像仪的预设数量的暗背景影像数据；
S2，通过单色仪的单色平行光照射获取所述待标定星载高光谱成像仪的相对光谱辐照度分布，采集所述待标定星载高光谱成像仪的定标影像；
S3，由所述待标定星载高光谱成像仪的起始波长和终止波长，得到多个定标影像；
S4，利用所述多个定标影像，确定各中心波长与定标影像中平均遥感影像像元亮度值的最大值所对应行号；
S5，基于所述相对光谱辐照度分布、所述多个定标影像和所述最大值所对应行号，计算所述待标定星载高光谱成像仪的相对光谱响应函数；
S6，根据所述相对光谱响应函数，确定所述待标定星载高光谱成像仪的各中心波长对应的半最大值全波；
其中，步骤S3包括：
从所述待标定星载高光谱成像仪的起始波长开始，以所述单色仪的预设中心波长变化步距重复步骤S2，直至达到所述待标定星载高光谱成像仪的终止波长，得到所述多个定标影像。
2.根据权利要求1所述的星载高光谱成像仪光谱定标方法，其特征在于，步骤S1之前还包括：
获取待标定星载高光谱成像仪的原始器件参数。
3.根据权利要求1所述的星载高光谱成像仪光谱定标方法，其特征在于，步骤S2包括：
所述单色仪接收宽谱照明光源之后输出单色光；
将所述单色光经过平行光管照射得到所述单色平行光；
由所述单色平行光照射所述待标定星载高光谱成像仪得到所述相对光谱辐照度分布；
以预设次数采集测量所述相对光谱辐照度分布，获得所述定标影像。
4.根据权利要求1所述的星载高光谱成像仪光谱定标方法，其特征在于，步骤S4包括：
S401，确定所述多个定标影像中任一中心波长定标影像，逐行计算所述任一中心波长定标影像的平均遥感影像像元亮度值，提取所述平均遥感影像像元亮度值中最大值所对应行号，基于所述最大值所对应行号确定任一中心波长输出行；
S402，逐个中心波长重复步骤S401，通过最小二乘拟合得到不同中心波长对应的中心波长输出行。
5.根据权利要求1所述的星载高光谱成像仪光谱定标方法，其特征在于，步骤S5包括：
S501，计算任一中心波长定标影像的平均遥感影像像元亮度值与对应单色光的相对光谱辐照度分布的比值，由所述比值除以所述比值在任一中心波长光谱通道波长下限与任一中心波长光谱通道波长上限之间的最大值，得到任一中心波长的相对光谱响应函数；
S502，采用最小二乘法对所述任一中心波长光谱通道的相对光谱响应函数进行拟合，得到所述待标定星载高光谱成像仪在任一中心波长对应光谱通道的相对光谱响应函数；
S503，逐个中心波长重复步骤S501和步骤S502，得到不同中心波长对应光谱通道的相对光谱响应函数。
6.根据权利要求1所述的星载高光谱成像仪光谱定标方法，其特征在于，步骤S6包括：
S601，确定任一中心波长对应光谱通道的相对光谱响应函数中的峰值响应二分之一且位于短波处的短波波长，以及位于长波处的长波波长，由所述长波波长减去所述短波波长，得到任一中心波长对应光谱通道的半最大值全波；
S602，逐个中心波长重复步骤S601，得到不同中心波长对应光谱通道的半最大值全波。
7.一种星载高光谱成像仪光谱定标***，其特征在于，包括：
采集模块，用于在暗背景成像模式下采集待标定星载高光谱成像仪的预设数量的暗背景影像数据；
获取模块，用于通过单色仪的单色平行光照射获取所述待标定星载高光谱成像仪的相对光谱辐照度分布，采集所述待标定星载高光谱成像仪的定标影像；
标定模块，用于由所述待标定星载高光谱成像仪的起始波长和终止波长，得到多个定标影像；
确定模块，用于利用所述多个定标影像，确定各中心波长与定标影像中平均遥感影像像元亮度值的最大值所对应行号；
计算模块，用于基于所述相对光谱辐照度分布、所述多个定标影像和所述最大值所对应行号，计算所述待标定星载高光谱成像仪的相对光谱响应函数；
综合模块，用于根据所述相对光谱响应函数，确定所述待标定星载高光谱成像仪的各中心波长对应的半最大值全波；
所述标定模块具体用于：
从所述待标定星载高光谱成像仪的起始波长开始，以所述单色仪的预设中心波长变化步距重复所述获取模块中的步骤，直至达到所述待标定星载高光谱成像仪的终止波长，得到所述多个定标影像。
8.一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至6任一项所述星载高光谱成像仪光谱定标方法。
9.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至6任一项所述星载高光谱成像仪光谱定标方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311443081.8A CN117168619B (zh) | 2023-11-02 | 2023-11-02 | 一种星载高光谱成像仪光谱定标方法及*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311443081.8A CN117168619B (zh) | 2023-11-02 | 2023-11-02 | 一种星载高光谱成像仪光谱定标方法及*** |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117168619A CN117168619A (zh) | 2023-12-05 |
| CN117168619B true CN117168619B (zh) | 2024-02-02 |
Family
ID=88947152
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311443081.8A Active CN117168619B (zh) | 2023-11-02 | 2023-11-02 | 一种星载高光谱成像仪光谱定标方法及*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117168619B (zh) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6002990A (en) * | 1997-10-16 | 1999-12-14 | Datex-Ohmeda, Inc. | Dynamic wavelength calibration for spectrographic analyzer |
| CN103698010A (zh) * | 2013-12-31 | 2014-04-02 | 中国科学院光电研究院 | 一种微型化线性渐变滤光片型成像光谱仪 |
| CN104457985A (zh) * | 2014-09-15 | 2015-03-25 | 中国资源卫星应用中心 | 一种色散型高光谱成像遥感器的在轨光谱定标方法 |
| CN105424186A (zh) * | 2015-11-04 | 2016-03-23 | 北京航空航天大学 | 一种光场成像光谱仪的光谱标定校正方法 |
| CN108007569A (zh) * | 2017-11-23 | 2018-05-08 | 首都师范大学 | 一种多光谱影像辐射定标方法 |
| CN109238465A (zh) * | 2018-10-10 | 2019-01-18 | 中国科学院上海技术物理研究所 | 一种适用于星载宽谱段相机的光谱定标*** |
| CN114739507A (zh) * | 2022-04-12 | 2022-07-12 | 中国科学院合肥物质科学研究院 | 一种高光谱辐亮度仪杂散光校正方法及测量装置 |
| WO2023159739A1 (zh) * | 2022-02-28 | 2023-08-31 | 中国科学院空天信息创新研究院 | 基于临近空间浮空器的光学卫星遥感传递定标方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002016905A2 (en) * | 2000-08-21 | 2002-02-28 | Euro-Celtique, S.A. | Near infrared blood glucose monitoring system |
| WO2021093676A1 (zh) * | 2019-11-14 | 2021-05-20 | 中国科学院上海技术物理研究所启东光电遥感中心 | 用于处理高光谱图像的方法 |
-
2023
- 2023-11-02 CN CN202311443081.8A patent/CN117168619B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6002990A (en) * | 1997-10-16 | 1999-12-14 | Datex-Ohmeda, Inc. | Dynamic wavelength calibration for spectrographic analyzer |
| CN103698010A (zh) * | 2013-12-31 | 2014-04-02 | 中国科学院光电研究院 | 一种微型化线性渐变滤光片型成像光谱仪 |
| CN104457985A (zh) * | 2014-09-15 | 2015-03-25 | 中国资源卫星应用中心 | 一种色散型高光谱成像遥感器的在轨光谱定标方法 |
| CN105424186A (zh) * | 2015-11-04 | 2016-03-23 | 北京航空航天大学 | 一种光场成像光谱仪的光谱标定校正方法 |
| CN108007569A (zh) * | 2017-11-23 | 2018-05-08 | 首都师范大学 | 一种多光谱影像辐射定标方法 |
| CN109238465A (zh) * | 2018-10-10 | 2019-01-18 | 中国科学院上海技术物理研究所 | 一种适用于星载宽谱段相机的光谱定标*** |
| WO2023159739A1 (zh) * | 2022-02-28 | 2023-08-31 | 中国科学院空天信息创新研究院 | 基于临近空间浮空器的光学卫星遥感传递定标方法 |
| CN114739507A (zh) * | 2022-04-12 | 2022-07-12 | 中国科学院合肥物质科学研究院 | 一种高光谱辐亮度仪杂散光校正方法及测量装置 |
Non-Patent Citations (1)
| Title |
|---|
| 顾及成像几何的国产遥感卫星交叉辐射定标方法研究;郭霏霏;《中国优秀硕士论文全文数据库 基础科学辑》;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117168619A (zh) | 2023-12-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11566941B2 (en) | Systems and methods for calibrating, configuring and validating an imaging device or system for multiplex tissue assays | |
| Kleipool et al. | Pre-launch calibration results of the TROPOMI payload on-board the Sentinel-5 Precursor satellite | |
| Mann et al. | Revised filter profiles and zero points for broadband photometry | |
| EP1303742B1 (en) | Spectral drift and correction technique for hyperspectral imaging systems | |
| Bertello et al. | Correlation between sunspot number and Ca II K emission index | |
| JP2016510408A5 (zh) | ||
| US20080118181A1 (en) | Method and Apparatus for Reconstructing Optical Spectra in a Static Multimode Multiplex Spectrometer | |
| US10323985B2 (en) | Signal processing for tunable Fabry-Perot interferometer based hyperspectral imaging | |
| JP2012242271A (ja) | 分光推定パラメーター生成装置およびその方法、分光推定装置、並びにコンピュータープログラム | |
| Katrašnik et al. | Radiometric calibration and noise estimation of acousto-optic tunable filter hyperspectral imaging systems | |
| Kobelski et al. | Calibrating data from the Hinode/X-Ray Telescope and associated uncertainties | |
| CN117168618B (zh) | 一种星载高光谱成像仪辐射定标方法及*** | |
| Zhang et al. | Spectral shift correction and adaptive band selection for multispectral imaging | |
| CN117168619B (zh) | 一种星载高光谱成像仪光谱定标方法及*** | |
| Jemec et al. | Deconvolution-based restoration of SWIR pushbroom imaging spectrometer images | |
| EP3327411A1 (en) | Optical system | |
| Lenhard et al. | Impact of improved calibration of a NEO HySpex VNIR-1600 sensor on remote sensing of water depth | |
| Jehle et al. | Detection and correction of radiance variations during spectral calibration in APEX | |
| Gilchrist et al. | Developing the IEEE P4001 standard for characterisation and calibration of hyperspectral imaging devices | |
| Rice et al. | Hyperspectral image compressive projection algorithm | |
| Statella et al. | Cross-calibration of the Rosetta Navigation Camera based on images of the 67P comet nucleus | |
| CN118129905B (zh) | 一种星载高光谱成像仪曲线型辐射定标方法及*** | |
| Ekpenyong | Calibration and Characterization of Hyperspectral Imaging Systems Used for Natural Scene Imagery | |
| Milton et al. | Estimating the spectral response function of the CASI-2 | |
| US20040218183A1 (en) | Method of determining an image recording device characterisation |
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 |