CN109990899B - 转换均匀或非均匀采样干涉图产生光谱数据的方法和设备 - Google Patents
转换均匀或非均匀采样干涉图产生光谱数据的方法和设备 Download PDFInfo
- Publication number
- CN109990899B CN109990899B CN201811476928.1A CN201811476928A CN109990899B CN 109990899 B CN109990899 B CN 109990899B CN 201811476928 A CN201811476928 A CN 201811476928A CN 109990899 B CN109990899 B CN 109990899B
- Authority
- CN
- China
- Prior art keywords
- matrix
- interferogram
- interferometer
- data cube
- wavelength
- 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
- 238000000034 method Methods 0.000 title claims description 29
- 230000003595 spectral effect Effects 0.000 title claims description 23
- 239000011159 matrix material Substances 0.000 claims abstract description 88
- 230000001934 delay Effects 0.000 claims abstract description 24
- 230000000737 periodic effect Effects 0.000 claims abstract description 16
- 230000006870 function Effects 0.000 claims description 36
- 230000003287 optical effect Effects 0.000 claims description 36
- 230000010287 polarization Effects 0.000 claims description 28
- 239000004973 liquid crystal related substance Substances 0.000 claims description 27
- 230000001419 dependent effect Effects 0.000 claims description 21
- 238000012545 processing Methods 0.000 claims description 13
- 230000004044 response Effects 0.000 claims description 8
- 230000010363 phase shift Effects 0.000 claims description 4
- 230000010355 oscillation Effects 0.000 claims description 2
- 230000007704 transition Effects 0.000 claims description 2
- 230000000638 stimulation Effects 0.000 claims 2
- 238000005314 correlation function Methods 0.000 claims 1
- 238000003491 array Methods 0.000 abstract description 2
- 230000000875 corresponding effect Effects 0.000 description 23
- 238000005070 sampling Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 238000001228 spectrum Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 239000000523 sample Substances 0.000 description 6
- 238000005259 measurement Methods 0.000 description 5
- 210000004027 cell Anatomy 0.000 description 4
- 238000000701 chemical imaging Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 210000002858 crystal cell Anatomy 0.000 description 3
- 239000006185 dispersion Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 238000000844 transformation 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/02—Details
- G01J3/0275—Details making use of sensor-related data, e.g. for identification of sensor parts or optical elements
-
- 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/0224—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows using polarising or depolarising elements
-
- 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/0256—Compact construction
-
- 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
-
- 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
-
- 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/45—Interferometric 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/45—Interferometric spectrometry
- G01J3/453—Interferometric spectrometry by correlation of the amplitudes
- G01J3/4531—Devices without moving parts
-
- 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/45—Interferometric spectrometry
- G01J3/453—Interferometric spectrometry by correlation of the amplitudes
- G01J3/4537—Devices with refractive scan
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
-
- 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/283—Investigating the spectrum computer-interfaced
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
用于计算高光谱数据立方体的重建矩阵包含周期函数行。所述重建矩阵的每一行对应于选定波长,并且每一列对应于干涉仪的选定延迟。所述周期函数具有所述对应行的作为参数的所述选定波长,并以所述对应列中的每一列的所述选定延迟进行采样。获得干涉图数据立方体并且所述干涉图数据立方体包含一个或多个同时测量的干涉图的阵列。所述干涉图数据立方体的每一行对应于所述选定延迟中的一个,并且每一列对应于与所述同时测量的干涉图不同的干涉图。通过将所述重建矩阵与所述干涉图数据立方体的列相乘来形成针对所述干涉图中的每一个的一组矩阵向量积,从而形成所述高光谱数据立方体。
Description
发明内容
本公开涉及一种用于转换均匀或非均匀采样的干涉图以产生预定义波长处的光谱数据的方法和设备。在一个实施例中,确定干涉仪的一组参考延迟,以及对应于高光谱数据立方体的光谱切片的一组波长。形成包含周期函数行的重建矩阵。所述重建矩阵的每一行对应于所述一组波长中的选定波长,并且所述重建矩阵的每一列对应于所述参考延迟中的选定延迟。所述周期函数具有所述对应行的作为参数的所述选定波长,并以所述对应列中的每一列的所述选定延迟进行采样。获得包含一个或多个同时测量的干涉图的阵列的干涉图数据立方体。所述干涉图数据立方体的每一行对应于所述选定延迟中的一个,并且所述干涉图数据立方体的每一列对应于与所述同时测量的干涉图不同的干涉图。形成针对所述干涉图中的每一个的一组矩阵向量积。所述矩阵向量积中的每一个包含所述重建矩阵与所述干涉图数据立方体的列的矩阵乘法。所述一组矩阵向量积形成所述高光谱数据立方体。
附图说明
图1是根据示例实施例的设备的图示;
图2是根据示例实施例的示出通过光学装置处理的干涉图的一组曲线;
图3是根据示例实施例的用于光学装置的延迟与时间的关系图;
图4是根据示例实施例的示出干涉仪的零路径延时的波长相关的测量结果的一组曲线;
图5是根据示例实施例的重建矩阵的图示;
图6是根据示例实施例的干涉图数据立方体的示意图;
图7是根据示例实施例的表示干涉图数据立方体的矩阵的图示;
图8是根据示例实施例的系统响应矩阵的图示;
图9是根据示例实施例的示出斜坡和窗口函数对重建矩阵的影响的一组曲线;以及
图10是根据示例实施例的方法的流程图。
具体实施方式
本公开涉及高光谱成像系统中的信号处理。本文描述的高光谱成像系统使用偏振干涉仪,所述偏振干涉仪被配置成在穿过干涉仪的光的分量中引入可变光路延时(或延迟)。将引起路径延时的装置(称为可变光学延迟器)放置在两个偏振器之间,使得在入射偏振方向上的第一光线与正交偏振方向上的第二光线之间引入可变路径延时(例如,普通光线和非常光线),两组光线追踪共同的路径。此路径延时导致第一光线与第二光线之间波长相关的相移。路径延时使离开偏振干涉仪的光形成通过光学传感器检测的干涉图,例如焦平面阵列。
偏振干涉仪可以使用一个或多个液晶(LC)单元作为可变光学延迟器。这种可变光学延迟器在本文中称为液晶可变延迟器(LCVR)。通常,液晶(LC) 材料是具有可以通过施加外部刺激,如电场或磁场而选择性地改变的一些结晶性质(例如,内部结构的朝向,如指示LC分子的局部平均配向的LC指向矢) 的液体。LC指向矢的朝向变化改变LC材料的光学性质,例如改变LC双折射的光轴。
LCVR在穿过液晶的两个正交偏振光之间生成可变光路延时或可变延迟。 LCVR内的一个或多个液晶单元用作电调谐双折射元件。通过改变跨液晶单元的电极的电压,单元分子改变它们的朝向,并且在一段时间内可控制地改变光路延时是可能的。
为了创建具有LCVR的偏振干涉仪,将LCVR放置在具有标称平行或垂直偏振轴的第一偏振器与第二偏振器之间。LCVR的慢轴(具有可变光路延时的偏振轴)相对于第一偏振器的偏振方向朝向标称45度。入射光被第一偏振器偏振到入射偏振方向。因为LCVR的慢轴相对于此入射偏振方向为45度,所以可以根据平行于LCVR的慢轴偏振的一部分光和垂直于此轴偏振的一部分光来描述偏振入射光。
当光通过LCVR时,它获得第一偏振与第二偏振之间的波长相关的相对相移,从而导致偏振态的波长相关的变化。与第一偏振器平行或垂直朝向的第二偏振器或分析器干扰平行于LCVR的慢轴偏振的那部分光,所述部分光垂直偏振,将LCVR的输出端处的波长相关的偏振态改变为可以由光学检测器或焦平面阵列感测的波长相关的强度图案。通过在改变LCVR的延迟的同时感测此强度,测量入射光的干涉图是可能的,其可以用于确定入射光的光谱性质。
如上所述,偏振干涉仪用于高光谱成像应用中,因为它能够将入射光的光谱信息编码成强度图案,用非光谱分辨检测器容易地测量所述强度图案。高光谱成像是指用于获取高光谱数据集或数据立方体的方法和装置,所述高光谱数据集或数据立方体可以包含在每个像素处提供了密集采样、精细分辨的光谱信息的图像。
由偏振干涉仪提供的波长相关的强度图案大致对应于入射光光谱的余弦变换。通过在偏振干涉仪的输出端处将每个像素的强度图案记录成LCVR的延迟的函数,可以同时对通过LCVR成像的场景的所有点生成的干涉图进行采样。由此,可以通过对记录的空间相关的干涉图沿着延迟轴应用如反余弦变换或傅里叶变换的变换,来名义上恢复高光谱数据立方体。
基于液晶偏振干涉仪的高光谱成像仪生成可能相对于光路延时是非均匀采样的干涉图。标准傅里叶变换对于将由这种装置生成的原始数据转换为高光谱图像可能不是最佳的。本公开描述了用于重建高光谱数据立方体的方法和设备,其具有比傅里叶变换更好的精度和更少的伪像,无论干涉图是否被均匀地采样。
为了增强对本文描述的高光谱图像处理的理解,图1中的框图示出了根据示例实施例的执行图像处理的设备100。设备100包含装置控制器102,所述装置控制器可以包含一个或多个处理器,如中央处理器、子处理器、图形处理器、数字信号处理器等。控制器102耦合到存储器104,所述存储器包含下文将更详细描述的功能模块。存储器104可以包含易失性和非易失性存储器的组合,并且可以存储本领域已知的指令和数据。
设备100包含光学部分106,所述光学部分具有从设备100外部接收光的外部光学接口108。外部光学接口108可以包含窗口、透镜、滤光片、光圈等,适合于将光109从设备100外部传递到内部光学部件。在此实例中,外部光学接口108被示出耦合到外部透镜110。
偏振干涉仪112位于设备100的光学部分106中。偏振干涉仪112例如通过电信号线耦合到控制器102。控制器102将信号施加到偏振干涉仪112,以在干涉仪112的一部分LCVR112a中引起时变光路延时或延迟。此时变光路延时在光109的不同偏振之间引起波长相关的相移,导致作为光路延时的函数变化的干涉图。由图像传感器114(例如,传感器像素阵列或焦平面阵列)检测干涉图,所述图像传感器也耦合到控制器102。
延迟控制器118指示装置控制器102将控制信号施加到LCVR 112a以实现时变延迟轨迹。图像处理器120使用此延迟轨迹作为时变路径延时以及在图像传感器114处检测到的干涉图的测量。可以通过计算作为LCVR 112a的对应位置处的路径延时的函数的变换来处理每个检测到的干涉图,并且一起处理的干涉图包括空间相关的光谱数据,例如高光谱数据立方体。
本文描述的所有处理可以发生在设备100内或通过由接口122或其某种组合耦合的外部计算机124。这由在计算机124上运行的图像处理器126指示。计算机124和/或设备100可以包含图形处理单元(GPU)128、130,其可以用于加速某些计算,如下所述。通常,GPU128、130将具有大量处理核(例如,数百或数千),其被配置用于可并行化的子计算,如矩阵或子矩阵乘法。假设子计算中的每一个是独立的,GPU可以大大减少图像处理时间。由于下面的计算涉及可并行化的矩阵乘法,因此可以使用GPU以提高处理速度。
本文描述的高光谱成像仪在每个像素处独立地生成相对于光路延时的干涉图的原始信号。当这些系统基于液晶偏振干涉仪时,干涉图的采样相对于光路延时可能不完全均匀。结果,典型的快速傅里叶变换(FFT)算法可能无法从干涉图计算高光谱数据立方体。更进一步,即使干涉图采样是均匀的,干涉图的FFT可能不会导致光谱数据的最佳估计,例如,高光谱数据立方体。可以有意地引入非均匀采样,例如通过在较短路径延时处的较密集采样,以增强某些波长处的灵敏度。
在将FFT应用于干涉图之前可以对非均匀采样进行补救的一种方法是通过插值以均匀的光路延时或延迟间隔对干涉图重新采样。可以假设每个记录的干涉图样本的光路延时是已知的,例如通过获得作为时间函数的已知波长的单色参考光源的干涉图的相位。如果通常情况下光路延时作为时间的函数单调地改变,则光路延时与时间之间将存在一对一的关系。将在已知延迟下获取各个干涉图样本,并且可以使用均匀间隔的延迟网格通过内插干涉图样本来计算内插干涉图。
然而,这种插值方法在计算上可能是昂贵的,尤其是如果使用高阶(例如,三次样条)方法或图像尺寸大的话。此外,在干涉图的可以在奈奎斯特极限附近采样的区域中,或者每个干涉条纹可能具有很少的样本(特别是对于接近零路径延时的较短波长),插值可能导致信号的显著损失和/或失真,将是不利的并且可能导致重建的光谱数据中的伪像。本文描述的高光谱成像仪可能已经对较短波长不太敏感,并且这种插值仅仅加剧了这个问题。
在图2中,一组曲线示出了根据示例实施例的不同采样方法如何影响干涉图的信号强度和失真。顶部曲线200是基于LCVR通过偏振干涉仪测量的蓝色 LED的干涉图。以均匀间隔的时间间隔对干涉图进行采样,同时以标称恒定速率扫描LCVR的延迟。然而,由于LCVR的扫描速率的小波动,顶部曲线200 中的干涉图以非均匀的延迟间隔被采样。中间曲线202示出了以均匀延迟增量在顶部曲线200中通过内插重新采样干涉图的结果。这种插值将干涉图信号强度(通过干涉图的标准偏差测量)降低了近20%。曲线204示出了相同蓝色LED的单独测量的干涉图,其相对于延迟有意地非均匀地采样(尽管在均匀间隔的时间间隔下,如图3所示),以增加接近零路径延时的样本数量206到208,如在曲线200、202、204中可以看到的那样,位于短的(蓝色)波长的大部分信号的位置。在此实例中,204中干涉图的有意的非均匀采样将信号相对于干涉图的标准偏差在200中提高了约30%。通常,来自跨可检测光谱的波长(例如,对于CMOS检测器为400nm到1000nm)的贡献将存在于测量的干涉图中。上述有意的非均匀延迟采样可以具有相对于较长波长在较短波长处增加信号强度的效果,这通常需要维持跨可检测光谱的动态范围。
也可以通过用于非均匀采样数据的傅立叶变换开发的许多方法来处理非均匀延迟采样,如非均匀离散傅立叶变换,可以用NuFFT软件库计算。然而,跨所有频率的傅里叶变换可能是过度的,因为许多计算的频率甚至可能不对应于可检测光谱内的波长。相反,图像处理器120可以通过仅利用感兴趣的(均匀或非均匀采样的)傅立叶分量的干涉图的点积来近似傅立叶变换,或者等效地,仅仅那些(常规或非均匀的)离散傅里叶变换(DFT)矩阵行,有助于感兴趣的波长。
这种方法的优点是可以独立处理每个波长,以并入不同的切趾和波长相关的零路径延时,并且傅里叶分量可以采用与干涉图相同的非均匀性进行采样,使得非均匀性不再是问题。任何额外的线性变换也可以与此矩阵组合以形成整体“重建”矩阵,用于直接重建来自原始干涉图的光谱数据。这种矩阵甚至可以具有空间索引以结合空间校准,如空间相关的延迟。
为了生成重建矩阵,从传感器获得参考延迟对时间的曲线,如图3中的示例曲线300所示。例如,位于液晶单元一侧的单色光源116(参见图1)可以用于获得此曲线300。来自光源116的光由传感器114检测,而控制器102将电压施加到LCVR 112a的一个或多个LC单元,导致LCVR的延迟从第一延迟转变到第二延迟。然后,控制器102分析检测到的光源116的强度振荡,以确定时间相关的延迟,或等效地,对于干涉图的每个样本,光源116的干涉图的相位。获得曲线300的数据的其它方式包括使用类似的过程分析一个或多个LC 单元的电容测量。
系统可以具有零路径延时的波长和温度相关变化,这可能例如通过用另一种双折射材料补偿残余液晶延迟来产生。期望在正路径延时和负路径延时两者处测量干涉图,并且在一些情况下,残余液晶延迟防止这种情况并且必须进行补偿。补偿残余延迟的一种方式是在偏振干涉仪内包含相反符号的固定延迟,例如,具有垂直于LCVR的慢轴的慢轴的多阶波片。如果此补偿波片具有以不同于LCVR的延迟的方式取决于波长和/或温度的延迟,则偏振干涉仪的零路径延时将取决于波长和温度两者。在这种情况下,系统还可以针对每个波长并因此针对重建矩阵的每一行确定并存储相对于参考延迟曲线的零路径延时点。图 4中的一组曲线中示出了根据示例实施例确定零路径延时的实例。
在图4中,曲线400是以10nm增量步进的窄带光源的干涉图的一组测量,示出的每个堆叠在另一个上。波长沿垂直轴增加。曲线410是曲线400的部分 402的展开视图,并且曲线420是曲线410的部分412的子集。曲线420中突出显示的点422指示参考延迟值,其表示每个波长的实际零路径延时,σ(λ) =0,其中σ是实际干涉仪路径延时而不是测量的参考延迟。这些点422形成曲线而不是垂直线,因为一些液晶延迟由不同色散的波片补偿,因此导致跨可检测光谱的残余色散。
此残余色散与液晶材料和补偿波片两者的折射指数之和有关。如曲线430 和等式432所示,可以使用柯西近似(例如,高达6阶)来拟合数据。注意,等式的系数通常取决于温度。使用等式432,可以在柯西近似有效的任何波长处计算对应于零路径延时的参考延迟。假设已经表征了相对于温度的系数变化,也可以针对不同温度求解此等式。系数可以存储在设备100的存储器中,如本领域中已知的,并且由图像处理器在下面描述的计算中使用。
可以由图像处理器120使用的另一个参数是感兴趣的给定波长的相位延时的变化率(例如,λi,对应于重建矩阵的ith行)相对的变化率在参考波长λ0处的参考延迟曲线的相位。可以获得比例因子,例如,通过液晶折射指数的知识,作为波长λ和温度T的函数,其中e-和o-下标表示非常光线和普通光线的相应指数。以下等式(1)可以用于在波长λi处确定缩放因子αi或校准可以用于经验地确定值。
鉴于上述的,重建矩阵的元件Aij然后可以如等式(2)表示如下,其中Φ是对应于干涉图样本的参考延迟曲线的相位。相位Φ是通过分析参考光源的干涉图经验地发现的,例如,使用武田的方法;它涉及到参考延迟Γj作为Φj= 2πΓj/λ0角度是对应于波长λ的零路径延时参考延迟曲线的相位,例如,如等式 (3)中通过图4中获得的拟合系数计算。从参考延迟曲线的相位减去相位φiΦj对行中采样的周期函数进行相移i使得它在对应于波长处的干涉仪的零延迟点具有0相位λi。
Aij=cos(αi(Φj-φi));i∈(1...M},j∈{1...N} (2)
注意,每个图像采集只需要计算等式(2)中的重建矩阵元素,因为每次采集都具有特定的参考延迟曲线和温度。在图5中,图表示出了根据示例实施例的重建矩阵500。矩阵500包含如等式(2)中的周期函数502(例如,余弦函数)的行。每一行对应于一组波长504的选定波长,每个波长对应于高光谱数据立方体的光谱切片。可以预先选择波长504以具有特定的范围和分布。注意,先前描述的参考波长0可以或可以不被包含在一组波长504中。
矩阵500的每一列对应于一组参考延迟506的选定延迟。一组参考延迟506 与特定时间点下的干涉仪的状态相关联。注意,一组延迟506可以以非均匀延迟间隔分布(例如,如图3所示),并且初始延迟Γ1可以是或可以不是最小的延迟。
周期函数502具有对应行的选定波长504(或选定波长504的指数)作为参数,并在对应列中的每一个的参考延迟506处进行采样。矩阵500可以与干涉图数据立方体相乘以获得光谱数据。在图6中,示意图示出了根据示例实施例的干涉图数据立方体。元件600表示检测器的各个传感器,例如,焦平面阵列的像素,跨xy平面布置。t轴表示从Γ1到FN的偏振干涉仪的延迟的变化,其中每个对应于获取数据立方体的t1与tN之间的不同时间点。曲线602表示在采集时间段内在每个传感器600处测量的各个干涉图。因此,干涉图数据立方体由阵列形成,每个元件具有对应于特定传感器位置和延迟值的样本值。点604 表示对应于传感器600的延迟Γj的各个样本。
虽然示出的传感器600处于二维空间布置中,但是出于计算的目的,它们可以由单个索引标记,例如,pk,k∈{1...P}。这样,干涉图数据立方体可以表示为二维矩阵X 700,如图7所示。矩阵700包含对应于同时测量的干涉图的样本702。矩阵700的每一行对应于参考延迟506中的选定延迟,并且每一列对应于来自具有索引k的同时测量的干涉图的组的不同干涉图测量。如此处所示,干涉图样本由传感器位置704标记,传感器阵列具有P个位置,例如,P=像素行数*像素列数。
光谱数据或高光谱数据立方体的计算涉及形成用于干涉图中的每一个的一组矩阵向量积,矩阵向量积中的每一个包括矩阵500与干涉图数据立方体700 的列k的矩阵乘法。因此,所得高光谱数据立方体H可以通过矩阵乘法H=AX 形成。H的每一行可以组成对应于波长504中的一个的图像,图像的每个像素对应于对应波长下的强度和用索引k标记的空间位置。
在一些实施例中,可以为干涉仪的至少两个空间区域构建不同的重建矩阵 500。因此,图7中的不同像素位置pk可以乘以不同的重建矩阵。这在图7中由区域705到707指示,其可以各自与不同的重建矩阵500相关联。这种区域可以涵盖一个或多个单独的像素,并且还可以涵盖数据立方体矩阵700的非连续部分。以这种方式,可以以任何期望的粒度对偏振干涉仪(或其它光学部件) 的空间相关特性进行补偿。
从相对于零路径延时不对称的干涉图重建光谱数据需要特殊处理以防止形成低分辨率的光谱伪像(参见例如C.D.Porter和D.B.Tanner,“傅里叶光谱中相位误差的校正(Correction of phase errors in Fourier spectroscopy)”,《国际红外和毫米波杂志(Int.J.Infrared Millimeter Waves)》4,273-298(1983))。将这种特殊处理结合到重建矩阵500中的一种方法是将矩阵500的每一行与线性“斜坡”函数相乘,所述函数在最接近零路径延时的末端开始等于零,在(波长相关的)零路径延时点处经过一半,并继续到具有相同斜率的一个,此时它保持不变。在不失一般性的情况下,可以假设参考延迟曲线的相位对于Φj,j=1开始最小为0,并且为了增加样本数j而单调地增加。然后,波长相关的斜坡可以如下面的等式(4)中定义。等式(4)定义了矩阵,所述矩阵可以与重建矩阵成元相乘,以适应不对称干涉图。
重建矩阵500的行也可以乘以窗口或切趾函数,只要它对于每个波长的零路径延时是对称的。逐行计算重建矩阵的一个优点是可以对每个波长应用不同尺寸的窗口函数。例如,由于跨LCVR的固有路径延时不均匀性可能限制较短波长的最佳可获得光谱分辨率,较短波长通常生成干涉图,其大部分信号能量位于零路径延时附近。因此,当从干涉图重建光谱的较短波长部分时,有利的是首先将干涉图乘以窗口函数,所述窗口函数仅包含存在来自较短波长的信号的那些部分。使窗口变大会增加噪声但不会增加信号,因为与信号不同,噪声通常跨所有路径延时分布。并且,为了避免在具有一个或多个相同窗口函数的不同位置处重复干涉图的乘法,窗口函数可以与重建矩阵500的行元素相乘。
如果是重建矩阵500和是与在由重建矩阵500的行表示各个波长的系统响应列的矩阵,然后重建矩阵500可以预先乘以对角化矩阵()-1,以“对角化”重建。也就是说,全重构矩阵可以被设置为()-1因此由系统响应矩阵乘法(或通过在给定波长下测量的响应)产生同一性。这类似于从干涉图重建光谱数据的最小二乘解决方案。另外,由于可以是病态的,可能有必要提供某种形式的正则化,以防止上述的矩阵求逆的发散。例如,可以在反转之前添加同一性的倍数:(+)-1,其中,I是M×M的单位矩阵,β是乘法因子。重建矩阵500还可以包含关于干涉仪的光谱灵敏度的信息,以便均衡跨波长的系统响应。图8 中的图表示出了系统响应矩阵B 800。矩阵800的列是由对应于A的行中的每一个的波长504处的干涉仪的单位波长刺激生成的预期原始干涉图。
在图9中,一对堆叠曲线示出了如上所述的两个重建矩阵,并且与图4中的测量类似地绘制。曲线900的每一行是如等式(2)中的重建矩阵A的行,考虑非均匀延迟间隔、波长相关零路径延时和液晶色散。曲线902示出了额外处理的效果,其包含将每行乘以等式(4)的斜坡函数以避免对干涉图的对称部分进行双重计数,以及波长相关的窗口以分别定义每个波长的分辨率仓大小。
在图10中,流程图示出了根据示例实施例的方法。方法涉及确定1000干涉仪的一组参考延迟。例如,基于高光谱数据立方体的期望光谱范围,还确定 1001对应于高光谱数据立方体的光谱切片的一组波长。形成重建矩阵1002。重建矩阵具有多行采样周期函数。所述重建矩阵的每一行对应于所述一组波长中的选定波长,并且所述重建矩阵的每一列对应于所述参考延迟中的选定延迟。周期函数具有对应行的作为参数的选定波长,并且在对应列中的每一个的参考延迟处进行采样。
例如通过光学检测器获得干涉图数据立方体1003。干涉图数据立方体包含一个或多个同时测量的干涉图的阵列。可以相对于光路延时对同时测量的干涉图进行非均匀地采样。干涉图数据立方体的每一行对应于选定的延迟中的一个,并且干涉图数据立方体的每一列对应于与同时测量的干涉图的组不同的干涉图。形成用于干涉图中的每一个的一组矩阵向量积1004。所述矩阵向量积中的每一个包含所述重建矩阵与所述干涉图数据立方体的列的矩阵乘法。所述一组矩阵向量积形成所述高光谱数据立方体。
Claims (19)
1.一种方法,包括:
确定干涉仪的一组参考延迟;
确定对应于高光谱数据立方体的光谱切片的一组波长;
形成包括周期函数的行的重建矩阵,其中所述重建矩阵的每一行对应于所述一组波长中的选定波长,并且所述重建矩阵的每一列对应于所述参考延迟中的选定延迟,其中所述周期函数具有对应行的所述选定波长作为参数,并以对应列中的每一列的所述选定延迟进行采样;
获得包括一个或多个同时测量的干涉图的阵列的干涉图数据立方体,其中所述干涉图数据立方体的每一行对应于所述选定延迟中的一个,并且所述干涉图数据立方体的每一列对应于与所述同时测量的干涉图不同的干涉图;以及
形成针对所述干涉图中的每一个的一组矩阵向量积,所述矩阵向量积中的每一个包括所述重建矩阵与所述干涉图数据立方体的列的矩阵乘法,所述一组矩阵向量积形成所述高光谱数据立方体。
2.根据权利要求1所述的方法,其中所述干涉图对应于所述干涉仪内的不同图像点。
3.根据权利要求1所述的方法,其中所述同时测量的干涉图是相对于光路延时不均匀地采样的。
4.根据权利要求1所述的方法,其中所述干涉仪包括偏振干涉仪。
5.根据权利要求4所述的方法,其中所述偏振干涉仪包括液晶可变延迟器,并且其中每一行中的所述周期函数的相对振荡速率基于所述液晶可变延迟器的作为波长和温度中的至少一个的函数的液晶双折射率。
6.根据权利要求1所述的方法,其中计算发生在图形处理单元上。
7.根据权利要求1所述的方法,其中为所述干涉仪的至少两个空间区域构建不同的矩阵。
8.根据权利要求1所述的方法,其中所述周期函数是余弦函数。
9.根据权利要求1所述的方法,其中所述周期函数中的每一个针对所述选定波长中的每一个是相移的,使得0相位大致对应于所述干涉仪在所述选定波长中的每一个处的零延迟点。
10.根据权利要求9所述的方法,其中所述相移取决于所述干涉仪的温度。
11.根据权利要求1所述的方法,还包括将所述重建矩阵的每一行乘以波长相关的函数,以减少由不对称测量的干涉图引起的伪像。
12.根据权利要求1所述的方法,其中所述重建矩阵的每一行乘以波长相关的窗口或切趾函数。
14.根据权利要求13所述的方法,其中应用正则化技术以防止所述对角化矩阵发散。
15.一种设备,包括:
干涉仪,所述干涉仪可操作用于通过一组参考延迟进行转换;
光学传感器,所述光学传感器可操作用于当所述干涉仪通过所述一组参考延迟进行转换时同时测量对应于所述光学传感器的不同位置的干涉图;以及
控制器,所述控制器耦合到所述干涉仪和所述光学传感器,所述控制器可操作用于:
确定对应于高光谱数据立方体的光谱切片的一组波长;
形成包括周期函数的行的重建矩阵,其中每一行对应于所述一组波长中的选定波长,并且每一列对应于所述参考延迟中的选定延迟,其中所述周期函数各自具有对应行的所述选定波长作为参数,并以对应列的所述选定延迟进行采样;
获得包括所述同时测量的干涉图中的一个或多个的阵列的干涉图数据立方体,其中所述干涉图数据立方体的每一行对应于所述选定延迟中的一个,并且所述干涉图数据立方体的每一列对应于与所述同时测量的干涉图不同的干涉图;并且
形成针对所述干涉图中的每一个的一组矩阵向量积,所述矩阵向量积中的每一个包括所述重建矩阵与所述干涉图数据立方体的列的矩阵乘法,所述一组矩阵向量积形成所述高光谱数据立方体。
16.根据权利要求15所述的设备,其中所述周期函数中的每一个针对所述选定波长中的每一个是相移的,使得0相位大致对应于所述干涉仪在所述选定波长中的每一个处的零延迟点。
17.根据权利要求15所述的设备,其中所述控制器进一步被配置成将所述重建矩阵的每一行乘以波长相关函数,以减少由不对称测量的干涉图引起的伪像。
18.根据权利要求15所述的设备,其中所述重建矩阵的每一行乘以波长相关窗口或切趾函数。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/858,338 US10663346B2 (en) | 2017-12-29 | 2017-12-29 | Method and apparatus for transforming uniformly or non-uniformly sampled interferograms to produce spectral data |
US15/858338 | 2017-12-29 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109990899A CN109990899A (zh) | 2019-07-09 |
CN109990899B true CN109990899B (zh) | 2022-01-14 |
Family
ID=64665224
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811476928.1A Active CN109990899B (zh) | 2017-12-29 | 2018-12-04 | 转换均匀或非均匀采样干涉图产生光谱数据的方法和设备 |
Country Status (5)
Country | Link |
---|---|
US (1) | US10663346B2 (zh) |
EP (1) | EP3505886B1 (zh) |
JP (1) | JP7082564B2 (zh) |
KR (1) | KR102439398B1 (zh) |
CN (1) | CN109990899B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11095390B2 (en) * | 2018-08-23 | 2021-08-17 | International Business Machines Corporation | Polarization-insensitive optical link |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1846114A (zh) * | 2003-07-18 | 2006-10-11 | 凯米映像公司 | 多波长成像光谱仪的方法和装置 |
CN1900741A (zh) * | 2005-11-18 | 2007-01-24 | 北京航空航天大学 | 高光谱全偏振三维成像集成探测系统 |
CN105571715A (zh) * | 2014-10-29 | 2016-05-11 | 帕洛阿尔托研究中心公司 | 液晶傅立叶变换成像光谱仪 |
Family Cites Families (85)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4342516A (en) | 1980-03-07 | 1982-08-03 | The Perkin-Elmer Corporation | Filter indexing for spectrophotometer system |
US4461543A (en) | 1982-03-26 | 1984-07-24 | Sperry Corporation | Electro optic switch |
JPS62103529A (ja) | 1985-10-31 | 1987-05-14 | Toshiba Corp | 複合管 |
US4848877A (en) | 1987-09-29 | 1989-07-18 | Cambridge Research And Instrumentation, Inc. | Liquid crystal devices and systems using such devices |
US4905169A (en) | 1988-06-02 | 1990-02-27 | The United States Of America As Represented By The United States Department Of Energy | Method and apparatus for simultaneously measuring a plurality of spectral wavelengths present in electromagnetic radiation |
US5126869A (en) | 1990-12-03 | 1992-06-30 | Raytheon Company | Two-dimensional, phased-array optical beam steerer |
US5784162A (en) | 1993-08-18 | 1998-07-21 | Applied Spectral Imaging Ltd. | Spectral bio-imaging methods for biological research, medical diagnostics and therapy |
DE69226998T2 (de) | 1991-07-19 | 1999-04-15 | Sharp Kk | Optisches Modulationselement und Vorrichtungen mit einem solchen Element |
US5347382A (en) | 1992-04-23 | 1994-09-13 | Rumbaugh Scott H | Liquid crystal cell retarder with driving beyond retardance value and two cells for high speed |
JP2925064B2 (ja) | 1993-12-02 | 1999-07-26 | 矢崎総業株式会社 | 可変波長フィルタ |
JPH07212639A (ja) | 1994-01-25 | 1995-08-11 | Sony Corp | テレビジョン・カメラ用電子シャッター装置 |
US5751384A (en) | 1995-05-23 | 1998-05-12 | The Board Of Regents Of The University Of Colorado | Color polarizers for polarizing an additive color spectrum along a first axis and it's compliment along a second axis |
US5856842A (en) | 1996-08-26 | 1999-01-05 | Kaiser Optical Systems Corporation | Apparatus facilitating eye-contact video communications |
US6100944A (en) * | 1997-10-10 | 2000-08-08 | Boulder Nonlinear Systems, Inc. | Polarizing interferometer using multiorder and zero order birefringence switches |
JP3732326B2 (ja) | 1998-02-02 | 2006-01-05 | 株式会社沖データ | 光プリントヘッド |
JP2002519716A (ja) | 1998-04-08 | 2002-07-02 | コーニング・アプライド・テクノロジーズ | 高速度電気光変換器 |
US6421131B1 (en) | 1999-07-02 | 2002-07-16 | Cambridge Research & Instrumentation Inc. | Birefringent interferometer |
GB2354833A (en) | 1999-09-28 | 2001-04-04 | Sharp Lab Of Europ Ltd | Liquid crystal device |
US7116370B1 (en) | 2000-03-31 | 2006-10-03 | Sharp Laboratories Of Ameria, Inc. | Image processing system optical shifting mechanism |
US6373614B1 (en) | 2000-08-31 | 2002-04-16 | Cambridge Research Instrumentation Inc. | High performance polarization controller and polarization sensor |
US7973936B2 (en) | 2001-01-30 | 2011-07-05 | Board Of Trustees Of Michigan State University | Control system and apparatus for use with ultra-fast laser |
US6576886B1 (en) | 2001-02-20 | 2003-06-10 | General Photonics Corporation | Dynamic control of polarization of an optical signal |
JP4382317B2 (ja) | 2001-12-06 | 2009-12-09 | シチズンホールディングス株式会社 | 液晶可変波長フィルタ装置 |
US6744509B2 (en) | 2002-08-20 | 2004-06-01 | Meadowlark Optics, Inc. | Retardance sweep polarimeter and method |
US7067795B1 (en) | 2002-10-03 | 2006-06-27 | General Photonics Corporation | Methods and systems for dynamic control of polarization of an optical signal |
EP1411333B1 (en) | 2002-10-15 | 2006-03-08 | Centre National De La Recherche Scientifique (Cnrs) | A liquid crystal based polarimetric system, a process for its calibration, and a polarimetric measurement process |
US7897406B2 (en) | 2002-12-20 | 2011-03-01 | Fiso Technologies Inc. | Method and sensor for detecting a chemical substance using an optically anisotropic material |
US7782377B2 (en) | 2003-02-26 | 2010-08-24 | Canon Kabushiki Kaisha | Image sensing apparatus, control method therefor, storage medium, and program to create correction data |
US7321791B2 (en) * | 2003-09-23 | 2008-01-22 | Cambridge Research And Instrumentation, Inc. | Spectral imaging of deep tissue |
JP4285350B2 (ja) | 2004-02-26 | 2009-06-24 | セイコーエプソン株式会社 | 視角制御素子およびその製造方法、液晶表示装置、電子機器 |
JP2005348140A (ja) | 2004-06-03 | 2005-12-15 | Matsushita Electric Ind Co Ltd | カメラ装置 |
JP4400448B2 (ja) | 2004-12-22 | 2010-01-20 | コニカミノルタセンシング株式会社 | 分光輝度計の校正方法、及び校正システムの動作プログラム |
US6992809B1 (en) | 2005-02-02 | 2006-01-31 | Chemimage Corporation | Multi-conjugate liquid crystal tunable filter |
US20080158550A1 (en) | 2005-03-29 | 2008-07-03 | Yoel Arieli | Spectral Imaging Camera and Applications |
US7420663B2 (en) | 2005-05-24 | 2008-09-02 | Bwt Property Inc. | Spectroscopic sensor on mobile phone |
JP4684771B2 (ja) | 2005-06-30 | 2011-05-18 | Hoya株式会社 | 像振れ補正装置 |
JP4205704B2 (ja) | 2005-08-02 | 2009-01-07 | 国立大学法人 北海道大学 | 撮像偏光計測方法 |
US7639363B2 (en) | 2005-09-07 | 2009-12-29 | California Institute Of Technology | Electro-optic imaging Fourier transform spectrometer |
US7630022B1 (en) | 2005-09-20 | 2009-12-08 | Thomas Baur | Advanced liquid crystal tunable filters and method |
WO2007038410A2 (en) | 2005-09-27 | 2007-04-05 | Chemimage Corporation | Liquid crystal filter with tunable rejection band |
US7525712B2 (en) | 2005-09-29 | 2009-04-28 | Teledyne Scientific & Imaging, Llc | Broad spectral range polarization rotator |
DE102005048240A1 (de) | 2005-10-07 | 2007-04-19 | Stefan Steib | Verfahren zur spektralen, integrierten Kalibrierung eines Bildsensors mittels monochromatischer Lichtquellen |
US8081311B2 (en) | 2005-11-04 | 2011-12-20 | General Hospital Corporation | System for multispectral imaging |
US8644911B1 (en) | 2006-06-30 | 2014-02-04 | Hypermed Imaging, Inc. | OxyVu-1 hyperspectral tissue oxygenation (HTO) measurement system |
US8447087B2 (en) | 2006-09-12 | 2013-05-21 | Carestream Health, Inc. | Apparatus and method for caries detection |
WO2008068753A2 (en) | 2006-12-04 | 2008-06-12 | Ben-Gurion University Of The Negev - Research And Development Authority | Polarization independent birefringent tunable filters |
GB0717967D0 (en) | 2007-09-14 | 2007-10-24 | Cascade Technologies Ltd | Polarimetric hyperspectral imager |
WO2009122317A2 (en) | 2008-04-01 | 2009-10-08 | Koninklijke Philips Electronics N.V. | Spectral detector calibration |
WO2010019515A2 (en) | 2008-08-10 | 2010-02-18 | Board Of Regents, The University Of Texas System | Digital light processing hyperspectral imaging apparatus |
FR2937732B1 (fr) | 2008-10-29 | 2010-12-31 | Horiba Jobin Yvon Sas | Dispositif et procede de mesures polarimetriques spectroscopiques dans le domaine visible et proche infrarouge |
BRPI1006077B8 (pt) | 2009-01-10 | 2023-02-28 | Goldfinch Solutions Llc | Sistema para analisar uma ou mais propriedades de carne |
US8525970B2 (en) | 2009-05-21 | 2013-09-03 | The Hong Kong University Of Science And Technology | Photo-aligned liquid-crystal micropolarimeter array and its manufacturing method |
US8982313B2 (en) | 2009-07-31 | 2015-03-17 | North Carolina State University | Beam steering devices including stacked liquid crystal polarization gratings and related methods of operation |
US7999933B2 (en) | 2009-08-14 | 2011-08-16 | Princeton Instruments | Method for calibrating imaging spectrographs |
US8400574B2 (en) | 2010-04-16 | 2013-03-19 | Chemimage Corporation | Short wave infrared multi-conjugate liquid crystal tunable filter |
WO2011132455A1 (ja) | 2010-04-19 | 2011-10-27 | シャープ株式会社 | 表示装置 |
US20110299089A1 (en) * | 2010-06-04 | 2011-12-08 | University Of Washington | Compact liquid crystal based fourier transform spectrometer system |
US20120013922A1 (en) | 2010-07-19 | 2012-01-19 | Lotz Michael G | Print Data Format Modification Mechanism |
US9074993B2 (en) * | 2010-09-03 | 2015-07-07 | The Arizona Board Of Regents On Behalf Of The University Of Arizona | White light achromatic grating imaging polarimeter |
US9442015B2 (en) * | 2010-09-03 | 2016-09-13 | The Arizona Board Of Regents On Behalf Of The University Of Arizona | Snapshot spatial heterodyne imaging polarimetry |
US9046422B2 (en) * | 2011-04-20 | 2015-06-02 | The Arizona Board Of Regents On Behalf Of The University Of Arizona | Ultra-compact snapshot imaging fourier transform spectrometer |
US8422119B1 (en) | 2010-09-20 | 2013-04-16 | Disco Corporation | Compensation of beam walkoff in nonlinear crystal using cylindrical lens |
US8736777B2 (en) | 2011-01-19 | 2014-05-27 | Chemimage Technologies Llc | VIS-SNIR multi-conjugate liquid crystal tunable filter |
EP2551635A1 (de) | 2011-07-26 | 2013-01-30 | Hexagon Technology Center GmbH | Optisches Messsystem mit Filtereinheit zur Extraktion elektromagnetischer Strahlung |
JP6073883B2 (ja) | 2011-07-28 | 2017-02-01 | マサチューセッツ インスティテュート オブ テクノロジー | 内部空間の3次元画像化のためのカメラのコンフィグレーション |
JP2013096883A (ja) | 2011-11-02 | 2013-05-20 | Seiko Epson Corp | 分光測定装置 |
US20140362331A1 (en) | 2013-03-15 | 2014-12-11 | Chemlmage Corporation | Short-Wavelength Infrared (SWIR) Multi-Conjugate Liquid Crystal Tunable Filter |
US20140354868A1 (en) | 2013-06-04 | 2014-12-04 | Corning Incorporated | Portable hyperspectral imager |
US9513166B2 (en) | 2013-07-16 | 2016-12-06 | Utah State University Research Foundation | Split field spectral imager |
WO2015015722A1 (ja) | 2013-07-29 | 2015-02-05 | パナソニックIpマネジメント株式会社 | 光学フィルタおよびそれを用いた偏光撮像装置 |
US9631973B2 (en) | 2013-12-13 | 2017-04-25 | Raytheon Company | Multifunction imager |
CN106462017B (zh) | 2014-04-02 | 2019-08-09 | 依视路国际公司 | 消色差相位调制器和光学器件 |
US11490037B2 (en) | 2014-10-29 | 2022-11-01 | Palo Alto Research Center Incorporated | Liquid crystal fourier transform imaging spectrometer |
US10469771B2 (en) | 2014-10-29 | 2019-11-05 | Palo Alto Research Center Incorporated | Liquid crystal fourier transform imaging spectrometer |
US10151634B2 (en) | 2014-11-16 | 2018-12-11 | Ibrahim Abdulhalim | Multi-spectral polarimetric variable optical device and imager |
US10302494B2 (en) * | 2014-12-18 | 2019-05-28 | Palo Alto Research Center Incorporated | Obtaining spectral information from a moving object |
EP3243103B1 (en) | 2015-01-09 | 2020-04-29 | Essilor International | Liquid-crystal achromatic phase modulator |
US10228517B2 (en) | 2015-03-03 | 2019-03-12 | Nistica, Inc. | Optical arrangement for managing diversity and isolation between ports in a wavelength selective switch |
JP6457122B2 (ja) * | 2016-02-04 | 2019-01-23 | 日本分光株式会社 | フーリエ変換型分光装置を用いたスペクトル測定方法 |
JP2017191071A (ja) * | 2016-04-15 | 2017-10-19 | キヤノン株式会社 | 分光データ処理装置、撮像装置、分光データ処理方法および分光データ処理プログラム |
US10425598B2 (en) | 2016-06-20 | 2019-09-24 | Massachusetts Institute Of Technology | Methods and systems for time-encoded multiplexed imaging |
US10768497B2 (en) | 2016-10-03 | 2020-09-08 | Xerox Corporation | Hyperspectral imaging system |
CN110178062A (zh) | 2017-01-06 | 2019-08-27 | 尼斯迪卡有限公司 | 用于抑制波长选择开关中的带外串扰的光装置 |
US10884278B2 (en) | 2017-10-24 | 2021-01-05 | Palo Alto Research Center Incorporated | Liquid crystal temperature control by resistive heating |
US10175116B1 (en) * | 2017-12-29 | 2019-01-08 | Palo Alto Research Center Incorporated | Color filter used with liquid-crystal polarization interferometer |
-
2017
- 2017-12-29 US US15/858,338 patent/US10663346B2/en active Active
-
2018
- 2018-12-04 CN CN201811476928.1A patent/CN109990899B/zh active Active
- 2018-12-06 JP JP2018228850A patent/JP7082564B2/ja active Active
- 2018-12-12 EP EP18212124.4A patent/EP3505886B1/en active Active
- 2018-12-13 KR KR1020180160553A patent/KR102439398B1/ko active IP Right Grant
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1846114A (zh) * | 2003-07-18 | 2006-10-11 | 凯米映像公司 | 多波长成像光谱仪的方法和装置 |
CN1900741A (zh) * | 2005-11-18 | 2007-01-24 | 北京航空航天大学 | 高光谱全偏振三维成像集成探测系统 |
CN105571715A (zh) * | 2014-10-29 | 2016-05-11 | 帕洛阿尔托研究中心公司 | 液晶傅立叶变换成像光谱仪 |
Non-Patent Citations (1)
Title |
---|
Compact snapshot birefringent imaging Fourier transform spectrometer;Michael W. Kudenov 等;《Proceedings of SPIE》;20100813;第7812卷;第781206页 * |
Also Published As
Publication number | Publication date |
---|---|
EP3505886A2 (en) | 2019-07-03 |
JP2019120681A (ja) | 2019-07-22 |
JP7082564B2 (ja) | 2022-06-08 |
EP3505886A3 (en) | 2019-07-31 |
US20190204153A1 (en) | 2019-07-04 |
KR102439398B1 (ko) | 2022-09-02 |
EP3505886B1 (en) | 2020-09-23 |
KR20190082090A (ko) | 2019-07-09 |
CN109990899A (zh) | 2019-07-09 |
US10663346B2 (en) | 2020-05-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9880053B2 (en) | Image pickup apparatus, spectroscopic system, and spectroscopic method | |
KR102430145B1 (ko) | 상이한 두께의 액정 셀들을 이용한 액정 가변 지연자 | |
US9097585B2 (en) | Spectrographic polarimeter and method of recording state of polarity | |
CN111443045B (zh) | 光谱成像分析系统及光谱成像分析方法 | |
Han et al. | Iterative reconstruction for snapshot intensity-modulated linear imaging spectropolarimetry without Fourier transform and phase calibration | |
CN109990899B (zh) | 转换均匀或非均匀采样干涉图产生光谱数据的方法和设备 | |
WO2017201548A1 (en) | Miniaturized and integrated photonic spectrometers | |
del Hoyo et al. | High precision calibration method for a four-axis Mueller matrix polarimeter | |
US20220113674A1 (en) | Differential holography | |
US11530953B2 (en) | Snapshot Mueller matrix polarimeter | |
KR102483017B1 (ko) | 불균일한 리타던스 간격에서 액정 가변 지연기를 통한 경로 지연 측정 | |
Piazzo | Least squares image estimation for large data in the presence of noise and irregular sampling | |
Samuel et al. | Invariants in Polarimetric Interferometry: A Non-Abelian Gauge Theory | |
CN112781728B (zh) | 精确求解结合压缩感知的全偏振高光谱成像方法 | |
Logofătu et al. | Snapshot interferometric multispectral imaging using deconvolution and colorimetric fit | |
Hagen et al. | Generating high-performance polarization measurements with low-performance polarizers: demonstration with a microgrid polarization camera | |
Han et al. | Deep image prior plus sparsity prior: toward single-shot full-Stokes spectropolarimetric imaging with a multiple-order retarder | |
Peinado et al. | Polarization imaging with enhanced spatial resolution | |
Miles et al. | Nonscanning computed tomography imaging spectropolarimeter (NS-CTISP): design and calibration | |
CN110118601B (zh) | 一种基于lctf的光谱成像装置及成像方法 | |
Murali | Analysis of an interferometric Stokes imaging polarimeter | |
Weinschenk et al. | Restoration of broadband imagery steered with a liquid-crystal optical phased array | |
Homolya et al. | Numerical determination of the density profile of an inhomogeneous membrane: solution of the inverse vibration problem | |
CN115096443A (zh) | 一种基于穆勒矩阵降维的椭偏参数计算方法 | |
Greivenkamp | Sub-nyquist interferometry |
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 |