留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

Hybrid PSO-AMLS-based method for data fitting in the calibration of the infrared radiometer

Yang Sen

杨森. 基于PSO-AMLS的数据拟合方法在红外辐射计校准中的应用[J]. 红外与激光工程, 2021, 50(8): 20200471. doi: 10.3788/IRLA20200471
引用本文: 杨森. 基于PSO-AMLS的数据拟合方法在红外辐射计校准中的应用[J]. 红外与激光工程, 2021, 50(8): 20200471. doi: 10.3788/IRLA20200471
Yang Sen. Hybrid PSO-AMLS-based method for data fitting in the calibration of the infrared radiometer[J]. Infrared and Laser Engineering, 2021, 50(8): 20200471. doi: 10.3788/IRLA20200471
Citation: Yang Sen. Hybrid PSO-AMLS-based method for data fitting in the calibration of the infrared radiometer[J]. Infrared and Laser Engineering, 2021, 50(8): 20200471. doi: 10.3788/IRLA20200471

基于PSO-AMLS的数据拟合方法在红外辐射计校准中的应用

doi: 10.3788/IRLA20200471
详细信息
  • 中图分类号: TN215

Hybrid PSO-AMLS-based method for data fitting in the calibration of the infrared radiometer

More Information
    Author Bio:

    杨森,男,讲师,硕士生导师,博士,主要研究方向为红外辐射测量、路面状态检测、食品安全检测

  • 摘要: 针对红外辐射计标定中的数据拟合问题,提出了一种基于PSO-AMLS算法的数据拟合方法。该方法结合了粒子群算法(PSO)和自适应移动最小二乘法(AMLS)。采用的PSO算法可以实现自适应移动最小二乘法中参数的最优化,提高了红外辐射计的校准精度。区别于传统的应用于红外辐射计校准中的数据拟合方法,文中提出的基于局部近似的PSO-AMLS数据拟合方法,不仅可以得到红外辐射计输出电压与准直辐射源辐射之间的对应关系,而且在数据处理中通过灵活的数据处理机制,可以降低因拟合数据不均匀分布对拟合结果产生的不利影响。实验结果表明,文中提出的PSO-AMLS数据拟合方法,优于红外辐射计校准中所用的传统数据拟合方法,可以提高数据拟合精度。
  • Figure  1.  Flowchart of the hybrid PSO-AMLS-based method

    Figure  2.  The fitting errors δ of the LS and PSO-AMLS method when the independent variables are uniform distribution: (a) nnodes=11; (b) nnodes=21

    Figure  3.  The fitting errors δ of the ALS and the PSO-AMLS method when the independent variables are non-uniform distribution

    Figure  4.  Diagram of the infrared radiometer

    Figure  5.  The fitting errors δ of the LS and ALS method: (a) LS method; (b) ALS method

    Figure  6.  The fitting errors of the PSO-AMLS and ALS method

    Figure  7.  The deviation δ* obtained for different fitting methods

    Table  1.   The data used in the first example

    xyxyxy
    0173.6458144.7417
    1283.8284154.873
    22.414294165
    32.7321104.1623175.1231
    43114.3166185.2426
    53.2361124.4641195.3589
    63.4495134.6056205.4721
    下载: 导出CSV

    Table  2.   The data used in the second example

    yxyxyx
    11611.51146.5
    21.5716.51256.5
    32.5822.51367.5
    44.5929.51479.5
    57.51037.51592.5
    下载: 导出CSV

    Table  3.   Main specifications of the infrared radiometer

    ComponentParameter
    Optical lensDiameter: 70 mm
    Field of view: ±0.6°
    Spectral channelsMW: 3.7-4.8 µm
    DetectorType: Insb
    A/DPrecision: 16
    Frequency: 100 kHz
    Preamplifier106×
    下载: 导出CSV

    Table  4.   Main specifications of the blackbody-collimator

    ComponentParameter
    BlackbodyModel: HFY-200C
    Temperature: 5-1000 ℃
    ApertureDiameter: 0.1-13.5 mm
    CollimatorModel: HGD-1
    Focal length: 650 mm
    下载: 导出CSV

    Table  5.   The calibration data of T and V in the calibration of the infrared radiometer

    V /VT /℃V /VT /℃
    0.021001.354450
    0.051501.828500
    0.1172002.383550
    0.2272502.994600
    0.3943003.75650
    0.6323504.555700
    0.949400
    下载: 导出CSV
  • [1] Ville A, Steven W B, Thomas C L, et al. Comparison of absolute spectral irradiance responsivity measurement techniques using wavelength-tunable lasers [J]. Applied Optics, 2007, 46(20): 4228-4236. doi:  10.1364/AO.46.004228
    [2] Hartmann J. High-temperature measurement techniques for the application in photometry, radiometry and thermometry [J]. Physics Reports, 2009, 469(5-6): 205-269. doi:  10.1016/j.physrep.2008.09.001
    [3] Wang X, Gao Z Y, Zhang J Y, et al. Research on calibration method of three band infrared integrated radiometer[C]//Proceedings of SPIE, 2002, 4927: 133-138.
    [4] Ibrahim R, John R H, Julian G, et al. Calibrating pyrgeometers outdoors independent from the reference value of the atmospheric longwave irradiance [J]. Journal of Atmospheric and Solar-Terrestrial Physics, 2006, 68(12): 1416-1424. doi:  10.1016/j.jastp.2006.05.013
    [5] Anevsky S, Krutikov V, Minaeva O, et al. Method for the calibration of the spectral irradiance of tungsten filament transfer standard sources traceable to synchrotron radiation [J]. Applied Optics, 2013, 52(21): 5152-5157. doi:  10.1364/AO.52.005152
    [6] Huang L K, Cebula R P, Hilsenrath E. New procedure for interpolating NIST FEL lamp irradiances [J]. Metrologia, 1998, 35(4): 381-386. doi:  10.1088/0026-1394/35/4/27
    [7] Samedov F, Durak M, Bazkir O. Filter-radiometer-based realization of candela and establishment of photometric scale at UME [J]. Optics and Lasers in Engineering, 2005, 43(11): 1252-1266. doi:  10.1016/j.optlaseng.2004.12.001
    [8] Durak M, Aslan M H. Optical characterization of the silicon photodiodes for the establishment of national radiometric standards [J]. Optics & Laser Technology, 2004, 36(3): 223-227. doi:  10.1016/j.optlastec.2003.08.009
    [9] Minato H, Ishido Y. Development of a five-band multi-spectral infrared radiometer for low temperature measurements [J]. Review of Scientific Instruments, 2003, 74(5): 2863-2870. doi:  10.1063/1.1569400
    [10] Cordero R R, Seckmeyer G, Pissulla D, et al. Uncertainty evaluation of spectral UV irradiance measurements [J]. Measurement Science and Technology, 2008, 19(4): 045104. doi:  10.1088/0957-0233/19/4/045104
    [11] Revercomb H E, Buijs H, Howell H B, et al. Radiometric calibration of IR Fourier transform spectrometers: solution to a problem with the High-Resolution Interferometer Sounder [J]. Applied Optics, 1988, 27(15): 3210-3218. doi:  10.1364/AO.27.003210
    [12] Flynn D S, Marlow S A, Sisko R B, et al. Accuracy of aperture irradiances from a resistor-array projection system[C]//Proceedings of SPIE, 2005, 5785: 112-123.
    [13] Zhang H Q, Guo C X, Su X F, et al. Measurement data fitting based on moving least squares method [J]. Mathematical Problems in Engineering, 2015, 2015: 195023. doi:  10.1155/2015/195023
    [14] Zhang L, Gu T Q, Zhao J, et al. An improved moving least squares method for curve and surface fitting [J]. Mathematical Problems in Engineering, 2013, 2013: 159694. doi:  10.1155/2013/159694
    [15] Ghosh S, Das S, Kundu D, et al. An inertia-adaptive particle swarm system with particle mobility factor for improved global optimization [J]. Neural Computing and Applications, 2012, 21(2): 237-250. doi:  10.1007/s00521-010-0356-x
    [16] Shepard D. A two-dimensional interpolation function for irregularly spaced points[C]//Proceedings-1968 ACM National Conference, 1968: 517-524.
    [17] Lancaster P, Salkauskas K. Surfaces generated by moving least squares methods [J]. Mathematics of Computation, 1981, 37(155): 141-158. doi:  10.1090/S0025-5718-1981-0616367-1
    [18] Lancaster P, Salkauskas K. Curve and surface fitting: An introduction [J]. SIAM Rev, 1989, 31(1): 155-157. doi:  10.1137/1031035
    [19] Clerc M, Kennedy J. The particle swarm-explosion, stability, and convergence in a multidimensional complex space [J]. IEEE Transactions on Evolutionary Computation, 2002, 6(1): 58-73. doi:  10.1109/4235.985692
    [20] García N P J, García G E, Alonso F J R, et al. A hybrid PSO optimized SVM-based model for predicting a successful growth cycle of the Spirulina platensis from raceway experiments data [J]. Journal of Computational and Applied Mathematics, 2016, 291: 293-303. doi:  10.1016/j.cam.2015.01.009
  • [1] 李小明, 朱国帅, 郭名航, 刘赢泽, 张崇.  基于光学自准直的旋转轴平行度测量与不确定度分析 . 红外与激光工程, 2023, 52(5): 20220794-1-20220794-12. doi: 10.3788/IRLA20220794
    [2] Fu Yanjun, Cai Xiaoqi, Zhong Kejun, Ma Baiheng, Yan Zhanjun.  Method for phase-height mapping calibration based on fringe projection profilometry . 红外与激光工程, 2022, 51(4): 20210403-1-20210403-9. doi: 10.3788/IRLA20210403
    [3] 李凯, 王玄玉, 高艳卿, 董文杰.  石墨烯红外波段复折射率及消光性能研究 . 红外与激光工程, 2021, 50(4): 20200246-1-20200246-7. doi: 10.3788/IRLA20200246
    [4] He Zixin, Shao Hui, Guo Hang, Chen Jie.  Classification of coal/rock based on Hyperspectral LiDAR calibration-free signals . 红外与激光工程, 2021, 50(10): 20200518-1-20200518-9. doi: 10.3788/IRLA20200518
    [5] Yang Guoqing, Li Zhou, Zhao Chen, Yu Yi, Qiao Yanfeng, He Fengyun.  Nonlinear atmospheric correction based on neural network for infrared target radiometry . 红外与激光工程, 2020, 49(5): 20190413-20190413-8. doi: 10.3788/IRLA20190413
    [6] 赵志军, 许方宇, 高玲, 郭杰, 徐世春, 刘忠.  3~5μm红外天空亮度测量的多元定标模型 . 红外与激光工程, 2017, 46(10): 1004004-1004004(8). doi: 10.3788/IRLA201783.1004004
    [7] 殷世民, 高丽伟, 梁永波, 朱健铭, 梁晋涛, 陈真诚.  基于FPGA的干涉式红外成像光谱仪实时光谱复原研究 . 红外与激光工程, 2017, 46(7): 720001-0720001(6). doi: 10.3788/IRLA201746.0720001
    [8] 徐文斌, 陈伟力, 李军伟, 王广平, 武敬力.  采用长波红外高光谱偏振技术的目标探测实验 . 红外与激光工程, 2017, 46(5): 504005-0504005(7). doi: 10.3788/IRLA201746.0504005
    [9] 徐超, 何利民, 王霞, 金伟其.  红外偏振成像系统高速处理模块设计 . 红外与激光工程, 2017, 46(2): 204002-0204002(8). doi: 10.3788/IRLA201746.0204002
    [10] 周星光, 贺宇, 王岭雪, 蔡毅, 刘福平, 李茂忠, 陈骥.  单视点双曲面折反射红外全景成像系统设计与分析 . 红外与激光工程, 2016, 45(9): 918004-0918004(9). doi: 10.3788/IRLA201645.0918004
    [11] 顾有林, 王成, 杨丽, 欧宗伟, 胡以华, 李乐, 赵义正, 陈卫, 王鹏.  黑曲霉孢子灭活前后红外消光特性 . 红外与激光工程, 2015, 44(1): 36-41.
    [12] 王济洲, 李宏, 熊玉卿, 董茂进, 张玲, 李晨.  一种具有抗静电反红外诱导滤光片的设计与制备 . 红外与激光工程, 2015, 44(10): 3005-3009.
    [13] 陈青山, 牛春晖, 吕勇, 徐岱, 魏元.  对镱铒共掺硫氧化钇红外上转换材料的远场探测与识别 . 红外与激光工程, 2015, 44(9): 2603-2608.
    [14] 乔立永, 徐立新, 高敏.  带宽自适应均值偏移红外目标跟踪 . 红外与激光工程, 2015, 44(1): 354-362.
    [15] 张金玉, 孟祥兵, 杨正伟, 王冬冬, 陶胜杰.  红外锁相法涂层测厚数值模拟与分析 . 红外与激光工程, 2015, 44(1): 6-11.
    [16] 邓潺, 刘文, 黄彪, 詹振环, 张欣.  基于多核DSP 的微光与红外图像彩色融合系统 . 红外与激光工程, 2014, 43(9): 3141-3145.
    [17] 杨少魁, 刘文.  一种微光与红外图像彩色融合方法 . 红外与激光工程, 2014, 43(5): 1654-1659.
    [18] 赵纪金, 李晓霞, 郭宇翔, 杨莉.  膨胀石墨体积膨胀率对红外遮蔽性能的影响 . 红外与激光工程, 2014, 43(2): 434-437.
    [19] 王孟军, 赵翠玲, 韩邦杰, 耿亚光, 马萄.  基于谐衍射的共轴双波段红外光学系统设计 . 红外与激光工程, 2013, 42(10): 2732-2736.
    [20] 于洋, 蹇毅, 潘兆鑫, 金亚平, 汤心溢.  红外二次成像无热化光学系统设计与实测 . 红外与激光工程, 2013, 42(12): 3180-3184.
  • 加载中
图(7) / 表(5)
计量
  • 文章访问数:  388
  • HTML全文浏览量:  128
  • PDF下载量:  21
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-12-07
  • 修回日期:  2021-03-24
  • 刊出日期:  2021-08-25

Hybrid PSO-AMLS-based method for data fitting in the calibration of the infrared radiometer

doi: 10.3788/IRLA20200471
    作者简介:

    杨森,男,讲师,硕士生导师,博士,主要研究方向为红外辐射测量、路面状态检测、食品安全检测

  • 中图分类号: TN215

摘要: 针对红外辐射计标定中的数据拟合问题,提出了一种基于PSO-AMLS算法的数据拟合方法。该方法结合了粒子群算法(PSO)和自适应移动最小二乘法(AMLS)。采用的PSO算法可以实现自适应移动最小二乘法中参数的最优化,提高了红外辐射计的校准精度。区别于传统的应用于红外辐射计校准中的数据拟合方法,文中提出的基于局部近似的PSO-AMLS数据拟合方法,不仅可以得到红外辐射计输出电压与准直辐射源辐射之间的对应关系,而且在数据处理中通过灵活的数据处理机制,可以降低因拟合数据不均匀分布对拟合结果产生的不利影响。实验结果表明,文中提出的PSO-AMLS数据拟合方法,优于红外辐射计校准中所用的传统数据拟合方法,可以提高数据拟合精度。

English Abstract

    • To describe the measured radiation of infrared targets continuously and accurately, data fitting is necessary for the calibration of the infrared radiometer. Data fitting is essentially a typical function approximation problem. The fitting accuracy will be different owing to the differences between fitting methods, which will influence the measurement accuracy of the infrared radiometer.

      In the calibration of the infrared radiometer, the least-squares (LS) method has been most widely used in data fitting[1-4]. The commonly used approximation function of the LS method is polynomials[5-10], as well as others such as sinusoidal function and log function[11-12]. However, the LS method is a global approximation scheme, which is not suitable for irregular or non-uniform distribution cases. For the LS method, the oscillation phenomenon occurs and the oscillation increases dramatically when the order of polynomial becomes larger, which reduces the fitting accuracy. For the calibration of the infrared radiometer, problems observed contain the non-uniform distribution of fitting data, the dramatic change in the order of magnitude, and the limitation on the number of calibration points, which resulting in sufficient fitting precision cannot be achieved by the LS method.

      To improve the fitting accuracy, The Adaptive Moving Least Squares (AMLS) method is used in this paper for data fitting in the calibration of the infrared radiometer, which combines Adaptive Processing (AP) and Moving Least Squares (MLS). AP procedures such as data transformation or data translation are performed before the calibration data fitting to reduce the adverse effect resulting from irregular or non-uniform data distribution. The MLS method combines the concept of moving window and compact support weighting functions, which can be regarded as a combination of Weighted Least Squares (WLS) and Segmented Least Square (SLS). As a local approximation method, The MLS method can not only acquire higher precision even with low order basis functions but also has good stability due to its local approximation scheme[13]. However, the parameters setting in the MLS method are often chosen through individual experience or simple parameters adjustment[14], by which optimal fitting accuracy is difficult to achieve. Thus, to carry out the optimization mechanism corresponding to the setting of weighting function parameters in the MLS method, the Particle Swarm Optimization (PSO) technique is used with the AMLS method. As a “speed-position” searching algorithm, the particle adjusts its speed vector value and direction following the optimal particle, and then the optimal penalty parameters of samples can be obtained[15].

      In this paper, the hybrid PSO-AMLS-based method is proposed for data fitting in the calibration of the infrared radiometer. The present study is structured as follows: firstly, proposed methods and principles are explained in Section 1; secondly, numerical examples of curve fitting is carried out based on hybrid PSO-AMLS-based method in Section 2, containing the examples when the independent variables are uniform distribution and non-uniform distribution; thirdly, the hybrid PSO-AMLS-based method is applied to the calibration of the infrared radiometer in Section 3; and finally, the main conclusions of this research work are drawn.

    • As a conventional data fitting method, the MLS method was proposed by Shepard in 1968[16]. The basic theory of the MLS can be seen in references [17] and [18], so this paper only gives a brief introduction. In MLS, an arbitrary function f(x) can be approximated by

      $$ f(x) = \sum\limits_{i = 1}^m {{p_i}} (x){a_i}(x) = {P^{\rm T}}(x)a(x) $$ (1)

      where pi(x) are basis functions, ai(x) are the coefficients, the coefficients a(x)=[a1(x), a2(x),···, am(x)]T varied with x. In order to determine the coefficients ai(x), a function is defined as

      $$ J = \sum\limits_{I = 1}^N {{w_I}} (x){[P({x_I})a(x) - f({x_I})]^2} $$ (2)

      where f(xI) are the given nodes, ω(x) is the weighting function with compact support.

      One commonly used weighting function is the Gaussian function, which is very flexible for the MLS and is adopted in the following discussion. The Gaussian weighting function is

      $$ w(r) = \left\{ {\begin{array}{*{20}{c}} {\dfrac{{{\rm exp}( - {r^2}{\beta ^2}) - {\rm exp}( - {\beta ^2})}}{{1 - {\rm exp}( - {\beta ^2})}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} r \leqslant 1} \\ {{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} 0{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} r > 1} \end{array}} \right. $$ (3)

      where r=|xxI|/dmi is the relative distance, dmi is the influencing radius, and β is the shape parameter.

      Then, the matrix form for equation (2) can be rewritten as

      $$ J = {[Pa(x) - f]^{\rm T}}W[Pa(x) - f] $$ (4)

      where f =[f (x1), f (x2),···, f (xN)]T, P and W are

      $$ P = {\left[ {\begin{array}{*{20}{c}} {{P^{\rm T}}({x_1})} \\ {{P^{\rm T}}({x_2})} \\ \vdots \\ {{P^{\rm T}}({x_N})} \end{array}} \right]_{N \times m}} $$ (5)
      $$ W = \left[ {\begin{array}{*{20}{c}} {\omega (x - {x_1})}&0& \cdots &0 \\ 0&{\omega (x - {x_2})}& \cdots &0 \\ \vdots & \vdots & \ddots & \vdots \\ 0&0& \cdots &{\omega (x - {x_N})} \end{array}} \right] $$ (6)

      Minimizing equation (4) with respect to coefficients a(x), the following expression can be obtained

      $$ \frac{{\partial J}}{{\partial a(x)}} = [{P^{\rm T}}WP]a(x) - [{P^{\rm T}}W]f = 0 $$ (7)

      Further, the coefficients a(x) and function f(x) are obtained as follows

      $$ a(x) = {[{P^{\rm T}}WP]^{{\rm{ - }}1}}[{P^{\rm T}}WP]f $$ (8)
      $$ f(x) = {P^{\rm T}}(x){[{P^{\rm T}}WP]^{{\rm{ - }}1}}[{P^{\rm T}}WP]f $$ (9)

      However, the calculation process of the MLS method is more complicated and requires a larger amount of computation.

    • The PSO technique is a population-based search algorithm based on the simulation of the bird flocking. The basic theory of the MLS can be seen in reference [19], so this paper only gives a brief introduction. In the PSO, the parameters or possible set of solutions are contained in a vector xi(k), which is called a particle of the swarm and represents its position in the search space of possible solutions. The particle dimension is the number of parameters. The particle initial position xi(0) and its velocity vi(0) are chosen randomly. The value of the fitness function is then calculated for each particle and the velocities, and then the positions are updated depending on these values. The algorithm updates the positions and the velocities of the particles as follows:

      $$ \begin{split} {v_i}\left( {k + 1} \right) =& \omega \cdot {v_i}\left( k \right) + {c_1} \cdot {r_1} \cdot \left( {{p_{ibest}} - {x_i}\left( k \right)} \right) + {c_2} \cdot \\ & {r_2} \cdot \left( {{g_{best}} - {x_i}\left( k \right)} \right) \\ \end{split} $$ (10)
      $$ {x_i}\left( {k + 1} \right) = {x_i}\left( k \right) + {v_i}\left( {k + 1} \right) $$ (11)

      where ω is the inertia factor, pibest is the local best position; gbest is the global best position; c1 and c2 are learning factors; rand() is a random number in the interval [0, 1].

      The current position xi(k) and velocity vi(k) of each particle are adjusted depending on its status in the previous step, the particle’s local best position pibest, and the global best position gbest. The particles of the swarm make up a cloud that covers the whole search space in the initial iteration and gradually contracts its size as iterations advance to perform the exploration. So in the initial stages, the algorithm performs an exploration searching for plausible zones, and in the last iterations the best solution is improved[20]. In this paper, c1 and c2 are set to 2, maxgen is set to 120, sizepop is set to 20, ω is set to 1, and vmax is set to 5.

    • The PSO-AMLS-based method can be regarded as a combination of AP procedures, PSO technique and MLS fitting method. For the implementation of the PSO-AMLS method, firstly, the calibration data is pretreated using AP procedures; secondly, the MLS method is used for the fitting of pretreated calibration data, in which the PSO method is used to carry out the optimization mechanism to get the optical weighting function parameters in the MLS method. The flowchart of the PSO-AMLS-based method is shown in Figure 1. Through adaptive data processing in procedures (1) and (2), possible extreme value, non-normal distribution, or the dramatic change in the order of magnitude can be reduced or eliminated.

      Figure 1.  Flowchart of the hybrid PSO-AMLS-based method

      For basic functions selection in the MLS, lower-order polynomials are preferred. For parameters setting in Gaussian weighting function, larger influencing radius dmi means better fitting smoothness; larger shape parameter β means the locality is enhanced while the smoothness is declined. To optimize the MLS parameters, the PSO searches for the best parameters by comparing the fitness factor. The main fitness factor used in this paper is the determination coefficient R2, which is defined as follows:

      $$ {R^2} = \sum\limits_{i = 1}^n {{{\left( {{y_i} - {y_{if}}} \right)}^2}} $$ (12)

      where yi is the dependent variable of calibration point, yif is the fitting value, n is the number of calibration points.

    • In this Section, numerical examples are carried out to investigate the fitting performance of the PSO-AMLS-based method. Functions similar to the calibration data of the infrared radiometer are used in this section. In the first example, the independent variables are uniform distribution. In the second example, the independent variables are non-uniform distribution. The data were processed in this paper with code written in MATLAB.

    • The formula used in the first example is as follows:

      $$ f(x) = 1 + \sqrt x $$ (13)

      A set of points (xi, yi) are chosen, where i=0-20. The independent variables are uniform distribution from 0 to 20 with interval h=1, while the interval of dependent variables is gradually increased and determined by equation (13). The data used in the first example is shown in Table 1.

      Table 1.  The data used in the first example

      xyxyxy
      0173.6458144.7417
      1283.8284154.873
      22.414294165
      32.7321104.1623175.1231
      43114.3166185.2426
      53.2361124.4641195.3589
      63.4495134.6056205.4721

      The Gaussian function is chosen to be the weighting function. Considering the fitting accuracy and oscillation problem, the polynomial order is 6 for the LS method. The fitting performance is characterized by the fitting error δ between real value yi and fitting value yif

      $$ \delta = {y_i} - {y_{if}} $$ (14)

      The PSO is used to optimize the MLS parameters dmi, β, and m, where m represents the type of the basic function. The basic function is set to be p1(x) = [1], p2(x) = [1, x] and p3(x) = [1, x, x2] when m=1, m=2 and m=3, respectively. Search space is organized in three dimensions, one for each parameter. Because the independent variables are uniform distribution, there is no AP performed in this example. The fitting errors δ of the LS and PSO-AMLS method when the independent variables are uniform distribution are shown in Figure 2.

      Figure 2.  The fitting errors δ of the LS and PSO-AMLS method when the independent variables are uniform distribution: (a) nnodes=11; (b) nnodes=21

      The parameter “nnodes” is the number of nodes used in the PSO-AMLS method. For the cases when the independent variables of fitting data are uniform distribution, the PSO-AMLS method can acquire better results than the LS method. The PSO-AMLS method can follow the changes of the original function even with low order basis function, while, for the global approximation scheme like the LS, the oscillation phenomenon occurs, which increases the approximation error. When nnodes is 11, an outlier exists where δPSO-AMLS is larger than δLS. The position and amplitude of the outlier are x=1 and δPSO-AMLS=0.19. According to the amplitude δ, the results when nnodes is 21 are much better than the results when nnodes is 11, fewer nnodes means larger dmi is needed. When nnodes is too little, the fitting accuracy may be affected; when nnodes is too much, the regional feature may not be obvious but can be enhanced by the adjustment of parameters in the AMLS.

    • In this subsection, a set of points (xi, yi) are used as the calibration points, where i=1-15. The dependent variables are uniformly distributed with interval h=1, while the independent variables are non-uniformly distributed and the intervals are gradually increased. The data used in the second example is shown in Table 2.

      Table 2.  The data used in the second example

      yxyxyx
      11611.51146.5
      21.5716.51256.5
      32.5822.51367.5
      44.5929.51479.5
      57.51037.51592.5

      The parameters setting in the PSO and the AMLS are the same as in Subsection 3.1. Because the independent variables are non-uniform distribution, logarithmic transformation is performed before the MLS in this example. For the observed fitting results of the LS method are too bad when the independent variables are non-uniformly distribution, the Adaptive Least Squares (ALS) method is used in this subsection to compare with the PSO-AMLS method, where logarithmic transformation is performed before the LS. The fitting errors δ of the ALS and the PSO-AMLS method when the independent variables are non-uniform distribution are shown in Figure 3, in which the x-coordinate is the serial number of points.

      Figure 3.  The fitting errors δ of the ALS and the PSO-AMLS method when the independent variables are non-uniform distribution

      It is shown in Figure 3 that, for the cases when the independent variables of fitting data are non-uniform distribution, the PSO-AMLS method can acquire much better results than the LS method. The oscillation phenomenon occurs both for the ALS and the PSO-AMLS method, but the oscillation amplitude of the PSO-AMLS method is much smaller than that of the ALS method. Combined with the results obtained in Subsection 3.1, it is concluded that the fitting performance of the PSO-AMLS method is superior compared to the LS and ALS method.

    • In this Section, based on the infrared radiometer and the data fitting method described above, the calibration of the infrared radiometer has been carried out. The comparison of different data fitting methods is performed to validate the effectiveness of the PSO-AMLS method in the real calibration.

      The response voltages V of the infrared radiometer were measured in the calibration, where a blackbody-collimator was used as the standard source. Figure 4 is a diagram of the infrared radiometer. The basic components of the infrared radiometer are the optical lens, chopper, detector, preamplifier, lock-in amplifier, and A/D.

      Figure 4.  Diagram of the infrared radiometer

      The main specifications of the infrared radiometer are shown in Table 3. The basic components of the blackbody-collimator are the blackbody, collimator, and aperture. The main specifications of the blackbody-collimator are shown in Table 4. The blackbody temperatures T are uniform distribution from 100 ℃ to 700 ℃ with interval h=50 ℃. The corresponding calibration data of T and V are shown in Table 5.

      Table 3.  Main specifications of the infrared radiometer

      ComponentParameter
      Optical lensDiameter: 70 mm
      Field of view: ±0.6°
      Spectral channelsMW: 3.7-4.8 µm
      DetectorType: Insb
      A/DPrecision: 16
      Frequency: 100 kHz
      Preamplifier106×

      Table 4.  Main specifications of the blackbody-collimator

      ComponentParameter
      BlackbodyModel: HFY-200C
      Temperature: 5-1000 ℃
      ApertureDiameter: 0.1-13.5 mm
      CollimatorModel: HGD-1
      Focal length: 650 mm

      Table 5.  The calibration data of T and V in the calibration of the infrared radiometer

      V /VT /℃V /VT /℃
      0.021001.354450
      0.051501.828500
      0.1172002.383550
      0.2272502.994600
      0.3943003.75650
      0.6323504.555700
      0.949400

      It is shown from Table 5 that the dependent variables are uniform distribution with equal interval, while the independent variables are non-uniform distribution and the intervals increase gradually. The reasons for this phenomenon are as follows: (1) the response voltage of the infrared radiometer is affected significantly by the standard radiation source; (2) the response of the infrared radiometer is mainly determined by the detection ability of the detector.

      For the data characteristic mentioned above and the number limitation of calibration data, conventional data fitting methods commonly used in the calibration of infrared radiometer cannot achieve enough fitting accuracy. The fitting errors δ using the LS method and ALS method are shown in Figure 5. In Figure 5(a), polynomial approximate based on the LS method is carried out with polynomial orders from 5 to 6. In Figure 5(b), logarithmic transformation is carried out before the polynomial approximate with the polynomial orders from 5 to 6.

      Figure 5.  The fitting errors δ of the LS and ALS method: (a) LS method; (b) ALS method

      It can be obviously seen that the ALS method can acquire much better results than the LS method, which means the AP is effective to improve the fitting accuracy and stability when the fitting points are non-uniform distribution. But for the LS and ALS methods, the oscillation phenomenon occurs and the approximation errors change dramatically, especially when the number of calibration points is small and the fitting data are non-uniform distribution.

      To improve the fitting accuracy, the PSO-AMLS-based method is applied in the calibration of the infrared radiometer. The fitting errors of the PSO-AMLS method are compared to the ALS method to verify the performance of the PSO-AMLS method. For the PSO-AMLS with LT, only logarithmic transformation is contained in the AP, but for the PSO-AMLS with LT and DT, both logarithmic transformation and data translation are contained in the AP. The fitting errors of the PSO-AMLS and ALS method are shown in Figure 6.

      Figure 6.  The fitting errors of the PSO-AMLS and ALS method

      In Figure 6, the R2 of the PSO-AMLS method is 1.8437 compared with 3.6962 of the ALS method, which means better fitting results are achieved by the PSO-AMLS method. The superior performance of the PSO-AMLS method is more obvious in the first half of calibration points, while in the second half, an outlier exists in point 10 where δPSO-AMLS is slightly larger than δALS, and δPSO-AMLS and δALS are almost the same for point 12 and 13. It is supposed that the reason of this phenomenon is the number limitation of calibration points in the second half of the calibration data. Fitting errors δ will be increased when fewer calibration points are used. In the shape function matrix of the PSO-AMLS method, the locality is strong at the beginning, but decline with the increase of voltage. At the end of the calibration points, the locality is weakened but the globality is enhanced. In addition, it is shown that the performance of the PSO-AMLS with LT is only slightly better than the PSO-AMLS with LT and DT, which means the effect of logarithmic transformation is more obvious than data translation in the AP. The PSO-AMLS method should provide a better fitting result but suffer from the number limitation of fitting data, especially in the calibration of the infrared radiometer. As a future line of research, to improve the fitting performance when the samples of fitting data are fewer, more AP methods will be used.

    • In this section, the infrared radiometer calibrated using different data fitting methods is applied to the measurement of an infrared imaging simulator in the medium waveband. The measurement parameter is the equivalent blackbody temperature. The measured targets include the target channel and the interference channels. The source of target channel and interference channels are high-pressure xenon lamps, and the output radiation is controlled by the adjustment of the working current. The collimator and apertures used in the experiment are the same as in the calibration of the infrared radiometer. The equivalent blackbody temperature Tt measured by the thermal imaging camera ImageIR-5300 is used as the reference value. Tm measured by the infrared radiometer using two fitting methods is compared to the reference value Tt to obtain the deviation δ*:

      $$ {\delta ^*} = {T_t} - {T_m} $$ (15)

      The deviation δ* obtained for different fitting methods are shown in Figure 7, where TC represents Target Channel, IC 1, 2, and 3 represents Interference Channel 1, 2, and 3.

      It is shown from Figure 7 that for different measured targets, the measured Tm for the PSO-AMLS method is closer to the reference value compared to the ALS method. Besides, compared with the measurement results in section 3.1, more error sources are introduced in this measurement, including the repeatability and reproducibility, which results in the unexpected changes of the measured voltage. However, the most important result is that the PSO-AMLS method is more superior compared to other conventional fitting methods, which is effective for the improvement of measurement accuracy of the infrared radiometer.

      Figure 7.  The deviation δ* obtained for different fitting methods

    • The hybrid PSO-AMLS-based method was proposed for data fitting in this paper. The superior performance of the PSO-AMLS method was demonstrated by numerical examples and application in engineering. The main findings of this research work can be summarized as follows:

      (1) The hybrid PSO-AMLS-based method is successfully applied to data fitting in the calibration of the infrared radiometer. The fitting performance is confirmed by numerical examples and experiments. Compared to the conventional data fitting methods, better fitting results are achieved by the PSO-AMLS-based method;

      (2) From the obtained fitting results using the PSO-AMLS method, it is suggested that the use of the MLS approximation in combination with the PSO technique and AP process. In summary, the hybrid PSO-AMLS-based method is a promising method for data fitting, which can be applied to other calibration of instruments with success.

参考文献 (20)

目录

    /

    返回文章
    返回