Low-light Image Enhancement Based on Gray Scale Transformation and Improved Retinex
-
摘要: 针对低光照条件下拍摄的图像受光和环境的影响,其重要信息丢失严重,出现对比度低、细节模糊等问题,提出了一种基于灰度变换与改进Retinex的图像增强方法。首先采用引力搜索算法(gravitational search algorithm, GSA)优化的全局灰度变换函数对图像的RGB各通道灰度图像进行灰度变换,增强图像光照强度,使其更接近均匀光照场景;然后将图像转为HSV色彩空间,对V通道(亮度通道)采用改进的多尺度Retinex(MSR)算法处理,将基于范围的自适应双边滤波和Gabor滤波作为Retinex算法的环绕函数,结合两种滤波的特性来增强图像的亮度和细节。最后采用伽马校正避免图像融合造成的图像色偏。实验结果显示,该方法处理过的增强图像在主观和客观评价上优于其他方法,图像颜色失真较小,细节更清晰,为图像的后续应用做了铺垫。Abstract: Aiming at the problems of low contrast and blurred details of images captured under low-light conditions due to the influence of light and the environment, important information is lost, and an image enhancement method based on grayscale transformation and improved Retinex is proposed. First, the global grayscale transformation function optimized by the gravity search algorithm(GSA) is used to perform grayscale transformation on the grayscale image of each RGB channel of the image to enhance the image illumination intensity and make it closer to the uniform illumination scene. The image is then converted to the HSV color space. The V channel (luminance channel) is processed by the improved multi-scale Retinex (MSR) algorithm; range-based adaptive bilateral filtering and Gabor filtering are used as the surround function of the Retinex algorithm, and the characteristics of the two filters are combined to enhance the brightness and detail of the image. Finally, gamma correction is used to avoid image color casts caused by image fusion. The experimental results show that the enhanced image processed by this method is better than that processed by other methods in subjective and objective evaluation, the color distortion of the image is smaller, and the details are clearer, which paves the way for the subsequent application of the image.
-
0. 引言
随着我国近年来综合国力的迅猛发展,各种作战武器更新迭代速度日益加快,尤其各种舰船雨后春笋般迅速投入国防作战系统中,而这也加剧了周边部分国家的紧张感,频繁派遣海军舰船进入我国周边海域,导致区域紧张局势不断加剧。故为了维护我国国家安全与领土完整,需要对敌方作战舰船目标的红外辐射特性进行深入研究,以提高我方红外制导导弹在潜在作战中的实际效能[1-2]。在舰船红外辐射特征研究方面,西安电子科技大学毛佳晨利用计算机图形渲染技术,实现了天基探测平台下舰船目标的红外成像仿真,并分析了不同探测器参数及环境干扰对舰船目标红外成像特征的影响[3],国防科学技术大学的张胜辉基于大量实测与仿真舰船目标红外图像,分析了舰船目标的红外辐射特性,提取了辐射特征[4]。目前国内在舰船目标的红外辐射特征研究主要还集中在计算方法研究和国内目标研究方面,缺乏对国外非合作军事作战舰船的相关研究。因此为支撑作战筹划、实施和作战使用研究对典型作战舰船目标红外目标特性的需求,本文以“阿利•伯克”级驱逐舰为研究对象,通过构建三维舰船模型,采用数值仿真手段分析驱逐舰的流场特性与中、长波红外辐射特性[5-6]。
1. 物理模型
1.1 几何模型
采用商用建模软件CATIA对“阿利•伯克”级驱逐舰进行建模。通过查阅相关文献和资料[7],获得“阿利•伯克”级驱逐舰的几何参数,如表 1所示。通过对比实体模型,对“阿利•伯克”级驱逐舰进行参数化建模,根据后续仿真需求,对“阿利•伯克”级驱逐舰进行了必要的简化,如图 1和图 2所示。
表 1 “阿利·伯克”级驱逐舰主要参数Table 1. Main parameters of the "Ali Burke" class destroyerModel "Ali Burke" class destroyer Length/m 155.29 Width/m 20.4 Draft/m 6 Tonnage/t 9238 Engine 4*LM-2500 gas turbine engine Highest Speed 31knots Endurance 4200 nautical miles(20 knots) 1.2 数学模型
1.2.1 网格划分与边界条件
采用ICEM CFD对“阿利•伯克”级驱逐舰模型进行网格划分。由于“阿利•伯克”级驱逐舰的舰身结构不规则,采用非结构网格进行划分,远场采用结构网格进行划分。
“阿利•伯克”级驱逐舰的舰身网格如图 3所示,整体网格如图 4所示。
通过查阅相关文献资料[8]获得“阿利•伯克”级驱逐舰流场计算边界条件,驱逐舰模型来流方向远场网格边界设置为速度入口,速度为舰船航速28节(52 km/h),在底部方向为水面,忽略气、液交界面,设置为固壁,在其他方向边界设置为压力远场;烟囱入口为压力入口,根据舰船使用燃油组分结构式,计算得到燃烧排放尾气中CO2质量分数WCO2=0.0528,H2O质量分数WH2O=0.0216。
1.2.2 流场计算数值模型
在离散化网格模型基础上采用商业CFD软件ANAYS FLUENT对驱逐舰的流场特性进行求解。在FLUENT数值计算中采用基于压力的隐式求解器,湍流模型选择SST k-ε模型,辐射换热选用DO模型,烟囱内工质燃烧采用组份输运模型,工质选择不可压理想气体,工作压力为标准大气压。
本文对“阿利•伯克”级驱逐舰巡航稳定状态进行仿真,在巡航状态下,驱逐舰周边空气流场是一种复杂的非稳态流动,流体的各种物理参数,如温度、压力、速度等都随时间与空间发生随机的变化,非稳态的Navier-Stokes方程适用于湍流的复杂瞬时运动。本文湍流模型选择工程中应用最为广泛的k-ε模型,在k-ε模型中,为了使流体方程组封闭,引入了湍动能k,以及其散度率ε[9]。ε被定义为:
$$\varepsilon=\frac{\mu}{\rho} \overline{\left(\frac{\partial \mu_i^{\prime}}{\partial x_k}\right)\left(\frac{\partial \mu_i^{\prime}}{\partial x_k}\right)}$$ (1) 式中:μ是动力粘度;ρ是密度。
湍动粘度μt可表示成k和ε的函数,即:
$$\mu_{\mathrm{t}}=\rho C_\mu \frac{k^2}{\varepsilon}$$ (2) 式中:Cμ为常数。
在标准k-ε模型中,与k和ε对应的输运方程为:
$$\left\{\begin{array}{c}\frac{\partial(\rho k)}{\partial t}+\frac{\partial\left(\rho k \mu_i\right)}{\partial x_i}=\frac{\partial}{\partial x_j}\left[\left(\mu+\frac{\mu_t}{\rho_k}\right) \frac{\partial k}{\partial x_j}\right]+ \\ G_k+G_{\mathrm{b}}-\rho \varepsilon-Y_{\mathrm{M}}+S_k \\ \frac{\partial(\rho \varepsilon)}{\partial t}+\frac{\partial\left(\rho \varepsilon \mu_i\right)}{\partial x_i}=\frac{\partial}{\partial x_j}\left[\left(\mu+\frac{\mu_t}{\rho_k}\right) \frac{\partial k}{\partial x_j}\right]+ \\ G_{1 \varepsilon} \frac{\varepsilon}{k}\left(G_k+G_{3 \varepsilon} G_{\mathrm{b}}\right)-G_{2 \varepsilon} \rho \frac{\varepsilon^2}{k}+S_{\varepsilon}\end{array}\right.$$ (3) 式中:Gk是由于平均速度梯度引起的湍动能k的产生项;Gb是由于浮力引起的湍动能k的产生项;YM代表可压湍流中脉动扩张的贡献;C1ε、C2ε、C3ε为经验常数;ρk和ρε分别是与湍动能k和耗散率ε对应的Prandtl数,Sk和Sε是用户定义的源项[10]。
考虑舰船燃油燃烧后排出尾气具有多种气体成分,且不同气体成分具有不同的辐射能力,因此在流场计算过程中,需要考虑气体组分的影响,计算过程中引入组分输运方程。
在流场中,对于某组分Yi,有如下的守恒方程:
$$\frac{\partial}{\partial t}\left(\rho Y_i\right)+\nabla \cdot\left(\rho V Y_i\right)=-\nabla \cdot J_i+R_i+S_i$$ (4) 式中:Ri是化学反应过程中组分Yi的净生成率;Si是离散相引入的生成率与人工定义的源项之和;Ji是组分Yi的扩散流率,在湍流条件下可写为:
$$J_i=-\left(\rho D_{i, m}+\frac{\mu_t}{s c_t}\right) \nabla Y_i-D_{T, i} \frac{\nabla T}{T}$$ (5) 1.2.3 红外辐射计算模型
本文舰船目标的红外辐射计算采用基于逆向蒙特卡洛法(Reverse Monte Carlo Method,RMCM)模型,RMCM对物理过程进行直接模拟,不需要对入射平面的空间立体角进行离散,可以比较容易地解决包含复杂几何形状和各向异性散射等问题,在处理非均匀介质的辐射传输问题时,没有繁琐的积分求和运算,计算时只与表面及介质的本身辐射相联系,避免了求解固体壁面的有效辐射,因此可以取得很大的模拟量,提高了计算的精度[11]。
目标在探测平面处的入射光谱辐射照度Hλ计算公式为:
$$H_\lambda=\int_{\varOmega_{\mathrm{t}}} L_\lambda(\vec{s}) \cos \theta \mathrm{d} \varOmega_{\mathrm{d}}$$ (6) 式中:Lλ是λ波段辐射亮度;θ是探测平面的法线与入射方向夹角;Ω是探测平面立体角。
按照离散传递法的观点,求解上述积分式需要将入射立体角Ωd离散为很多微元立体角,然后在微元立体角中心求解辐射传输方程,得到特征射线在入射方向上到达探测点处的有效辐射亮度,从而将公式(6)转化为数值求和运算。而RMCM方法则不用将入射立体角Ωd离散,它会在Ωd立体角空间内随机地发射N条射线,利用概率密度函数判断射线的吸收点,以这N个具有统计意义的辐射源来代替实际的辐射源[12],因此可将公式(6)写为:
$$ {H_\lambda } = \sum\limits_{i = 1}^N {\left( {{L_{{\text{b}}\lambda }}\left( i \right)\cos {\theta _i}/N} \right)} \cdot {\varOmega _{\text{d}}} $$ (7) 式中:Lbλ(i)为第i条随机射线吸收点处的黑体光谱辐射亮度;θi为第i条随机射线与探测面处法线夹角,当探测点距离目标距离d很远时,cosθi≈1,将公式(7)代入辐射强度计算方程Iλ=Hλ×d2以及Ωλ=A/d2,可得到探测点处光谱辐射强度表达式:
$$ {I_\lambda } = {A_{\text{d}}}\sum\limits_{i = 1}^N {{L_{{\text{b}}\lambda }}\left( i \right)/N} $$ (8) 式中:Ad为目标区域在垂直于探测方位平面上的投影面积。
2. 流场特性计算结果
2.1 温度场
“阿利•伯克”级驱逐舰烟囱尾流温度场与壁面温度场的仿真计算结果如图 5、图 6所示,从图中可以看出,在烟囱出口处尾流温度最高,最高温度达到688.5 K,随着高温尾流与周围空气掺混,尾流温度沿着轴向与径向迅速降低,在一个舰身长度位置尾流温度与空气温度一致。
舰船壁面温度高温位置位于烟囱出口附近,因为靠近烟囱出口位置,受高温尾流热对流与热辐射影响,靠近烟囱出口的桅杆壁面温度较舰体其他壁面温度高,最高达360 K。
2.2 速度场
“阿利•伯克”级驱逐舰流场速度场的仿真计算结果如图 7所示,高温尾流经轮机燃烧从烟囱出口加速排出,尾流最大流速达45 m/s;受舰体遮挡,舰体壁面附近出现滞止区域,流速趋近于0。
2.3 CO2质量浓度场
“阿利•伯克”级驱逐舰流场CO2浓度质量分数的计算结果如图 8所示,驱逐舰燃料为重油,燃烧后尾流中主要气体成分为CO2等,在烟囱出口处CO2质量分数最高,达0.119,随着与空气逐渐掺混,尾流中CO2质量分数逐渐降低。
3. 红外特性计算结果
3.1 参数设置
基于ANSYS FLUENT的仿真计算结果,利用红外辐射理论,计算“阿利·伯克”级驱逐舰红外辐射特性。
红外仿真计算参数设置如表 2所示。
表 2 红外仿真计算参数设置Table 2. Infrared simulation calculation parameters settingInfrared calculation setting parameters Value Detection range/km 5 Number of detection points 60 Gas components H2O, CO2 Gas spectral band model Narrow spectral band Solid emissivity 0.8 Calculation method for atmospheric
transmittanceModtran Spectral calculation range/μm Medium-wave: 3-5
Long-wave: 8-12红外仿真探测点分布如图 9所示,共60个探测点。
3.2 中波红外计算结果
中波波段(3~5 μm)红外计算仿真计算结果如图 10所示,选择正上方、正侧方、前侧方、后侧方4个典型探测点进行展示,可以看出,在中波波段探测下,烟囱高温壁面及排出的高温尾流有着很强的红外辐射特征,而温度较低的舰身在中波探测下几乎不可见。
驱逐舰在中波波段(3~5 μm)典型探测点光谱曲线如图 11所示,在正侧方、前侧方、后侧方探测角度下的红外辐射主要来自于高温尾流,而高温尾流由于对称性,在这3个探测角度上的可视投影面积基本一致,因此3个探测角度下的中波波段红外光谱辐射强度基本一致,由于H2O和CO2分子的吸收与发散作用出现波峰与波谷,在4.6 μm和4.9 μm处为波峰,光谱辐射强度最大值在60000 W/(sr·μm)以上。
在正上方探测角度下的红外辐射由烟囱高温壁面和高温尾流贡献,因此不同于其他3个探测角度在强吸收波段内(3.4~4.2 μm、4.6~4.8 μm)光谱辐射强度值很低,由于高温壁面的贡献,在正上方探测角度下光谱辐射强度在这些波段内依然表现出很强的辐射特征,光谱辐射强度均在6000 W/(sr·μm)以上。但同时相比于其他探测角度,由于正上方探测角度下高温尾流可视面积小,导致在4.6 μm和4.9 μm两处波峰处的光谱辐射强度较低,最高为28000 W/(sr·μm)。
3.3 长波红外计算结果
长波波段(8~12 μm)仿真计算结果如图 12所示,选择正上方、正侧方、前侧方、后侧方4个典型探测点,可以看出,相较于中波波段,除高温尾流外,在长波红外下舰身辐射特征也非常明显,舰体轮廓清晰可见。
驱逐舰在长波波段(8~12 μm)典型探测点光谱曲线如图 13所示,舰体蒙皮是长波波段主要的红外辐射源,从正上方、正侧方、前侧方、后侧方视角,舰体投影面积逐渐减小,因此在舰体温度各方向基本一致的情况下,目标的光谱辐射强度也逐渐降低。
从仿真计算结果可看出,由于在正上方视角下驱逐舰目标可视投影面积最大,在后侧方视角下可视投影面积最小,因此长波波段下红外辐射强度在正上方视角下辐射强度最高,最大值达18200 W/(sr·μm),在后侧方长波辐射强度值最小。
3.4 仿真模型准确性验证
公开资料无法获取外场实测“阿利·伯克”级驱逐舰及相似军舰目标的红外辐射特征分布。因此为验证仿真模型与结果准确性,以民用渔船为验证对象(图 14),对比分析了中波和长波波段渔船外场实测数据与仿真数据,结果如图 15所示,从对比结果可以看出,仿真结果与外场试验实测结果整体趋势相同,整体吻合度较好,满足仿真准确性的要求。
4. 结论
针对红外导弹对驱逐舰目标的红外辐射特性的研究需求,本文搭建了“阿利•伯克”级驱逐舰目标红外辐射特性模型,首先建立了“阿利•伯克”级驱逐舰的三维物理模型,通过仿真模拟计算得到了模型在巡航状态下的流场特征与红外辐射特征,并采用民用渔船实测数据与仿真数据对比,验证了“阿利•伯克”级驱逐舰模型计算准确度,可用于反舰红外制导导弹打击舰船目标作战评估及战法推演。所得到的主要结论有:
1)在巡航状态下,“阿利•伯克”级驱逐舰高温部件主要是烟囱壁面和羽流,温度最高达688.5 K,尾流的最大流速达45 m/s,且尾流中CO2质量分数最高达0.119;
2)舰船目标在中波和长波波段下红外辐射特征差异明显,在中波波段探测下,仅可见高温尾流和高温壁面,在长波波段探测下,低温舰体和高温尾流均清晰可见;
3)舰船目标的中波波段光谱辐射主要由高温尾流贡献,长波波段光谱辐射主要由舰体壁面贡献,且中波波段光谱辐射强度最大值(65000 W/(sr·μm))明显大于长波波段光谱辐射强度最大值(18200 W/(sr·μm));
4)舰船目标的光谱辐射强度在不同探测角度下差异明显,主要是因为在不同探测角度下目标的可视投影面积区别较大。
综上所述,针对作战实际情况可对红外导弹探测波段选择提供支撑,在明确重点打击目标部位后,依据中长波探测差异综合判断选择,支撑单波段或者多波段红外导弹探测能力论证。
-
表 1 各算法参数
Table 1 Parameters of each algorithm
Method Parameter BSSR sigma _s={15, 70, 110};
sigma _s ={0.05, 0.1, 0.3};
Filter size N=2;Ref. [4] sigma ={64, 128, 256};
Color recovery factor C=1/3V-SSR sigma=256 Ours sigma _s={15, 70, 110};
sigma _d={0.05, 0.1, 0.3};
ρ=3; filter size N=2;
Weight factor m=0.4表 2 图像A、B、C客观质量评价
Table 2 Objective quality evaluation of images A, B and C
Method Image A Image B Image C IE MG PSNR IE MG PSNR IE MG PSNR Original image 6.8780 0.0141 - 6.6150 0.0255 - 7.1937 0.0375 - BSSR 7.5857 0.0380 52.6133 7.5051 0.0564 41.3939 7.6856 0.0982 36.2253 Ref.[4] 7.5335 0.0465 51.6841 7.4900 0.0574 37.5988 7.6190 0.0961 35.2424 V-SSR 7.4974 0.0349 49.9789 7.4368 0.0540 39.4727 7.4950 0.0904 34.0436 Ours1 7.6121 0.0344 68.0142 7.5252 0.0598 54.4110 7.8207 0.0651 55.6753 Ours 7.8976 0.0608 75.2760 7.5871 0.0499 64.4282 7.9459 0.1326 79.0920 -
[1] Dhal K G, Ray S, Das A, et al. A survey on nature-inspired optimization algorithms and their application in image enhancement domain[J]. Archives of Computational Methods in Engineering, 2019, 26(5): 1607-1638. DOI: 10.1007/s11831-018-9289-9
[2] GAO H, ZENG W, CHEN J. An improved gray-scale transformation method for pseudo-color image enhancement[J]. Computer Optics, 2019, 43(1): 78-82. http://www.researchgate.net/publication/332056355_An_improved_gray-scale_transformation_method_for_pseudo-color_image_enhancement
[3] Hassan N, Ullah S, Bhatti N, et al. The Retinex based improved underwater image enhancement[J]. Multimedia Tools and Applications, 2020(1): 1-19.
[4] SHAO Wenbo, LIU Lei, JIANG Jiawei, et al. Low-light-level image enhancement based on fusion and Retinex[J]. Journal of Modern Optics, 2020, 67(13): 1190-1196. DOI: 10.1080/09500340.2020.1823502
[5] 牟琦, 魏妍妍, 李姣, 等. 改进的Retinex低照度图像增强算法研究[J]. 哈尔滨工程大学学报, 2018, 39(12): 2001-2010. https://www.cnki.com.cn/Article/CJFDTOTAL-HEBG201812019.htm MU Q, WEI Y Y, LI J, et al. Research on the improved Retinex algorithm for low illumination image enhancement[J]. Journal of Harbin Engineering University, 2018, 39(12): 2001-2010. https://www.cnki.com.cn/Article/CJFDTOTAL-HEBG201812019.htm
[6] 许凤麟, 苗玉彬, 张铭. 基于彩色加权引导滤波-Retinex算法的导航图像增强[J]. 上海交通大学学报, 2019, 53(8): 921-927. https://www.cnki.com.cn/Article/CJFDTOTAL-SHJT201908007.htm XU F L, MIAO Y B, ZHANG M. Navigation image enhancement based on color weighted guided image filtering Retinex algorithm[J]. Journal of Shanghai Jiaotong University, 2019, 53(8): 921-927. https://www.cnki.com.cn/Article/CJFDTOTAL-SHJT201908007.htm
[7] Oh J G, Hong M C. Adaptive image rendering using a nonlinear mapping-function-based Retinex model[J]. Sensors, 2019, 19(4): 969-985. DOI: 10.3390/s19040969
[8] GAO H, ZENG W, CHEN J. An improved gray-scale transformation method for pseudo-color image enhancement[J]. Computer Optics, 2019, 43(1): 78-82.
[9] SUN X, XU Q, ZHU L. An effective Gaussian fitting approach for image contrast enhancement[J]. IEEE Access, 2019, 7: 31946- 31958.
[10] Verma H K, Pal S. Modified sigmoid function based gray scale image contrast enhancement using particle swarm optimization[J]. Journal of the Institution of Engineers, 2016, 97(2): 243-251.
[11] 刘宾, 赵鹏翔, 赵霞, 等. 变能量X射线融合图像的增强算法研究[J]. 光学学报, 2020, 40(18): 216-221. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB202018025.htm LIU B, ZHAO P X, ZHAO X, et al. Enhancement algorithm of variable energy X-Ray fusion images[J]. Acta Optica Sinica, 2020, 40(18): 216-221. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB202018025.htm
[12] Katirciolu F, Cngz Z. A novel gray image enhancement using the regional similarity transformation function and dragonfly algorithm[J]. El-Cezeri Fen ve Mühendislik Dergisi, 2020, 7(3): 1201-1219.
[13] JIANG S, JI Z, SHEN Y. A novel hybrid particle swarm optimization and gravitational search algorithm for solving economic emission load dispatch problems with various practical constraints[J]. International Journal of Electrical Power & Energy Systems, 2014, 55: 628-644.
[14] 魏赟, 欧阳鹏. 基于快速亮通双边滤波器的Retinex图像增强算法[J]. 小型微型计算机系统, 2021, 42(9): 1944-1949. https://www.cnki.com.cn/Article/CJFDTOTAL-XXWX202109026.htm WEI Y, OU Y P. Retinex image enhancement algorithm based on fast bright-pass bilateral filtering[J]. Journal of Chinese Computer Systems, 2021, 42(9): 1944-1949. https://www.cnki.com.cn/Article/CJFDTOTAL-XXWX202109026.htm
[15] WANG Z, WANG Z. A generic approach for cell segmentation based on gabor filtering and area-constrained ultimate erosion[J]. Artificial Intelligence in Medicine, 2020, 107: 101929.
[16] 王奎, 黄福珍. 基于光照补偿的HSV空间多尺度Retinex图像增强[J/OL]. 激光与光电子学进展[2023-01-29]. http://kns.cnki.net/kcms/detail/31.1690.TN.20210713.1635.050.html. WANG K, Huang F Z. Multi-scale Retinex Image Enhancement in HSV Space Based on Illumination Compensation[J/OL]. Laser & Optoelectronics Progress [2023-01-29]. http://kns.cnki.net/kcms/detail/31.1690.TN.20210713.1635.050.html
[17] ZHANG Z, DAI Y, MA C. Design of selectively multilayered periodic gratings by PSO algorithm for radiative cooling[J]. Optics Communications, 2021, 500: 127323.
[18] Bakhshipour M, Namdari F, Samadinasab S. Optimal coordination of overcurrent relays with constraining communication links using DE–GA algorithm[J]. Electrical Engineering, 2021: 1-15.
[19] 吕晓宁, 刘扬阳, 谭政, 等. 一种偏振普适性多尺度实时的图像去雾算法[J]. 光子学报, 2019, 48(8): 117-127. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB201908014.htm LU Xiaoning, LIU Yangyang, TAN Zheng, et al. A polarizing universal multi-scale and real-time image defogging algorithm[J]. Acta Photonica Sinica, 2019, 48(8): 117-127. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB201908014.htm
[20] 吴鹏飞, 王辉亮, 雷思琛, 等. 大气湍流环境下光斑中心定位算法[J/OL]. 光子学报[2022-03-15]. http://kns.cnki.net/kcms/detail/61.1235.O4.20211022.1807.008.html. WU Pengfei, WANG Huiliang, LEI Sichen, et al. Spot center localization algorithm in atmospheric turbulence environment[J/OL]. Acta Photonica Sinica [2022-03-15]. http://kns.cnki.net/kcms/detail/61.1235.O4.20211022.1807.008.html.