CN218157520U - 一种穆勒矩阵测量系统 - Google Patents
一种穆勒矩阵测量系统 Download PDFInfo
- Publication number
- CN218157520U CN218157520U CN202221625027.6U CN202221625027U CN218157520U CN 218157520 U CN218157520 U CN 218157520U CN 202221625027 U CN202221625027 U CN 202221625027U CN 218157520 U CN218157520 U CN 218157520U
- Authority
- CN
- China
- Prior art keywords
- lcvr
- phase retarder
- light
- polarized light
- mueller matrix
- 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
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
一种穆勒矩阵测量系统，属于偏振探测技术领域，解决如何快速且完整的测量样品的穆勒矩阵问题，通过光源竖直向上发射调制光，第一垂直线偏振片用于将入射的调制光转换为线偏振光；第一相位延迟片、第二相位延迟片调制线偏振光的偏振态，数字微反射镜加载有条纹图片，带有条纹图片的线偏振光在待测生物组织中发生散射，散射光依次经过第三相位延迟片、第一相位延迟片、第二垂直线偏振片后由相机捕捉；通过第一相位延迟片、第二相位延迟片、第三相位延迟片、第四相位延迟片控制调制线偏振光的偏振态，来控制偏振光对待测生物组织的测量深度，共调制16偏振光的偏振态，能够测量出样品的完整穆勒矩阵。
Description
技术领域
本实用新型属于偏振探测技术领域，涉及一种穆勒矩阵测量系统。
背景技术
生物组织携带丰富的光学信息(如透射、反射、衰减信息、延迟信息、退偏信息)，能够反映生物组织相关的生理特征。例如，一般病理组织(例如癌症)因为高度组织化的显微结构被破环，组织形态固有的异质性丧失，从而更好的保留了探测光的偏振状态，具有较低的退偏性。疤痕组织含有密集的胶原纤维，探测光束在组织样本中传播时发生更多的正交偏振态相移，从而具有高线性延迟性。生物组织基质的异常改变(例如胶原蛋白的重塑)导致其组织散射和吸收特性的改变，从而使组织的去偏振特性改变。
使用斯托克斯矢量来定义光的偏振状态，并使用穆勒矩阵来描述光通过设计的系统时偏振的变化。斯托克斯矢量表达式如下公式：S＝[S0,S1,S2,S3]T，S0表示总光强，即水平偏振分量和垂直偏振分量之和，S1表示水平偏振分量和垂直偏振分量之差，S2表示45°和135°偏振分量之差，S3表示右旋圆偏振和左旋圆偏振分量之差。
穆勒矩阵的形式如下：
Mueller矩阵作为入射Stokes矢量和出射Stokes矢量之间的变换矩阵，满足下面的公式：Sout＝MSin，式中，Sin表示入射光的Stokes矢量，Sout表示出射光的Stokes矢量。
为了对病理组织进行相关分析，一种方法是通过偏振光成像技术检测病理组织的光学信息。并且这也是目前广泛采用的方法。例如，通过偏振成像为射频消融治疗心率失常症状提供帮助，通过成像技术精密确定射频损伤的范围。通过偏振光成像用以确定结肠和宫颈中健康组织和癌变组织的范围。通过偏振光成像测量浑浊介质中的葡萄糖水平。除了提供结构和形态信息，偏振光成像还提供了胶原组织的偏振信息，包括去偏振，线性延迟和旋光等。因此，使用偏振光成像测量胶原组织的穆勒矩阵(Mueller)信息成为可行的手段。
目前，结合空间频率成像和偏振光成像技术的偏振空间频率成像正在被深入研究。例如，Bin Yang等人设计了一套偏振空间频率系统装置用以检测牛肌腱组织的纤维方向，并取得了良好的实验结果。WILL GOTH设计了相似的装置并对猪心脏瓣膜瓣膜中胶原纤维排列进行检测。Samuel V.Jett等人更进一步对二尖瓣叶组织胶原在机械负载下的纤维排列作出了定量分析，也测出了部分的穆勒矩阵信息。但是他们的实验装置存在着些许不足之处：在实验测量中，起偏系统和检偏系统共用一组偏振片，提供的偏振光偏振态范围较小，仅仅考虑线偏振光的偏振态，没有考虑圆偏振光和其他偏振态的影响，对测量待测生物组织穆勒矩阵的完整性有所不足。
实用新型内容
本实用新型的目的在于提出一种新的穆勒矩阵测量系统，以解决现有技术中待测生物组织穆勒矩阵测量的完整性不足的问题。
本实用新型是通过以下技术方案解决上述技术问题的：
一种穆勒矩阵测量系统，包括：光源系统(11)、起偏子系统(12)、待测样品台(13)、检偏子系统(14)；所述的光源系统(11)包括：光源(111)、透镜(112)、带通滤波器(113)；所述的起偏子系统(12)包括：数字微反射镜(121)、第一垂直线偏振片(P1)、第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)、第一成像透镜(L1)、第一反射镜(122)；所述的检偏子系统(14)包括：第二反射镜(141)、第三相位延迟片(LCVR3)、第四相位延迟片(LCVR4)、第二垂直线偏振片(P2)、第二成像透镜(L2)、相机(142)；
所述的光源(111)竖直向上发射调制光，经过透镜(112)、带通滤波器(113)后产生单一波长的调制光，所述的数字微反射镜(121)将单一波长的调制光水平反射出去，沿着水平调制光的光路、从左到右依次设置第一垂直线偏振片(P1)、第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)、第一成像透镜(L1)、第一反射镜(122)，所述的数字微反射镜(121)加载有条纹图片，所述的第一垂直线偏振片(P1)用于将入射的调制光转换为线偏振光，所述的第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)用于调制线偏振光的偏振态；水平调制光经过第一垂直线偏振片(P1)、第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)、第一成像透镜(L1)后，再由第一反射镜(122)将带有条纹图片的线偏振光投射到待测样品台(13)上放置的待测生物组织表面；带有条纹图片的线偏振光在待测生物组织中发生散射，散射光通过第二反射镜(141)反射后水平射出，再依次经过第三相位延迟片(LCVR3)、第四相位延迟片(LCVR4)、第二垂直线偏振片(P2)、第二成像透镜(L2)后，由相机(142)捕捉。
本实用新型的穆勒矩阵测量系统的光源系统(11)竖直向上发射调制光，第一垂直线偏振片(P1)用于将入射的调制光转换为线偏振光；第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)调制线偏振光的偏振态，数字微反射镜(121)加载有条纹图片，带有条纹图片的线偏振光在待测生物组织中发生散射，散射光依次经过第三相位延迟片(LCVR3)、第一相位延迟片(143)、第二垂直线偏振片(P2)后，由相机(142)捕捉；通过第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)、第三相位延迟片(LCVR3)、第四相位延迟片(LCVR4)控制调制线偏振光的偏振态，来控制偏振光对待测生物组织的测量深度，共调制16偏振光的偏振态，能够测量出样品的完整穆勒矩阵。
在一种实施例中，所述的第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)的快轴与水平方向夹角分别为45°和90°。
在一种实施例中，所述的第三相位延迟片(LCVR3)、第四相位延迟片(LCVR4)的快轴与水平方向夹角分别为90°和45°。
在一种实施例中，所述的光源系统(11)采用复合光源。
在一种实施例中，所述的相机(142)采用工业级CCD相机。
本实用新型的优点在于：
本实用新型的一种穆勒矩阵测量系统，光源系统(11)竖直向上发射调制光，第一垂直线偏振片(P1)用于将入射的调制光转换为线偏振光；第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)调制线偏振光的偏振态，数字微反射镜(121)加载有条纹图片，带有条纹图片的线偏振光在待测生物组织中发生散射，散射光依次经过第三相位延迟片(LCVR3)、第一相位延迟片(143)、第二垂直线偏振片(P2)后，由相机(142)捕捉；通过第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)、第三相位延迟片(LCVR3)、第四相位延迟片(LCVR4)控制调制线偏振光的偏振态，来控制偏振光对待测生物组织的测量深度，共调制16偏振光的偏振态，能够测量出样品的完整穆勒矩阵。
附图说明
图1是本实用新型实施例的一种穆勒矩阵测量系统的结构图。
具体实施方式
为使本实用新型实施例的目的、技术方案和优点更加清楚，下面将结合本实用新型实施例，对本实用新型实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本实用新型一部分实施例，而不是全部的实施例。基于本实用新型中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本实用新型保护的范围。
下面结合说明书附图以及具体的实施例对本实用新型的技术方案作进一步描述：
实施例一
如图1所示，一种穆勒矩阵测量系统，包括：光源系统(11)、起偏子系统(12)、待测样品台(13)、检偏子系统(14)；所述的光源系统(11)包括：光源(111)、透镜(112)、带通滤波器(113)；所述的起偏子系统(12)包括：数字微反射镜(121)、第一垂直线偏振片(P1)、第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)、第一成像透镜(L1)、第一反射镜(122)；所述的检偏子系统(14)包括：第二反射镜(141)、第三相位延迟片(LCVR3)、第四相位延迟片(LCVR4)、第二垂直线偏振片(P2)、第二成像透镜(L2)、相机(142)。
所述的光源(111)竖直向上发射调制光，经过透镜(112)、带通滤波器(113)后产生单一波长的调制光，所述的光源系统(11)可以采用复合光源。所述的数字微反射镜(121)将单一波长的调制光水平反射出去，沿着水平调制光的光路、从左到右依次设置第一垂直线偏振片(P1)、第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)、第一成像透镜(L1)、第一反射镜(122)，所述的数字微反射镜(121)加载有条纹图片，所述的第一垂直线偏振片(P1)用于将入射的调制光转换为线偏振光，所述的第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)用于调制线偏振光的偏振态，第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)的快轴与水平方向夹角分别为45°和90°；水平调制光经过第一垂直线偏振片(P1)、第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)、第一成像透镜(L1)后，再由第一反射镜(122)将带有条纹图片的线偏振光投射到待测样品台(13)上放置的待测生物组织表面；带有条纹图片的线偏振光在待测生物组织中发生散射，散射光通过第二反射镜(141)反射后水平射出，再依次经过第三相位延迟片(LCVR3)、第四相位延迟片(LCVR4)、第二垂直线偏振片(P2)、第二成像透镜(L2)后，由相机(142)捕捉，所述的第三相位延迟片(LCVR3)、第四相位延迟片(LCVR4)的快轴与水平方向夹角分别为90°和45°；所述的相机(142)采用工业级CCD(电荷耦合器件，Charge Coupled Device)相机。
测量方法如下：
1、确定条纹图片空间频率
在实验测量中，需要对生物组织的表皮层进行测量，因此，需要确定在表皮层成像深度下，条纹图片的空间频率fx是多少。运用空间频率成像方法(空间频率成像方法是现有的方法，具体参见文献：Gioux S,Mazhar A,Cuccia DJ.Spatial frequency domainimaging in 2019:principles,applications,and perspectives)进行三次测量，在每一次测量中投影一张不同的条纹图片，图案的明暗条纹按三角函数形式周期变换，条纹图片的相位分别为0°，120°，240°。
设入射在待测样品的入射光强度的函数表达式为：
其中，s0、M0、fx、x、
分别是光源的光照强度、空间调制深度、条纹图片的空间频率、空间坐标以及空间相位。
在测量过程中，结构光在样品组织中散射并由探测相机捕获到散射图像。散射光光照强度I(x,fx)可以被分解为直流DC分量IDC(x)和交流AC分量IAC(x,fx)。
I(x,fx)＝IDC(x)+IAC(x,fx) (2)
DC分量IDC只与空间位置x相关，AC分量IAC是空间位置x和空间频率fx的函数。具有如下表达式：
MAC代表反射光子密度的振幅。一般利用三相位移法进行求解。在测量过程中投影了三张条纹图片，三张图像光强值表示I0(x,fx),I120(x,fx),I240(x,fx)。有以下表达式：
同时也有DC分量IDC(x,fx)表达式：
求出以上分量后，在用相同的步骤测量漫反射板的相关分量。则样品组织的漫反射率Rd可以表示为：
其中分母为漫反射板的AC分量，Rd,ref(x,fx)为漫反射板的漫反射率(0.99)。
又根据漫反射近似理论，漫反射率，空间频率，吸收系数，散射系数存在以下表达式：
其中，A为一个常数，
utr＝ua+us'。ua是吸收系数，us'是散射系数。所以整个表达式空间频率fx是自变量，Rd是因变量。利用最小二乘法将不同空间频率和对应的漫反射率进行拟合，即可求出吸收系数和散射系数。
求得吸收系数和散射系数后，成像深度与吸收系数，散射系数，空间频率具有如下表达式。
设定成像深度，即可确定条纹图片的空间频率fx，继而进行步骤2的测量。
2、使用穆勒矩阵测量系统，进行16组偏振态组合的测量
由于穆勒矩阵有16个元素，在整个测量中需要对偏振光进行16次不同的偏振态调制。所以，需要对起偏系统进行4次偏振态调制(通过对液晶相位延迟器加载不同电压)，对检偏系统进行4次调制，两两组合，可以得到16次不同的偏振态组合。而又因为使用了空间频率成像方法，需要投影三张条纹图片，所以，在一种偏振态调制下，需要对待测生物组织成像三次才能得到最终的光强值结果。即起偏系统确认一种偏振态，检偏系统确认一种偏振态，在这种偏振态组合条件下，依次投影三张条纹图片得到三张成像图片。重复以上步骤，进行剩下15组偏振态组合的测量，共计得到16组偏振态组合下的成像结果，总计48张成像图片。
3、进行穆勒矩阵计算的相关准备
根据步骤1所示，计算16组偏振态下散射光的光强值。
Iout(x,fx)＝IDC(x,fx)+IAC(x,fx) (9)
根据图1所示，经目标反射后CCD探测到的出射光斯托克斯矢量Iout(x,fx)可以表示为：
IOut＝MP2*MLCVR4*MLCVR3*MSample*MLCVR2*MLCVR1*MP1*Iin (10)
其中，Iin为入射光的斯托克斯矢量，MP1为起偏系统垂直偏振片的穆勒矩阵，MLCVR1和MLCVR2分别为液晶相位延迟器LCVR1和LCVR2的穆勒矩阵，MSample为待测生物组织的穆勒矩阵，MLCVR4和MLCVR3分别为液晶相位延迟器LCVR4和LCVR3的穆勒矩阵，MP2为检偏系统垂直偏振片的穆勒矩阵。
在步骤2的测量过程中，LCVR1和LCVR2快轴与水平方向夹角分别为45°和90°(LCVR4和LCVR3也为45°和90°)，由于入射光被一个垂直线偏振片调制，因此入射斯托克斯矢量Iin＝[1,-1,0,0]T。
因此从LCVR2出射光的斯托克斯矢量ILCVR2可以表示为：
其中，δ1和δ2分别是两个液晶相位延迟器的相位延迟，其通过设置在液晶相位延迟器上的加载电压决定。
将(11)式代入(10)式，并令MPi＝MP2*MLCVR4*MLCVR3可得：
在步骤2中设置检偏系统液晶相位延迟器加载电压，分别使检偏系统进行四次的偏振态调制。
MLCVR4，MLCVR3与起偏系统液晶相位延迟器LCVR1,LCVR2穆勒矩阵相同，由加载电压确定。垂直偏振片的穆勒矩阵同样可求，所以是MPi可知的。
其中，δ1′和δ2′分别是两个液晶相位延迟器的相位延迟；
单独取IOut的第一个数(总光强值)，即只计算MPi的第一行，则有：
其中，a为由δ1，δ2，δ1′和δ2′决定的系数，对于每一组测量，该值不同。总计测量16组数据，则可得到16个等式方程组：I＝A*M，通过解矩阵方程可求出最终的穆勒矩阵。
以上实施例仅用以说明本实用新型的技术方案，而非对其限制；尽管参照前述实施例对本实用新型进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本实用新型各实施例技术方案的精神和范围。
Claims (5)
1.一种穆勒矩阵测量系统，其特征在于，包括：光源系统(11)、起偏子系统(12)、待测样品台(13)、检偏子系统(14)；所述的光源系统(11)包括：光源(111)、透镜(112)、带通滤波器(113)；所述的起偏子系统(12)包括：数字微反射镜(121)、第一垂直线偏振片(P1)、第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)、第一成像透镜(L1)、第一反射镜(122)；所述的检偏子系统(14)包括：第二反射镜(141)、第三相位延迟片(LCVR3)、第四相位延迟片(LCVR4)、第二垂直线偏振片(P2)、第二成像透镜(L2)、相机(142)；
所述的光源(111)竖直向上发射调制光，经过透镜(112)、带通滤波器(113)后产生单一波长的调制光，所述的数字微反射镜(121)将单一波长的调制光水平反射出去，沿着水平调制光的光路、从左到右依次设置第一垂直线偏振片(P1)、第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)、第一成像透镜(L1)、第一反射镜(122)，所述的数字微反射镜(121)加载有条纹图片，所述的第一垂直线偏振片(P1)用于将入射的调制光转换为线偏振光，所述的第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)用于调制线偏振光的偏振态；水平调制光经过第一垂直线偏振片(P1)、第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)、第一成像透镜(L1)后，再由第一反射镜(122)将带有条纹图片的线偏振光投射到待测样品台(13)上放置的待测生物组织表面；带有条纹图片的线偏振光在待测生物组织中发生散射，散射光通过第二反射镜(141)反射后水平射出，再依次经过第三相位延迟片(LCVR3)、第四相位延迟片(LCVR4)、第二垂直线偏振片(P2)、第二成像透镜(L2)后，由相机(142)捕捉。
2.根据权利要求1所述的穆勒矩阵测量系统，其特征在于，所述的第一相位延迟片(LCVR1)、第二相位延迟片(LCVR2)的快轴与水平方向夹角分别为45°和90°。
3.根据权利要求2所述的穆勒矩阵测量系统，其特征在于，所述的第三相位延迟片(LCVR3)、第四相位延迟片(LCVR4)的快轴与水平方向夹角分别为90°和45°。
4.根据权利要求1所述的穆勒矩阵测量系统，其特征在于，所述的光源系统(11)采用复合光源。
5.根据权利要求1所述的穆勒矩阵测量系统，其特征在于，所述的相机(142)采用工业级CCD相机。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202221625027.6U CN218157520U (zh) | 2022-06-24 | 2022-06-24 | 一种穆勒矩阵测量系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202221625027.6U CN218157520U (zh) | 2022-06-24 | 2022-06-24 | 一种穆勒矩阵测量系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN218157520U true CN218157520U (zh) | 2022-12-27 |
Family
ID=84576139
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202221625027.6U Active CN218157520U (zh) | 2022-06-24 | 2022-06-24 | 一种穆勒矩阵测量系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN218157520U (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI831641B (zh) * | 2023-03-10 | 2024-02-01 | 國立聯合大學 | 雙液晶穆勒偏振檢測系統及檢測非酵素性糖化膠原蛋白的方法 |
-
2022
- 2022-06-24 CN CN202221625027.6U patent/CN218157520U/zh active Active
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI831641B (zh) * | 2023-03-10 | 2024-02-01 | 國立聯合大學 | 雙液晶穆勒偏振檢測系統及檢測非酵素性糖化膠原蛋白的方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10634560B2 (en) | Multi-spectral polarimetric variable optical device and imager | |
| US10126111B2 (en) | Associating optical coherence tomography (OCT) data with visual imagery of a sample | |
| US6262818B1 (en) | Method for simultaneous amplitude and quantitative phase contrast imaging by numerical reconstruction of digital holograms | |
| WO2010021148A1 (ja) | 形状・傾斜検知及び／又は計測光学装置及び方法並びにその関連装置 | |
| US6208416B1 (en) | Method and apparatus for measuring shape of objects | |
| US7826059B2 (en) | Method and apparatus for polarization-sensitive optical coherence tomography | |
| Angelsky et al. | Polarization visualization and selection of biotissue image two-layer scattering medium | |
| Ushenko et al. | Laser diagnostics of anisotropy in birefringent networks of biological tissues in different physiological conditions | |
| US20210240134A1 (en) | Holographic reconstruction apparatus and method | |
| CN218157520U (zh) | 一种穆勒矩阵测量系统 | |
| WO2020108388A1 (zh) | 基于非瑞利散斑场的关联成像光谱相机及其成像方法 | |
| CN108931478A (zh) | 单次采集无色散移相全场光学相干层析成像装置及方法 | |
| Aspden et al. | Determination of the direction of preferred orientation and the orientation distribution function of collagen fibrils in connective tissues from high-angle X-ray diffraction patterns | |
| Peyvasteh et al. | Evolution of raw meat polarization‐based properties by means of Mueller matrix imaging | |
| Baroni et al. | Extending quantitative phase imaging to polarization-sensitive materials | |
| Zabolotna et al. | Methods and systems of polarization reproduction and analysis of the biological layers structure in the diagnosis of pathologies | |
| Sakamoto et al. | Exploring the influence of wavelength and polarization of illuminating light on polarization imaging for birefringent objects concealed behind scattering structure | |
| CN115200712A (zh) | 一种穆勒矩阵测量系统及测量方法 | |
| Liu et al. | Depth-resolved imaging of intercellular structures of large intestine with polarization-sensitive FFOCT | |
| CN113008427B (zh) | 一种液晶光弹性应力测量仪及测量方法 | |
| Zabolotna et al. | Orientational tomography of optical axes directions distributions of multilayer biological tissues birefringent polycrystalline networks | |
| CN118175277A (zh) | 基于空间光调制器的舌象形质三维结构检测仪 | |
| Šćepanović | Light scattering spectroscopy clinical imaging device implementation | |
| Fediv et al. | Polarization-phase diagnostics of latent course of cholelithiasis in patients with chronic cholecystitis combined with diabetes mellitus type 2 | |
| Abookasis et al. | Towards 3D mapping and correction of optical properties in turbid media based on spatially modulated illumination |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| GR01 | Patent grant | ||
| GR01 | Patent grant |