一种光场成像光谱仪的光谱标定校正方法
【技术领域】
[0001] 本发明涉及光学成像领域，具体涉及一种光场成像光谱仪的光谱标定校正方法。
【背景技术】
[0002] 成像光谱仪是作为遥感探测技术之一，可以获取的景物或目标光谱信息，在目标 识别、特征提取和精确分类等领域具有广泛的应用价值。具有图谱合一的优势，被搭载在航 空、航天器上对陆地、大气、海洋进行观测。根据光谱三维数据立方体的获取方式不同，可以 分为扫描型和快照。扫描型成像光谱技术在空间或光谱维进行扫描获取完整数据立方体， 而快照型成像光谱技术一次曝光即获取目标完整数据立方体，在动态目标监测和追踪方面 快照型成像光谱技术具有优势和广泛的应用前景。另一方面，成像光谱仪的应用以遥感信 息的定量化为基础，通过实验或物理模型将采集到的遥感信息和观察目标参量联系起来， 并反演或推算出地学、生物学和大气等观测目标参量。遥感信息定量化要求对成像光谱仪 进行精确的光谱通道定标，因此成像光谱仪的光谱性能参数的准确定标是高光谱遥感器数 据提取和定量应用的基本前提。
[0003] 基于光场成像技术的光场成像光谱仪是通过在光场相机的前置成像系统（主镜） 的瞳面出放置由线性渐变滤光片组成的滤光片阵列，实现快照式获取景物目标的三维图谱 数据立方体信息，其主要由主镜、滤波片阵列、微透镜阵列和CCD阵列构成，如图1所示[1]。 对目标成像时，从目标发出的光线经主镜头和微透镜后投影到CCD阵列上形成宏像素，宏 像素中的一个像元对应主镜头孔径的一个采样（子孔径）。滤波片阵列由具有不同波长范 围的线性渐变滤光片组成，滤波片不同位置对应子孔径不同单位位置，当目标发出的光线 经过主镜时，在每个子孔径处分别得到目标不同的波长的光，再经过微透镜成像后被相应 的像元接收。此时，每个微透镜对应一个宏像素，每个宏像素对应目标场景的一个空间采 样，宏像素中每个像元对应该空间位置目标的一个光谱采样，因此一个微透镜对应的宏像 素内的像元包含了空间目标的不同谱段信息。光场成像光谱仪需提取目标景物通过滤光片 阵列的各光谱通道以获取对应的光谱图像。由于系统结构误差、光学衍射和安装误差等引 入了偏移，同一滤光片单元对应的子孔径在不同微透镜后对应的位置不同，造成各个微透 镜后的光谱信息分布存在差异，即任意微透镜后的宏像素对应一个对立的光谱仪具有各自 独立的光谱中心波长和带宽。因此，以各个微透镜后宏像素的固定位置像元代所提取图像 并不是真实的某一特定光谱通道图像。
[0004] [1]发明专利：苏丽娟，袁艳，胡亮，一种快照式高通量的光谱成像方法和光谱成 像仪，授权号ZL201210451761.X。
【发明内容】
[0005] 本发明旨在至少解决上述技术问题之一。
[0006] 为此，本发明的目的在于提出一种成像光谱仪的调试方法。
[0007] 为了实现上述目的，本发明的第一方面的实施例公开了一种成像光谱仪的调试方 法，包括以下步骤：S1 :利用普通光源、单色仪和积分球组成谱段可调节单色均匀面光源， 其中，所述普通光源、所述单色仪和所述积分球依次连接；S2 :利用所述单色均匀面光源照 射成像光谱仪，获取所述成像光谱仪在各个谱段下的系统响应；S3 :根据所述成像光谱仪 在各个谱段下的系统响应，确定每个微透镜下各个像元对应的中心波长和光谱带宽；S4 : 根据每个微透镜下各个像元对应的所述中心波长和所述光谱带宽，确定每个微透镜的光谱 采样波长位置；S5 :根据所述单色均匀面光源的光谱数据、所述中心波长和所述光谱带宽 计算每个微透镜在各个谱段的相对响应非均匀系数；S6 :针对实际场景图像，提取每个微 透镜覆盖像元的灰度值，并根据每个微透镜的所述光谱采样波长位置和所述相对响应非均 匀系数计算得到校正后的光谱数据；以及S7 :根据所述校正后的光谱数据，按照目标空间 位置依次排列重构实际场景光谱数据立方体。
[0008] 根据本发明实施例的一种光场成像光谱仪的光谱标定校正方法，光谱定标装置可 实现对光场成像光谱仪全视场的定标，无需局部扫描；提出的确立各个微透镜各自的波长 采样位置机理，更符合光场成像光谱仪自身机理，有利于保证信息提取的准确性；本发明所 提出数据立方体重构方法在光谱维度上进行数据重构，减少了混叠等因素的影响。
[0009] 另外，根据本发明上述实施例的一种光场成像光谱仪的光谱标定校正方法，还可 以具有如下附加的技术特征：
[0010] 进一步地，所述步骤S1进一步包括，通过调节所述单色仪输出光谱段调节所述单 色面均匀光源对应的谱段λi，利用辐射度计监视并记录不同谱段输出的光源强度I&U。
[0011] 进一步地，所述步骤S2进一步包括，利用计算机控制所述单色仪在响应波长范围 内，以预设的步长Αλ进行扫描，使所述积分球输出不同谱段的单色面光源，由所述成像 光谱仪采集相应谱段的系统响应图像，根据所述系统响应图像获取所述成像光谱仪在各个 谱段下的系统响应。
[0012] 进一步地，所述步骤S3进一步包括，对某一微透镜后覆盖的任一像元，提取所述 像元在所述谱段λ灰度值Iλι，所述灰度值Ιλι满足下列公式：
[0013]
[0014] 其中，I。是最大灰度响应值，λ。是该像元对应的中心波长，Δλ是该像元对应的 光谱带宽，根据步骤S2中采样获取的各个谱段的灰度和上述灰度值Iλι的公式进行拟合， 计算获得所述像元对应的中心波长和光谱带宽。。
[0015] 进一步地，所述步骤S4进一步包括，对所有微透镜根据步骤S3确定各个微透镜各 自覆盖的的像元对应中心波长和光谱带宽，从而确定各个微透镜对应的所述光谱采样波长 位置。
[0016] 进一步地，所述步骤S5进一步包括，获取单色均匀面光源的光谱数据
[ΚλΟ,ΚλΙΚλΙ…IUJ]，并根据所述采样中心波长[λ1;λ2,λ3，…和所述 光谱带宽[ΑλΔλ2,Δλ3，…ΔλJ进行重采样获取的参考光谱数据[IJλD，IJλ2)， IJλ3)，…IJλJ]，根据所述光谱数据和所述参考光谱数据并根据以下公式所述每个微透 镜在各个谱段的相对响应非均匀系数E(Ai) ==ΚλΟ/Ι^λ^，其中，m为微透 镜的编号。
[0017] 进一步地，步骤S6进一步包括，针对实际场景图像，提取每个微透镜覆盖像元的 灰度值，并根据所述光谱采样波长位置获取每个微透镜对应目标A的初始光谱数据[ΙΑ (λ i)，ΙΑ (λ2)，ΙΑ (λ3)，…ΙΑ (λJ]，并根据所述相对响应非均匀系数根据以下公式计算实际光 谱响应ζΙ^λΟ/Ε^λ》，其中，所述目标A的光谱响应曲线为[UaUJ，〗'A( λ2)，Ι' Α(λ3)，…I' Α(λ?)]。
[0018] 本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变 得明显，或通过本发明的实践了解到。
【附图说明】
[0019] 本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变 得明显和容易理解，其中：
[0020] 图1为本发明针对的光谱仪的原理示意图；
[0021] 图2为本发明一个实施例的标定和数据重构的步骤流程图；
[0022] 图3为本发明一个实施例的光谱定标装置示意图；
[0023] 图4为本发明一个实施例的某一波长标定结果局部放大图；
[0024] 图5为本发明一个实施例的某一像元实测响应及拟合的曲线图；
[0025] 图6为本发明一个实施例的三个不同微透镜后各个像元对应的中心波长示意图；
[0026] 图7为本发明一个实施例的响应非均匀性标定装置示意图；
[0027] 图8为本发明一个实施例的初始提取的光谱数据立方体；
[0028] 图9为本发明一个实施例的插值重建的光谱数据立方体。
【具体实施方式】
[0029] 下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终 相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附 图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
[0030] 在本发明的描述中，需要理解的是，术语"中心"、"纵向"、"横向"、"上"、"下"、"前"、 "后"、"左"、"右"、"竖直"、"水平"、"顶"、"底"、"内"、"外"等指示的方位或位置关系为基于 附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所 指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发 明的限制。此外，术语"第一"、"第二"仅用于描述目的，而不能理解为指示或暗示相对重要 性。
[0031] 在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语"安装"、"相 连"、"连接"应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可 以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是 两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本 发明中的具体含义。
[0032]参照下面的描述和附图，将清楚本发明的实施例的这些和其他方面。在这些描述 和附图中，具体公开了本发明的实施例中的一些特定实施方式，来表示实施本发明的实施 例的原理的一些方式，但是应当理解，本发明的实施例的范围不受此限制。相反，本发明的 实施例包括落入所附加权利要求书的精神和内涵范围内的所有变化、修改和等同物。
[0033] 以下结合附图描述根据本发明实施例的一种光场成像光谱仪的光谱标定校正方 法。
[0034] 请参考图2,本发明的实施例的一种成像光谱仪的调试方法，包括以下步骤：
[0035] S1 :利用普通光源、单色仪和积分球组成谱段可调节单色均匀面光源可调节单色 均匀面光源输出，并保证波长可调。其中，所述普通光源、所述单色仪和所述积分球依次连 接，
[0036] 如图3所示，在本发明的一个实施例中，普通光源发出的光经过单色仪分光后成 为单色光进入积分球，经过积分球后形成均匀分布的单色面光源，充满被测系统孔径。通过 调节单色仪输出光谱段调节均匀单色面光源对应的谱段λi，利用辐射度计监视并记录不 同谱段输出的光源强度ku。
[0037]