CN104463961B - 几何光学光线追迹仿真微粒衍射图方法 - Google Patents

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

Info

Publication number
CN104463961B
CN104463961B CN201410727208.3A CN201410727208A CN104463961B CN 104463961 B CN104463961 B CN 104463961B CN 201410727208 A CN201410727208 A CN 201410727208A CN 104463961 B CN104463961 B CN 104463961B
Authority
CN
China
Prior art keywords
light
diffraction
imaging
ray tracing
particulate
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
CN201410727208.3A
Other languages
English (en)
Other versions
CN104463961A (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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN201410727208.3A priority Critical patent/CN104463961B/zh
Publication of CN104463961A publication Critical patent/CN104463961A/zh
Application granted granted Critical
Publication of CN104463961B publication Critical patent/CN104463961B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/12—Symbolic schematics

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

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

Description

几何光学光线追迹仿真微粒衍射图方法
技术领域
本发明涉及一种微粒衍射成像的仿真方法。特别是涉及一种针对微米量级微粒，在相干光照明条件下采用对显微成像系统进行光线追迹获取衍射图像的几何光学光线追迹仿真微粒衍射图方法。
背景技术
X射线衍射测量是测定物质微观结构的一种常用方法，同样利用紫外、可见和红外波段的光电磁波也可对微米量级的微粒如细胞结构进行测定。单个细胞的尺寸恰好位于微米量级，处于这些波段的电磁波衍射测量范围内。因此出现了以相干光激发、以衍射成像方法测量相干散射光、最后确定细胞结构的方法，其中利用流式方法得到细胞队列，并对队列中的单个微粒进行衍射成像测量的方法称为衍射成像流式细胞术(diffraction imagingflow 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 microscopeobjective”,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 byflow 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 imagingand 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)获得远场散射光场：
采用时域有限差分法或离散偶极子近似方法计算获得的微粒散射米勒矩阵，由远场散射光强斯托克斯向量、入射光强斯托克斯向量及米勒矩阵之间的关系得到远场散射光强斯托克斯向量：
其中，为远场散射光强斯托克斯向量，为相干入射光强斯托克斯向量，为米勒矩阵；光强斯托克斯向量中I表示总光强，Q表示x轴方向线偏振光分量，U表示45度线偏振光分量，V表示右旋圆偏振光分量，r是位置矢量绝对值；k＝2π/λ，λ为光波长；
(2)将矩阵元S11投影到显微系统入射面x＝x0，所述显微系统入射面x＝x0上的点(y,z)的光强表示为：
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)获得像面衍射图
完成所有光线的追迹计算后，得到光线和成像面的相交点，相交点的分布密度图即是所仿真的衍射图。
米勒矩阵表示了微粒的衍射特性，在每个方位(θs,φs)，其中θs为仰角，φs为方位角，米勒矩阵中的S11可以表示为S11(θs,φs)，代表了单色平面波照明情况下(θs,φs)方向上的散射光强度。
步骤2)和3)还能够利用商用光学软件实现。
本发明的几何光学光线追迹仿真微粒衍射图方法，具有如下效果：
1、把微粒的远场散射场处理成几何光学所需的光源，利用光线追迹计算显微光学系统衍射成像，可以避免对结构复杂的显微光学成像系统的相干衍射光场计算，从而大大降低计算成本。
2、本发明采用商用光线追迹软件，可以对各种复杂的光学成像系统实现衍射图像计算，图形显示界面极为方便，易学易用。
附图说明
图1是本发明方法所使用的衍射成像流式细胞仪显微成像系统的结构示意图；
图2是微粒与衍射成像流式细胞仪显微成像系统入射面位置图：
图3是菲涅耳定律光线示意图；
图4是标准微球粒子在显微成像系统入射面的散射光强效果图；
图5是生成的光源光线效果图；
图6是微球光线追迹示意图；
图7a是仿真得到的CCD处微球衍射图
图7b是实测衍射图。
图中
1：微粒的球面波散射光 2：显微成像系统入射面位置
3：无穷远显微物镜 4：筒镜
5：CCD传感器 6：显微成像系统
7：流体室 8：流体(水)
9：微粒(细胞) 10：入射平面波激光
具体实施方式
下面结合实施例和附图对本发明的几何光学光线追迹仿真微粒衍射图方法做出详细说明。
本发明的几何光学光线追迹仿真微粒衍射图方法，是在衍射成像流式细胞仪显微成像系统上运行的。所述衍射成像流式细胞仪显微成像系统如图1所示，包括有：置于流体室(玻璃材质)7内的流体(水)8，置于流体室7内的流体(水)8内的微粒(细胞)9，流体室7外侧对应所述的微粒9设置有衍射成像流式细胞仪显微成像系统6，所述的衍射成像流式细胞仪显微成像系统6是由依次设置的无穷远显微物镜3、筒镜4和CCD传感器构成。其中的微粒或细胞9在流体室7中排列成单列，逐个经过观测区域，在平面波激光10(如图2)照射下发生散射，经衍射成像流式细胞仪显微成像系统6在CCD传感器5上形成衍射图像。
本发明的几何光学光线追迹仿真微粒衍射图方法，仿真系统的输入为衍射成像流式细胞仪显微成像系统入射平面的散射光场，输出为衍射成像流式细胞仪显微成像系统CCD传感器上所接收的衍射图像。具体包括如下步骤：
1)通过远场散射光场生成光源光线，包括：
(1)获得远场散射光场：
采用时域有限差分法(FDTD)或离散偶极子近似方法(DDA)计算获得的微粒散射米勒矩阵，由远场散射光强斯托克斯向量，入射光强斯托克斯向量及微粒米勒矩阵的关系得到远场散射光强斯托克斯向量：
其中，为远场散射光强斯托克斯向量，为入射光强斯托克斯向量，光强斯托克斯向量中I表示总光强，Q表示x轴方向线偏振光分量，U表示45度线偏振光分量，V表示右旋圆偏振光分量，r是位置矢量绝对值；k＝2π/λ，λ为光波长；
米勒矩阵表示了微粒的散射特性，在每个散射方向(θs,φs)(其中θs为仰角，φs为方位角)，米勒矩阵中的S11可以表示为S11(θs,φs)，代表了单色平面波照明情况下(θs,φs)方向上的散射光强度。
(2)将矩阵元S11投影到显微系统入射面x＝x0，所述显微系统入射面x＝x0上的点(y,z)的光强表示为：
x＝x0面上所有(y,z)点的光强I(y,z)就组成了显微系统入射面上的散射光场光强分布图；
(3)将衍射成像流式细胞仪显微成像系统入射面x＝x0上所有点的光强I(y,z)及方向(θs,φs)组成光线追迹的光源，点(y,z)发出的光线的数量N(y,z)正比于光强I(y,z)，方向为(θs,φs)；
(4)对所有通过衍射成像流式细胞仪显微成像系统成像的点(y,z)建立光线，即完成了光线追迹仿真系统光源光线的生成。
上述通过远场散射光场生成光源光线，可参考如下文献：K.Dong,Y.Feng,K.M.Jacobs,J.Q.Lu,R.S.Brock,L.V.Yang,F.E.Bertrand,M.A.Farwell,X.H.Hu,“Label-free classification of cultured cells through diffraction imaging”,BiomedicalOptics Express,2,1717-1726(2011)
R.S.Brock,X.H.Hu,P.Yang,J.Q.Lu,“Evaluation of a parallel FDTD codeand application to modeling of light scattering by deformed red blood cells”,Optics Express,13,5279-5292(2005)
2)追迹计算
对每一条生成的光线，依据菲涅耳折射定律：
n1sinI＝n2sinI' (3)
按照显微光学系统各光学镜片的折射特性，都获得光线的传播轨迹，当一条光线到达衍射成像流式细胞仪显微成像系统的成像平面，或者逃逸出衍射成像流式细胞仪显微成像系统，完成对这条光线的计算，直至所有光线计算完成，如图3所示，其中PO为入射光，OP'和OP″折射光和反射光，I'为对应的折射角，I为对应的入射角，n1为介质1的折射率，n2介质2的折射率；
本发明还能够利用商用光学软件进行追迹计算。
3)获得像面衍射图
完成所有光线的追迹计算后，得到光线和成像面的相交点，相交点的分布密度图即是所仿真的衍射图。本发明还能够利用商用光学软件获得像面衍射图。
下面以聚苯乙烯微球粒子为例说明本发明的几何光学光线追迹仿真微粒衍射图方法，仿真过程分为3个步骤：
1)通过远场散射光场生成光源光线
对于聚苯乙烯微球，具有简单的球对称结构，其远场的散射光强I(y,z)如图4所示，亦即显微成像系统入射面2上的散射光强图，把散射光强I(y,z)做为光线追迹的光源，对所有可能通过显微系统成像的点(y,z)建立光线如图5所示。
2)利用通用光学追迹软件计算追迹光线
这一步可通过商用软件Zemax进行，根据实际情况，选取适合的50X平场无穷远显微物镜和75mm筒镜，选定步骤1中生成的光源，进行光线追迹计算，光线轨迹如图6所示。
3)获得像面衍射图
这一步仍可通过商用软件Zemax进行，Zemax中有像面显示(Detector Viewer)功能，可显示获得的衍射图像，见附图7a，图7b是实测的衍射图，对比可见仿真结果非常相似。
以上结合附图对本发明的具体实施方式作了说明，但本发明的保护范围在于“把微粒的散射光强处理为几何光学可以处理的光线，通过光线追迹方法仿真显微光学系统对微粒衍射成像”。具体可表述为：
将微粒远场的散射光强I(y,z)，作为显微光学系统的入射面，将其强度分布转化为光线分布，作为光学追迹方法的光源，利用通用光学追迹方法仿真获得像面衍射图像。

Claims (2)

1.一种几何光学光线追迹仿真微粒衍射图方法，其特征在于，是在衍射成像流式细胞仪显微成像系统上运行，包括如下步骤：
1)通过远场散射光场生成光源光线，包括：
(1)获得远场散射光场：
采用时域有限差分法或离散偶极子近似方法计算获得的微粒散射米勒矩阵，由远场散射光强斯托克斯向量、入射光强斯托克斯向量及米勒矩阵之间的关系得到远场散射光强斯托克斯向量：
其中，为远场散射光强斯托克斯向量，为相干入射光强斯托克斯向量，为米勒矩阵；光强斯托克斯向量中I表示总光强，Q表示x轴方向线偏振光分量，U表示45度线偏振光分量，V表示右旋圆偏振光分量，r是位置矢量绝对值；
k＝2π/λ，λ为光波长；
米勒矩阵表示了微粒的衍射特性，在每个方位(θs,φs)，其中θs为仰角，φs为方位角，米勒矩阵中的S11可以表示为S11(θs,φs)，代表了单色平面波照明情况下(θs,φs)方向上的散射光强度；
(2)将矩阵元S11投影到显微系统入射面x＝x0，所述显微系统入射面x＝x0上的点(y,z)的光强表示为：
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)追迹计算
对每一条生成的光线，依据菲涅耳折射定律：
n1sin I＝n2sin I' (3)
按照显微光学系统各光学镜片的折射特性，都获得光线的传播轨迹，当一条光线到达衍射成像流式细胞仪显微成像系统的成像平面，或者逃逸出衍射成像流式细胞仪显微成像系统，完成对这条光线的计算，直至所有光线计算完成，式中，I′为对应的折射角，I为对应的入射角，n1为介质1的折射率，n2介质2的折射率；
3)获得像面衍射图
完成所有光线的追迹计算后，得到光线和成像面的相交点，相交点的分布密度图即是所仿真的衍射图。
2.根据权利要求1所述的几何光学光线追迹仿真微粒衍射图方法，其特征在于，步骤2)和3)还能够利用商用光学软件实现。
CN201410727208.3A 2014-12-03 2014-12-03 几何光学光线追迹仿真微粒衍射图方法 Active CN104463961B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410727208.3A CN104463961B (zh) 2014-12-03 2014-12-03 几何光学光线追迹仿真微粒衍射图方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410727208.3A CN104463961B (zh) 2014-12-03 2014-12-03 几何光学光线追迹仿真微粒衍射图方法

Publications (2)

Publication Number Publication Date
CN104463961A CN104463961A (zh) 2015-03-25
CN104463961B true CN104463961B (zh) 2017-06-16

Family

ID=52909943

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410727208.3A Active CN104463961B (zh) 2014-12-03 2014-12-03 几何光学光线追迹仿真微粒衍射图方法

Country Status (1)

Country Link
CN (1) CN104463961B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109615651B (zh) * 2019-01-29 2022-05-20 清华大学 基于光场显微系统的三维显微成像方法及系统
CN111965837B (zh) * 2020-07-30 2022-07-29 北京应用物理与计算数学研究所 用于光线衍射传输问题的求解方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101667136A (zh) * 2009-09-27 2010-03-10 北京航空航天大学 一种基于前向光线跟踪技术的星图模拟方法
CN103927421A (zh) * 2014-04-23 2014-07-16 清华大学 三维光学系统的实现方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8570531B2 (en) * 2011-12-11 2013-10-29 Tokyo Electron Limited Method of regenerating diffraction signals for optical metrology systems

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101667136A (zh) * 2009-09-27 2010-03-10 北京航空航天大学 一种基于前向光线跟踪技术的星图模拟方法
CN103927421A (zh) * 2014-04-23 2014-07-16 清华大学 三维光学系统的实现方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Analysis of cellular objects through diffraction images acquired by flow cytometry;Jun Zhang 等;《OPTICS EXPRESS》;20131009;第21卷(第21期);24819-24828 *
Label-free classification of cultured cells through diffraction imaging;Ke Dong 等;《BIOMEDICAL OPTICS EXPRESS》;20110601;第2卷(第6期);第2节，第3节 *
Simulations of light scattering from a biconcave red blood cell using the finite-difference time-domain method;Jun Q.Lu 等;《Journal of Biomedical Optics》;20050404;第10卷(第2期);024022-1-024022-10 *
一种基于光线追迹的遥感卫星光学系统成像仿真方法;张璇;《全国空间探测学术交流会》;20101231;摘要 *

Also Published As

Publication number Publication date
CN104463961A (zh) 2015-03-25

Similar Documents

Publication Publication Date Title
CN104374677B (zh) 一种粉尘浓度测量装置及方法
Bianco et al. Identification of microplastics based on the fractal properties of their holographic fingerprint
Wu et al. Holography and micro-holography of particle fields: a numerical standard
Skupsch et al. Multiple-plane particle image velocimetry using a light-field camera
Cerbino et al. Dark field differential dynamic microscopy enables accurate characterization of the roto-translational dynamics of bacteria and colloidal clusters
Tien et al. Color-coded three-dimensional micro particle tracking velocimetry and application to micro backward-facing step flows
Hannel et al. Holographic characterization of imperfect colloidal spheres
CN104463961B (zh) 几何光学光线追迹仿真微粒衍射图方法
Shanmugam et al. Lensless fluorescence imaging with height calculation
Gurtner et al. Twin-beam real-time position estimation of micro-objects in 3D
CN112666061B (zh) 一种基于无透镜成像系统光强模型的准球细胞测量方法
CN204177707U (zh) 一种大粒径范围的粉尘浓度测量装置
CN102735188B (zh) 一种测量球面曲率半径的方法
Miller et al. Observation of elliptically polarized light from total internal reflection in bubbles
Hillman et al. Microscopic particle discrimination using spatially-resolved Fourier-holographic light scattering angular spectroscopy
Li et al. Impact of aerosols on the polarization patterns of full-sky background radiation
CN108872152B (zh) 颗粒折射率测量方法、计算机装置及计算机可读存储介质
Gómez et al. Nonparaxial Mie theory of image formation in optical microscopes and characterization of colloidal particles
Wang et al. Measuring the three-dimensional volume scattering functions of microsphere suspension: Design and laboratory experiments
Li et al. Holographic imaging platform for particle discrimination based on simultaneous mass density and refractive index measurements
Zhan et al. Surface parameter based image estimation from application of a scattering model to polarized light measurements
Jin et al. An improved calculation model of weight coefficient for three-dimensional flame chemiluminescence tomography based on lens imaging theory
Yoo et al. Multimodal imaging Mueller polarimetric microscope to study polarimetric properties of spheroidal microparticles
Heinson Light scattering studies of irregularly shaped particles
Zhang et al. Characteristics analysis of infrared polarization for several typical artificial objects

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant