CN105486408B - 基于彩色CCD的三波长Stokes矢量偏振测量方法 - Google Patents
基于彩色CCD的三波长Stokes矢量偏振测量方法 Download PDFInfo
- Publication number
- CN105486408B CN105486408B CN201510837226.1A CN201510837226A CN105486408B CN 105486408 B CN105486408 B CN 105486408B CN 201510837226 A CN201510837226 A CN 201510837226A CN 105486408 B CN105486408 B CN 105486408B
- Authority
- CN
- China
- Prior art keywords
- mtd
- mrow
- msub
- mtr
- 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.)
- Expired - Fee Related
Links
- 239000013598 vector Substances 0.000 title claims abstract description 37
- 230000010287 polarization Effects 0.000 title claims abstract description 34
- 238000000691 measurement method Methods 0.000 title claims description 10
- 238000000034 method Methods 0.000 claims abstract description 16
- 238000006243 chemical reaction Methods 0.000 claims abstract description 10
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 238000001228 spectrum Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 64
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 5
- 230000015556 catabolic process Effects 0.000 claims description 4
- 238000005259 measurement Methods 0.000 abstract description 25
- 238000000711 polarimetry Methods 0.000 abstract description 4
- 238000001514 detection method Methods 0.000 abstract description 3
- 238000003384 imaging method Methods 0.000 abstract description 3
- 230000003595 spectral effect Effects 0.000 abstract 1
- 230000010354 integration Effects 0.000 description 4
- 241001670157 Gymnura Species 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009466 transformation Effects 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
- G01J4/00—Measuring polarisation of light
- G01J4/04—Polarimeters using electric detection means
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Color Television Image Signal Generators (AREA)
- Polarising Elements (AREA)
Abstract
本发明公开了基于彩色CCD的三波长Stokes矢量偏振测量系统及方法,用于同步测量三个波长下Stokes矢量的值,该方法在Stokes偏振测量装置中采用红、绿、蓝对应的三个波长的激光器作为光源,并采用彩色CCD作为光强探测器件;考虑彩色CCD红绿蓝(RGB)三个通道量子转换效率的光谱特性,补偿三个通道量子转换效率光谱重叠引起的串扰,从而通过RGB三个通道下测得的光强值计算三个波长下的真实光强值;同时考虑偏振分析器(PSA)在三个波长下的特征偏振态,并基于此由三个波长下的真实光强值最终计算出三个波长下的Stokes矢量。该方法可仅通过四次光强测量即获得三个波长下的Stokes矢量值,从而大大增加了获取的偏振信息量。该方法对于偏振成像和偏振测量探测具有重要意义。
Description
技术领域
本发明涉及偏振测量领域,特别是涉及一种基于彩色CCD的三波长Stokes矢量偏振测量方法。
背景技术
偏振信息作为光波的基本物理信息之一,可以提供其它光波信息所不能提供的被测物信息,因此基于偏振信息的测量在许多领域有着十分广泛的应用。Stokes矢量描述了光波的偏振态,包含了最基本的偏振信息。Stokes矢量的测量也由此成为偏振测量领域的主要方向之一,因此提高Stokes矢量的测量精度对于提高偏振测量技术的水平具有重要意义。测量数据的方差是影响测量精度的关键因素。通常的Stokes矢量测量是通过光强的测量来实现的:首先测量偏振分析器(PSA)四个不同状态下的透射光强信息,其中PSA四个不同状态下对应的四次光强测量的积分时间相同;然后根据PSA测量矩阵和光强的测量值计算Stokes矢量。之前的均分光强积分时间的方案并没有考虑光强积分时间对于Stokes矢量测量方差的影响,因此在某些情况下,不能实现最小化的Stokes矢量测量方差和最优化的测量精度。
发明内容
基于上述现有技术存在的问题,本发明提出了一种基于彩色CCD的三波长Stokes矢量偏振测量方法,利用彩色CCD相机RGB模式三通道下光强串扰的直接影响,并通过测量RGB三通道中共十二个维度下的光强信息,得到被测物体在不同波长光照射下分别的偏振信息,进而得到物体在多波长光照射下各自的Stokes矢量的测量值。
本发明提出了一种基于光强测量积分时间优化的Stokes矢量测量方法,用于同步测量三个波长下Stokes矢量的值,该方法包括以下步骤:
步骤一、由偏振片和四分之一波片构成的偏振分析器,其Mueller矩阵M表示为:
M=MPMR
式中,MP为与水平相夹角度为θP的偏振片的Mueller矩阵;MR为相位延迟量为α、其快轴方向与水平相夹角度为θR的波片的Mueller矩阵;取其第一行行矩阵并设其为T1,当偏振片与四分之一波片处于四种不同状态时,获得四个不同的行矩阵求出三个波长下相应的测量矩阵W:
利用彩色CCD相机采集偏振分析器四个不同状态下的下的四幅RGB三维彩色图像;
步骤二、将采集得到的彩色图像的光强信息进行图像分解和处理,获得彩色光强信息在R、G、B三个独立色彩空间下各自的单色光光强;
步骤三、利用已知的彩色CCD相机RGB三通道下的量子转换效率光谱曲线对单色光光强进行修正,光强修正系数矩阵K为:
其转置矩阵KT为:
进而得到CCD相机入射光在三个不同波长下的真实光强;实际入射光对应的光强矩阵I:
I=KTIs
Is为彩色CCD相机在PSA处于四种不同偏振状态下直接获得的光强矩阵;
步骤四、根据修正后单色光光强矩阵I以及每个波长下的测量矩阵W,分别测得三波长光照射下,物体对应于每一波长下的Stokes参量,具体包括以下计算:
求解矩阵I,并分离其行矩阵即可将获得在每一波长下的光强矩阵
i1,i2,i3,i4为偏振分析器四种不同偏振状态下的光强;
每一波长下的Stokes矢量Sλ表示为:
Sλ=Wλ -1Iλ
其中λ=450nm,532nm,633nm,Wλ -1其为每一波长下的测量矩阵W的转置矩阵。
其中所述步骤二的对彩色图像的光强信息进行图像分解和处理的方法是:将彩色图像的光强信息分解为RGB空间下三色单色光光强信息,并将其以同一维度图像的单色光强按R、G、B三色顺序排成一列的规律制成一3×4的光强矩阵。
其中所述步骤三的对单色光光强信息进行修正的方法是:将已知CCD相机的RGB串扰曲线在入射光对应波长下的RGB串扰系数编制成一3×3矩阵,将该矩阵作为系数矩阵与光强矩阵相乘,所获得的新的3×4矩阵即为原光强矩阵的修正矩阵。
本发明的有益效果及优点在于:
1、本发明方法利用彩色CCD三通道同步光强测量的特点,将Stokes矢量测量方法推广至多光谱领域,并且不增加光强测量次数。测量多波长条件下下的Stokes矢量对彩色偏振成像研究具有重要意义。
2、本发明方法只需使用CCD相机采集4次图像即可得到三个波长下共12个维度的物体偏振信息,大大增加了获取的偏振信息量。
附图说明
图1为基于彩色CCD的三波长Stokes矢量测量系统结构示意图;
图2为本实验室采用的彩色CCD相机波长与RGB三个通道对于不同波长的转换效率曲线。
图3为实验过程中采集的彩色图像在RGB空间分解得到的三色灰度(光强)图像;
图4为图3中经消串扰处理后得到的RGB各个通道归一化真实光强图像;
图5为最终处理得到的三波长下Stokes参量各分量图像信息。
附图标记:1、450nm、532nm、633nm激光光源组(三波长激光光源组),2、反射镜,3、半透半反镜组,4、扩束装置组,5、反射板(表面贴有被测物体),6、四分之一波片,7、偏振片,8、彩色CCD相机(Stingray F-033C型)。
具体实施方式
如图1所示,450nm、532nm、633nm激光光源组1发出的光分别经过反射镜2和半透半反镜组3入射到扩束装置4上,经反射板5反射后形成固定的待测Stokes矢量,再经过由四分之一波片6和偏振片7组成的偏振态分析器(PSA)后进入彩色CCD相机8。通过调节四分之一波片6和偏振片7的角度可以实现测量矩阵W中的各个偏振状态。本系统所选用的光强探测器件是Stingray F-033C型彩色CCD相机。
基于彩色CCD的三波长Stokes矢量测量方法具体步骤如下:
对于一个任意波片,设其相位延迟量为α,其快轴方向与水平相夹角度为θR。则其Mueller矩阵为:
而对于一个与水平相夹角度为θP的线偏振片而言,其Mueller矩阵为:
本发明中的PSA是由一个线偏振片和一个四分之一波片构成的。此时其Mueller矩阵M可以表示为:
M=MPMR (11)
取其第一行行矩阵并设其为T1。当偏振片与波片处于四种不同状态时,可获得四个不同的行矩阵则此时W矩阵为:
需要说明的是,波片的相位延迟量为α是随波长变化的。通过获得实验中所用波片在不同波长的下的相位延迟量,并将偏振片和波片的角度一同代入公式(12)中的W矩阵即可得到矩阵W在不同波长下的具体参数。
在本实施例中,PSA四个状态下偏振片和波片在的角度分别为(θP,θR)=(0°,0°);(90°,0°);(45°,45°);(45°,0°)。以PSA偏振器件组在450nm、532nm、633nm波长下的测量矩阵的逆矩阵为例,该三个波长对应的测量矩阵分别W为:
本实施案例所用的Stingray F-033C型彩色CCD相机RGB三个通道在不同波长下的转换效率如图2所示。根据图2,R通道在450nm、532nm、633nm波长下的转换效率分别为0.0061,0.0142,0.2678;G通道在450nm、532nm、633nm波长下的转换效率分别为0.0497,0.3010,0.0081;B通道在450nm、532nm、633nm波长下的转换效率分别为0.2330,0.0477,0。则其光强修正系数矩阵K为:
其转置矩阵KT为:
首先记录在PSA(θP,θR)=(0°,0°);(90°,0°);(45°,45°);(45°,0°)四个状态下对应的四副彩色CCD图像。根据记录得到的彩色图像,可以获得RGB三个通道对应的灰度图像,如图3所示。
随后,根据光强修正系数矩阵K,利用公式(6)便可消除串扰,得到经消串扰处理后RGB各个通道真实光强图像,即450nm、532nm、633nm三个波长的真实的光强二维分布,如图4所示。
根据图4、公式(7)以及公式(8),便可计算可以得到450nm、532nm、633nm三个波长下Stokes矢量的值,如图5所示。其中,中心矩形区域三波长下的Stokes参量分别为:
S450=[0.1682,0.0747,0.0943,-0.0602]
S532=[0.2120,0.1501,0.0311,-0.0943]
S633=[0.1582,0.0770,0.0445,-0.0014]
可见,利用彩色CCD,并基于上述方法,可以仅通过四次光强测量(图像采集),即计算出三个波长下的Stokes矢量。
图3中每一列图像均为PSA处于某一特定状态下采集的RGB三色灰度图像。
图4中每行图像由上而下依次表示PSA处于实施例所述四种偏振状态下450nm、532nm、633nm波长的灰度图像。
图5中每行图像由上而下依次表示450nm、532nm、633nm波长下被测物体Stokes矢量S0、S1、S2、S3参量的测量值。
本发明的理论依据:
在测量单色光照射下的条件下,考虑以下基于光强测量实现Stokes矢量测量的方法:
I=WS (1)
其中I=[i1,i2,i3,i4]T一般为四阶光强列向量,其每一个分量表示经过PSA过滤后被CCD记录的光强;W表示的是4×4的测量矩阵,每一行都可以认为是一组Stokes矢量,表征着PSA的状态;S=[S0,S1,S2,S3]T表示待测的Stokes矢量。
在不考虑彩色CCD各通道串扰的情况下,当物体同时被多个波长的光束照射时,利用彩色CCD相机可以得到目标物体的彩色图像,但此时彩色图像并不只具备灰度这一概念。因此本发明将CCD相机拍摄的彩色图像在RGB空间进行分解,从而可以从每一幅彩色图像中获得R、G、B三个独立色彩空间的三幅灰度图像,即光强信息。利用分解得到的各个波长下的光强信息,可以依照上述测量单色光Stokes参量方法测量RGB三色空间下对应的Stokes参量。
在考虑分别对应于R、G、B三个通道的633nm、532nm和450nm的三色激光同时入射被测物体的条件下,设其对应的四阶光强向量分别为I1,I2,I3;其4x4测量矩阵分别为W1,W2,W3;被测物体在三个波长下的Stokes参量分别为S1,S2,S3,则有:
实际上,由于彩色CCD相机的工作原理是利用CCD芯片前的Bayer滤光片将物体的灰度信息模拟成为彩色图像。在成像过程中,彩色CCD相机的RGB三通道对应的滤光光谱的部分重叠将导致RGB三个通道间将存在光强的串扰。举例而言,当一束以633nm波长入射的红色激光进入彩色CCD相机时;彩色CCD相机会将入射的一部分红光误判为蓝光和绿光,从而导致RGB三个通道所测的光强信息与真实值产生偏差。
为了排除该串扰对Stokes参量测量的影响,本发明利用彩色CCD波长与RGB三色灰度转换效率曲线,求得了在特定入射波长组合下的光强修正系数矩阵K。该矩阵在入射光波长分别为450nm、532nm、633nm的条件下可以表示为:
在该矩阵中参数的定义为:波长为λ的入射光在R、G、B空间下各自的光强转化效率系数。
设彩色CCD相机在PSA处于四种不同偏振状态下直接获得的光强矩阵为:
矩阵第一行至第三行行矩阵分别为在R、G、B空间下测得的四维光强矩阵。
而实际入射光对应的光强矩阵也可表示为:
矩阵第一行至第三行行矩阵分别为在450nm、532nm、633nm波长下的四维光强矩阵。
此时I、Is、K之间应满足关系:
I=KTIs (6)
求解I矩阵并分离其行矩阵即可将获得在每一波长下的光强矩阵:
其中λ=450nm,532nm,633nm。此时每一波长下的Stokes矢量Sλ可表示为:
Sλ=Wλ -1Iλ (8)。
Claims (3)
1.一种基于彩色CCD的三波长Stokes矢量偏振测量方法,用于同步测量三个波长下Stokes矢量的值,其特征在于,该方法包括以下步骤:
步骤一、由偏振片和四分之一波片构成的偏振分析器,其Mueller矩阵M表示为:
M=MPMR
<mrow>
<msub>
<mi>M</mi>
<mi>P</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mi>cos</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>P</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>sin</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>P</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>cos</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>P</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>P</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>sin</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>P</mi>
</msub>
<mi>cos</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>P</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>sin</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>P</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>sin</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>P</mi>
</msub>
<mi>cos</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>P</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>P</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mi>R</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>R</mi>
</msub>
<mo>+</mo>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>R</mi>
</msub>
<mi>cos</mi>
<mi>&alpha;</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>sin</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>R</mi>
</msub>
<mi>cos</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>cos</mi>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>sin</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>R</mi>
</msub>
<mi>sin</mi>
<mi>&alpha;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>sin</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>R</mi>
</msub>
<mi>cos</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>cos</mi>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>R</mi>
</msub>
<mo>+</mo>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>R</mi>
</msub>
<mi>cos</mi>
<mi>&alpha;</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>cos</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>R</mi>
</msub>
<mi>sin</mi>
<mi>&alpha;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>sin</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>R</mi>
</msub>
<mi>sin</mi>
<mi>&alpha;</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>cos</mi>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>R</mi>
</msub>
<mi>sin</mi>
<mi>&alpha;</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>cos</mi>
<mi>&alpha;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中,MP为与水平相夹角度为θP的偏振片的Mueller矩阵;MR为相位延迟量为α、其快轴方向与水平相夹角度为θR的波片的Mueller矩阵;取其第一行行矩阵并设其为T1,当偏振片与四分之一波片处于四种不同状态时,获得四个不同的行矩阵求出三个波长下相应的测量矩阵W:
<mrow>
<mi>W</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>T</mi>
<mn>1</mn>
<mn>1</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>T</mi>
<mn>2</mn>
<mn>1</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>T</mi>
<mn>3</mn>
<mn>1</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>T</mi>
<mn>4</mn>
<mn>1</mn>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
利用彩色CCD相机采集偏振分析器四个不同状态下的下的四幅RGB三维彩色图像;
步骤二、将采集得到的彩色图像的光强信息进行图像分解和处理,获得彩色光强信息在R、G、B三个独立色彩空间下各自的单色光光强;
步骤三、利用已知的彩色CCD相机RGB三通道下的量子转换效率光谱曲线对单色光光强进行修正,光强修正系数矩阵K为:
<mrow>
<mi>K</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0.0061</mn>
</mtd>
<mtd>
<mn>0.0142</mn>
</mtd>
<mtd>
<mn>0.2678</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0.0497</mn>
</mtd>
<mtd>
<mn>0.3010</mn>
</mtd>
<mtd>
<mn>0.0081</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0.2330</mn>
</mtd>
<mtd>
<mn>0.0477</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其转置矩阵KT为:
<mrow>
<msup>
<mi>K</mi>
<mi>T</mi>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0.0213</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.7049</mn>
</mrow>
</mtd>
<mtd>
<mn>4.4416</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.1041</mn>
</mrow>
</mtd>
<mtd>
<mn>3.4431</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.7317</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>3.7392</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.1665</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.0624</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
进而得到CCD相机入射光在三个不同波长下的真实光强;实际入射光对应的光强矩阵I:
I=KTIs
Is为彩色CCD相机在PSA处于四种不同偏振状态下直接获得的光强矩阵;
步骤四、根据修正后单色光光强矩阵I以及每个波长下的测量矩阵W,分别测得三波长光照射下,物体对应于每一波长下的Stokes参量,具体包括以下计算:
求解矩阵I,并分离其行矩阵即可将获得在每一波长下的光强矩阵
<mrow>
<msub>
<mi>I</mi>
<mi>&lambda;</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>i</mi>
<mn>1</mn>
<mi>&lambda;</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>i</mi>
<mn>2</mn>
<mi>&lambda;</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>i</mi>
<mn>3</mn>
<mi>&lambda;</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>i</mi>
<mn>4</mn>
<mi>&lambda;</mi>
</msubsup>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
i1,i2,i3,i4为偏振分析器四种不同偏振状态下的光强;
每一波长下的Stokes矢量Sλ表示为:
Sλ=Wλ -1Iλ
其中λ=450nm,532nm,633nm,Wλ -1其为每一波长下的测量矩阵W的转置矩阵。
2.如权利要求1所述的一种基于彩色CCD的三波长Stokes矢量偏振测量方法,其特征在于,其中所述步骤二的对彩色图像的光强信息进行图像分解和处理的方法是:将彩色图像的光强信息分解为RGB空间下三色单色光光强信息,并将其以同一维度图像的单色光强按R、G、B三色顺序排成一列的规律制成一3×4的光强矩阵。
3.如权利要求1所述的一种基于彩色CCD的三波长Stokes矢量偏振测量方法,其特征在于,其中所述步骤三的对单色光光强信息进行修正的方法是:将已知CCD相机的RGB串扰曲线在入射光对应波长下的RGB串扰系数编制成一3×3矩阵,将该矩阵作为系数矩阵与光强矩阵相乘,所获得的新的3×4矩阵即为原光强矩阵的修正矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510837226.1A CN105486408B (zh) | 2015-11-26 | 2015-11-26 | 基于彩色CCD的三波长Stokes矢量偏振测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510837226.1A CN105486408B (zh) | 2015-11-26 | 2015-11-26 | 基于彩色CCD的三波长Stokes矢量偏振测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105486408A CN105486408A (zh) | 2016-04-13 |
CN105486408B true CN105486408B (zh) | 2017-10-27 |
Family
ID=55673539
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510837226.1A Expired - Fee Related CN105486408B (zh) | 2015-11-26 | 2015-11-26 | 基于彩色CCD的三波长Stokes矢量偏振测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105486408B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106504291A (zh) * | 2016-10-26 | 2017-03-15 | 天津大学 | 基于场景Stokes矢量测量及其RGB图像优化的多物体区分方法 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106569402A (zh) * | 2016-11-15 | 2017-04-19 | 西北工业大学 | 一种多波长数字全息图的色彩分离与数字处理方法 |
CN107817050B (zh) * | 2017-09-07 | 2019-12-17 | 天津大学 | 基于彩色相机的多光谱彩色偏振成像对比度优化的方法 |
CN110458769B (zh) * | 2019-07-22 | 2023-04-07 | 天津大学 | 一种基于消除通道间串扰的彩色偏振图像复原方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0622624A1 (en) * | 1993-04-23 | 1994-11-02 | Research Development Corporation Of Japan | A method for observing film thickness and/or refractive index |
JP2000266666A (ja) * | 1999-03-19 | 2000-09-29 | Hitachi Koki Co Ltd | エリプソメータ |
WO2004097383A1 (en) * | 2003-04-29 | 2004-11-11 | Surfoptic Limited | Measuring a surface characteristic |
-
2015
- 2015-11-26 CN CN201510837226.1A patent/CN105486408B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0622624A1 (en) * | 1993-04-23 | 1994-11-02 | Research Development Corporation Of Japan | A method for observing film thickness and/or refractive index |
JP2000266666A (ja) * | 1999-03-19 | 2000-09-29 | Hitachi Koki Co Ltd | エリプソメータ |
WO2004097383A1 (en) * | 2003-04-29 | 2004-11-11 | Surfoptic Limited | Measuring a surface characteristic |
Non-Patent Citations (1)
Title |
---|
Optimal distribution of integration time for intensity measurements in Stokes polarimetry;Xiaobo Li et al;《OPTICS EXPRESS》;20151013;全文 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106504291A (zh) * | 2016-10-26 | 2017-03-15 | 天津大学 | 基于场景Stokes矢量测量及其RGB图像优化的多物体区分方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105486408A (zh) | 2016-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105486408B (zh) | 基于彩色CCD的三波长Stokes矢量偏振测量方法 | |
CN101124462B (zh) | 多光谱和超光谱成像系统 | |
KR102069935B1 (ko) | 디스플레이를 테스트하기 위한 색채측정기 시스템 | |
US9030580B2 (en) | Color filter modules for plenoptic XYZ imaging systems | |
US7768641B2 (en) | Spatial image modulation to improve performance of computed tomography imaging spectrometer | |
Shrestha et al. | Multispectral imaging using LED illumination and an RGB camera | |
CN107615022A (zh) | 具有用于创建空间编码图像的图像分散性的成像设备 | |
CN104181546B (zh) | 彩色三维扫描激光雷达颜色信息采集显示方法 | |
DE112018001633T5 (de) | Mehrkamera-Bilderfassung zum Messen der Beleuchtung | |
Bongiorno et al. | Spectral characterization of COTS RGB cameras using a linear variable edge filter | |
CN103884427A (zh) | 高光谱高空间分辨率的图像获取方法及装置 | |
CN106153193B (zh) | 一种利用多光谱相机双光源响应值获取光谱反射比的方法 | |
US11622084B2 (en) | Hyperspectral camera based on continuously variable film filter and coating method thereof | |
CN106153192A (zh) | 一种利用多光谱相机虚拟响应值获取光谱反射比的方法 | |
US7876434B2 (en) | Color camera computed tomography imaging spectrometer for improved spatial-spectral image accuracy | |
JP2018040976A (ja) | フィルタ、画像撮影装置および画像撮影システム | |
Picollo et al. | Application of hyper-spectral imaging technique for colorimetric analysis of paintings | |
Lee et al. | Building a two-way hyperspectral imaging system with liquid crystal tunable filters | |
Tisserand | Custom Bayer filter multispectral imaging: emerging integrated technology | |
CN110332991B (zh) | 一种光谱的重构方法、装置和电子设备 | |
Gilchrist et al. | Developing the IEEE P4001 standard for characterisation and calibration of hyperspectral imaging devices | |
Lin et al. | Spectral image reconstruction by a tunable led illumination | |
Salesin et al. | Diy hyperspectral imaging via polarization-induced spectral filters | |
CN111323123B (zh) | 一种基于渐变薄膜滤光片的高光谱相机及其镀膜方法 | |
Shin et al. | Dispersed Structured Light for Hyperspectral 3D Imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171027 Termination date: 20181126 |
|
CF01 | Termination of patent right due to non-payment of annual fee |