几何光学光线追迹仿真微粒衍射图方法技术

一种几何光学光线追迹仿真微粒衍射图方法：1)通过远场散射光场生成光源光线，具体有获得远场散射光场，将矩阵元S11投影到显微系统入射面x＝x0，得到散射光场光强分布图，组成光线追迹的光源；对所有通过衍射成像流式细胞仪显微成像系统成像的点(y,z)建立光线；2)依据菲涅耳折射定律对每一条生成的光线进行追迹计算；3)完成所有光线的追迹计算后，得到光线和成像面的相交点，相交点的分布密度图即是所仿真的衍射图。本发明专利技术把微粒的远场散射场处理成几何光学所需的光源，利用光线追迹计算显微光学系统衍射成像，可以避免对结构复杂的显微光学成像系统的相干衍射光场计算，可以对各种复杂的光学成像系统实现衍射图像计算，图形显示界面极为方便。

【技术实现步骤摘要】

本专利技术涉及一种微粒衍射成像的仿真方法。特别是涉及一种针对微米量级微粒，在相干光照明条件下采用对显微成像系统进行光线追迹获取衍射图像的几何光学光线追迹仿真微粒衍射图方法。
技术介绍
X射线衍射测量是测定物质微观结构的一种常用方法，同样利用紫外、可见和红外波段的光电磁波也可对微米量级的微粒如细胞结构进行测定。单个细胞的尺寸恰好位于微米量级，处于这些波段的电磁波衍射测量范围内。因此出现了以相干光激发、以衍射成像方法测量相干散射光、最后确定细胞结构的方法，其中利用流式方法得到细胞队列，并对队列中的单个微粒进行衍射成像测量的方法称为衍射成像流式细胞术(diffraction imaging flow cytometry,简写为DIFC)。该方法由液流系统形成微粒队列，由离焦显微光学系统拍摄微粒的衍射图像。该方法最近还增加了测量相干散射光的偏振状态的功能，因此也称为偏振衍射成像流式细胞术(polarized diffraction imaging flow cytometry,简写为p-DIFC)。参考文献：K.M.Jacobs,L.V.Yang,J.Ding,A.E.Ekpenyong,R.Castellone,J.Q.Lu,X.H.Hu,“Diffraction imaging of spheres and melanoma cells with a microscope objective”,Journal of Biophotonics,2,521-527(2009)；K.M.Jacobs,J.Q.Lu,X.H.Hu,“Development of a diffraction imaging flow cytometer”,Optics Letters,34,2985-2987(2009)；J.Zhang,Y.Feng,M.S.Moran,J.Q.Lu,L.V.Yang,Y.Sa,N.Zhang,L.Dong,X.H.Hu,\Analysis of cellular objects through diffraction images acquired by flow cytometry”,Optics Express,21,24819–24828(2013)；Y.Feng,N.Zhang,K.M.Jacobs,W.Jiang,L.V.Yang,Z.Li,J.Zhang,J.Q.Lu,X.H.Hu,\Polarization imaging and classification of Jurkat T and Ramos B cells using a flow cytometer”,Cytometry Part A,85,817-826(2014)。在DIFC技术中，获得的微粒衍射像多为斑点状，无法直观地获取信息，必须通过一定的算法提取细胞的形态特征。所提取的图像特征参数和真正的细胞结构的对应关系是一个非常关键的问题。通常的做法是：对已知结构的微粒如细胞进行数字建模，采用时域有限差分方法(Finite Difference Time Domain,简写为FDTD)或离散偶极子近似方法(Discrete Dipole Approximation，简写为DDA)等获得其相干散射光的空间分布和衍射图像，然后建立计算衍射图像和已知微粒结构之间的关联。然而，DIFC技术通常采用显微光学系统获取微粒衍射像，若直接应用FDTD或DDA等方法计算包含有显微光学元件的衍射光场，则会对系统建模带来很大挑战，可以预见计算成本极大提高。
技术实现思路
本专利技术所要解决的技术问题是，提供一种可以对各种复杂的光学成像系统实现衍射图像计算的几何光学光线追迹仿真微粒衍射图方法。本专利技术所采用的技术方案是：一种几何光学光线追迹仿真微粒衍射图方法，是在衍射成像流式细胞仪显微成像系统上运行，包括如下步骤：1)通过远场散射光场生成光源光线，包括：(1)获得远场散射光场：采用时域有限差分法或离散偶极子近似方法计算获得的微粒散射米勒矩阵，由远场散射光强斯托克斯向量、入射光强斯托克斯向量及米勒矩阵之间的关系得到远场散射光强斯托克斯向量：IsQsUsVs=1k2r2S11S12S13S14S21S22S23S24S31S32S33S34S41S42S43S44IiQiUiVi---(1)]]>其中，IsQsUsVs]]>为远场散射光强斯托克斯向量，IiQiUiVi]]>为相干入射光强斯托克斯向量，S11S12S13S14S21S22S23S24S31S32S33S34S41S42S43S44]]>为米勒矩阵；光强斯托克斯向量中I表示总光强，Q表示x轴方向线偏振光分量，U表示45度线偏振光分量，V表示右旋圆偏振光分量，r是位置矢量绝对值；k＝2π/λ，λ为光波长；(2)将矩阵元S11投影到显微系统入射面x＝x0，所述显微系统入射面x＝x0上的点(y,z)的光强表示为：I(y,z)=|cosφssinθs|x02(1+tan2φs+1+tan2φstan2θs)S11(θs,φs)---(2)]]>I(y,z)即显微系统入射面上的散射光场光强分布图；(3)将衍射成像流式细胞仪显微成像系统入射面x＝x0上所有点的光强I(y,z)及方向(θs,φs)组成光线追迹的光源，点(y,z)发出的光线的数量N(y,z)正比于光强I(y,z)，方向为(θs,φs)；(4)对所有通过衍射成像流式细胞仪显微成像系统成像的点(y,z)建立光线，即完成了光线追迹仿真系统光源光线的生成。2)追迹计算对每一条生成的光线，依据菲涅耳折射定律：n1sinI＝n2sinI'          (3)按照显微光学系统各光学镜片的折射特性，都获得光线的传播轨迹，当一条光线到达衍射成像流式细胞仪显微成像系统的成像平面，或者逃逸出衍射成像流式细胞仪显微成像系统，完成对这条光线的计算，直至所有光线计算完成，式中，I为对应的折射角，I为对应的入射角，n1为介质1的折射率，n2介质2的折射率；3)获得像面衍射图完成所有光线的追迹计算后，得到光线和成像面的相交点，相交点的分布密度图即是所仿真的衍射图。米勒矩阵S11S12S13S14S21S22S23S24S31S32S33S34S41S42S43S44]]>表示了微粒的衍射特性，在每个方位(θs,φs)，其中θs为仰角，φs为方位角，米勒矩阵中的S11可以表示为S11(θs,φs)，代表了单色平面波照明情况下(θs,φs)方向上的散射光强度。步骤2)和3)还能够利用本文档来自技高网...

【技术保护点】
一种几何光学光线追迹仿真微粒衍射图方法，其特征在于，是在衍射成像流式细胞仪显微成像系统上运行，包括如下步骤：1)通过远场散射光场生成光源光线，包括：(1)获得远场散射光场：采用时域有限差分法或离散偶极子近似方法计算获得的微粒散射米勒矩阵，由远场散射光强斯托克斯向量、入射光强斯托克斯向量及米勒矩阵之间的关系得到远场散射光强斯托克斯向量：IsQsUsVs=1k2r2S11S12S13S14S21S22S23S24S31S32S33S34S41S42S43S44IiQiUiVi--(1)]]>其中，IsQsUsVs]]>为远场散射光强斯托克斯向量，IiQiUiVi]]>为相干入射光强斯托克斯向量，S11S12S13S14S21S22S23S24S31S32S33S34S41S42S43S44]]>为米勒矩阵；光强斯托克斯向量中I表示总光强，Q表示x轴方向线偏振光分量，U表示45度线偏振光分量，V表示右旋圆偏振光分量，r是位置矢量绝对值；k＝2π/λ，λ为光波长；(2)将矩阵元S11投影到显微系统入射面x＝x0，所述显微系统入射面x＝x0上的点(y,z)的光强表示为：I(y,z)=|cosφssinθs|x02(1+tan2φs+1+tan2φstan2θs)S11(θs,φs)---(2)]]>I(y,z)即显微系统入射面上的散射光场光强分布图；(3)将衍射成像流式细胞仪显微成像系统入射面x＝x0上所有点的光强I(y,z)及方向(θs,φs)组成光线追迹的光源，点(y,z)发出的光线的数量N(y,z)正比于光强I(y,z)，方向为(θs,φs)；(4)对所有通过衍射成像流式细胞仪显微成像系统成像的点(y,z)建立光线，即完成了光线追迹仿真系统光源光线的生成。2)追迹计算对每一条生成的光线，依据菲涅耳折射定律：n1sinI＝n2sinI'    (3)按照显微光学系统各光学镜片的折射特性，都获得光线的传播轨迹，当一条光线到达衍射成像流式细胞仪显微成像系统的成像平面，或者逃逸出衍射成像流式细胞仪显微成像系统，完成对这条光线的计算，直至所有光线计算完成，式中，I为对应的折射角，I为对应的入射角，n1为介质1的折射率，n2介质2的折射率；3)获得像面衍射图完成所有光线的追迹计算后，得到光线和成像面的相交点，相交点的分布密度图即是所仿真的衍射图。...

【技术特征摘要】
1.一种几何光学光线追迹仿真微粒衍射图方法，其特征在于，是在衍射成像流式细胞
仪显微成像系统上运行，包括如下步骤：
1)通过远场散射光场生成光源光线，包括：
(1)获得远场散射光场：
采用时域有限差分法或离散偶极子近似方法计算获得的微粒散射米勒矩阵，由远场散射
光强斯托克斯向量、入射光强斯托克斯向量及米勒矩阵之间的关系得到远场散射光强斯托克
斯向量：
IsQsUsVs=1k2r2S11S12S13S14S21S22S23S24S31S32S33S34S41S42S43S44IiQiUiVi--(1)]]>其中，IsQsUsVs]]>为远场散射光强斯托克斯向量，IiQiUiVi]]>为相干入射光强斯托克斯向量，
S11S12S13S14S21S22S23S24S31S32S33S34S41S42S43S44]]>为米勒矩阵；光强斯托克斯向量中I表示总光强，Q表示x轴方向线偏
振光分量，U表示45度线偏振光分量，V表示右旋圆偏振光分量，r是位置矢量绝对值；
k＝2π/λ，λ为光波长；
(2)将矩阵元S11投影到显微系统入射面x＝x0，所述显微系统入射面x＝x0上的点(y,z)
的光强表示为：
I(y,z)=|cosφssinθs|x02(1+tan2φs+1+tan2φstan2θs)S11(θs,φs)---(2)]]>I(y,z)即显微系统入射面上的散射光场光强分布图；
(3)将衍射成像流式细胞仪显微成像系统入射面x＝x0上所有点的...

【专利技术属性】
技术研发人员：撒昱，冯远明，胡新华，
申请(专利权)人：天津大学，
类型：发明
国别省市：天津;12