发明名称 几何光学光线追迹仿真微粒衍射图方法
摘要 一种几何光学光线追迹仿真微粒衍射图方法:1)通过远场散射光场生成光源光线,具体有获得远场散射光场,将矩阵元S<sub>11</sub>投影到显微系统入射面x=x<sub>0</sub>,得到散射光场光强分布图,组成光线追迹的光源;对所有通过衍射成像流式细胞仪显微成像系统成像的点(y,z)建立光线;2)依据菲涅耳折射定律对每一条生成的光线进行追迹计算;3)完成所有光线的追迹计算后,得到光线和成像面的相交点,相交点的分布密度图即是所仿真的衍射图。本发明把微粒的远场散射场处理成几何光学所需的光源,利用光线追迹计算显微光学系统衍射成像,可以避免对结构复杂的显微光学成像系统的相干衍射光场计算,可以对各种复杂的光学成像系统实现衍射图像计算,图形显示界面极为方便。
申请公布号 CN104463961A 申请公布日期 2015.03.25
申请号 CN201410727208.3 申请日期 2014.12.03
申请人 天津大学 发明人 撒昱;冯远明;胡新华
分类号 G06T17/00(2006.01)I;G06F17/50(2006.01)I 主分类号 G06T17/00(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 杜文茹
主权项 一种几何光学光线追迹仿真微粒衍射图方法,其特征在于,是在衍射成像流式细胞仪显微成像系统上运行,包括如下步骤:1)通过远场散射光场生成光源光线,包括:(1)获得远场散射光场:采用时域有限差分法或离散偶极子近似方法计算获得的微粒散射米勒矩阵,由远场散射光强斯托克斯向量、入射光强斯托克斯向量及米勒矩阵之间的关系得到远场散射光强斯托克斯向量:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>I</mi><mi>s</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mi>s</mi></msub></mtd></mtr><mtr><mtd><msub><mi>U</mi><mi>s</mi></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mi>s</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfrac><mn>1</mn><mrow><msup><mi>k</mi><mn>2</mn></msup><msup><mi>r</mi><mn>2</mn></msup></mrow></mfrac><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>S</mi><mn>11</mn></msub></mtd><mtd><msub><mi>S</mi><mn>12</mn></msub></mtd><mtd><msub><mi>S</mi><mn>13</mn></msub></mtd><mtd><msub><mi>S</mi><mn>14</mn></msub></mtd></mtr><mtr><mtd><msub><mi>S</mi><mn>21</mn></msub></mtd><mtd><msub><mi>S</mi><mn>22</mn></msub></mtd><mtd><msub><mi>S</mi><mn>23</mn></msub></mtd><mtd><msub><mi>S</mi><mn>24</mn></msub></mtd></mtr><mtr><mtd><msub><mi>S</mi><mn>31</mn></msub></mtd><mtd><msub><mi>S</mi><mn>32</mn></msub></mtd><mtd><msub><mi>S</mi><mn>33</mn></msub></mtd><mtd><msub><mi>S</mi><mn>34</mn></msub></mtd></mtr><mtr><mtd><msub><mi>S</mi><mn>41</mn></msub></mtd><mtd><msub><mi>S</mi><mn>42</mn></msub></mtd><mtd><msub><mi>S</mi><mn>43</mn></msub></mtd><mtd><msub><mi>S</mi><mn>44</mn></msub></mtd></mtr></mtable></mfenced><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>I</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>U</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000624039250000011.GIF" wi="1317" he="349" /></maths>其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>I</mi><mi>s</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mi>s</mi></msub></mtd></mtr><mtr><mtd><msub><mi>U</mi><mi>s</mi></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mi>s</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000624039250000012.GIF" wi="126" he="312" /></maths>为远场散射光强斯托克斯向量,<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>I</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>U</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000624039250000013.GIF" wi="122" he="315" /></maths>为相干入射光强斯托克斯向量,<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>S</mi><mn>11</mn></msub></mtd><mtd><msub><mi>S</mi><mn>12</mn></msub></mtd><mtd><msub><mi>S</mi><mn>13</mn></msub></mtd><mtd><msub><mi>S</mi><mn>14</mn></msub></mtd></mtr><mtr><mtd><msub><mi>S</mi><mn>21</mn></msub></mtd><mtd><msub><mi>S</mi><mn>22</mn></msub></mtd><mtd><msub><mi>S</mi><mn>23</mn></msub></mtd><mtd><msub><mi>S</mi><mn>24</mn></msub></mtd></mtr><mtr><mtd><msub><mi>S</mi><mn>31</mn></msub></mtd><mtd><msub><mi>S</mi><mn>32</mn></msub></mtd><mtd><msub><mi>S</mi><mn>33</mn></msub></mtd><mtd><msub><mi>S</mi><mn>34</mn></msub></mtd></mtr><mtr><mtd><msub><mi>S</mi><mn>41</mn></msub></mtd><mtd><msub><mi>S</mi><mn>42</mn></msub></mtd><mtd><msub><mi>S</mi><mn>43</mn></msub></mtd><mtd><msub><mi>S</mi><mn>44</mn></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000624039250000014.GIF" wi="454" he="309" /></maths>为米勒矩阵;光强斯托克斯向量中I表示总光强,Q表示x轴方向线偏振光分量,U表示45度线偏振光分量,V表示右旋圆偏振光分量,r是位置矢量绝对值;k=2π/λ,λ为光波长;(2)将矩阵元S<sub>11</sub>投影到显微系统入射面x=x<sub>0</sub>,所述显微系统入射面x=x<sub>0</sub>上的点(y,z)的光强表示为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>I</mi><mrow><mo>(</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>|</mo><mi>cos</mi><msub><mi>&phi;</mi><mi>s</mi></msub><mi>sin</mi><msub><mi>&theta;</mi><mi>s</mi></msub><mo>|</mo></mrow><mrow><msup><msub><mi>x</mi><mn>0</mn></msub><mn>2</mn></msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>tan</mi><mn>2</mn></msup><msub><mi>&phi;</mi><mi>s</mi></msub><mo>+</mo><mfrac><mrow><mn>1</mn><mo>+</mo><msup><mi>tan</mi><mn>2</mn></msup><msub><mi>&phi;</mi><mi>s</mi></msub></mrow><mrow><msup><mi>tan</mi><mn>2</mn></msup><msub><mi>&theta;</mi><mi>s</mi></msub></mrow></mfrac><mo>)</mo></mrow></mrow></mfrac><msub><mi>S</mi><mn>11</mn></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>s</mi></msub><mo>,</mo><msub><mi>&phi;</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000624039250000015.GIF" wi="1144" he="193" /></maths>I(y,z)即显微系统入射面上的散射光场光强分布图;(3)将衍射成像流式细胞仪显微成像系统入射面x=x<sub>0</sub>上所有点的光强I(y,z)及方向(θ<sub>s</sub>,φ<sub>s</sub>)组成光线追迹的光源,点(y,z)发出的光线的数量N(y,z)正比于光强I(y,z),方向为(θ<sub>s</sub>,φ<sub>s</sub>);(4)对所有通过衍射成像流式细胞仪显微成像系统成像的点(y,z)建立光线,即完成了光线追迹仿真系统光源光线的生成。2)追迹计算对每一条生成的光线,依据菲涅耳折射定律:n<sub>1</sub>sinI=n<sub>2</sub>sinI'    (3)按照显微光学系统各光学镜片的折射特性,都获得光线的传播轨迹,当一条光线到达衍射成像流式细胞仪显微成像系统的成像平面,或者逃逸出衍射成像流式细胞仪显微成像系统,完成对这条光线的计算,直至所有光线计算完成,式中,I为对应的折射角,I为对应的入射角,n<sub>1</sub>为介质1的折射率,n<sub>2</sub>介质2的折射率;3)获得像面衍射图完成所有光线的追迹计算后,得到光线和成像面的相交点,相交点的分布密度图即是所仿真的衍射图。
地址 300072 天津市南开区卫津路92号