CN117091700B - 一种使用非阵列探测器实现高光谱成像的系统 - Google Patents
一种使用非阵列探测器实现高光谱成像的系统 Download PDFInfo
- Publication number
- CN117091700B CN117091700B CN202311107878.0A CN202311107878A CN117091700B CN 117091700 B CN117091700 B CN 117091700B CN 202311107878 A CN202311107878 A CN 202311107878A CN 117091700 B CN117091700 B CN 117091700B
- Authority
- CN
- China
- Prior art keywords
- imaging
- spectrum
- matrix
- target object
- modulation
- 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
- 238000000701 chemical imaging Methods 0.000 title claims abstract description 36
- 238000003384 imaging method Methods 0.000 claims abstract description 77
- 239000000463 material Substances 0.000 claims abstract description 75
- 239000002159 nanocrystal Substances 0.000 claims abstract description 73
- 238000001228 spectrum Methods 0.000 claims abstract description 73
- 230000001105 regulatory effect Effects 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 80
- 230000003595 spectral effect Effects 0.000 claims description 25
- 239000011521 glass Substances 0.000 claims description 22
- 238000005259 measurement Methods 0.000 claims description 18
- 239000002707 nanocrystalline material Substances 0.000 claims description 17
- 238000000034 method Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000003287 optical effect Effects 0.000 claims description 10
- KFZMGEQAYNKOFK-UHFFFAOYSA-N Isopropanol Chemical compound CC(C)O KFZMGEQAYNKOFK-UHFFFAOYSA-N 0.000 claims description 7
- 238000004140 cleaning Methods 0.000 claims description 7
- 239000003960 organic solvent Substances 0.000 claims description 6
- 238000000411 transmission spectrum Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 4
- 238000004506 ultrasonic cleaning Methods 0.000 claims description 4
- 238000002360 preparation method Methods 0.000 claims description 3
- 239000002904 solvent Substances 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 230000006798 recombination Effects 0.000 claims description 2
- 238000005215 recombination Methods 0.000 claims description 2
- 239000010408 film Substances 0.000 claims 1
- 239000010409 thin film Substances 0.000 claims 1
- 238000005516 engineering process Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000033001 locomotion Effects 0.000 description 4
- VLKZOEOYAKHREP-UHFFFAOYSA-N n-Hexane Chemical group CCCCCC VLKZOEOYAKHREP-UHFFFAOYSA-N 0.000 description 4
- 239000000428 dust Substances 0.000 description 3
- 229910052981 lead sulfide Inorganic materials 0.000 description 3
- 229940056932 lead sulfide Drugs 0.000 description 3
- WUPHOULIZUERAE-UHFFFAOYSA-N 3-(oxolan-2-yl)propanoic acid Chemical compound OC(=O)CCC1CCCO1 WUPHOULIZUERAE-UHFFFAOYSA-N 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 229910052980 cadmium sulfide Inorganic materials 0.000 description 2
- UHYPYGJEEGLRJD-UHFFFAOYSA-N cadmium(2+);selenium(2-) Chemical compound [Se-2].[Cd+2] UHYPYGJEEGLRJD-UHFFFAOYSA-N 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 238000005137 deposition process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- GGYFMLJDMAMTAB-UHFFFAOYSA-N selanylidenelead Chemical compound [Pb]=[Se] GGYFMLJDMAMTAB-UHFFFAOYSA-N 0.000 description 2
- 238000000862 absorption spectrum Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
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
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B82—NANOTECHNOLOGY
- B82Y—SPECIFIC USES OR APPLICATIONS OF NANOSTRUCTURES; MEASUREMENT OR ANALYSIS OF NANOSTRUCTURES; MANUFACTURE OR TREATMENT OF NANOSTRUCTURES
- B82Y30/00—Nanotechnology for materials or surface science, e.g. nanocomposites
-
- 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/02—Details
- G01J3/0205—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows
- G01J3/0208—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows using focussing or collimating elements, e.g. lenses or mirrors; performing aberration correction
-
- 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/2803—Investigating the spectrum using photoelectric array detector
- G01J2003/2813—2D-array
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Nanotechnology (AREA)
- Composite Materials (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Materials Engineering (AREA)
- Crystallography & Structural Chemistry (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及高光谱成像技术领域,提供一种使用非阵列探测器实现高光谱成像的系统,包括成像模块、振幅调制模块、光谱调制模块、接收模块和计算机;成像模块对目标物体的反射光进行收集,成像到数字微镜器件的靶面上;振幅调制模块实现数字微镜器件对其靶面上的成像进行振幅调制,由反射镜调整出射方向;光谱调制模块使用纳米晶体材料滤波片对出射信号光进行光谱调制;接收模块对信号光聚焦并将模拟电压信号转化为数字信号,实现加载到数字微镜器件上的调制图案与光电探测器输出的电压值信号对应的关系;计算机重构出目标物体的三维数据立方体。上述系统解除了高光谱成像系统对阵列探测器的依赖,并将光谱成像的范围拓展到特殊波段。
Description
技术领域
本发明涉及一种用于实现高光谱成像的系统,尤其是涉及应用单像素计算成像技术与纳米晶体材料,使用非阵列探测器实现高光谱成像的系统,属于高光谱成像技术领域。
背景技术
高光谱成像技术是认知世界的一种重要手段,能够同时捕获空间信息和光谱信息,常常用来检测目标物体的物理结构、化学成分等信息。高光谱成像技术已广泛应用在遥感领域、食品安全领域、生物医药领域以及军事等领域。
目前已有的高光谱成像技术主要有以下几种:点扫描、线阵扫描、光谱阵列扫描以及快照式成像。点扫描使用光栅或棱镜等分光器件,与线阵探测器配合,测量目标物体中单个像素点处的所有光谱信息。之后通过两个自由度的机械运动控制,在x,y空间中的每一个像素点进行点扫描。这种成像方式由于引入了复杂的机械运动装置,所以成像速度慢,空间分辨率低。同时,这种成像方式需要线阵探测器对分光信号进行探测,对这种器件存在依赖,限制了使用场景。
线阵扫描能够在一次测量中测量目标物体中一条线上所有像素点的光谱信息,之后通过使用机械运动装置,在另一个维度上进行扫描。线阵扫描的工作原理与点扫描类似,不同点在于点扫描获取的是单个像素点的光谱信息,线阵扫描是获取一条线上所有像素点的光谱信息。由于这种成像方式同样引入了机械运动装置,存在和点扫描形同的成像速度慢,依赖探测器阵列的不足。
光谱阵列扫描无需在空间维度进行扫描,能够使用窄带滤波器和面阵探测器,每次测量中都获取一个窄带波段内的2D空间信息。这种方式与点扫描、线阵扫描技术相似、能量利用率也较低。快照式光谱成像是指在探测器单个积分周期内测量3D数据立方体。而成像探测器是2D探测器,单个积分周期只能获取2D信息,因此快照式光谱成像需要将数据立方体划分为多个2D数据以同时测量,再通过后处理将2D数据重新组合为3D立方体。最简易的快照式光谱成像仪为通道式光谱成像仪。
快照式光谱成像解决了2D成像系统不能一次采集3D数据的问题。每个通道由滤光片、光学系统和探测器构成,每个通道配置不同中心波长的窄带滤光片,一个通道获取一个窄带光谱图像,形成3D数据立方体。这种成像方式具有原理简单、易于实现、相机可靠的优点。但由于多通道并行采集,各通道存在视差,各通道采集的窄带光谱图像需要配准;同时,一个成像通道对应一个窄带波段,光谱通道数取决于光学系统与探测器的数目,因此光谱通道数量较少,一般为3、4;采用多个探测器,体积、质量均比较大。
发明内容
本发明针对现有高光谱成像技术存在的问题,提供一种使用非阵列探测器实现高光谱成像的系统,该系统应用单像素成像技术和纳米晶体材料,能够在不使用阵列探测器的情况下,实现对目标物体空间信息和光谱信息的采集,获取三维数据立方体。
本发明的使用非阵列探测器实现高光谱成像的系统,采用以下技术方案。
该系统,包括成像模块、振幅调制模块、光谱调制模块、接收模块和计算机;
所述成像模块:通过红外光源照明目标物体之后,使用成像镜头对目标物体的反射光进行收集,成像到数字微镜器件的靶面上;
所述振幅调制模块:数字微镜器件作为空间光调制器,对其靶面上的成像进行振幅调制;数字微镜器件上投影的调制图案是经过随机化的哈达玛矩阵;经过振幅调制之后的信号光由反射镜调整出射方向,出射的信号光进入光谱调制模块;
所述光谱调制模块:使用纳米晶体材料滤波片对从振幅调制模块出射的信号光进行光谱调制;经过振幅调制和光谱调制之后的信号光直接出射后进入接收模块。
所述接收模块:包括接收端的收集透镜、光电探测器以及数据采集器件;收集透镜使用凸透镜,对信号光进行聚焦,汇聚到光电探测器的靶面上;光电探测器是单像素探测器(带有放大电路的光电二极管,又称为桶探测器),将经过调制之后的光信号转化为模拟电压信号,调制过程中不同强度的信号光会转化为不同的电压数值;数据采集器件(数据采集卡)将模拟电压信号转化为计算机能够识别的数字信号,同时具有同步触发功能,实现加载到数字微镜器件上的调制图案与光电探测器输出的电压值信号一一对应的关系;
所述计算机:接收数据采集器件输送的电压数据,通过将光电探测器输出的一维强度信号作为测量值矩阵,将随机化的哈达玛矩阵和纳米晶体材料滤波片的透射光谱作为测量矩阵,进行压缩感知数据处理和运算,重构出目标物体的三维数据立方体。
所述纳米晶体材料滤波片的制备过程如下所述:
(1)将纳米晶体材料分散到有机溶剂中,得到纳米晶体材料溶液;
所述纳米晶体材料如硫化铅、硒化铅、硫化镉、硒化镉等。
所述有机溶剂为正己烷。
(2)玻璃片清洗处理;
将玻璃片放入装有异丙醇的容器中超声清洗三十分钟,烘干后放入真空等离子体清洗机中进行清洗,这样可以提高纳米晶体材料在玻璃片表面的附着效果。
(3)将纳米晶体材料溶液滴加到清洗后的玻璃片上,布满整个玻璃片,罩住玻璃片,静置使溶剂完全挥发,在玻璃片上形成纳米晶体材料薄膜,得到纳米晶体材料滤波片。
通过控制纳米晶体材料薄膜的厚度调控吸收特性,可以通过调整(1)中纳米晶体材料与有机溶剂的比例,得到不同厚度的纳米晶体材料薄膜。本发明中对纳米晶体材料薄膜的厚度不做要求,任何厚度都能达到本发明的目的,较佳的纳米晶体材料薄膜厚度为50nm-200nm。
所述数字微镜器件上投影的调制图案的生成过程如下所述:
调制图案的分辨率与目标物体的分辨率直接相关(如:当目标物体的分辨率为8×8时,生成哈达玛矩阵的维度为64×64);在目标物体的成像生成哈达玛矩阵之后进行随机化处理,即将生成的哈达玛矩阵与一个随机生成的相同大小的二值矩阵进行逻辑异或,之后采用差分测量的方式对目标物体进行振幅调制,差分测量是指在加载调制图案时,除了投影一张已生成的调制图案之外,另外投影一张将该调制图案取反之后得到的调制图案;生成调制图案时将得到的矩阵的所有行都单独提取,并转化行矩阵为一个二维矩阵,并以此作为加载到数字微镜器件上的调制图案。
所述使用纳米晶体材料滤波片对从振幅调制模块出射的信号光进行光谱调制的过程是;
每在光路中放置一个纳米晶体材料滤波片,都进行一次单像素成像,获得一张单像素成像重构结果图像,这张图像包含了在该纳米晶体材料滤波片调制下的图像信息;接着切换到下一个纳米晶体材料滤波片,重复进行单像素成像;这样一直重复,直至切换完所有的k个纳米晶体材料滤波片,获得k张单像素重构结果图像;根据这一系列图像进行光谱重构,具体是提取每一张重构结果图像中相同空间位置的像素点的强度值,再以所有光谱调制作为测量矩阵进行压缩感知计算重构,获得该像素点的光谱信息,接着选择下一个像素点,重复进行光谱重构,这样一直重复,直至完成所有像素点的光谱重构;对得到的所有光谱数据进行排列组合,得到重构结果。
所述计算机重构出目标物体的三维数据立方体的过程如下所述:
设目标物体的空间分辨率为n×n,设目标物体光谱的离散化程度为t,即目标物体T是一个t×n×n的三维数据立方体,这个矩阵经过重组表示为T, T∈R^(t×n2);使用哈达玛矩阵对目标物体进行全采样(即使用与目标物体分辨率相同的采样次数);设投影在数字微镜器件上的调制图案为Pi, Pi∈R^( n2×n2), i=1,2,3……n2,其中i代表投影调制图案的索引;设纳米晶体材料滤波片F共有k个,则Fj∈R^(k×t), j=1,2,3……k,其中j代表使用的纳米晶体材料滤波片的索引,设系统对光谱的响应为D(其中既有单像素探测器的响应率曲线,也有系统中的镜头、汇聚透镜等),这是一个一维向量,即D∈R^t;由此,计算过程表示为:
Y=F×T×P,
其中矩阵Y代表单像素探测器输出的测量值,这是一个k行n2列的二维矩阵;矩阵P代表成像系统对目标物体的空间振幅调制;矩阵F代表成像系统对目标物体的光谱调制,T代表目标物体;
测量值矩阵Y有n2列数据表示数字微镜器件投影了n2次,这个矩阵有k行表示一共使用了k个纳米晶体材料滤波片进行光谱调制,这个矩阵中的元素代表在数字微镜器件第i次投影,在第j个纳米晶体材料滤波片调制下,单像素探测器输出的测量值;测量值矩阵Y矩阵同时包含目标物体的空间信息和光谱信息,通过关联运算对目标物体的信息进行求解;
使用空间光调制器对目标物体经过n2次采样之后,通过关联运算的方式求解出目标物体的空间分布;测量值矩阵Y的每一列都应用单像素成像方式,重构出在第j个纳米晶体材料滤波片调制下的目标图像;测量值矩阵Y的每一行都应用光谱计算重构方式,重构出该像素点的光谱曲线;上述两种方式在获取数据时分步进行,在计算重构时同时进行。
上述系统,为了解除高光谱成像系统对阵列探测器的依赖,并将光谱成像的范围拓展到特殊波段,使用单像素成像技术获取目标物体的成像信息,使用纳米晶体材料滤波片获取目标物体在特殊波段的光谱信息,通过使用单像素探测器输出的一维强度信号进行高维重构,实现高光谱单像素成像。
本发明具有以下特点:
1.应用纳米晶体材料滤波片和单像素探测器,使用非阵列探测器实现特殊波段高光谱成像,通过使用不同数量的自制滤波器(纳米晶体材料滤波片)实现特殊波段内的高光谱成像;如使用50个自制滤波器(纳米晶体材料滤波片)实现了在1050nm-1630nm范围内,平均光谱分辨率为8.59nm的重构质量,图像分辨率达128×128。
2.通过使用滤波器的柔性制备工艺和微机电系统,将高光谱成像系统小型化。
3.解除了高光谱成像技术对阵列探测器的依赖,将光谱成像范围拓展到特殊波段。
附图说明
图1为本发明使用非阵列探测器实现高光谱成像的系统的结构原理框图。
图2为本发明使用非阵列探测器实现高光谱成像的系统的结构示意图。
图3为本发明提供的使用的振幅调制图案的设计示意图。
图4为本发明使用的重构算法的图形化表示示意图。
图中:1.光源,2.目标物体,3.成像镜头,4.反射镜,5.纳米晶体材料滤波片,6.收集透镜,7.单像素探测器,8.数字微镜器件,9.数据采集卡,10.计算机。
具体实施方式
本发明提供了一种应用单像素成像技术和纳米晶体材料(如硫化铅、硒化铅、硫化镉、硒化镉等)、使用非阵列探测器实现高光谱成像的系统,能够在不使用阵列探测器的情况下,实现对目标物体空间信息和光谱信息的采集,获取三维数据立方体,解决了以下问题:
1.使用非阵列探测器实现高光谱成像。
2.将光谱成像的范围拓展到特殊波段。
3.缩小高光谱成像系统的尺寸。
为使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图对本发明作进一步的详细说明,但不仅限于此,本发明未详尽说明的,均按本领域常规技术。另外还需要说明的是,为了方便描述,附图中仅示出了与本发明相关的部分而并非全部结构。
图1给出了本发明使用非阵列探测器实现高光谱成像的系统的研究思路框图。为了实现在红外等特殊波段,不使用阵列探测器实现高光谱成像,本发明主要从用滤波片为调制器件获取目标物体光谱信息为思路,寻找解决方案,以使用数字微镜器件进行振幅调制,使用纳米晶体材料滤波片进行光谱调制为中心思想。
本发明目的的实现过程可简要概括为:
1.制备纳米晶体材料滤波片;
2.搭建包含振幅调制和光谱调制的成像系统;
3.在成像系统中添加纳米晶体材料滤波片;
4.系统对目标探测获取数据,得到一维强度信息;
5.利用压缩感知算法对目标物体进行高维重构。
纳米晶体材料滤波片5的透射光谱与纳米晶体材料的吸收光谱直接相关,通过使用这种滤波片作为系统光谱滤波的方式,能够获取目标物体在特殊波段的光谱信息。
以下详细说明纳米晶体材料滤波片5的制备过程。
将纳米晶体材料(如硫化铅)分散到有机溶剂正己烷中,得到纳米晶体材料溶液。将玻璃片放入装有异丙醇的烧杯中,放入超声波清洗机中超声清洗三十分钟。超声清洗之后将玻璃片烘干,放入真空等离子体清洗机中对表面进行清洗,提高纳米晶体材料在玻璃片表面的附着效果。在经过处理之后,使用移液枪将部分纳米晶体材料溶液滴加到玻璃片的中心,让溶液自行延展直至布满整个玻璃片。在这之后,为了防止纳米晶体材料在沉积过程中,空气中的灰尘与材料混杂,影响制备效果,需要通过防尘罩罩住玻璃片,这一方面能够防止灰尘等杂质污染,也可以减缓纳米晶体材料在玻璃片上的沉积过程,能够增加材料与玻璃片的粘连性。经过一至两天的放置,可以使溶剂完全挥发,完成纳米晶体材料的沉积,在玻璃片上形成一层薄膜,得到纳米晶体材料滤波片。
本发明中对纳米晶体材料薄膜的厚度不做要求,任何厚度都能达到本发明的目的,通过控制纳米晶体材料薄膜的厚度调控吸收特性,可以通过调整纳米晶体材料分散到有机溶剂的比例,得到不同厚度的纳米晶体材料薄膜。较佳的纳米晶体材料薄膜厚度为50nm-200nm。
以下对搭建的成像系统详细说明。
本发明从单像素成像技术为出发点,搭建的成像系统包括成像模块、振幅调制模块、光谱调制模块、接收模块和计算机。以下结合图2对该系统各模块的构成及功能详细说明。
成像模块通过红外光源1照明目标物体2之后,使用成像镜头3对目标物体2的反射光进行收集,成像到数字微镜器件8的靶面上。
振幅调制模块使用数字微镜器件8作为空间光调制器,对其靶面上的像进行振幅调制。数字微镜器件8上投影的调制图案是经过随机化的哈达玛矩阵,经过振幅调制之后的信号光由反射镜4调整出射方向,进入光谱调制模块。
光谱调制模块使用纳米晶体材料滤波片5,对从振幅调制模块(数字微镜器件8)出射的信号光进行光谱调制;经过振幅调制和光谱调制之后的信号光直接出射后,进入接收模块。
接收模块包括接收端的收集透镜6、光电探测器以及数据采集器件。收集透镜6使用凸透镜,对信号光进行聚焦,汇聚到光电探测器的靶面上。
数据采集器件采用数据采集卡9,将模拟电压信号转化为计算机能够识别的数字信号。同时,数据采集卡还具有同步触发功能,实现了加载到数字微镜器件上的调制图案,与单像素探测器7输出的电压值信号一一对应的关系。
数据采集卡7将电压数据上传至计算机10,通过将光电探测器(单像素探测器7)输出的一维强度信号作为测量值矩阵,将随机化的哈达玛矩阵和纳米晶体材料滤波片的透射光谱作为测量矩阵,进行压缩感知数据处理和运算,重构出目标物体的三维数据立方体。
上述系统,通过将成像系统获取的原始数据,利用单像素成像算法和压缩感知算法进行数据处理和运算,就能够得到目标物体的三维数据立方体。其中包括目标物体的二维空间信息,以及目标物体的一维光谱信息。为了更直观的展示系统的计算结果,对数据立方体进行伪彩色赋值。通过将目标物体的光谱曲线,映射到CIE1931色彩空间中的光谱三刺激值曲线中,能够获取目标物体的伪彩色图。
图2给出的系统中,首先使用台灯或太阳光等恒稳光源1照亮实验的目标物体2,其反射光透过成像镜头3成像到数字微镜器件8上。为了得到清晰的像,需要让目标物体2和数字微镜器件8上的像满足成像关系。通过在数字微镜器件8上加载不同的调制图案,实现对所成的像的振幅调制。使用反射镜4调整经过振幅调制的光场的出射方向,使其能够通过纳米晶体材料滤波片5。纳米晶体材料滤波片5根据在玻璃片载体上沉积的纳米晶体材料的不同,完成对入射光场光谱的调制。光场经过调制之后,经由收集透镜6完成聚焦,被单像素探测器7接收。单像素探测器7将两次调制的光场进行光电转换,输出模拟电压值信号,该信号通过数据采集卡9实现模数转换,将模拟电压值信号转为数字信号,传输给计算机10。在计算机10中就可以通过将电压值与加载到数字微镜器件8上的调制图案、纳米晶体材料滤波片5的透射光谱进行关联运算,求解出目标物体2的空间信息和光谱信息,得到目标物体2的三维数据立方体。
每在光路中放置一个纳米晶体材料滤波片5,都进行一次单像素成像,获得一张单像素成像重构结果图像。这张图像就包含了在该纳米晶体材料滤波片调制下的图像信息。接着切换到下一个纳米晶体材料滤波片,重复进行单像素成像。这样一直重复,直至切换完所有的k个纳米晶体材料滤波片,获得k张单像素重构结果图像。下一步是根据这一系列图像进行光谱重构,具体是提取每一张重构结果图像中,相同空间位置的像素点的强度值,再以成像系统中的所有光谱调制作为测量矩阵进行压缩感知计算重构,获得该像素点的光谱信息。接着选择下一个像素点,重复进行光谱重构。这样一直重复,直至完成所有像素点的光谱重构。最后是对得到的所有光谱数据进行排列组合,得到重构结果。
通过光谱重构能够获得每一个像素点的光谱曲线,能够按照像素点的空间位置关系重组出数据立方体,但并不能获得像素点与像素点之间的强度关系。为了获得这种强度关系,需要再进行一次单像素成像实验,这一次的单像素成像实验的目的是为了获取像素点与像素点之间的强度关系,所以不需要在光路中放置纳米晶体材料滤波片。在得到这样的重构结果之后,取其中每个像素点的强度值,作为该像素点的光谱曲线的权重,对直接重组的数据立方体进行加权,就能够得到包含目标物体的空间分布,光谱信息和像素间强度关联的数据立方体,它的每一个沿光谱的切片都是该波长下的单像素图像重构结果。
图3给出了使用的振幅调制图案。由于数字微镜器件8上加载的调制图案是对目标物体2的成像进行的,所以调制图案的分辨率与目标物体2的分辨率直接相关。在图4中,为了方便表示,以目标物体2的分辨率为8×8为例。当目标物体2的分辨率为8×8时,生成哈达玛矩阵的维度为64×64。在生成矩阵之后,对其进行随机化处理,即将生成的哈达玛矩阵与一个随机生成的相同大小的二值矩阵进行逻辑异或。通过这样的方式,能够使得生成的调制图案中,亮点与暗点的数量大致相同,即能够使得单像素探测器7输出的电压值大致相同,接近均值。在这之后,为了提高信噪比,采用差分测量的方式对目标物体2进行振幅调制。差分测量是指在加载调制图案时,除了投影一张已生成的调制图案之外,另外投影一张将该调制图案取反之后得到的调制图案。通过这样的形式能够减少成像系统中始终存在的噪声,如环境噪声或暗噪声等。生成调制图案时,是将得到的矩阵的所有行都单独提取,并转化行矩阵为一个二维矩阵,并以此作为加载到数字微镜器件8上的调制图案。
当目标物体2为8×8时,得到的随机化后的哈达玛矩阵维度为64×64。经过差分处理,能够得到与其维度相同的大小为64×64的调制矩阵。通过提取行矩阵并进行重组,共可生成128张调制图案。
图4是本发明使用的重构算法的图形化表示。计算中,设目标物体T的空间分辨率为n×n,设目标物体光谱的离散化程度为t,即目标物体T是一个t×n×n的三维数据立方体,这个矩阵经过重组可以表示为T, T∈R^(t×n2)。正如前文介绍的那样,在单像素成像系统中使用正交哈达玛矩阵在理论上能够得到目标图像的完美解。所以在成像系统中,为了得到质量最高的图像重构结果,使用哈达玛矩阵对目标物体进行全采样,即使用和目标物体分辨率相同的采样次数。设投影在数字微镜器件上的调制图案为Pi, Pi∈R^( n2×n2), i=1,2,3……n2,其中i代表投影调制图案的索引。计算时还需要考虑成像系统对目标物体的影响,这包括使用的用于对入射光场进行调制的纳米晶体材料滤波片的透射光谱,也包括成像系统本身对光谱的响应。设纳米晶体材料滤波片F共有k个,则Fj∈R^(k×t), j=1,2,3……k,其中j代表使用的纳米晶体材料滤波片的索引。设系统对光谱的响应为D,其中既有桶探测器的响应率曲线,也有系统中的镜头,汇聚透镜等,这是一个一维向量,即D∈R^t。由此,成像系统的计算过程可以表示为:
Y=F×T×P
其中矩阵Y代表单像素探测器7输出的测量值,这是一个k行n2列的二维矩阵。矩阵P代表成像系统对目标物体2的空间振幅调制,矩阵F代表成像系统对目标物体2的光谱调制,T代表目标物体2。
接下来对测量值矩阵Y进行分析。测量值矩阵Y有n2列数据表示数字微镜器件8投影了n2次,这个矩阵有k行表示一共使用了k个纳米晶体材料滤波片5进行光谱调制,这个矩阵中的元素代表在数字微镜器件第i次投影,在第j个纳米晶体材料滤波片调制下,单像素探测器7输出的测量值。这个矩阵同时包含目标物体的空间信息和光谱信息,下一步就是需要通过关联运算对目标物体的信息进行求解。
由单像素关联成像理论可以了解到,使用空间光调制器对目标物体经过n2次采样之后,就可以通过关联运算的方式求解出目标物体的空间分布。测量值矩阵Y的每一列都可以应用单像素成像方式,重构出在第j个纳米晶体材料滤波片调制下的目标图像;测量值矩阵Y的每一行都可以应用光谱计算重构方式,重构出该像素点的光谱曲线。上述两种方式在获取数据时分步进行,在计算重构时同时进行。
使用50个纳米晶体材料滤波片5(纳米晶体材料薄膜厚度为50nm-200nm)实现了在1050nm-1630nm范围内,平均光谱分辨率为8.59nm的重构质量,图像分辨率达128×128。
Claims (6)
1.一种使用非阵列探测器实现高光谱成像的系统,其特征是:包括成像模块、振幅调制模块、光谱调制模块、接收模块和计算机;
所述成像模块:通过红外光源照明目标物体之后,使用成像镜头对目标物体的反射光进行收集,成像到数字微镜器件的靶面上;
所述振幅调制模块:数字微镜器件作为空间光调制器,对其靶面上的成像进行振幅调制;数字微镜器件上投影的调制图案是经过随机化的哈达玛矩阵;经过振幅调制之后的信号光由反射镜调整出射方向,出射的信号光进入光谱调制模块;
所述光谱调制模块:使用纳米晶体材料滤波片对从振幅调制模块出射的信号光进行光谱调制;经过振幅调制和光谱调制之后的信号光直接出射后进入接收模块;
所述接收模块:包括接收端的收集透镜、光电探测器以及数据采集器件;收集透镜使用凸透镜,对信号光进行聚焦,汇聚到光电探测器的靶面上;光电探测器是单像素探测器,将经过调制之后的光信号转化为模拟电压信号,调制过程中不同强度的信号光会转化为不同的电压数值;数据采集器件将模拟电压信号转化为计算机能够识别的数字信号,同时具有同步触发功能,实现加载到数字微镜器件上的调制图案与光电探测器输出的电压值信号一一对应的关系;
所述计算机:接收数据采集器件输送的电压数据,通过将光电探测器输出的一维强度信号作为测量值矩阵,将随机化的哈达玛矩阵和纳米晶体材料滤波片的透射光谱作为测量矩阵,进行压缩感知数据处理和运算,重构出目标物体的三维数据立方体;
所述光谱调制模块中使用纳米晶体材料滤波片对从振幅调制模块出射的信号光进行光谱调制的过程是;
每在光路中放置一个纳米晶体材料滤波片,都进行一次单像素成像,获得一张单像素成像重构结果图像,这张图像包含了在该纳米晶体材料滤波片调制下的图像信息;接着切换到下一个纳米晶体材料滤波片,重复进行单像素成像;这样一直重复,直至切换完所有的k个纳米晶体材料滤波片,获得k张单像素重构结果图像;根据这一系列图像进行光谱重构,具体是提取每一张重构结果图像中相同空间位置的像素点的强度值,再以所有光谱调制作为测量矩阵进行压缩感知计算重构,获得该像素点的光谱信息,接着选择下一个像素点,重复进行光谱重构,这样一直重复,直至完成所有像素点的光谱重构;对得到的所有光谱数据进行排列组合,得到重构结果。
2.根据权利要求1所述的使用非阵列探测器实现高光谱成像的系统,其特征是,所述振幅调制模块中数字微镜器件上投影的调制图案的生成过程如下所述:
调制图案的分辨率与目标物体的分辨率直接相关;在目标物体的成像生成哈达玛矩阵之后进行随机化处理,即将生成的哈达玛矩阵与一个随机生成的相同大小的二值矩阵进行逻辑异或,之后采用差分测量的方式对目标物体进行振幅调制,差分测量是指在加载调制图案时,除了投影一张已生成的调制图案之外,另外投影一张将该调制图案取反之后得到的调制图案;生成调制图案时将得到的矩阵的所有行都单独提取,并转化行矩阵为一个二维矩阵,并以此作为加载到数字微镜器件上的调制图案。
3.根据权利要求1所述的使用非阵列探测器实现高光谱成像的系统,其特征是,所述光谱调制模块中纳米晶体材料滤波片的制备过程如下所述:
(1)将纳米晶体材料分散到有机溶剂中,得到纳米晶体材料溶液;
(2)玻璃片清洗处理;
(3)将纳米晶体材料溶液滴加到清洗后的玻璃片上,布满整个玻璃片,罩住玻璃片,静置使溶剂完全挥发,在玻璃片上形成纳米晶体材料薄膜,得到纳米晶体材料滤波片。
4.根据权利要求3所述的使用非阵列探测器实现高光谱成像的系统,其特征是,所述(2)玻璃片清洗处理的过程是:将玻璃片放入装有异丙醇的容器中超声清洗三十分钟,烘干后放入真空等离子体清洗机中进行清洗。
5.根据权利要求3所述的使用非阵列探测器实现高光谱成像的系统,其特征是,所述(3)中纳米晶体材料薄膜的厚度为50nm-200nm。
6.根据权利要求1所述的使用非阵列探测器实现高光谱成像的系统,其特征是,所述计算机重构出目标物体的三维数据立方体的过程如下所述:
设目标物体的空间分辨率为n×n,设目标物体光谱的离散化程度为t,即目标物体T是一个t×n×n的三维数据立方体,这个矩阵经过重组表示为T, T∈R^(t×n2 );使用哈达玛矩阵对目标物体进行全采样;设投影在数字微镜器件上的调制图案为Pi, Pi∈R^( n2×n2), i=1,2,3……n2,其中i代表投影调制图案的索引;设纳米晶体材料滤波片F共有k个,则Fj∈R^(k×t), j=1,2,3……k,其中j代表使用的纳米晶体材料滤波片的索引,设系统对光谱的响应为D,这是一个一维向量,即D∈R^t;由此,计算过程表示为:
Y=F×T×P,
其中矩阵Y代表单像素探测器输出的测量值,这是一个k行n2列的二维矩阵;矩阵P代表成像系统对目标物体的空间振幅调制;矩阵F代表成像系统对目标物体的光谱调制,T代表目标物体;
测量值矩阵Y有n2列数据表示数字微镜器件投影了n2次,这个矩阵有k行表示一共使用了k个纳米晶体材料滤波片进行光谱调制,这个矩阵中的元素代表在数字微镜器件第i次投影,在第j个纳米晶体材料滤波片调制下,单像素探测器输出的测量值;测量值矩阵Y同时包含目标物体的空间信息和光谱信息,通过关联运算对目标物体的信息进行求解;
使用数字微镜器件对目标物体经过n2次采样之后,通过关联运算的方式求解出目标物体的空间分布;测量值矩阵Y的每一列都应用单像素成像方式,重构出在第j个纳米晶体材料滤波片调制下的目标图像;测量值矩阵Y的每一行都应用光谱计算重构方式,重构出该像素点的光谱曲线;上述两种方式在获取数据时分步进行,在计算重构时同时进行。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311107878.0A CN117091700B (zh) | 2023-08-31 | 2023-08-31 | 一种使用非阵列探测器实现高光谱成像的系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311107878.0A CN117091700B (zh) | 2023-08-31 | 2023-08-31 | 一种使用非阵列探测器实现高光谱成像的系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117091700A CN117091700A (zh) | 2023-11-21 |
CN117091700B true CN117091700B (zh) | 2024-02-06 |
Family
ID=88781139
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311107878.0A Active CN117091700B (zh) | 2023-08-31 | 2023-08-31 | 一种使用非阵列探测器实现高光谱成像的系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117091700B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109923401A (zh) * | 2016-11-08 | 2019-06-21 | 南加州大学 | 高光谱成像系统 |
CN111856497A (zh) * | 2020-07-29 | 2020-10-30 | 北京深测科技有限公司 | 一种单光子成像方法和系统 |
CN113709325A (zh) * | 2021-07-20 | 2021-11-26 | 辽宁大学 | 基于Hadamard频域变换矩阵阈值滤波的单像素成像方法 |
CN216852177U (zh) * | 2022-03-14 | 2022-06-28 | 河北大学 | 一种基于树莓派的彩色单像素成像系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3742135B1 (fr) * | 2019-05-20 | 2022-01-19 | Centre National de la Recherche Scientifique | Imagerie mono-pixel hyper-spectrale résolue en temps |
-
2023
- 2023-08-31 CN CN202311107878.0A patent/CN117091700B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109923401A (zh) * | 2016-11-08 | 2019-06-21 | 南加州大学 | 高光谱成像系统 |
CN111856497A (zh) * | 2020-07-29 | 2020-10-30 | 北京深测科技有限公司 | 一种单光子成像方法和系统 |
CN113709325A (zh) * | 2021-07-20 | 2021-11-26 | 辽宁大学 | 基于Hadamard频域变换矩阵阈值滤波的单像素成像方法 |
CN216852177U (zh) * | 2022-03-14 | 2022-06-28 | 河北大学 | 一种基于树莓派的彩色单像素成像系统 |
Non-Patent Citations (2)
Title |
---|
Quantitative imaging for optical field via a single-pixel detector;Xianye Li 等;Signal Processing;第188卷;1-9 * |
单像素成像在特殊波段及三维成像的应用发展;孙宝清 等;红外与激光工程;第49卷(第3期);0303016-1—0303016-13 * |
Also Published As
Publication number | Publication date |
---|---|
CN117091700A (zh) | 2023-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cao et al. | A prism-mask system for multispectral video acquisition | |
US11175489B2 (en) | Smart coded access optical sensor | |
EP2603828B1 (en) | Decreasing image acquisition time for compressive imaging devices | |
US8970740B2 (en) | Overlap patterns and image stitching for multiple-detector compressive-sensing camera | |
US8885073B2 (en) | Dedicated power meter to measure background light level in compressive imaging system | |
US20130002715A1 (en) | Image Sequence Reconstruction based on Overlapping Measurement Subsets | |
US8717484B2 (en) | TIR prism to separate incident light and modulated light in compressive imaging device | |
CN106993121B (zh) | 一种基于压缩感知的高光谱图像采集成像系统与控制方法 | |
CN109828285B (zh) | 一种双波段时域压缩感知高速成像方法及装置 | |
CN107942338B (zh) | 一种基于数字微镜器件的多波长关联成像系统 | |
US11231323B2 (en) | Time-resolved hyper-spectral single-pixel imaging | |
Osorio Quero et al. | Single-pixel imaging: An overview of different methods to be used for 3D space reconstruction in harsh environments | |
CN114279564B (zh) | 一种并行压缩感知计算层析成像光谱仪及其成像光谱重建方法 | |
CN114264370B (zh) | 一种压缩感知计算层析成像光谱仪系统和成像方法 | |
CN114659634A (zh) | 微型快照式压缩光谱成像探测装置及探测方法 | |
Murguia et al. | Compact visible/near-infrared hyperspectral imager | |
CN117091700B (zh) | 一种使用非阵列探测器实现高光谱成像的系统 | |
CN113884234A (zh) | 一种互补单像素质心探测系统及方法 | |
CN114689174A (zh) | 一种芯片级多光谱相机系统及其操作方法 | |
CN210862922U (zh) | 一种基于可调谐滤光片的单像素光谱成像系统 | |
CN112556849B (zh) | 高光谱成像装置 | |
US20030133109A1 (en) | Real time LASER and LED detection system using a hyperspectral imager | |
CN114895449B (zh) | 一种基于压缩感知的四维高速荧光显微成像装置 | |
CN114923569A (zh) | 一种基于计算成像的多光谱相机 | |
Martins et al. | Fast hyperspectral single-pixel imaging |
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 |