CN111443045B - 光谱成像分析系统及光谱成像分析方法 - Google Patents
光谱成像分析系统及光谱成像分析方法 Download PDFInfo
- Publication number
- CN111443045B CN111443045B CN202010338391.3A CN202010338391A CN111443045B CN 111443045 B CN111443045 B CN 111443045B CN 202010338391 A CN202010338391 A CN 202010338391A CN 111443045 B CN111443045 B CN 111443045B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- liquid crystal
- incident
- spectral
- array
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 56
- 238000000701 chemical imaging Methods 0.000 title claims description 53
- 238000001228 spectrum Methods 0.000 claims abstract description 157
- 239000004973 liquid crystal related substance Substances 0.000 claims abstract description 72
- 230000003287 optical effect Effects 0.000 claims abstract description 51
- 230000010287 polarization Effects 0.000 claims abstract description 49
- 238000003384 imaging method Methods 0.000 claims abstract description 37
- 239000000126 substance Substances 0.000 claims abstract description 36
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000010183 spectrum analysis Methods 0.000 claims abstract description 12
- 230000006870 function Effects 0.000 claims description 59
- 230000003595 spectral effect Effects 0.000 claims description 55
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- 230000004044 response Effects 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 4
- 230000000979 retarding effect Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 15
- 238000005516 engineering process Methods 0.000 abstract description 4
- 210000002858 crystal cell Anatomy 0.000 description 10
- 238000010276 construction Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 210000004027 cell Anatomy 0.000 description 3
- 238000007405 data analysis Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012634 optical imaging Methods 0.000 description 3
- 238000000611 regression analysis Methods 0.000 description 3
- 102000005717 Myeloma Proteins Human genes 0.000 description 1
- 108010045503 Myeloma Proteins Proteins 0.000 description 1
- CLOMYZFHNHFSIQ-UHFFFAOYSA-N clonixin Chemical compound CC1=C(Cl)C=CC=C1NC1=NC=CC=C1C(O)=O CLOMYZFHNHFSIQ-UHFFFAOYSA-N 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/255—Details, e.g. use of specially adapted sources, lighting or optical systems
-
- 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
- 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/447—Polarisation spectrometry
-
- 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
- G01J2003/2826—Multispectral imaging, e.g. filter imaging
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Spectrometry And Color Measurement (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供了一种光谱成像分析系统及其成像分析方法,光谱成像分析系统包括成像透镜阵列、光计算芯片和图像传感器;成像透镜阵列获取待测目标的原始图像;光计算芯片包括起偏振器、液晶相位延迟器阵列和检偏振器,光计算芯片先通过起偏振器偏振处理,后通过液晶相位延迟器阵干涉处理,以获得入射光谱和产生一组基频与谐波,液晶相位延迟器阵列根据基频和谐波构建该入射光谱的光谱回归函数,通过检偏振器检测入射光谱的偏振干涉态信息;图像传感器生成一组含有空间信息和偏振干涉态信息的光谱图像序列,并根据光谱回归函数计算出入射光谱的物质物理化学量信息。与相关技术相比,本发明实现了对动态目标进行实时全光谱成像和光谱分析。
Description
【技术领域】
本发明涉及光谱成像技术领域,尤其涉及一种光谱成像分析系统及光谱成像分析方法。
【背景技术】
光谱成像技术通常用于物质成分及组分分析,光谱成像技术可获得物质轮廓纹理信息,光谱分析和成像相结合的光谱成像分析系统可测得图像每个像素的对应物质的化学成分,可在复杂背景下检测和识别高价值目标。
在相关技术中,基于液晶可调谐滤波器、声光调谐滤波器分光的凝视式光谱成像分析系统采用可编程调谐的滤光片进行光谱分光,一次采集一个窄带光谱图像,分时多次采集完成图像每像素的光谱拼接,且对各波段图像需要像素级的配准以获得光谱图像数据;基于光栅分光的推扫式光谱成像分析系统通过光栅的缝隙将目标切割成多条线目标,通过推扫的方式获取多条线目标的光谱数据,且同一线目标上各点的光谱信息是同时采集的,最后将多条线目标的光谱信息进行拼接以获得观察目标光谱图像数据;快照式光谱成像分析系统通常把多个像素分为一组,组合成一个“光谱像素”,在“光谱像素”的每个像素镀上特定波长的窄带滤光材料,获取多个波段的光谱信息。
但凝视式、推扫式或快照式的光谱成像分析系统均需要利用光谱图像数据处理软件对光谱数据进行分析,提取各像素光谱曲线,建立目标光谱与物质物理化学量信息的回归模型,才能进对图像中各像素对应的物质进行光谱分类或识别。
然而,相关技术中,凝视式和推扫式的光谱成像分析系统需要进行光谱扫描,光谱扫描耗时,且光谱数据分析所需的大容量数据传输和大规模计算也需要较长时间,另外还要求在光谱扫描过程中目标保持固定的空间位置,难以处理空间位置发生变化的运动目标,适用性差,难以进行动态目标实时光谱成像分析;而快照式光谱成像分析系统对光谱数据分析的过程需要耗费的时间较长,亦难以进行实时光谱成像分析。
因此,实有必要提供一种新的光谱成像分析系统及光谱成像分析方法解决上述技术问题。
【发明内容】
本发明的目的在于提供一种光谱成像分析系统及光谱成像分析方法,以解决难以对动态目标进行实时分析的问题。
为达到上述目的,本发明提供一种光谱成像分析系统,所述光谱成像分析系统沿光入射方向依次包括:
成像透镜阵列,所述成像透镜阵列获取待测目标的原始图像;
光计算芯片,所述光计算芯片的焦面在所述成像透镜阵列的焦面上,所述光计算芯片包括沿所述光入射方向依次设置的起偏振器、液晶相位延迟器阵列和检偏振器,所述起偏振器位于所述液晶相位延迟器阵列靠近所述成像透镜阵列的一侧,所述检偏振器位于所述液晶相位延迟器阵列远离所述成像透镜阵列的一侧,所述光计算芯片接收所述原始图像,并通过所述起偏振器对所述原始图像进行偏振处理,再通过所述液晶相位延迟器阵列对偏振后的所述原始图像进行干涉处理以获得入射光谱且同时产生一组与所述入射光谱对应的基频与谐波,所述液晶相位延迟器阵列根据所述基频和所述谐波构建该入射光谱的光谱回归函数,通过所述检偏振器检测所述入射光谱的偏振干涉态信息;以及,
图像传感器,所述图像传感器接收所述入射光谱,利用压缩光学传感算法生成一组含有空间信息和偏振干涉态信息的光谱图像序列,并根据所述光谱回归函数对所述入射光谱进行光谱分析处理以计算出与该光谱图像序列对应的物质物理化学量信息。
优选的,所述起偏振器的偏振方向与所述液晶相位延迟器阵列的光轴方向之间的夹角为45°。
优选的,所述检偏振器包括相互连接的第一区域和第二区域,所述第一区域的偏振方向与所述液晶相位延迟器的光轴方向之间的夹角为45°,所述第二区域的偏振方向与的所述液晶相位延迟器的光轴方向之间的夹角为135°。
优选的,所述成像透镜阵列包括多个呈阵列分布的透镜单元,所述液晶相位延迟器阵列包括多个呈阵列分布的液晶单元,所述液晶单元与所述透镜单元一一对应设置;一个所述透镜单元获取一个所述原始图像,一个所述液晶单元对与其对应的所述透镜单元所获取的所述原始图像进行干涉处理以获得一个所述入射光谱和一组与该入射光谱对应的所述基频与所述谐波。
本发明提供一种光谱成像分析方法,该方法包括以下步骤:
步骤S1,获取待测目标的原始图像;
步骤S2,对所述原始图像进行偏振处理,对偏振后的所述原始图像进行可调相位延迟量的干涉处理以获得入射光谱和一组与该入射光谱对应的基频与谐波,根据所述基频和所述谐波构建该入射光谱的光谱回归函数;其中,可通过调节驱动电压的大小调节所述可调相位延迟量的大小,进而调整所述入射光谱的偏振干涉态并重构所述光谱回归函数;
步骤S3,检测所述入射光谱的偏振干涉态信息;
步骤S4,利用压缩光学传感算法对所述入射光谱进行计算处理以生成一组含有空间信息和偏振干涉态信息的光谱图像序列,并根据所述光谱回归函数对所述入射光谱进行光谱分析处理以计算出与该光谱图像序列对应的物质物理化学量信息;
步骤S5,判断环境光与光谱图像序列之间的响应状态是否发生变化,若是,则根据判断结果实时调节所述驱动电压的大小以调节所述可调相位延迟量的大小,并返回所述步骤S2。
优选的,所述原始图像包括多个,多个所述原始图像含有相同的空间信息,一个所述原始图像经过偏振干涉处理后生成一个所述入射光谱。
本发明提供一种光谱成像分析系统,其包括处理器以及存储器,所述存储器中存储有用于所述处理器执行的控制程序,其中,所述控制程序被所述处理器执行时实现本发明所述的光谱成像分析方法的步骤。
本发明提供一种计算机可读存储介质,其存储有计算机程序,所述计算机程序被处理器执行时实现本发明所述的光谱成像分析系统的步骤。
与相关技术相比,本发明的光谱成像分析系统和光谱成像分析方法中,通过光计算芯片对待测目标的原始图像进行偏振干涉处理以获得入射光谱和一组与该入射光谱对应的基频与谐波,同时根据干涉处理后获得的一组与入射光谱对应的基频与谐波构建该入射光谱的光谱回归函数,图像传感器利用压缩光学传感算法对入射光谱进行计算处理以生成一组含有空间信息和偏振干涉态信息的光谱图像序列,并根据光谱回归函数对入射光谱进行光谱分析处理以计算出与该光谱图像序列对应的物质物理化学量信息;上述结构和方法中,在图像传感器中同时通过压缩光学传感算法计算出目标的光谱图像序列和通过可重构的光谱回归函数对入射光谱进行光谱数据分析以获得目标的物质物理化学量信息,此处,即把光谱成像和光谱回归分析相融合,使目标在通过光谱成像分析系统的成像过程中,使得光谱图像信息的采集和分析能够同时完成,实现了对动态目标实时的高灵敏度的全光谱成像分析,另外,光谱回归函数可以根据实际应用的需求进行实时的调整,并利用可编程光计算器件,根据环境光及系统光谱响应实时构建光计算所需的透过谱线,具有优良的光环境适应能力,适用性强。
【附图说明】
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图,其中:
图1为本发明光谱成像分析系统的立体结构示意图;
图2为本发明光计算芯片的结构示意图;
图3为本发明光谱成像分析系统的流程示意图;
图4为本发明的构建的正、余弦函数示意图;
图5为本发明的光谱回归函数的光学构建示意图。
【具体实施方式】
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
请同时参阅图1-5所示,本发明提供一种光谱成像分析系统100,其沿光入射方向依次包括成像透镜阵列1、光计算芯片2以及图像传感器3。
所述成像透镜阵列1用于获取待测目标4的原始图像。
在本实施方式中,所述成像透镜阵列1包括多个呈阵列分布的透镜单元11,所述透镜单元11的具体数量可以根据实际应用的需要进行具体的设置,其中,一个透镜单元11能够获取一个原始图像。
所述光计算芯片2的焦面在所述成像透镜阵列1的焦面上,用于获取所述原始图像,并对所述原始图像进行计算处理以获得入射光谱以及构建与该入射光谱对应的光谱回归函数。
在本实施方式中,所述光计算芯片2包括沿所述光入射方向依次设置的起偏振器21、液晶相位延迟器阵列22和检偏振器23;更优的,在其他的实施方式中,所述光计算芯片2还包括设置于所述液晶相位延迟器阵列22与所述检偏振器23之间的固定相位延迟片。
所述起偏振器21和所述检偏振器23均为偏振片,所述起偏振器21位于所述液晶相位延迟器阵列22靠近所述成像透镜阵列1的一侧,其用于对所述原始图像进行偏振处理;所述液晶相位延迟器阵列22用于对偏振后的所述原始图像进行干涉处理以获得入射光谱且同时产生一组与所述入射光谱对应的基频与谐波,且所述液晶相位延迟器阵列22根据所述基频和所述谐波构建该入射光谱的光谱回归函数;所述检偏振器23位于所述液晶相位延迟器阵列22远离所述成像透镜阵列1的一侧,其用于检测所述入射光谱的偏振干涉态信息;在设置固定相位延迟片的实施例当中,该固定相位延迟片用于对所述入射光谱进行相位延迟量为的干涉处理。
进一步的,所述起偏振器21的偏振方向D1与所述液晶相位延迟器阵列22的光轴方向D2之间的夹角为45°,所述检偏振器23包括相互连接的第一区域231和第二区域232,所述第一区域231的偏振方向D3与所述液晶相位延迟器22的光轴方向D2之间的夹角为45°,所述第二区域232的偏振方向D4与的所述液晶相位延迟器22的光轴方向D2之间的夹角为135°。
值得一提的是,所述液晶相位延迟器阵列22包括多个呈阵列分布的液晶单元221,所述液晶单元221与所述透镜单元11一一对应设置,即所述液晶单元221的数量及其空间位置均与所述透镜单元11的数量及空间位置对应设置;其中,一个所述液晶单元221相当于构建了一个“光谱像素”,而各所述液晶单元221用于对与其对应的所述透镜单元11所获取的所述原始图像进行干涉处理以获得一个入射光谱和一组与该入射光谱对应的所述基频与所述谐波。
所述图像传感器3接收所述入射光谱,利用压缩光学传感算法生成一组含有空间信息和偏振干涉态信息的光谱图像序列,并根据所述光谱回归函数对所述入射光谱进行光谱分析处理以计算出与该光谱图像序列对应的物质物理化学量信息,该物质物理化学量信息至少包括轮廓纹理信息、物质成分信息以及各组分信息中的任意一种。
本发明还提供一种应用于上述光谱成像系统100的光谱成像分析方法,其包括以下步骤:
步骤S1,获取待测目标的原始图像。
在所述步骤S1中,通过所述透镜成像阵列1获取所述待测目标4的所述原始图像,更具体的,通过所述透镜成像阵列1的多个透镜单元11获取多个原始图像,多个所述原始图像含有相同的空间信息。
步骤S2,对所述原始图像进行偏振处理,对偏振后的所述原始图像进行可调相位延迟量的干涉处理以获得入射光谱和一组与该入射光谱对应的基频与谐波,根据所述基频和所述谐波构建该入射光谱的光谱回归函数;其中,可通过调节驱动电压的大小调节所述可调相位延迟量的大小,进而调整所述入射光谱的偏振干涉态并重构所述光谱回归函数。
在所述步骤S2中,所述透镜成像阵列1获取的所述原始图像成像于所述光计算芯片2的焦面上,所述光计算芯片2首先通过所述起偏振器21对所述原始图像进行偏振处理,然后通过所述液晶相位延迟器阵列22对所述镜成像阵列1所获取的多个所述原始图像进行可调相位延迟量的干涉处理以获得多个所述入射光谱和一组与该入射光谱对应的所述基频与所述谐波,更具体的,一个所述液晶单元221处理一个与其对应的透镜单元11所获取的所述原始图像进行干涉处理以获得一个入射光谱和一组与该入射光谱对应的所述基频与所述谐波,同时,通过所述液晶相位延迟器阵列22根据所述基频和所述谐波构建该入射光谱的光谱回归函数。在实际应用中,为了能够使该光谱成像系统100能够具有优良的光环境适应能力,提高其适用性,可以根据实际的待测目标4实际所处的光环境状态对所述液晶相位延迟器阵列22的驱动电压的大小进行调节,从而调节所述可调相位延迟量的大小,进而调整所述入射光谱的偏振干涉态并重构所述光谱回归函数。
更优的,在其他增设了固定相位延迟片的实施方式中,在所述步骤S2中还包括以下步骤:通过所述固定相位延迟片对所述入射光谱进行相位延迟量为的干涉处理,该步骤发生在所述检偏振器23检测所述入射光谱的偏振干涉态信息的步骤之前;由于所述固定相位延迟片的引入,使得所述检偏振器23检测到入射光谱的相位延迟量增大了/>
步骤S3,检测所述入射光谱的偏振干涉态信息。
在所述步骤S3中,所述光计算芯片2通过所述检偏振器23检测所述入射光谱的偏振干涉态信息。
步骤S4,利用压缩光学传感算法对所述入射光谱进行计算处理以生成一组含有空间信息和偏振干涉态信息的光谱图像序列,并根据所述光谱回归函数对所述入射光谱进行光谱分析处理以计算出与该光谱图像序列对应的物质物理化学量信息。
在所述步骤S4中,所述图像传感器3利用压缩光学传感算法对多个所述入射光谱进行计算处理以生成一组含有空间信息和偏振干涉态信息的光谱图像序列,同时,所述图像传感器3根据所述光谱回归函数分别对所述入射光谱进行光谱分析处理以计算出与该光谱图像序列对应的物质物理化学量信息。
步骤S5,判断环境光与光谱图像序列之间的响应状态是否发生变化,若是,则根据判断结果实时调节所述驱动电压的大小以调节所述可调相位延迟量的大小,并返回所述步骤S2。
在上述过程中,还需要说明的是光谱回归函数的构建和重建过程,具体如下:
(1)、余弦函数、正弦函数光学构建与重构
则两个单元对应的空间位置透射光强的差为:
在液晶相位延迟器阵列22上增加固定相位延迟片,即引入固定相位延迟对于相位延迟量为/>的液晶单元221,对应起偏振器21与检偏振器23偏振方向垂直、起偏振器21与检偏振器23偏振方向平行的偏振干涉透射光强分别为:
则两个区域的差为:
(2)、谐波函数光学构建与重构
液晶相位延迟器阵列22各单元相位延迟量由液晶厚度d和液晶双折射率差Δr之积确定,使液晶相位延迟器阵列22各液晶单元221的相位延迟量为该液晶相位延迟量随电压增大而减小,调节电压,使各单元透过光强分别为:
cos(δv)、sin(δv)、cos(2δv)、sin(2δv)、cos(3δv)、sin(3δv)、......、cos(kδν)、sin(mδv),其中n为自然数且n<m。即可以调节液晶相位延迟器阵列22的驱动电压可以得到正弦和余弦正交的基频和谐波函数序列,该函数序列可作为回归函数的傅里叶级数的正交频率分量,线性组合成所需的光谱回归函数。
由于液晶相位延迟器阵列22的相位延迟量可以由驱动电压来调节,基频和谐波函数序列cos(δv)、sin(δv)、cos(2δν)、sin(2δv)、cos(3δv)、sin(3δν)、......、cos(kδv)、sin(mδv)在液晶单元221的工作范围内任意调节,即可以通过改变驱动电压来重构所需的光谱回归函数。
(3)、光谱回归函数光学构建
若所述待测目标4含有M个组分,各组分的物理化学参量为c1、c2、……、cM,M为自然数,测得的对应的光谱为f1(v)→[a11,a12,...,a1N]
f2(v)→[a21,a22,...,a2N]
...
fM(v)→[aM1,aM2,...,aMN]
其中,N为自然数,与光谱采样分辨率相关。建立光谱与物质物理化学参量关系:
a11x1+a12x2+...+a1NxN=c1
a21x1+a22x2+...+a2NxN=c2
......
aM1x1+aM2x2+...+aMNxN=cN
当基频cos(δv)、sin(δv)在光谱分析波长范围内有一个周期时,若n值较大时,例如n=20,由基频和谐波组成的近似的光谱回归函数为:
可以充分逼近X(v),即用液晶相位延迟器阵列22的有限的n个液晶单元221偏振干涉产生的基频与谐波光强函数,详细如图4所示,可以以高精度构建光谱回归函数,分别用4×4、6×4、8×4、10×4、12×4数量的液晶单元221偏振干涉产生的基频与谐波光强函数,当液晶单元221的数量达到48时,基频和有限的谐波组合的函数可以以高精度表达光谱回归函数。
(4)、光计算
在成像光路中加入光计算芯片2,则待测目标4中各组分fi(v)的成像可以表示为:
则光学成像过程即是回归分析的积分过程,积分值为图像传感器3获得的图像灰度值。把经过光计算芯片2的光学成像过程称为光计算光谱成像。可根据回归模型、光源、光学成像镜头及传感器光谱响应来设计。
(5)、光计算芯片重构
环境光或成像系统发生变化,或待测目标4的光谱性质改变,驱动光计算芯片2构建新的光谱回归函数:
光计算芯片2的滤波函数也发生变换。根据感兴趣的光谱范围,选择基频,驱动液晶,并计算出系数后,通过调整图像传感器3的各区域图像灰度值,组合成新的滤波函数。
本发明提供一种光谱成像分析系统,其包括处理器以及存储器,所述存储器中存储有用于所述处理器执行的控制程序,其中,所述控制程序被所述处理器执行时实现本发明所述的光谱成像分析方法的步骤。
本发明提供一种计算机可读存储介质,其存储有计算机程序,所述计算机程序被处理器执行时实现本发明所述的光谱成像分析系统的步骤。
与相关技术相比,本发明的光谱成像分析系统和光谱成像分析方法中,通过光计算芯片对待测目标的原始图像进行偏振干涉处理以获得入射光谱和一组与该入射光谱对应的基频与谐波,同时根据干涉处理后获得的一组与入射光谱对应的基频与谐波构建该入射光谱的光谱回归函数,图像传感器利用压缩光学传感算法对入射光谱进行计算处理以生成一组含有空间信息和偏振干涉态信息的光谱图像序列,并根据光谱回归函数对入射光谱进行光谱分析处理以计算出与该光谱图像序列对应的物质物理化学量信息;上述结构和方法中,在图像传感器中同时通过压缩光学传感算法计算出目标的光谱图像序列和通过可重构的光谱回归函数对入射光谱进行光谱数据分析以获得目标的物质物理化学量信息,此处,即把光谱成像和光谱回归分析相融合,使目标在通过光谱成像分析系统的成像过程中,使得光谱图像信息的采集和分析能够同时完成,实现了对动态目标实时的高灵敏度的全光谱成像分析,另外,光谱回归函数可以根据实际应用的需求进行实时的调整,并利用可编程光计算器件,根据环境光及系统光谱响应实时构建光计算所需的透过谱线,具有优良的光环境适应能力,适用性强。
以上所述的仅是本发明的实施方式,在此应当指出,对于本领域的普通技术人员来说,在不脱离本发明创造构思的前提下,还可以做出改进,但这些均属于本发明的保护范围。
Claims (10)
1.一种光谱成像分析系统,其特征在于,其沿光入射方向依次包括:
成像透镜阵列,所述成像透镜阵列获取待测目标的原始图像;
光计算芯片,所述光计算芯片的焦面在所述成像透镜阵列的焦面上,所述光计算芯片包括沿所述光入射方向依次设置的起偏振器、液晶相位延迟器阵列和检偏振器,所述起偏振器位于所述液晶相位延迟器阵列靠近所述成像透镜阵列的一侧,所述检偏振器位于所述液晶相位延迟器阵列远离所述成像透镜阵列的一侧,所述光计算芯片接收所述原始图像,并通过所述起偏振器对所述原始图像进行偏振处理,再通过所述液晶相位延迟器阵列对偏振后的所述原始图像进行干涉处理以获得入射光谱且同时产生一组与所述入射光谱对应的基频与谐波,所述液晶相位延迟器阵列根据所述基频和所述谐波构建该入射光谱的光谱回归函数,通过所述检偏振器检测所述入射光谱的偏振干涉态信息;以及,
图像传感器,所述图像传感器接收所述入射光谱,利用压缩光学传感算法生成一组含有空间信息和偏振干涉态信息的光谱图像序列,并根据所述光谱回归函数对所述入射光谱进行光谱分析处理以计算出与该光谱图像序列对应的物质物理化学量信息;
其中,所述光谱回归函数通过以下步骤构建:
定义待测目标含有M个组分,每一所述组分的物理化学参量为c1、c2、……、cM,M为自然数,所述待测目标对应的所述入射光谱f(v)满足:
其中,N为自然数;
建立所述光谱与所述物理化学参量的关系为:
2.根据权利要求1所述的光谱成像分析系统,其特征在于,所述起偏振器的偏振方向与所述液晶相位延迟器阵列的光轴方向之间的夹角为45°。
3.根据权利要求1所述的光谱成像分析系统,其特征在于,所述检偏振器包括相互连接的第一区域和第二区域,所述第一区域的偏振方向与所述液晶相位延迟器的光轴方向之间的夹角为45°,所述第二区域的偏振方向与的所述液晶相位延迟器的光轴方向之间的夹角为135°。
5.根据权利要求1所述的光谱成像分析系统,其特征在于,所述成像透镜阵列包括多个呈阵列分布的透镜单元,所述液晶相位延迟器阵列包括多个呈阵列分布的液晶单元,所述液晶单元与所述透镜单元一一对应设置;一个所述透镜单元获取一个所述原始图像,一个所述液晶单元对与其对应的所述透镜单元所获取的所述原始图像进行干涉处理以获得一个所述入射光谱和一组与该入射光谱对应的所述基频与所述谐波。
6.一种光谱成像分析方法,其特征在于,所述光谱成像分析方法包括以下步骤:
步骤S1,获取待测目标的原始图像;
步骤S2,对所述原始图像进行偏振处理,对偏振后的所述原始图像进行可调相位延迟量的干涉处理以获得入射光谱和一组与该入射光谱对应的基频与谐波,根据所述基频和所述谐波构建该入射光谱的光谱回归函数;其中,通过调节驱动电压的大小调节所述可调相位延迟量的大小,进而调整所述入射光谱的偏振干涉态并重构所述光谱回归函数,其中,所述光谱回归函数通过以下步骤构建:
定义待测目标含有M个组分,每一所述组分的物理化学参量为c1、c2、……、cM,M为自然数,所述待测目标对应的所述入射光谱f(v)满足:
其中,N为自然数;
建立所述光谱与所述物理化学参量的关系为:
步骤S3,检测所述入射光谱的偏振干涉态信息;
步骤S4,利用压缩光学传感算法对所述入射光谱进行计算处理以生成一组含有空间信息和偏振干涉态信息的光谱图像序列,并根据所述光谱回归函数对所述入射光谱进行光谱分析处理以计算出与该光谱图像序列对应的物质物理化学量信息;
步骤S5,判断环境光与光谱图像序列之间的响应状态是否发生变化,若是,则根据判断结果实时调节所述驱动电压的大小以调节所述可调相位延迟量的大小,并返回所述步骤S2。
7.根据权利要求6所述的光谱成像分析方法,其特征在于,所述原始图像包括多个,多个所述原始图像含有相同的空间信息,一个所述原始图像经过偏振干涉处理后生成一个所述入射光谱。
9.一种光谱成像分析系统,其特征在于,所述光谱成像分析系统包括处理器以及存储器,所述存储器中存储有用于所述处理器执行的控制程序,其中,所述控制程序被所述处理器执行时实现如上权利要求6至8中任一项所述的光谱成像分析方法的步骤。
10.一种计算机可读存储介质,其存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如上权利要求6至8中任一项所述的光谱成像分析系统的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010338391.3A CN111443045B (zh) | 2020-04-26 | 2020-04-26 | 光谱成像分析系统及光谱成像分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010338391.3A CN111443045B (zh) | 2020-04-26 | 2020-04-26 | 光谱成像分析系统及光谱成像分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111443045A CN111443045A (zh) | 2020-07-24 |
CN111443045B true CN111443045B (zh) | 2023-07-04 |
Family
ID=71656202
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010338391.3A Active CN111443045B (zh) | 2020-04-26 | 2020-04-26 | 光谱成像分析系统及光谱成像分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111443045B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112163627B (zh) * | 2020-10-09 | 2024-01-23 | 北京环境特性研究所 | 目标物体的融合图像生成方法、装置及系统 |
CN113359069B (zh) * | 2021-06-04 | 2022-06-07 | 中国科学院国家天文台南京天文光学技术研究所 | 一种高效率全Stokes分量偏振测量方法 |
CN116448689B (zh) * | 2023-06-16 | 2023-08-29 | 天津博霆光电技术有限公司 | 基于像素微分的单光频率物体全构造检测方法及设备 |
CN117011505B (zh) * | 2023-10-07 | 2024-03-15 | 深圳市中达瑞和科技有限公司 | 基于光谱数据的识别方法、系统及相关设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106872037A (zh) * | 2017-02-22 | 2017-06-20 | 西安交通大学 | 快照式紧凑小型化光场成像全偏振光谱探测装置及方法 |
CN108007574A (zh) * | 2017-11-17 | 2018-05-08 | 西安电子科技大学 | 分辨率可调型快照式图像光谱线偏振探测装置及方法 |
CN108627247A (zh) * | 2017-03-23 | 2018-10-09 | 罗伯特·博世有限公司 | 傅里叶变换光谱仪和用于运行傅里叶变换光谱仪的方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008042766A1 (en) * | 2006-09-29 | 2008-04-10 | Chemimage Corporation | Spectral imaging system |
US8125641B2 (en) * | 2009-03-27 | 2012-02-28 | N&K Technology, Inc. | Method and apparatus for phase-compensated sensitivity-enhanced spectroscopy (PCSES) |
US10768497B2 (en) * | 2016-10-03 | 2020-09-08 | Xerox Corporation | Hyperspectral imaging system |
US10642056B2 (en) * | 2016-10-19 | 2020-05-05 | CSEM Centre Suisse d'Electronique et de Microtechnique SA—Recherche et Développement | Multispectral or hyperspectral imaging and imaging system based on birefringent subwavelength resonating structure |
-
2020
- 2020-04-26 CN CN202010338391.3A patent/CN111443045B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106872037A (zh) * | 2017-02-22 | 2017-06-20 | 西安交通大学 | 快照式紧凑小型化光场成像全偏振光谱探测装置及方法 |
CN108627247A (zh) * | 2017-03-23 | 2018-10-09 | 罗伯特·博世有限公司 | 傅里叶变换光谱仪和用于运行傅里叶变换光谱仪的方法 |
CN108007574A (zh) * | 2017-11-17 | 2018-05-08 | 西安电子科技大学 | 分辨率可调型快照式图像光谱线偏振探测装置及方法 |
Non-Patent Citations (2)
Title |
---|
A snapshot light field imaging spectrometer;Lijuan Sua et al.;《optik》;20151231;第126卷;877-881页 * |
光谱法鉴别朱墨时序的实验研究;陈淑萍等;《法制博览》;20171231(第6期);160-161页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111443045A (zh) | 2020-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111443045B (zh) | 光谱成像分析系统及光谱成像分析方法 | |
Zhao et al. | Multi-band polarization imaging | |
Gupta et al. | Acousto-optic tunable filter based visible-to near-infrared spectropolarimetric imager | |
CN108007574B (zh) | 分辨率可调型快照式图像光谱线偏振探测装置及方法 | |
Hegyi et al. | Hyperspectral imaging with a liquid crystal polarization interferometer | |
US9097585B2 (en) | Spectrographic polarimeter and method of recording state of polarity | |
CN102879097B (zh) | 圆偏振高光谱成像探测系统 | |
Sattar et al. | Review of spectral and polarization imaging systems | |
Martínez-Ponce et al. | Hybrid complete Mueller polarimeter based on phase modulators | |
EP3877734B1 (en) | Method and system for polarimetry using static geometric polarization manipulation | |
CN103528688B (zh) | 一种全偏振高光谱干涉成像装置及方法 | |
Han et al. | Iterative reconstruction for snapshot intensity-modulated linear imaging spectropolarimetry without Fourier transform and phase calibration | |
CN102721470A (zh) | 一种静态光谱偏振成像仪 | |
Vaughn et al. | A portable imaging mueller matrix polarimeter based on a spatio-temporal modulation approach: theory and implementation | |
Ning et al. | Single-shot measurement of polarization state at low light field using Mueller-mapping star test polarimetry | |
US7023546B1 (en) | Real-time imaging spectropolarimeter based on an optical modulator | |
Vila-Francés et al. | Improving the performance of acousto-optic tunable filters in imaging applications | |
O’holleran et al. | Methodology for imaging the 3D structure of singularities in scalar and vector optical fields | |
Arteaga et al. | Mueller Matrix Imaging | |
Zhou et al. | Fast compressed channeled spectropolarimeter for full Stokes vector measurement | |
EP4407287A1 (en) | Polarimetric phase imaging | |
D’Acquisto et al. | Full field automated evaluation of the quarter wave plate retardation by phase stepping technique | |
Vila-Francés et al. | Configurable bandwidth imaging spectrometer based on acousto-optic tunable filter | |
Pozhar et al. | Acousto-optical imaging spectroscopy | |
US20190204153A1 (en) | Method and apparatus for transforming uniformly or non-uniformly sampled interferograms to produce spectral data |
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 |