CN113984676A - 一种多次散射介质穆勒矩阵分解方法 - Google Patents
一种多次散射介质穆勒矩阵分解方法 Download PDFInfo
- Publication number
- CN113984676A CN113984676A CN202111218316.4A CN202111218316A CN113984676A CN 113984676 A CN113984676 A CN 113984676A CN 202111218316 A CN202111218316 A CN 202111218316A CN 113984676 A CN113984676 A CN 113984676A
- Authority
- CN
- China
- Prior art keywords
- matrix
- mueller
- mueller matrix
- multiple scattering
- follows
- 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.)
- Granted
Links
- 239000011159 matrix material Substances 0.000 title claims abstract description 126
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 22
- 230000010287 polarization Effects 0.000 claims abstract description 32
- 230000028161 membrane depolarization Effects 0.000 claims abstract description 21
- 239000000463 material Substances 0.000 claims abstract description 20
- 239000000126 substance Substances 0.000 claims abstract description 16
- 230000001427 coherent effect Effects 0.000 claims abstract description 10
- 238000012360 testing method Methods 0.000 claims abstract description 8
- 238000005259 measurement Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 9
- 239000013598 vector Substances 0.000 claims description 9
- 230000002999 depolarising effect Effects 0.000 claims description 5
- 239000002689 soil Substances 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241000219000 Populus Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229920002678 cellulose Polymers 0.000 description 1
- 239000001913 cellulose Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- CPBQJMYROZQQJC-UHFFFAOYSA-N helium neon Chemical group [He].[Ne] CPBQJMYROZQQJC-UHFFFAOYSA-N 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000004088 simulation 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/21—Polarisation-affecting properties
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Pathology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种多次散射介质穆勒矩阵分解方法,该方法包括:测试多次散射物质的穆勒矩阵;计算多次散射物质穆勒矩阵对应的相干矩阵的特征值并分析特征值;计算多次散射物质的退偏振矩阵及穆勒琼斯矩阵;将穆勒琼斯矩阵和退偏振矩阵建立分解物理方程,采用最小二乘法进行穆勒琼斯矩阵和散射矩阵的最优化估计,并对偏振噪声进行分析。本发明可以处理表面粗糙介质各项异性退偏的问题,对具有粗糙表面的多次散射物质的穆勒矩阵进行穆勒矩阵分解,进行穆勒‑琼斯矩阵的最优化估计,对于多次散射物质具有更强的适用性。
Description
技术领域
本发明属于偏振技术领域,特别涉及一种多次散射介质穆勒矩阵分解方法。
背景技术
穆勒矩阵测量技术是一种重要的光学测量技术,近年来,穆勒矩阵测量技术在目标识别,物质分类,生物分析,医学诊断等领域的应用更加广泛。随着穆勒矩阵测量技术的快速发展和应用,穆勒矩阵分解方法对于物质退偏振研究和实验数据解释分析是一种有效的研究方法。
近年来,研究人员提出了几种典型的方法,包括Cloude分解法, Christian 分解法, Roy-Brehonnet 方法,Lu-Chipman 方法等。Cloude分解方法用于进行穆勒琼斯矩阵的最优化估计,并对于偏振噪声进行分析。Christian 分解法用于计算退偏振系数,计算偏振残差。Roy-Brehonnet 方法是基于Cloude分解法和Christian 分解法,可以对退偏振物质进行偏振分解的最优化估计及噪声误差分析。Lu-Chipman 方法可以对光学系统中的不同偏振特征包括,衰减,延迟,退偏振进行子矩阵提取。
常见光学器件的穆勒矩阵的不同偏振态的偏振度基本一致,而大多数散射物质的退偏振由于多次散射作用,其不同的偏振态具有不同的偏振度,例如粗糙的岩石,泥土,沙子等。其退偏振机理机理复杂,多次散射在穆勒矩阵的测量中主要会导致:
a)多次散射引起附加的退偏振效应;b)多次散射改变了入射光的偏振态,导致部分信息丢失;c)多次散射强度弱,动态范围大。
因而,这些分解方法在进行散射物质的穆勒矩阵分解时就遇到了挑战,其应用也受到了一定的限制。
发明内容
本发明的目的在于提供一种多次散射介质穆勒矩阵分解方法。
实现本发明目的的技术解决方案为:一种多次散射介质穆勒矩阵分解方法,包括以下步骤:
步骤1、测试多次散射物质的穆勒矩阵;
步骤2、计算多次散射物质穆勒矩阵对应的相干矩阵的特征值并分析特征值;
步骤3、计算多次散射物质的退偏振矩阵:
步骤4、将穆勒矩阵进行分解,根据最小二乘法进行穆勒琼斯矩阵的最优化估计,并对偏振噪声进行分析。
进一步的,在步骤2中,穆勒矩阵对应的相干矩阵计算式如下:
进一步的,穆勒矩阵元素计算式如下:
进一步的,在步骤2中,对特征值的分析包括物理可靠性分析,判断穆勒矩阵对应的相干矩阵的特征值是否准确,如果有特征值小于0,则表示穆勒矩阵的测量具有噪声和误差。
进一步的,在步骤2中,所述相干矩阵分解公式如下:
进一步的,在步骤4中,实验残差2-范数的计算式为
进一步的,采用穆勒偏振成像仪作为测试系统,所述穆勒矩阵的表达式如下:
与现有技术相比,本发明具有如下有益效果:本发明可以分解复杂散射物质的穆勒矩阵,同时可以提取穆勒琼斯矩阵;本发明可以处理表面粗糙介质各项异性退偏的问题,对具有粗糙表面的多次散射物质的穆勒矩阵进行穆勒矩阵分解,进行穆勒-琼斯矩阵的最优化估计,对于多次散射物质具有更强的适用性。比传统方法实用,有效,物理意义明确。
附图说明
图1是本发明的多次散射介质穆勒矩阵分解方法流程图。
图2是本发明的实验装置原理图。
图3是本发明方法(EAMMD)与Roy-Brehonnet 方法(RBM)的比较示意图,其中(a)、(b)分别为土壤样品和叶片样品。
具体实施方式
以下将结合附图,对本发明的技术方案进行详细说明。
本发明提出了一种多次散射介质穆勒矩阵分解方法,如图1所示,步骤如下:
步骤1、测试粗糙表面物质的穆勒矩阵。
在步骤1中,优选地,采用穆勒偏振成像仪作为测试系统。
穆勒矩阵的表达式如下:
步骤2、 计算粗糙表面物质穆勒矩阵对应的相干矩阵的特征值并分析特征值。
在步骤2中,优选地,特征值分析包括物理可靠性分析,该过程是通过计算相干矩阵的特征值,判断穆勒矩阵的测量是否准确,即特征值,由于相干矩阵是厄米特半正定矩阵,因此满足,如果有特征值,则穆勒矩阵的测量数据具有噪声和误差,如果是一个较大的负值,说明粗糙物质的穆勒矩阵不准确,测量结果不可靠,因此需要进行重复测量以得到更加准确可靠的结果。
步骤3、计算粗糙表面物质多次散射物质穆勒矩阵的退偏振穆勒矩阵:
得到:
穆勒矩阵的退偏振系数d的计算方法如下:
步骤4、将穆勒矩阵进行分解,建立分解的物理方程,根据最小二乘法进行穆勒琼斯矩阵的最优化估计,并对偏振噪声进行分析。
分解物理方程如下:
多次散射物质的退偏振穆勒矩阵可以分解为偏振部分和退偏振部分,其偏振部分对应物质的穆勒琼斯矩阵。退偏振矩阵对应多次散射的散射噪声、实验误差、背景等。对于多次散射物质,上述分解物理方程表明穆勒矩阵和物理意义上是相同的。由此,可以得到穆勒琼斯矩阵的最佳估计,退偏振矩阵以及偏振噪声。如果实验误差和噪声超过系统一般误差,一般为±2%,则在计算当中误差不能忽略,可以利用上述分解物理方程予以消除。
下面结合实施例对本发明进行详细说明。
实施例
本发明中所利用的实验装置原理图如图2所示,其中,激光器为氦氖激光器,波长为632.8nm, W1、W2、W3为石英四分之一波片,P1、P2、P3为线偏光片。
检测器是 CCD 相机,图像尺寸:800×600 像素,4096 个灰度级。本发明中的Mueller 矩阵的获得是基于传统的穆勒矩阵测量方法,使用偏振态发生器(PSG、P2 和 W2)和偏振态分析仪(PSA、P3 和 W3)。测量样品的入射角从 40°到 70°,包括六个入射斯托克斯矢量,分别为线性水平、垂直、45°、-45°、右圆偏振光和左圆偏振光用于产生不同偏振态的偏振光。需要采集 36 张图像来计算样本的穆勒矩阵。实验在光学暗室中进行。
实验仿真部分:选取两种样品进行算法实验,(a)土壤是一种混合物,由各种矿物质组成。土壤样品是一种常见的具有微观结构和成分的散射材料,从室外的地面采集,样品的粒径从 0.2 mm 到 0.4 mm 不等。(b)杨树叶样品是心形扁平形状,其为厚度0.5 mm。叶片样品的主要成分是具有多层结构的纤维素。 叶子样品的粗糙表面具有网状结构。偏振测试的入射角度为40°到75°。
图3为本发明的方法(EAMMD)与Roy-Brehonnet 方法(RBM)的比较,其中(a)、(b)分别对应土壤样品和叶片样品。从图中可以发现本发明方法的实验残差小于Roy-Brehonnet方法。如图3(a)所示,对于泥土样品,当入射角度为65°时,本发明方法的实验残差不到RBM的1/3。
实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (9)
5.根据权利要求1所述的多次散射介质穆勒矩阵分解方法,其特征在于,在步骤2中,对特征值的分析包括物理可靠性分析,判断穆勒矩阵对应的相干矩阵的特征值是否准确,如果有特征值小于0,则表示穆勒矩阵的测量具有噪声和误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111218316.4A CN113984676B (zh) | 2021-10-20 | 2021-10-20 | 一种多次散射介质穆勒矩阵分解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111218316.4A CN113984676B (zh) | 2021-10-20 | 2021-10-20 | 一种多次散射介质穆勒矩阵分解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113984676A true CN113984676A (zh) | 2022-01-28 |
CN113984676B CN113984676B (zh) | 2023-03-24 |
Family
ID=79739496
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111218316.4A Active CN113984676B (zh) | 2021-10-20 | 2021-10-20 | 一种多次散射介质穆勒矩阵分解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113984676B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116482123A (zh) * | 2023-06-20 | 2023-07-25 | 泉州装备制造研究所 | 一种光学元件表面微观缺陷检测方法、系统、设备及介质 |
WO2023236406A1 (zh) * | 2022-06-09 | 2023-12-14 | 温州医科大学 | 背向散射Mueller矩阵的可逆极分解方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102262223A (zh) * | 2010-05-26 | 2011-11-30 | 中国科学院对地观测与数字地球科学中心 | 一种基于散射矩阵的相干目标分解方法 |
JP2017194340A (ja) * | 2016-04-20 | 2017-10-26 | 国立研究開発法人産業技術総合研究所 | 散乱試料の偏光特性の測定方法 |
CN111122452A (zh) * | 2019-12-28 | 2020-05-08 | 天津大学 | 基于穆勒矩阵的去散射成像方法 |
CN112147080A (zh) * | 2020-09-28 | 2020-12-29 | 福建师范大学 | 全自动时分探测的自由空间型 Mueller OCT成像系统及成像方法 |
CN112417370A (zh) * | 2020-11-12 | 2021-02-26 | 南京航空航天大学 | 粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法 |
WO2021087345A1 (en) * | 2019-11-01 | 2021-05-06 | Unm Rainforest Innovations | In-line angular optical multi-point scatterometry for nanomanufacturing systems |
-
2021
- 2021-10-20 CN CN202111218316.4A patent/CN113984676B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102262223A (zh) * | 2010-05-26 | 2011-11-30 | 中国科学院对地观测与数字地球科学中心 | 一种基于散射矩阵的相干目标分解方法 |
JP2017194340A (ja) * | 2016-04-20 | 2017-10-26 | 国立研究開発法人産業技術総合研究所 | 散乱試料の偏光特性の測定方法 |
WO2021087345A1 (en) * | 2019-11-01 | 2021-05-06 | Unm Rainforest Innovations | In-line angular optical multi-point scatterometry for nanomanufacturing systems |
CN111122452A (zh) * | 2019-12-28 | 2020-05-08 | 天津大学 | 基于穆勒矩阵的去散射成像方法 |
CN112147080A (zh) * | 2020-09-28 | 2020-12-29 | 福建师范大学 | 全自动时分探测的自由空间型 Mueller OCT成像系统及成像方法 |
CN112417370A (zh) * | 2020-11-12 | 2021-02-26 | 南京航空航天大学 | 粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法 |
Non-Patent Citations (4)
Title |
---|
GHOSH, WOOD, AND VITKIN: "Mueller matrix decomposition for extraction of individual polarization parameters from complex turbid media exhibiting multiple scattering, optical activity, and linear birefringence", 《JOURNAL OF BIOMEDICAL OPTICS》 * |
M. IQBAL ET AL.: "Comparative study of Mueller matrix transformation and polar decomposition for optical characterization of turbid media", 《OPTIK》 * |
权亚楠等: "基于相干矩阵完全分解的PolSAR图像分解", 《北京测绘》 * |
章延隽,王霞,贺思: "基于偏振双向反射分布函数的粗糙表面偏振特性", 《光学学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023236406A1 (zh) * | 2022-06-09 | 2023-12-14 | 温州医科大学 | 背向散射Mueller矩阵的可逆极分解方法 |
CN116482123A (zh) * | 2023-06-20 | 2023-07-25 | 泉州装备制造研究所 | 一种光学元件表面微观缺陷检测方法、系统、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113984676B (zh) | 2023-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113984676B (zh) | 一种多次散射介质穆勒矩阵分解方法 | |
Ourmazd et al. | Quantitative chemical lattice imaging: theory and practice | |
Dunkley et al. | Prospects for polarized foreground removal | |
Xue et al. | Blind image quality assessment using joint statistics of gradient magnitude and Laplacian features | |
Simoncelli | Modeling the joint statistics of images in the wavelet domain | |
Bobin et al. | Sparsity and adaptivity for the blind separation of partially correlated sources | |
Tu | Unsupervised signature extraction and separation in hyperspectral images: a noise-adjusted fast independent component analysis | |
Santos et al. | Estimate of the cosmological bispectrum from the maxima-1 cosmic microwave background map | |
Idier et al. | On the superresolution capacity of imagers using unknown speckle illuminations | |
CN112417370B (zh) | 粗糙表面物质的穆勒琼斯矩阵估计及偏振噪声分析方法 | |
Iordache et al. | Recent developments in sparse hyperspectral unmixing | |
CN105894018B (zh) | 基于深度多示例学习的极化sar图像分类方法 | |
CN106056070A (zh) | 基于低秩矩阵恢复和稀疏表示的sar目标识别方法 | |
WO2014065694A1 (ru) | Способ и устройство для оптического измерения распределения параметров дисперсных частиц в жидкостях и газах | |
CN103093431B (zh) | 基于pca字典和结构先验信息的压缩感知重构方法 | |
Veeraraghavan et al. | Increasing depth resolution of electron microscopy of neural circuits using sparse tomographic reconstruction | |
CN112837294A (zh) | 基于卷积自编码器图像扩增的热成像缺陷检测方法 | |
Pépy et al. | Cylindrical nanochannels in ion-track polycarbonate membranes studied by small-angle X-ray scattering | |
US10215642B2 (en) | System and method for polarimetric wavelet fractal detection and imaging | |
Liu et al. | The application of principal component analysis and non-negative matrix factorization to analyze time-resolved optical waveguide absorption spectroscopy data | |
Marcos-Caballero et al. | Scale-dependent dipolar modulation and the quadrupole-octopole alignment in the CMB temperature | |
US11227403B2 (en) | Anisotropic twicing for single particle reconstruction using autocorrelation analysis | |
Feng et al. | High-dimensional rank tests for sphericity | |
CN116429990A (zh) | 一种基于纤维素酶的燕麦种抗氧化活性检测方法和装置 | |
CN112307795A (zh) | 基于聚类分析提取统计特征量的物质甄别设备及其方法 |
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 |