CN112417370B - 粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法 - Google Patents

粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法 Download PDF

Info

Publication number
CN112417370B
CN112417370B CN202011263093.9A CN202011263093A CN112417370B CN 112417370 B CN112417370 B CN 112417370B CN 202011263093 A CN202011263093 A CN 202011263093A CN 112417370 B CN112417370 B CN 112417370B
Authority
CN
China
Prior art keywords
matrix
polarization
mueller
stokes
rough surface
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
Application number
CN202011263093.9A
Other languages
English (en)
Other versions
CN112417370A (zh
Inventor
李宏哲
赵志敏
俞晓磊
顾慈勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202011263093.9A priority Critical patent/CN112417370B/zh
Publication of CN112417370A publication Critical patent/CN112417370A/zh
Application granted granted Critical
Publication of CN112417370B publication Critical patent/CN112417370B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/21Polarisation-affecting properties
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/21Polarisation-affecting properties
    • G01N2021/217Measuring depolarisation or comparing polarised and depolarised parts of light

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Chemical & Material Sciences (AREA)
  • Operations Research (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biochemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Computing Systems (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法,步骤:测试粗糙表面物质的穆勒矩阵;计算粗糙表面物质穆勒矩阵对应的相干矩阵的特征值并分析特征值;计算粗糙表面物质穆勒矩阵的退偏振系数;计算粗糙表面物质出射光斯托克斯矢量的平均偏振度;将穆勒矩阵和斯托克斯矢量同时分解,建立分解物理方程,进行穆勒琼斯矩阵的最优化估计,并对偏振噪声进行分析。

Description

粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法
技术领域
本发明属于偏振技术领域,特别涉及了粗糙表面物质的穆勒矩阵分解方法。
背景技术
近年来,偏振技术在目标识别,物质分类,遥感,雷达技术,医学检测等领域的应用更加广泛。粗糙表面物质的退偏振研究是偏振研究的重点之一,退偏振与出射光的空间或时间相干性的缺失有关。很多光学元件的退偏振数值比较小,而金属,薄膜物质,自然界的物质其退偏振数值比较大,对于光学系统具有重要的影响。
粗糙表面物质退偏振的影响因素复杂,如穆勒矩阵的内部结构,表面散射,光学元件的加工缺陷,仪器的光谱分辨率,穆勒琼斯矩阵的非相干叠加等。当物质的穆勒矩阵为退偏振矩阵时,琼斯矩阵和穆勒矩阵就不能完全对应,求解琼斯矩阵就需要进行最佳估计或者近似。穆勒矩阵分解方法对于物质退偏振研究和实验数据解释分析是一种有效的研究方法。
穆勒矩阵分解方法有三种典型的方法,包括乘积分解法,微分分解法,求和分解法。当物质的退偏振效应较小时,几种分解方法的结果基本一致。当粗糙表面物质的退偏振效应较大时,不同分解方法的结果会有较大的区别,原因是因为计算偏振效应的方法不同。在常见的分解方法中缺少一种通用的分解方法,能够同时对于斯托克斯矢量和穆勒矩阵进行分解,以进行穆勒琼斯矩阵的最优化估计,并对于偏振噪声进行分析。
发明内容
为了解决上述背景技术提到的技术问题,本发明提出了粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法。
为了实现上述技术目的,本发明的技术方案为:
一种粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法,包括以下步骤:
(1)测试粗糙表面物质的穆勒矩阵;
(2)计算粗糙表面物质穆勒矩阵对应的相干矩阵的特征值并分析特征值;
(3)计算粗糙表面物质穆勒矩阵的退偏振系数;
(4)计算粗糙表面物质出射光斯托克斯矢量的平均偏振度;
(5)将穆勒矩阵和斯托克斯矢量同时分解,建立分解物理方程,进行穆勒琼斯矩阵的最优化估计,并对偏振噪声进行分析。
进一步地,在步骤(1)中,采用穆勒偏振成像仪作为测试系统。
进一步地,在步骤(1)中,所述穆勒矩阵的表达式如下:
S'=MS(ε,θ)
其中,S'为出射光的斯托克斯矢量,S(ε,θ)为入射光的斯托克斯矢量,ε和θ为椭圆偏振参数,M为粗糙表面物质的穆勒矩阵。
进一步地,在步骤(2)中,对特征值的分析包括物理可靠性分析,判断穆勒矩阵对应的相干矩阵的特征值是否准确,如果有特征值小于0,则表示穆勒矩阵的测量具有噪声和误差。
进一步地,在步骤(4)中,出射光斯托克斯矢量的偏振度计算式如下:
其中,p为出射光斯托克斯矢量的偏振度,Sj'为出射光斯托克斯矢量的第j个元素,j=0,1,2,3;
出射光斯托克斯矢量的平均偏振度计算式如下:
其中,为出射光斯托克斯矢量的平均偏振度,pn为第n个出射光斯托克斯矢量的偏振度,N为斯托克斯矢量的数目。
进一步地,在邦加球坐标内取N组不同的椭圆偏振参数(ε,θ)进行数值方法处理,得到N个入射光斯托克斯矢量,从而得到N个出射光斯托克斯矢量。
进一步地,在步骤(5)中,所述分解物理方程如下:
其中,为入射光斯托克斯矢量的平均值,/>为出射光斯托克斯矢量的平均值,/>为/>中的第j个元素,j=0,1,2,3;d为步骤(3)计算的退偏振系数,/>λi为步骤(2)求出的穆勒矩阵对应的相干矩阵的第i个特征值,i=1,2,3,4;M1,1为穆勒矩阵的第一个元素;MJ为粗糙表面物质的穆勒琼斯矩阵估计,MD为完全退偏振矩阵,用MD表征偏振噪声。
采用上述技术方案带来的有益效果:
本发明能够同时对于斯托克斯矢量和穆勒矩阵进行分解,可以对粗糙表面物质进行穆勒-琼斯矩阵的准确估计,并对于偏振噪声进行分析。本发明简单,实用,有效,物理意义明确。
附图说明
图1是本发明的方法流程图。
具体实施方式
以下将结合附图,对本发明的技术方案进行详细说明。
本发明设计了一种粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法,如图1所示,步骤如下:
步骤1、测试粗糙表面物质的穆勒矩阵。
在步骤1中,优选地,采用穆勒偏振成像仪作为测试系统。
穆勒矩阵的表达式如下:
S'=MS(ε,θ)
其中,S'为出射光的斯托克斯矢量,S(ε,θ)为入射光的斯托克斯矢量,ε和θ为椭圆偏振参数,M为粗糙表面物质的穆勒矩阵。
步骤2、计算粗糙表面物质穆勒矩阵对应的相干矩阵的特征值并分析特征值。
在步骤2中,优选地,特征值分析包括物理可靠性分析,该过程是通过计算相干矩阵的特征值,判断穆勒矩阵的测量时否准确,即特征值λi≥0(i=1,2,3,4),如果有特征值λi<0,则穆勒矩阵的测量具有噪声和误差,如果λi是一个较大的负值,说明粗糙物质的穆勒矩阵不准确,测量结果不可靠,因此需要进行重复测量以得到更加准确可靠的结果。
步骤3、计算粗糙表面物质穆勒矩阵的退偏振系数。
在步骤3中,优选地,穆勒矩阵的退偏振系数的计算方法如下:
其中,d为退偏振系数。
步骤4、计算粗糙表面物质出射光斯托克斯矢量的平均偏振度。
在步骤4中,优选地,在邦加球坐标内取N组不同的(ε,θ)进行数值方法处理:
εn=-π/4+nπ/2000,θn=nπ/1000
然后根据下式计算相应的N个入射光斯托克斯矢量:
其中,Sj(j=0,1,2,3)为入射光的斯托克斯矢量元素。
根据入射光的斯托克斯矢量得到出射光的斯托克斯矢量,根据下式计算出射光斯托克斯矢量的平均偏振度:
其中,为出射光斯托克斯矢量的平均偏振度,pn为第n个出射光斯托克斯矢量的偏振度,S'j(j=0,1,2,3)为出射光的斯托克斯矢量元素。
步骤5、将穆勒矩阵和斯托克斯矢量同时分解,建立分解物理方程,进行穆勒琼斯矩阵的最优化估计,并对偏振噪声进行分析。
在步骤5中,优选地,分解物理方程如下:
其中,为入射光斯托克斯矢量的平均值,/>为出射光斯托克斯矢量的平均值,为/>中的第j个元素;d为退偏振系数,M1,1为穆勒矩阵的第一个元素,MJ为粗糙表面物质的穆勒琼斯矩阵估计,MD为完全退偏振矩阵,用MD表征偏振噪声。
粗糙物质的退偏振穆勒矩阵可以分解为偏振部分和退偏振部分,其偏振部分对应物质的穆勒琼斯矩阵。退偏振部分对应穆勒矩阵测量中的散射噪声、实验误差、背景等,其中退偏振部分的强度系数可以利用粗糙物质穆勒矩阵相对应的相干矩阵的特征值的函数表示。因为各项同性物质出射光的偏振度相等,因此可以将出射光的斯托克斯矢量分解为偏振部分,退偏振部分,其中退偏振部分可以看做偏振噪声。对于粗糙物质,上述分解物理方程表明穆勒矩阵和斯托克斯矢量两种不同的分解过程物理意义上是相同的。由此,可以得到穆勒琼斯矩阵的最佳估计MJ,以及偏振噪声MD。如果偏振噪声超过系统噪声,一般为2%,则在计算当中噪声不能忽略,可以利用上述分解物理方程予以消除。
实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。

Claims (4)

1.一种粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法,其特征在于,包括以下步骤:
(1)测试粗糙表面物质的穆勒矩阵;
(2)计算粗糙表面物质穆勒矩阵对应的相干矩阵的特征值并分析特征值;
(3)计算粗糙表面物质穆勒矩阵的退偏振系数;
(4)计算粗糙表面物质出射光斯托克斯矢量的平均偏振度;
(5)将穆勒矩阵和斯托克斯矢量同时分解,建立分解物理方程,进行穆勒琼斯矩阵的最优化估计,并对偏振噪声进行分析;
在步骤(4)中,出射光斯托克斯矢量的偏振度计算式如下:
其中,p为出射光斯托克斯矢量的偏振度,S'j为出射光斯托克斯矢量的第j个元素,j=0,1,2,3;
出射光斯托克斯矢量的平均偏振度计算式如下:
其中,为出射光斯托克斯矢量的平均偏振度,pn为第n个出射光斯托克斯矢量的偏振度,N为斯托克斯矢量的数目;
在邦加球坐标内取N组不同的椭圆偏振参数(ε,θ)进行数值方法处理,得到N个入射光斯托克斯矢量,从而得到N个出射光斯托克斯矢量;
在步骤(5)中,所述分解物理方程如下:
其中,为入射光斯托克斯矢量的平均值,/>为出射光斯托克斯矢量的平均值,/>为/>中的第j个元素,j=0,1,2,3;d为步骤(3)计算的退偏振系数,/>λi为步骤(2)求出的穆勒矩阵对应的相干矩阵的第i个特征值,i=1,2,3,4;M1,1为穆勒矩阵的第一个元素;MJ为粗糙表面物质的穆勒琼斯矩阵估计,MD为完全退偏振矩阵,用MD表征偏振噪声。
2.根据权利要求1所述粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法,其特征在于,在步骤(1)中,采用穆勒偏振成像仪作为测试系统。
3.根据权利要求1所述粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法,其特征在于,在步骤(1)中,所述穆勒矩阵的表达式如下:
S'=MS(ε,θ)
其中,S'为出射光的斯托克斯矢量,S(ε,θ)为入射光的斯托克斯矢量,ε和θ为椭圆偏振参数,M为粗糙表面物质的穆勒矩阵。
4.根据权利要求1所述粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法,其特征在于,在步骤(2)中,对特征值的分析包括物理可靠性分析,判断穆勒矩阵对应的相干矩阵的特征值是否准确,如果有特征值小于0,则表示穆勒矩阵的测量具有噪声和误差。
CN202011263093.9A 2020-11-12 2020-11-12 粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法 Active CN112417370B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011263093.9A CN112417370B (zh) 2020-11-12 2020-11-12 粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011263093.9A CN112417370B (zh) 2020-11-12 2020-11-12 粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法

Publications (2)

Publication Number Publication Date
CN112417370A CN112417370A (zh) 2021-02-26
CN112417370B true CN112417370B (zh) 2024-04-30

Family

ID=74831867

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011263093.9A Active CN112417370B (zh) 2020-11-12 2020-11-12 粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法

Country Status (1)

Country Link
CN (1) CN112417370B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113984676B (zh) * 2021-10-20 2023-03-24 南京航空航天大学 一种多次散射介质穆勒矩阵分解方法
CN114003195B (zh) * 2021-12-29 2022-04-19 苏州大学 快速计算三维偏振维度的方法、设备及存储介质
WO2023171058A1 (ja) * 2022-03-08 2023-09-14 東レ株式会社 偏光評価装置
CN116482123A (zh) * 2023-06-20 2023-07-25 泉州装备制造研究所 一种光学元件表面微观缺陷检测方法、系统、设备及介质
CN117805852B (zh) * 2024-02-26 2024-05-10 天津市天开海洋科技有限公司 一种基于虚拟照明调制实现水下散射抑制的偏振成像方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111060711A (zh) * 2019-12-31 2020-04-24 华中科技大学 一种基于斯托克斯矢量的光学转速测量系统及测量方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111060711A (zh) * 2019-12-31 2020-04-24 华中科技大学 一种基于斯托克斯矢量的光学转速测量系统及测量方法

Also Published As

Publication number Publication date
CN112417370A (zh) 2021-02-26

Similar Documents

Publication Publication Date Title
CN112417370B (zh) 粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法
Oyaizu et al. Photometric redshift error estimators
CN105486655B (zh) 基于红外光谱智能鉴定模型的土壤有机质快速检测方法
CN110751087B (zh) 一种基于eof的无人机信号识别系统和方法
CN106056070A (zh) 基于低秩矩阵恢复和稀疏表示的sar目标识别方法
Brombacher et al. Calibration of the repeatability of foraminiferal test size and shape measures with recommendations for future use
US10215642B2 (en) System and method for polarimetric wavelet fractal detection and imaging
CN111241991B (zh) 一种基于遥感影像的水质分类方法
CN107316009B (zh) 基于张量线性判别分析降维的高光谱图像目标检测方法
CN111879709A (zh) 湖泊水体光谱反射率检验方法及装置
CN117889965B (zh) 一种中短波双色红外探测器的性能测试方法
CN115420707A (zh) 一种污水近红外光谱的化学需氧量评估方法及系统
CN114611582A (zh) 一种基于近红外光谱技术分析物质浓度的方法及系统
CN117113008B (zh) 基于x射线的板厚测量方法
Li et al. Spectral–spatial depth-based framework for hyperspectral underwater target detection
Liu et al. Detection of Apple Taste Information Using Model Based on Hyperspectral Imaging and Electronic Tongue Data.
CN110793920B (zh) 一种化学成像与高光谱联用的大气遥测方法
CN111366573B (zh) 基于libs光谱成分分析结果的评价方法
CN113984676A (zh) 一种多次散射介质穆勒矩阵分解方法
Głowacz et al. Comparison of various data analysis techniques applied for the classification of oligopeptides and amino acids by voltammetric electronic tongue
CN109738392A (zh) 面向tdlas在线检测瓶内氧气浓度的压缩感知重构方法及系统
Li et al. Voltammetric electronic tongue for the qualitative analysis of milk adulterated with urea combined with multi-way data analysis
CN115015131A (zh) 红外光谱训练集样本筛选方法
CN111624194B (zh) 基于激光诱导击穿光谱的叶片元素无损测量方法
Krishnamurthy et al. Target detection performance bounds in compressive 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