高精度激光定位散斑机理分析及抑制 Mechanism analysis and suppression of speckle in high-precision laser location 1 3 first-author 徐 晶（1977-），女，黑龙江宾县人，博士，副教授，2003年、2009年于长春理工大学分别获得硕士和博士学位，主要从事数字图像处理、机器视觉等方面的研究。E-mail： xujing@cust.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714762&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714781&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714778&type= 22.01332855 32.00401306 徐 晶 （1977-），女，黑龙江宾县人，博士，副教授，2003年、2009年于长春理工大学分别获得硕士和博士学位，主要从事数字图像处理、机器视觉等方面的研究。E-mail： xujing@cust.edu.cn xujing@cust.edu.cn 1 corresp E-mail： wqcaaa2022@163.com E-mail： wqcaaa2022@163.com 王钎潮（1998-），男，湖北黄冈人，硕士研究生，主要从事机器视觉方面的研究。E-mail：wqcaaa2022@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714793&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714796&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714795&type= 22.01332855 32.00399780 王钎潮 （1998-），男，湖北黄冈人，硕士研究生，主要从事机器视觉方面的研究。E-mail： wqcaaa2022@163.com wqcaaa2022@163.com 2 2 2 长春理工大学 计算机科学技术学院，吉林 长春 130022 长春理工大学 计算机科学技术学院，吉林 长春 130022 School of Computer Science and Technology， Changchun University of Science and Technology， Changchun 130022， China School of Computer Science and Technology， Changchun University of Science and Technology， Changchun 130022， China 长春理工大学 人工智能学院，吉林 长春 130022 长春理工大学 人工智能学院，吉林 长春 130022 School of Artificial Intelligence， Changchun University of Science and Technology， Changchun 130022， China School of Artificial Intelligence， Changchun University of Science and Technology， Changchun 130022， China 长春理工大学 中山研究院，广东 中山 528403 长春理工大学 中山研究院，广东 中山 528403 Zhongshan Institute of Changchun University of Science and Technology， Zhongshan 528403， China Zhongshan Institute of Changchun University of Science and Technology， Zhongshan 528403， China 为有效解决高精度视觉定位测量系统中激光随机散斑导致的光斑质心定位精度降低问题，对散斑机理及抑制方法进行研究。从激光光斑的成像过程分析了激光随机散斑出现的原因，并推导出空间中任意一点总波扰表达式。建立了成像介质吸光性、粗糙度和表面出射光信噪比与随机散斑强度的关系模型，得出增加成像介质粗糙度和吸光性来提高出射光信噪比，进而降低随机散斑强度的结论。通过改变成像介质的物理特性进行对比实验，获得光斑质心位置稳定性随成像介质物理特性变化的趋势，验证了关系模型的正确性。实验结果表明，当激光功率为50 mW时，选取吸光性强、粗糙度与厚度适中的成像介质作为成像靶面，激光光斑质心的提取精度稳定在0.04个pixel之内，实现了对激光散斑的抑制，提高了光斑质心的定位精度。 To solve the problem of the reduction in the position accuracy of the laser spot centroid caused by laser random speckle in the high-precision visual positioning measurement system， the speckle mechanism and suppression methods were studied. First， the causes of laser random speckle in the laser spot imaging process were analyzed， and the expression for the total wave disturbance at any point in space was derived. Then， a model for the relationships among the light absorption， roughness of the imaging medium， surface outgoing light signal-to-noise ratio， and random speckle intensity was established. It is concluded that increasing the roughness and light absorption of the imaging medium can increase the signal-to-noise ratio of the outgoing light and thus reduce the random speckle intensity. Finally， comparison experiments were performed to obtain the trend chart of the stability of the spot centroid position with changes in the physical characteristics of the imaging medium， and the results validated the relationship model. The experimental results indicate that when the laser power is 50 mW， the imaging medium had strong light absorption， and the imaging target surface had moderate roughness and thickness， the stability of laser spot centroid extraction is within 0.04 pixels. Laser speckle is suppressed， the stability of imaging spot coordinate extraction is improved， and the accuracy of spot centroid positioning is increased. 激光光斑 动态散斑 光的散射 分子运动论 微表面模型 laser spot dynamic speckle scattering of light molecular motion theory microfacet model 1 引 言 1 激光具有良好的单一方向性、高能量性和高相干性，常作为高精度视觉定位测量系统中的成像载体 ［ 1 1 - 2 2 1-2 ］ 。激光照射成像介质时，成像介质内分子接收到激光能量后会发生受激运动，导致激光投影成像时存在随机复合散斑 ［ 3 3 - 4 4 3-4 ］ ，这会影响成像光斑的稳定性，进而增加光斑质心的定位误差。因此，抑制随机散斑是从根源上提高视觉定位测量系统中光斑质心提取精度的关键。 近年来，许多学者针对散斑抑制方法开展了研究。王锐 ［ 5 5 ］ 等提出建立多光束分光照明系统对激光散斑效应进行抑制。张圣涛 ［ 6 6 ］ 使用多个互不相干光源，利用半导体激光二极管阵列的合理排布，引进角度多样性降低激光散斑。这类方法需要改进系统光源的功率或数量，会导致光斑过饱和，增加光斑质心的提取难度。梁传样 ［ 7 7 ］ 等使用衍射光学元件实现激光整形和散斑抑制。该方法改变了激光成像的高斯型，很难完成光斑质心的高精度提取。Deng ［ 8 8 ］ 等采用旋转球透镜来抑制激光投影散斑。付瀚毅 ［ 9 9 ］ 提出转动随机透镜阵列抑制散斑的结构。孙鸣捷 ［ 10 10 ］ 等提出对激光照明系统中混光棒施加连续作用力，使混光棒连续形变从而实现激光散斑抑制。这种方法虽抑制了散斑，但也会让成像光斑位置发生物理偏移，增大成像光斑质心的提取误差。上述方法虽能抑制散斑，但在视觉定位测量系统中，光源的光强与波长、成像靶面的位置都需动态改变。这种动态改变会降低光斑质心的提取精度，所以不适用于以光斑质心定位为基础的视觉定位测量系统。鉴于此，本文主要针对视觉定位测量中系统保持相对静止时的随机散斑进行抑制，提高成像光斑的稳定性及光斑质心的定位精度。 本文理论分析了视觉定位测量系统中随机散斑产生的原因，推导出空间中任意一点总波扰表达式。定量分析成像介质密度、粗糙度、吸光性、激光功率与成像光斑图像信噪比之间的关系，并在此基础上设计验证实验。提取不同物理特性成像介质上的光斑质心，获得光斑质心位置稳定性随成像介质物理特性变化的趋势，验证理论分析的正确性。 2 原 理 1 2.1　随机激光散斑的成因 2 图1 图1 是以激光为载体的高精度视觉定位测量系统。理论上，当相机和成像靶面相对静止时，成像系统记录的激光光斑应处于较稳定的状态。但实际使用这类视觉定位测量系统时，在相机、成像介质和激光器位置保持不变时，相机获得的光斑成像图稳定性却较差，严重影响光斑质心位置的提取精度。 稳定激光器输出的激光模式为高斯光束，其光场分布可用麦克斯韦方程组来表示 ［ 11 11 ］ 。结合麦克斯韦方程组与介质方程，得到自由空间稳态传输的电磁场所满足的亥姆霍兹方程： <math display="block" id="M1"><mfenced separators="|"><mrow><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mo>∂</mo><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mo>∂</mo><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mo>∂</mo><msup><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>+</mo><msup><mrow><mi>k</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced><mi>E</mi><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714805&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714804&type= 51.81600189 10.41399956 ， 式中： <math id="M2"><mfenced separators="|"><mrow><msub><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>,</mo><msub><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714826&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714823&type= 16.25600052 4.82600021 分别为 <math id="M3"><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714817&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714829&type= 9.31333351 3.72533321 坐标轴的单位矢量， <math id="M4"><mi>k</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714835&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714834&type= 1.69333339 2.79399991 是波数， <math id="M5"><mi>E</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714854&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714839&type= 2.70933342 2.79399991 是振幅。在近轴条件下，激光器产生的高斯光束是方程（1）的一个特解，可表示为： <math display="block" id="M6"><mi>E</mi><mo>=</mo><msub><mrow><mi>E</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>Ψ</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi></mrow></mfenced><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced separators="|"><mrow><mo>-</mo><mi mathvariant="normal">i</mi><mi>k</mi><mi>z</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714859&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714858&type= 45.21200180 4.82600021 ， 式中： <math id="M7"><msub><mrow><mi>E</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714847&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714845&type= 3.80999994 3.72533321 为最大幅值， <math id="M8"><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714867&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714849&type= 12.53066730 4.82600021 为空间中任意一点， <math id="M9"><mi>Ψ</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714883&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714869&type= 16.59466553 4.82600021 可以看成是一个沿 <math id="M10"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714887&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714872&type= 1.77800000 2.79399991 轴缓慢变化的复函数， <math id="M11"><mi mathvariant="normal">i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714876&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714891&type= 0.93133330 2.79399991 为虚部。将 式（2） 式（2） 代入 式（1） 式（1） ，忽略 <math id="M12"><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>Ψ</mi></mrow><mrow><mo>∂</mo><msup><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714900&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714897&type= 6.68866634 8.38199997 项得： <math display="block" id="M13"><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>Ψ</mi></mrow><mrow><mo>∂</mo><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>Ψ</mi></mrow><mrow><mo>∂</mo><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>-</mo><mi mathvariant="normal">i</mi><mi>k</mi><mi>z</mi><mfrac><mrow><mo>∂</mo><mi>Ψ</mi></mrow><mrow><mo>∂</mo><mi>z</mi></mrow></mfrac><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714914&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714881&type= 40.13199997 9.31333351 . 式（3） 式（3） 为抛物线型微分方程，用试解法求解，最后得到高斯光束的表达式： <math display="block" id="M14"><mtable columnalign="left"><mtr><mtd><mi>E</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi></mrow></mfenced><mo>=</mo><mfrac><mrow><msub><mrow><mi>E</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mi>w</mi><mfenced separators="|"><mrow><mi>z</mi></mrow></mfenced></mrow></mfrac><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced open="[" close="]" separators="|"><mrow><mo>-</mo><mfrac><mrow><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><msup><mrow><mi>w</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mfenced separators="|"><mrow><mi>z</mi></mrow></mfenced></mrow></mfrac></mrow></mfenced><mo>⋅</mo></mtd></mtr><mtr><mtd><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced open="[" close="]" separators="|"><mrow><mo>-</mo><mi mathvariant="normal">i</mi><mi>k</mi><mfenced separators="|"><mrow><mi>z</mi><mo>+</mo><mfrac><mrow><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mn mathvariant="normal">2</mn><mi>R</mi><mo stretchy="false">(</mo><mi>z</mi><mo stretchy="false">)</mo></mrow></mfrac><mo>+</mo><mi mathvariant="normal">i</mi><mi>φ</mi><mfenced separators="|"><mrow><mi>z</mi></mrow></mfenced></mrow></mfenced></mrow></mfenced><mo>,</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714907&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714916&type= 56.30332947 22.77533340 式中： <math id="M15"><mi>R</mi><mo stretchy="false">(</mo><mi>z</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714911&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714909&type= 8.12800026 3.55599999 表示高斯光束的等相面曲率半径， <math id="M16"><mi>w</mi><mfenced separators="|"><mrow><mi>z</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714948&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714928&type= 8.55133343 3.97933316 表示光斑半径， <math id="M17"><mi>φ</mi><mfenced separators="|"><mrow><mi>z</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714952&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714930&type= 7.87400007 4.06400013 表示附加相移。 式（4） 式（4） 含振幅和相位两部分，它反映高斯光束场分布及光束在传播过程中的变化规律。 将 式（4） 式（4） 代入惠更斯-菲涅尔公式中，则空间中位置矢量为 <math id="M18"><msup><mrow><mi mathvariant="bold-italic">r</mi></mrow><mrow><mi>'</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714989&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714980&type= 2.28600001 3.04800010 的散射微元所产生的散射子波扰为： <math display="block" id="M19"><mi>μ</mi><mfenced separators="|"><mrow><msup><mrow><mi mathvariant="bold-italic">r</mi></mrow><mrow><mi>'</mi></mrow></msup></mrow></mfenced><mo>=</mo><mi>E</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi></mrow></mfenced><mfrac><mrow><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced separators="|"><mrow><mi mathvariant="normal">i</mi><mi>k</mi><msup><mrow><mi mathvariant="bold-italic">r</mi></mrow><mrow><mi>'</mi></mrow></msup></mrow></mfenced></mrow><mrow><msup><mrow><mi mathvariant="bold-italic">r</mi></mrow><mrow><mi>'</mi></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714957&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714956&type= 44.19599915 8.29733276 . 当空间中子波叠加时，空间中任意一点 <math id="M20"><mi>P</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715006&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715000&type= 12.19200039 3.97933316 总波扰可表示为： <math display="block" id="M21"><mi>ζ</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced><mo>=</mo><mfrac><mrow><mi mathvariant="normal">i</mi><mi>μ</mi><mfenced separators="|"><mrow><msup><mrow><mi mathvariant="bold-italic">r</mi></mrow><mrow><mi>'</mi></mrow></msup></mrow></mfenced></mrow><mrow><mi>λ</mi></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∬</mo><mrow><mi>S</mi></mrow><mrow/></munderover></mstyle><mrow><mfrac><mrow><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced separators="|"><mrow><mi mathvariant="normal">i</mi><mi>k</mi><mi mathvariant="bold-italic">R</mi></mrow></mfenced></mrow><mrow><mi mathvariant="bold-italic">R</mi></mrow></mfrac><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mi>K</mi><mfenced separators="|"><mrow><mi>χ</mi></mrow></mfenced></mrow><mi mathvariant="normal">d</mi><mi>S</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714965&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714962&type= 62.56866837 9.65200043 ， 式中： <math id="M22"><mi>K</mi><mfenced separators="|"><mrow><mi>χ</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714970&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714968&type= 8.97466660 4.74133301 为倾斜因子，主要用来调整子波朝着各个方向传播波幅； <math id="M23"><msup><mrow><mi mathvariant="bold-italic">r</mi></mrow><mrow><mi>'</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714989&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714980&type= 2.28600001 3.04800010 为空间中某位置矢量； <math id="M24"><mi mathvariant="bold-italic">R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714985&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45714983&type= 2.87866688 2.79399991 为子波源 <math id="M25"><mi>μ</mi><mfenced separators="|"><mrow><msup><mrow><mi mathvariant="bold-italic">r</mi></mrow><mrow><mi>'</mi></mrow></msup></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715043&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715053&type= 8.21266651 3.97933316 到 <math id="M26"><mi>P</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715006&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715000&type= 12.19200039 3.97933316 点的位置矢量； <math id="M27"><mi>S</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715050&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715030&type= 2.20133328 2.79399991 是散射积分曲面； <math id="M28"><mi>ζ</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715060&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715045&type= 11.59933376 4.82600021 为某位置矢量总波扰。 图2 图2 为折射光在成像介质内部碰撞的过程。由分子运动论、光的折射定律和光的波粒二象性定理可知，用激光照射成像介质，折射光会进入到成像介质内部，内部部分分子吸收其他分子传递的能量后，不只在平衡位置振动，从而导致成像介质内部分子始终在随机运动 ［ 12 12 - 13 13 12-13 ］ 。因此，散射面 <math id="M29"><mi>S</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715050&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715030&type= 2.20133328 2.79399991 所产生的散射子波在方向和强度上随机。由 式（6） 式（6） 可知，如果散射面 <math id="M30"><mi>S</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715050&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715030&type= 2.20133328 2.79399991 射出的子波 <math id="M31"><mi>μ</mi><mfenced separators="|"><mrow><msup><mrow><mi mathvariant="bold-italic">r</mi></mrow><mrow><mi>'</mi></mrow></msup></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715043&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715053&type= 8.21266651 3.97933316 在强度和方向上发生随机变化，那么总波扰 <math id="M32"><mi>ζ</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715060&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715045&type= 11.59933376 4.82600021 的强度也会发生变化。这些随机子波在空间中相互干涉，经过光学系统后，在像面上形成许多随机分布的明暗动态变化图样，出现随机激光散斑现象。 2.2　成像介质表面出射光信噪比分析 2 由2.1节分析可知，成像介质表面随机散射光导致随机激光散斑的出现。因反射光不会进入成像介质内部，所以在不考虑空气扰动的情况下，反射光本身是稳定的。以激光为光源，从成像介质表面进入到空间中的反射光 <math id="M33"><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715096&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715093&type= 3.30200005 3.72533321 和散射光 <math id="M34"><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715266&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715252&type= 3.47133350 3.72533321 可表示为： <math display="block" id="M35"><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub><mo>=</mo><msub><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><msub><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715105&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715112&type= 17.69533348 3.80999994 ， <math display="block" id="M36"><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub><mo>=</mo><msub><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>+</mo><msub><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715118&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715116&type= 17.78000069 3.80999994 ， 式中： <math id="M37"><msub><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715303&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715293&type= 2.37066650 3.72533321 表示纯反射光， <math id="M38"><msub><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715161&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715157&type= 2.37066650 3.72533321 表示反射光中与散射光发生干涉的部分， <math id="M39"><msub><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715153&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715150&type= 2.37066650 3.72533321 表示纯散射光， <math id="M40"><msub><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715161&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715157&type= 2.37066650 3.72533321 和 <math id="M41"><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715612&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715611&type= 3.47133350 3.72533321 成正比。可见光斑的不稳定性主要来自散射光自干涉，以及散射光与反射光的交叉干涉。 图3 图3 表示了3种光的分布情况，3部分扇形内线条越凌乱说明其稳定性越差。 因相机视场范围有限，并非所有从成像介质表面出来的光都能进入到相机中。根据 式（7） 式（7） 和 式（8） 式（8） ，像面上记录的能量可表示为： <math display="block" id="M42"><msub><mrow><mi>M</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">U</mi></mrow></msub><mo>=</mo><munder><mo movablelimits="false">∑</mo><mrow><mi>Ω</mi></mrow></munder><mrow><mfenced separators="|"><mrow><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub><mo>+</mo><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub></mrow></mfenced></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715195&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715206&type= 30.56466866 6.60400009 ， 式中： <math id="M43"><mi>Ω</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715230&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715213&type= 2.70933342 2.79399991 为视场范围， <math id="M44"><msub><mrow><mi>M</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">U</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715626&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715624&type= 6.51933336 3.72533321 是相机像面中记录的激光总信号。 <math id="M45"><msub><mrow><mi>M</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">U</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715626&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715624&type= 6.51933336 3.72533321 决定最终像面中光斑的整体亮度。若将 <math id="M46"><msub><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715303&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715293&type= 2.37066650 3.72533321 定义为稳定的信号，将 <math id="M47"><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715266&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715252&type= 3.47133350 3.72533321 定义为噪声，则信噪比 <math id="M48"><msub><mrow><mi>R</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">N</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715275&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715270&type= 5.58799982 3.72533321 可表示为： <math display="block" id="M49"><msub><mrow><mi>R</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">N</mi></mrow></msub><mo>=</mo><munder><mo movablelimits="false">∑</mo><mrow><mi>Ω</mi></mrow></munder><mrow><mfrac><mrow><msub><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub></mrow></mfrac></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715291&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715289&type= 19.13466644 8.63599968 . 由 式（10） 式（10） 可知，为提高信噪比，可以增大 <math id="M50"><msub><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715303&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715293&type= 2.37066650 3.72533321 ，但这会使 <math id="M51"><msub><mrow><mi>M</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">U</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715626&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715624&type= 6.51933336 3.72533321 值过大，成像光斑过亮，不利于提取质心。因此，只能在减小 <math id="M52"><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715674&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715684&type= 3.47133350 3.72533321 的基础上提高信噪比。 2.3　成像介质粗糙度对光照的影响 2 根据微平面理论，在微观尺度下，任何表面都是由若干个细小平面组成 ［ 14 14 ］ 。微表面粗糙程度 <math id="M53"><mi>α</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715548&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715524&type= 1.94733346 2.79399991 可表示为： <math display="block" id="M54"><mi>α</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>m</mi></mrow></mfrac><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>m</mi></mrow></munderover></mstyle><mrow><mfenced open="|" close="|" separators="|"><mrow><msub><mrow><mi>y</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfenced></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715356&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715352&type= 20.82799911 8.63599968 ， 式中： <math id="M55"><mi>m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715364&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715340&type= 2.87866688 2.79399991 为采样点个数， <math id="M56"><msub><mrow><mi>y</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715379&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715345&type= 2.87866688 4.57200003 为采样点与微平面平均高度的差值。 由双向反射分布函数（Bidirectional Reflected Distribution Function， BRDF）的定义可知，只有当微平面法线 <math id="M57"><mi mathvariant="bold-italic">n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715384&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715371&type= 2.11666679 2.79399991 与 <math id="M58"><mi mathvariant="bold-italic">l</mi><mo>⋅</mo><mi mathvariant="bold-italic">v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715390&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715376&type= 5.67266655 2.79399991 的半程向量 <math id="M59"><mi mathvariant="bold-italic">h</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715401&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715409&type= 2.20133328 2.79399991 方向一致时 ［ 15 15 ］ ，相机才能接收到微平面反射光。如 图4 图4 所示（彩图见期刊电子版），虽然 A ， B 两点都能被相机和光照到，但只有 B 点的光能反射到相机中。此外，微表面法线分布的不规则性也会导致散射光从成像介质内部出来的方向更加随机。 Trowbridge等在 式（11） 式（11） 基础上，使用高斯分布来近似统计微表面的法线分布与半程向量 <math id="M60"><mi mathvariant="bold-italic">h</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715482&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715493&type= 2.20133328 2.79399991 重合的比率 ［ 16 16 ］ ，即： <math display="block" id="M61"><mi>D</mi><mfenced separators="|"><mrow><mi mathvariant="bold-italic">n</mi><mo>,</mo><mi mathvariant="bold-italic">h</mi><mo>,</mo><mi>α</mi></mrow></mfenced><mo>=</mo><mfrac><mrow><msup><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mi>π</mi><msup><mrow><mfenced separators="|"><mrow><msup><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">n</mi><mo>⋅</mo><mi mathvariant="bold-italic">h</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mfenced separators="|"><mrow><msup><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced><mo>+</mo><mn mathvariant="normal">1</mn></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715431&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715429&type= 59.94400024 12.19200039 ， <math display="block" id="M62"><mi mathvariant="bold-italic">h</mi><mo>=</mo><mfrac><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">l</mi><mo>+</mo><mi mathvariant="bold-italic">v</mi></mrow></mfenced></mrow><mrow><mfenced open="|" close="|" separators="|"><mrow><mi mathvariant="bold-italic">l</mi><mo>+</mo><mi mathvariant="bold-italic">v</mi></mrow></mfenced></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715436&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715450&type= 18.96533203 9.56733322 ， 式中： <math id="M63"><mi>α</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715548&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715524&type= 1.94733346 2.79399991 为粗糙度， <math id="M64"><mi mathvariant="bold-italic">n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715538&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715535&type= 2.11666679 2.79399991 为微平面法线， <math id="M65"><mi mathvariant="bold-italic">l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715557&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715556&type= 1.35466671 2.79399991 为光位置矢量， <math id="M66"><mi mathvariant="bold-italic">v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715577&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715575&type= 2.03200006 2.79399991 为视线位置矢量， <math id="M67"><mi mathvariant="bold-italic">h</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715482&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715493&type= 2.20133328 2.79399991 为半程向量。 式（12） 式（12） 是通过统计微平面法线分布与半程向量重合度来估计进入到相机像面上的能量，但这种关系并非绝对成立。因为粗糙微表面会导致微平面间相互遮蔽，这种遮蔽会损耗光线能量 ［ 17 17 ］ 。 图4 图4 蓝色入射光展示了微平面对入射光和出射光的遮蔽现象，从 图4 图4 可以看出，微平面遮挡关系依赖于相机视线方向和入射光方向。 为量化微平面这种不规则性排列对光能量的影响，Smith利用统计学理论，近似求得微平面之间未被遮挡部分的比例，公式 ［ 18 18 ］ 如下： <math display="block" id="M68"><msub><mrow><mi>G</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub><mfenced separators="|"><mrow><mi mathvariant="bold-italic">n</mi><mo>,</mo><mi mathvariant="bold-italic">d</mi><mo>,</mo><mi>k</mi></mrow></mfenced><mo>=</mo><mfrac><mrow><mi mathvariant="bold-italic">n</mi><mo>⋅</mo><mi mathvariant="bold-italic">d</mi></mrow><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">n</mi><mo>⋅</mo><mi mathvariant="bold-italic">d</mi></mrow></mfenced><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>-</mo><mi>k</mi></mrow></mfenced><mo>+</mo><mi>k</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715502&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715513&type= 54.86399841 9.22866726 ， <math display="block" id="M69"><mi>k</mi><mo>=</mo><mfrac><mrow><msup><mrow><mfenced separators="|"><mrow><mi>α</mi><mo>+</mo><mn mathvariant="normal">1</mn></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mn mathvariant="normal">8</mn></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715508&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715517&type= 20.23533440 9.05933285 ， 式中： <math id="M70"><mi mathvariant="bold-italic">d</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715512&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715522&type= 2.45533323 2.87866688 为射线位置矢量， <math id="M71"><mi>α</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715548&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715524&type= 1.94733346 2.79399991 为粗糙度， <math id="M72"><mi mathvariant="bold-italic">n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715538&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715535&type= 2.11666679 2.79399991 为微平面法向量。 考虑到粗糙度会影响反射和入射光线能量，由 式（13） 式（13） 可知，任意微平面未被遮挡的概率可表示为： <math display="block" id="M73"><mi>G</mi><mfenced separators="|"><mrow><mi mathvariant="bold-italic">l</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo>,</mo><mi mathvariant="bold-italic">n</mi></mrow></mfenced><mo>=</mo><msub><mrow><mi>G</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub><mfenced separators="|"><mrow><mi mathvariant="bold-italic">n</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo>,</mo><mi>k</mi></mrow></mfenced><msub><mrow><mi>G</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub><mfenced separators="|"><mrow><mi mathvariant="bold-italic">n</mi><mo>,</mo><mi mathvariant="bold-italic">l</mi><mo>,</mo><mi>k</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715540&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715550&type= 61.46799850 3.97933316 ， 式中： <math id="M74"><mi mathvariant="bold-italic">l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715557&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715556&type= 1.35466671 2.79399991 为入射光位置矢量， <math id="M75"><mi mathvariant="bold-italic">v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715577&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715575&type= 2.03200006 2.79399991 为视线位置矢量。 结合 式（12） 式（12） 和 式（15） 式（15） 可知，成像介质粗糙度对入射光、散射光和反射光的影响可表示为： <math display="block" id="M76"><mi>R</mi><mo>=</mo><mi>D</mi><mfenced separators="|"><mrow><mi mathvariant="bold-italic">n</mi><mo>,</mo><mi mathvariant="bold-italic">h</mi><mo>,</mo><mi>α</mi></mrow></mfenced><mi>G</mi><mfenced separators="|"><mrow><mi mathvariant="bold-italic">l</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo>,</mo><mi mathvariant="bold-italic">n</mi></mrow></mfenced><msub><mrow><mi>L</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715584&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715580&type= 44.02666473 3.97933316 ， 式中： <math id="M77"><msub><mrow><mi>L</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715618&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715607&type= 3.21733332 3.72533321 为单束入射光或散射光能量， <math id="M78"><mi>R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715597&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715571&type= 2.79399991 2.79399991 为 <math id="M79"><msub><mrow><mi>L</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715618&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715607&type= 3.21733332 3.72533321 最后进入到相机像面的能量。 综上所述，增加粗糙度会导致散射光和镜面反射光的自损耗，并提高超出相机视场的概率，引起 <math id="M80"><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715612&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715611&type= 3.47133350 3.72533321 减少， <math id="M81"><msub><mrow><mi>M</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">U</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715626&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715624&type= 6.51933336 3.72533321 降低，进而提高光斑稳定性，减弱像面上记录的光斑亮度。虽然增加粗糙度能抑制散射光在空间中的传播，但反射光也会在一定程度上被抑制。此外，过大的粗糙度导致像面上的能量过小，增加光斑质心提取的难度，因此粗糙度要适中。 2.4　成像介质吸光性对光照的影响 2 根据朗伯比尔定律，物质本身可以吸收光子，物质的颜色是对不同波长可见光选择性吸收的结果。某物质吸收了大部分自然光，那么该物质呈现黑色，反之则表现为白色 ［ 19 19 ］ 。物质对光的吸收度 <math id="M82"><mi>A</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715648&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715629&type= 3.04800010 2.79399991 可表示为： <math display="block" id="M83"><mi>A</mi><mo>=</mo><mi>K</mi><mi>b</mi><mi>c</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715639&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715636&type= 13.97000027 2.79399991 ， 式中： <math id="M84"><mi>K</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715659&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715641&type= 3.04800010 2.79399991 为摩尔吸光系数， <math id="M85"><mi>b</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715667&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715663&type= 1.86266661 2.79399991 为吸收层厚度， <math id="M86"><mi>c</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715680&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715670&type= 1.52400005 2.79399991 为吸光物质浓度。 由 式（17） 式（17） 可知，相比于同浓度、同厚度的白色物质，暗色物质有更强的吸光性，所以用加厚的暗色物质作为成像介质，能让更多光子在成像介质内部消耗，使成像介质表面射出的散射子波数量更少、能量更小，即抑制 <math id="M87"><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715674&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715684&type= 3.47133350 3.72533321 的产生，提高光斑稳定性，减小成像平面的光斑亮度。结合2.2节的分析，粗糙度和吸光性均会一定程度地抑制散射光，选择合适粗糙度和吸光性的成像介质，能提高成像介质表面出射光的信噪比，抑制视觉定位测量系统中随机激光散斑的产生。 3 实 验 1 本文从成像介质密度、表面粗糙度、吸光性和激光输出功率4个方面开展实验。在成像介质实验中，为排除背景光对漏孔靶面光斑提取的影响，实验在暗室进行。其他实验在自然光下完成。实验装置包括相机、全站仪和成像靶面。选取具有不同物理特性的成像介质作为激光照射的靶面，来验证靶面与散射光强度之间的关系。以全站仪的激光器作为激光光源射向靶面，相机仰视看向靶面并将采集的数据实时显示在电脑屏幕上。实验过程中保持相机、激光器和成像靶面的位置不变，激光照射方向可变。其中，相机分辨率为1 280 <math id="M88"><mo>×</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715713&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715709&type= 2.53999996 2.79399991 1 024，帧率为250 frame/s，激光器功率在5~200 mW之间，激光器距离靶面10 m。本文使用一种改进的高斯拟合法 ［ 20 20 ］ 对靶面中的激光光斑进行质心提取，以坐标稳定性来评判不同材质的成像靶面对激光散斑的抑制能力。为了防止激光光斑能量过低时光斑提取受到外界环境干扰，或光斑能量过高出现饱和现象，本文将激光器的输出功率设定为50 mW，并采集连续4 s内的1 000张照片进行质心提取。为评价光斑稳定性，使用贝塞尔公式计算光斑质心标准差，然后统计光斑质心的均值（Mean）、极差值（Range）和标准差（STD）。 3.1　材料密度数据对比实验 2 由2.1节分析可知，激光照射处靶面内部分子的随机运动导致激光随机散斑现象。为验证该结论，本实验选取材质相同、漏孔尺寸分别为0.1，1，2 mm的白色尼龙布（WNC）作为靶面。 图5 图5 中方框内是成像靶面上的激光光斑。 图5 图5 左下角为靶面材料。为展示成像光斑的细节，将实验中靶面上的光斑放大12倍， 图5 图5 右上角为靶面成像光斑放大图。后续的靶面材料粗糙度实验和吸光性实验中，对靶面成像光斑进行相同的放大处理。随着靶面密度的减小，会有更多的光透过靶面，即单位面积内产生的散射子波更少，在成像靶面上记录的光斑稳定性增强，亮度降低。具体实验数据如 表1 表1 所示。 2 Target surface type 2 Mean 2 Range 2 STD X Y X Y X Y 0.1 mm WNC 895.627 644.556 0.121 0.098 0.094 0.090 1 mm WNC 887.694 634.385 0.055 0.093 0.047 0.061 2 mm WNC 891.228 658.119 0.012 0.009 0.009 0.007 从 表1 表1 可知，漏孔越大，靶面密度越小，光斑质心坐标的极差值和标准差越小，光斑越稳定。2 mmWNC靶面上光斑质心坐标的极差值和标准差最小，光斑质心坐标提取的标准差在0.009 pixel内。靶面密度与光斑稳定性成反比，与理论分析结果一致。但漏孔靶面不适用于实际测量，因为大多数实验在多光源环境下进行，背景光会影响靶面光斑提取。 3.2　粗糙度数据对比实验 2 当光照环境、相机和激光器位置以及激光功率不变时，粗糙度增加导致散射光和镜面反射光自损耗的增大，超出相机视场的概率也会相应增加。为验证粗糙度对激光散斑的影响，本实验选取黑色羊绒布（BCC）、黑色玻纤幕布（BFC）、黑色光面亚克力板（BAC）作为靶面，如 图6 图6 所示，3种材料的表面粗糙度依次下降。 图6 图6 中方形框内是成像靶面上的激光光斑，右上角为靶面成像光斑放大图，右下角为靶面材料。 由2.3节的理论分析可知，在一定范围内，粗糙度越大，光斑稳定性越高，亮度越弱。具体数据如 表2 表2 所示。 2 Target surface type 2 Mean 2 Range 2 STD X Y X Y X Y BCC 625.256 345.185 0.087 0.071 0.084 0.073 BFC 624.328 346.726 0.301 0.276 0.281 0.260 BAS 627.871 347.443 0.671 0.588 0.392 0.496 由 表2 表2 可知，粗糙度最大的BCC靶面的光斑坐标极差值和标准差最小，光斑稳定性最高，光斑质心坐标的标准差在0.084 pixel内，所得到的实验结果和理论分析一致。 3.3　吸光性数据对比实验 2 因物质对光的吸收能力与物质本身的摩尔吸光系数以及厚度成正比。为验证吸光性对散射光的影响，本实验选取白色羊绒布（WCC）、黑色羊绒布（BCC）和加厚版黑色羊绒布（TBCC）作为靶面。3种靶面均为同种材质，且密度相同。如 图7 图7 所示，方形框内是成像靶面上的激光光斑，右上角为靶面成像光斑放大图，右下角为靶面材料。 根据 式（16） 式（16） 可知，3种靶面吸光性从左至右逐渐增强，所以加厚版黑色羊绒布上的光斑应最稳定且光斑亮度最低，实验数据如 表3 表3 所示。 2 Target surface type 2 Mean 2 Range 2 STD X Y X Y X Y WCC 624.413 344.166 0.512 0.404 0.297 0.312 BCC 625.256 345.185 0.087 0.071 0.084 0.073 TBCC 625.513 345.416 0.032 0.044 0.023 0.031 由 表3 表3 可知，TBCC靶面的光斑稳定性最高，光斑质心坐标的标准差在0.031 pixel内，这个结果与理论分析的结果一致。 3.4　激光功率数据对比实验 2 由 表2 表2 和 表3 表3 中的实验数据可知，加厚黑色羊绒布的标准差最小。为测试该材质对动态散斑的抑制范围，本实验将激光器的输出功率分别调为10，50，70，90，110，130，170和200 mW，并计算了质心坐标的均值、极差值和标准差，实验数据如 表4 表4 所示。 图8 图8 为8种功率下靶面成像光斑放大13倍后的细节图。 图9 图9 为8种功率下光斑质心STD的变化曲线。 2 Laser power/mW 2 Mean 2 Range 2 STD X Y X Y X Y 10 684.339 516.197 0.031 0.027 0.021 0.011 50 683.946 516.556 0.037 0.072 0.029 0.026 70 684.778 515.727 0.061 0.067 0.049 0.037 90 683.315 516.831 0.394 0.411 0.197 0.201 110 686.324 517.960 0.667 0.842 0.373 0.433 130 686.871 514.713 1.330 1.522 0.811 0.861 170 683.556 513.177 1.897 2.131 1.426 1.821 200 687.334 518.463 2.529 2.614 1.756 1.937 由 表4 表4 和 图9 图9 可知，随着功率的增加，光斑稳定性逐渐降低。激光功率在10~70 mW时，光斑质心的STD较为稳定，而且功率的增加并未明显影响其稳定性。然而，当激光功率大于70 mW时，光斑逐渐呈现饱和状态，光斑质心计算呈现较大的随机性。因此，实际应用时应在激光光斑不出现饱和的状态下选择最佳材质。随着入射光能量的增加，随机散射光的数量和能量也会增加。由2.2节对信噪比的分析可知， <math id="M89"><msub><mrow><mi>L</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715837&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715822&type= 3.47133350 3.72533321 增大，光斑不稳定性会增强，并且 <math id="M90"><msub><mrow><mi>M</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">U</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715841&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45715839&type= 6.51933336 3.72533321 也会提升，使光斑出现饱和状态，理论分析结论和实验结果一致。 4 结 论 1 本文针对高精度视觉定位测量系统中激光随机散斑所引起的光斑不稳定性，进而导致定位误差的问题开展研究，结合激光发光原理和分子运动论，分析出随机散斑是成像介质内分子受激后进行无规则运动，进而导致散射子波在空间中随机干涉所产生的现象，这种现象与散射子波的强度以及方向有关。从图像信噪比角度分析散射光和反射光，得出只有抑制散射光才能提高光斑稳定性和质心提取准确性的结论。在此基础上，使用粗糙度适当和吸光性的成像介质能有效抑制散射光的出现。为验证理论分析的正确性，从成像介质表面密度、粗糙度、吸光性和激光器输出功率4个方面展开实验，得到的实验结果与理论分析的结论完全一致。实验数据表明，激光功率在50 mW左右时，选择STD最小的加厚黑色羊绒布材质作为成像靶面，光斑稳定性在0.04 pixel内，能较好地抑制随机散斑，提升光斑质心的定位精度。 参考文献 朱凡 ， 李颖先 ， 谭久彬 . 高精度激光光束准直系统设计 ［J］. 光学 精密工程 ， 2020 ， 28 （ 4 ）： 817 - 826 . 高精度激光光束准直系统设计 ZHU F ， LI Y X ， TAN J B . Design of high-accuracy laser beam collimation system ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 4 ）： 817 - 826 . （in Chinese） Design of high-accuracy laser beam collimation system 刘伟 ， 周建行 ， 杨延西 . 基于平行激光的便携式视觉裂缝测量系统 ［J］. 仪器仪表学报 ， 2020 ， 41 （ 11 ）： 235 - 243 . doi: 10.19650/j.cnki.cjsi.J2006451 http://dx.doi.org/10.19650/j.cnki.cjsi.J2006451 基于平行激光的便携式视觉裂缝测量系统 LIU W ， ZHOU J H ， YANG Y X . Portable vision crack measurement system based on parallel laser ［J］. Chinese Journal of Scientific Instrument ， 2020 ， 41 （ 11 ）： 235 - 243 . （in Chinese） . doi: 10.19650/j.cnki.cjsi.J2006451 http://dx.doi.org/10.19650/j.cnki.cjsi.J2006451 Portable vision crack measurement system based on parallel laser 贺培鑫 ， 赵彤 ， 沈振民 . 混沌激光对激光显示中的散斑抑制研究 ［J］. 激光与光电子学进展 ， 2022 ， 59 （ 7 ）： 0714001 . doi: 10.3788/LOP202259.0714001 http://dx.doi.org/10.3788/LOP202259.0714001 混沌激光对激光显示中的散斑抑制研究 HE P X ， ZHAO T ， SHEN ZH M . Study on speckle suppression of chaotic laser in laser display ［J］. Laser & Optoelectronics Progress ， 2022 ， 59 （ 7 ）： 0714001 . （in Chinese） . doi: 10.3788/LOP202259.0714001 http://dx.doi.org/10.3788/LOP202259.0714001 Study on speckle suppression of chaotic laser in laser display 邬丹丹 ， 姚康 ， 管凯捷 ， 等 . 基于激光散斑衬比成像技术的脑血流分析 ［J］. 光学 精密工程 ， 2020 ， 28 （ 11 ）： 2411 - 2420 . doi: 10.37188/OPE.20202811.2411 http://dx.doi.org/10.37188/OPE.20202811.2411 基于激光散斑衬比成像技术的脑血流分析 WU D D ， YAO K ， GUAN K J ， et al . Cerebral blood flow analysis based on laser speckle contrast imaging technology ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 11 ）： 2411 - 2420 . （in Chinese） . doi: 10.37188/OPE.20202811.2411 http://dx.doi.org/10.37188/OPE.20202811.2411 Cerebral blood flow analysis based on laser speckle contrast imaging technology 王锐 ， 史瑞新 . 基于多光束照明的回波光场散斑抑制机理 ［J］. 光学 精密工程 ， 2017 ， 25 （ 9 ）： 2333 - 2338 . doi: 10.3788/ope.20172509.2333 http://dx.doi.org/10.3788/ope.20172509.2333 基于多光束照明的回波光场散斑抑制机理 WANG R ， SHI R X . Suppression mechanics of returning wave speckle with multibeams illumination ［J］. Opt. Precision Eng. ， 2017 ， 25 （ 9 ）： 2333 - 2338 . （in Chinese） . doi: 10.3788/ope.20172509.2333 http://dx.doi.org/10.3788/ope.20172509.2333 Suppression mechanics of returning wave speckle with multibeams illumination 张圣涛 ， 高文宏 ， 赵鹏飞 ， 等 . 激光阵列光源角度多样性抑制散斑方法 ［J］. 中国激光 ， 2013 ， 40 （ 1 ）： 48 - 53 . doi: 10.3788/cjl201340.0102005 http://dx.doi.org/10.3788/cjl201340.0102005 激光阵列光源角度多样性抑制散斑方法 ZHANG SH T ， GAO W H ， ZHAO P F ， et al . Laser speckle reduction due to angular diversity introduced by 2D laser diode array ［J］. Chinese Journal of Lasers ， 2013 ， 40 （ 1 ）： 48 - 53 . （in Chinese） . doi: 10.3788/cjl201340.0102005 http://dx.doi.org/10.3788/cjl201340.0102005 Laser speckle reduction due to angular diversity introduced by 2D laser diode array 梁传样 ， 张巍 ， 芮大为 ， 等 . 基于硅基液晶实现激光投影显示系统匀光整形和散斑抑制 ［J］. 中国激光 ， 2016 ， 43 （ 12 ）： 1201001 . doi: 10.3788/cjl201643.1201001 http://dx.doi.org/10.3788/cjl201643.1201001 基于硅基液晶实现激光投影显示系统匀光整形和散斑抑制 LIANG CH Y ， ZHANG W ， RUI D W ， et al . Uniform light shaping and speckle suppression of laser projection display system based on liquid crystal on silicon ［J］. Chinese Journal of Lasers ， 2016 ， 43 （ 12 ）： 1201001 . （in Chinese） . doi: 10.3788/cjl201643.1201001 http://dx.doi.org/10.3788/cjl201643.1201001 Uniform light shaping and speckle suppression of laser projection display system based on liquid crystal on silicon DENG L X ， DONG T H ， FANG Y W ， et al . Speckle reduction in laser projection based on a rotating ball lens ［J］. Optics & Laser Technology ， 2021 ， 135 ： 106686 . doi: 10.1016/j.optlastec.2020.106686 http://dx.doi.org/10.1016/j.optlastec.2020.106686 Speckle reduction in laser projection based on a rotating ball lens 付瀚毅 ， 刘伟奇 ， 魏忠伦 ， 等 . 转动随机微透镜阵列对激光显示中散斑的抑制 ［J］. 液晶与显示 ， 2015 ， 30 （ 2 ）： 359 - 364 . doi: 10.3788/yjyxs20153002.0359 http://dx.doi.org/10.3788/yjyxs20153002.0359 转动随机微透镜阵列对激光显示中散斑的抑制 FU H Y ， LIU W Q ， WEI ZH L ， et al . Speckle reduction using rotating random micro-lens array in laser display ［J］. Chinese Journal of Liquid Crystals and Displays ， 2015 ， 30 （ 2 ）： 359 - 364 . （in Chinese） . doi: 10.3788/yjyxs20153002.0359 http://dx.doi.org/10.3788/yjyxs20153002.0359 Speckle reduction using rotating random micro-lens array in laser display 孙鸣捷 ， 于康龙 ， 孟照魁 . 一种采用振动混光棒的散斑抑制技术 ［J］. 中国激光 ， 2013 ， 40 （ 10 ）： 64 - 69 . doi: 10.3788/cjl201340.1002013 http://dx.doi.org/10.3788/cjl201340.1002013 一种采用振动混光棒的散斑抑制技术 SUN M J ， YU K L ， MENG ZH K . Speckle suppression method with vibrating light pipe ［J］. Chinese Journal of Lasers ， 2013 ， 40 （ 10 ）： 64 - 69 . （in Chinese） . doi: 10.3788/cjl201340.1002013 http://dx.doi.org/10.3788/cjl201340.1002013 Speckle suppression method with vibrating light pipe 张亦斌 ， 方爱平 ， 齐泽坤 ， 等 . 激光散斑形成及漂移现象探究 ［J］. 物理与工程 ， 2021 ， 31（ 6 ） 129 - 132 . doi: 10.3969/j.issn.1009-7104.2021.06.024 http://dx.doi.org/10.3969/j.issn.1009-7104.2021.06.024 激光散斑形成及漂移现象探究 ZHANG Y B ， FANG A P ， QI Z K ， et al . Research on the formation and drift of laser speckle ［J］. Physics and Engineering ， 2021 ， 31（ 6 ） 129 - 132 （in Chinese） . doi: 10.3969/j.issn.1009-7104.2021.06.024 http://dx.doi.org/10.3969/j.issn.1009-7104.2021.06.024 Research on the formation and drift of laser speckle PAUL B ， KAMAL S ， QURESHI T . Multipath wave-particle duality in classical optics ［J］. Optics Letters ， 2020 ， 45 （ 12 ）： 3204 . doi: 10.1364/ol.392762 http://dx.doi.org/10.1364/ol.392762 Multipath wave-particle duality in classical optics 陈明 ， 丁丹 . 聚集态分子内运动触发的生物光热诊疗 ［J］. 科学通报 ， 2021 ， 66 （ 18 ）： 2276 - 2287 . doi: 10.1360/tb-2020-1102 http://dx.doi.org/10.1360/tb-2020-1102 聚集态分子内运动触发的生物光热诊疗 CHEN M ， DING D . Biological photothermal diagnosis and treatment triggered by intramolecular motion in aggregated state ［J］. Chinese Science Bulletin ， 2021 ， 66 （ 18 ）： 2276 - 2287 . （in Chinese） . doi: 10.1360/tb-2020-1102 http://dx.doi.org/10.1360/tb-2020-1102 Biological photothermal diagnosis and treatment triggered by intramolecular motion in aggregated state CHENG J ， CHENG X L ， MENG X C ， et al . A Monte Carlo emissivity model for wind-roughened sea surface ［J］. Sensors ， 2019 ， 19 （ 9 ）： 2166 . doi: 10.3390/s19092166 http://dx.doi.org/10.3390/s19092166 A Monte Carlo emissivity model for wind-roughened sea surface SUN L L ， ZHAO F . Bidirectional reflectance distribution function algorithm based on the Poynting vector analysis ［J］. Optical Engineering ， 2021 ， 60 （ 6 ）： 063104 . doi: 10.1117/1.oe.60.6.063104 http://dx.doi.org/10.1117/1.oe.60.6.063104 Bidirectional reflectance distribution function algorithm based on the Poynting vector analysis TROWBRIDGE T S ， REITZ K P . Average irregularity representation of a rough surface for ray reflection ［J］. Journal of the Optical Society of America ， 1975 ， 65 （ 5 ）： 531 - 536 . doi: 10.1364/josa.65.000531 http://dx.doi.org/10.1364/josa.65.000531 Average irregularity representation of a rough surface for ray reflection BUTLER S D ， NAUYOKS S E ， MARCINIAK M A . Comparison of microfacet BRDF model to modified Beckmann-Kirchhoff BRDF model for rough and smooth surfaces ［J］. Optics Express ， 2015 ， 23 （ 22 ）： 29100 . doi: 10.1364/oe.23.029100 http://dx.doi.org/10.1364/oe.23.029100 Comparison of microfacet BRDF model to modified Beckmann-Kirchhoff BRDF model for rough and smooth surfaces IGOUCHKINE O ， ZHANG Y B ， MA K L . Multi-material volume rendering with a physically-based surface reflection model ［J］. IEEE Transactions on Visualization and Computer Graphics ， 2018 ， 24 （ 12 ）： 3147 - 3159 . doi: 10.1109/tvcg.2017.2784830 http://dx.doi.org/10.1109/tvcg.2017.2784830 Multi-material volume rendering with a physically-based surface reflection model MAYERHÖFER T G ， PAHLOW S ， POPP J . The bouguer-beer-lambert law： shining light on the obscure ［J］. Chemphyschem： a European Journal of Chemical Physics and Physical Chemistry ， 2020 ， 21 （ 18 ）： 2029 - 2046 . doi: 10.1002/cphc.202000464 http://dx.doi.org/10.1002/cphc.202000464 The bouguer-beer-lambert law： shining light on the obscure 任守峰 ， 唐新明 ， 祝小勇 ， 等 . 一种高分七号卫星激光足印光斑质心的高精度提取算法 ［J］. 光学学报 ， 2021 ， 41 （ 10 ）： 1012001 . doi: 10.3788/AOS202141.1012001 http://dx.doi.org/10.3788/AOS202141.1012001 一种高分七号卫星激光足印光斑质心的高精度提取算法 REN SH F ， TANG X M ， ZHU X Y ， et al . A high-precision algorithm for extracting the centroid of laser footprint of high-scoring satellite No.7 ［J］. Acta Optica Sinica ， 2021 ， 41 （ 10 ）： 1012001 . （in Chinese） . doi: 10.3788/AOS202141.1012001 http://dx.doi.org/10.3788/AOS202141.1012001 A high-precision algorithm for extracting the centroid of laser footprint of high-scoring satellite No.7 10.37188/OPE.20233113.1900 TN24 TN249 A 1004-924X（2023）13-1900-09 徐晶,王钎潮,任正玮等.高精度激光定位散斑机理分析及抑制[J].光学精密工程,2023,31(13):1900-1908. 徐晶,王钎潮,任正玮等.高精度激光定位散斑机理分析及抑制[J].光学精密工程,2023,31(13):1900-1908. XU Jing,WANG Qianchao,REN Zhengwei,et al.Mechanism analysis and suppression of speckle in high-precision laser location[J].Optics and Precision Engineering,2023,31(13):1900-1908. XU Jing,WANG Qianchao,REN Zhengwei,et al.Mechanism analysis and suppression of speckle in high-precision laser location[J].Optics and Precision Engineering,2023,31(13):1900-1908. 光电信息控制和安全技术重点实验室项目（No. 2021JCJQLB055011） 2022-11-22 2023-01-13 2023-07-10 2023-07-17 39565280 光学精密工程 13 31