Inertial stabilization technology in optical-electric tracking system
-
摘要: 在侦查探测、激光通讯等领域,光电跟踪系统的闭环精度是其重要技术指标之一。为了提高闭环精度,一般可使用图像稳定技术,惯性稳定技术或整体自稳定技术。惯性稳定技术因其良好的稳定效果,已在光电跟踪系统中得到广泛应用。采用对比分析的方法对光电跟踪系统中的机架惯性稳定、反射镜惯性稳定以及惯性基准光稳定技术进行了原理分析,优势比较以及发展展望,总结出多种惯性稳定技术交叉使用的复合轴惯性稳定仍是未来一段时间的发展趋势。Abstract: The closed-loop accuracy of the optical-electric tracking system is one of the important technical index in the fields of reconnaissance and detection, laser communication, etc. Researchers usually use image stabilization techniques, inertial stabilization techniques, or overall self-stabilization techniques to improve the closed-loop accuracy. Inertial stabilization techniques has been widely used in photoelectric tracking system for its good stabilization effect. This paper adopts the method of comparative analysis to analyze the principle, compare the advantages and forecast the development prospect of the frame inertial stabilization, mirror inertial stabilization and inertial reference light stabilization techniques in the photoelectric tracking system. It is proposed that the composite axis inertial stabilization using multiple inertial stabilization techniques is still a development tendency in the near future.
-
工业CT技术[1-3]不受被检测物体材料、形状、表面状况等限制,能够给出被检测物体二维、三维图像,成像直观,分辨率高。因此,工业CT被广泛应用在我国航空、航天、兵器、汽车制造、铁路、考古等领域,应用范围涵盖缺陷检测、材料密度表征、尺寸测量、装配结构分析、逆向工程等[4-10]。随着工业CT技术在高端装备制造业中的应用和发展,精密复杂零部件的内部缺陷检测需求日益增加,对缺陷的检出尺度及测量精度要求极高,在特定的CT设备和工艺条件下,获得的CT图像中,小缺陷往往被噪声所掩盖,难以进行有效识别和表征[11-13],直接影响工业CT缺陷测量的精度和准确性。因此,开展小缺陷工业CT无损检测方法研究是至关重要的,但当前相关技术研究主要集中在设备、算法和检测工艺等方面。缺陷尺度接近CT系统检出极限的内部小缺陷,在CT图像上表现出容积效应,或由于缺陷表面与CT扫描平面不垂直等因素,造成CT图像中缺陷边缘有较宽的灰度过渡区,形成了渐变边缘或称弱边缘,这时缺陷测量难度较高[14-20]。与此同时,缺陷类型、结构特点及伪影特性等多方面因素都对缺陷定量精度产生较大影响。
对于金属材料内部不规则形貌缺陷的工业CT无损检测,当前普遍采用半高宽缺陷尺寸计算法[21-23],CT检测人员通过观察及经验判定缺陷位置,人为分割缺陷边缘。该方法适用于尺度远大于CT图像像素尺寸的缺陷评价,但当缺陷尺寸接近或小于像素尺寸时,CT检测人员判定的缺陷测量及评价结果存在很大误差,由此会产生“超标”误判,造成不必要的浪费。因此,在金属材料工业CT检测过程中,有必要研究一种适用不同CT设备及工艺参数的缺陷精确测量方法。
本文针对金属材料内部缺陷CT检测的特点,通过引入基于CT图像的点扩散函数(PSF)计算模型[24-26],建立适用于金属材料内部缺陷的线阵工业CT测量模型,结合缺陷表面垂直度、缺陷大小等影响因素,对CT图像中缺陷边缘退化进行预测,利用当前图像PSF曲线进行缺陷的定量检测,以有效解决金属材料检测过程中小缺陷测量精度低的问题。
1. 缺陷评价方法
1.1 CT图像PSF测量模型
采用工业CT系统穿透被检测金属试样横截面,X射线透过衰减吸收被另一侧的线阵探测器接收,利用图像重建技术将内部结构及缺陷表现出来。工业CT图像中缺陷弱边缘是一种图像退化过程,该过程可以用原始图像与点扩散函数的卷积来描述,求解点扩散函数就是图像反卷积,将图像转换到频域上求解逆问题,考虑噪声影响,通过获取垂直于X射线平面的试样边缘图像灰度分布得到成像系统的近似点扩散函数(PSF)。
P(x,y)=FFT(g(x,y))−FFT(n(x,y))FFT(f(x,y))⋅eiφ (1) 式中:P(x, y)为系统PSF;f(x, y)为试样横截面图像;g(x, y)为CT扫描获得的退化图像;n(x, y)为加性噪声。
在特定工业CT系统下,电压、电流、切片厚度与切片位置等扫描参数确定后即可实现该函数的计算。将金属试样置于工业CT检测系统中,根据工业CT系统获取断层扫描图像直接计算得到系统近似点扩散函数。
在检测过程中,检测横截面应位于试样垂直界面中段,如图 1所示。若通过公式(1)直接计算系统调制传递函数方法,需要对噪声进行预测,该过程相当复杂。本文通过分段拟合统计一维边界函数方法可以较好地解决此问题。选取CT图像中一个长方形计算区域,使得试样边缘图像包含在该长方形区域中,对该区域进行二值化分割,在标准试块边缘上均匀提取N条等距离垂直于边缘的轨迹lN,各条轨迹上点的灰度值用一维数组表示,长度为n,则该一维数组lN(i), i∈(0, n)。将N条轨迹进行合并求平均,得到最初边缘响应函数(ERF):
E=N∑i=1lN(i)N (2) 式中:E为最初边缘响应函数(ERF)。
为了进一步降低噪声对计算的干扰,采用分段最小二乘近似方法。假设由公式(2)得到的E数组长度为n,设e[i, i+t]为E的索引从i到i+t子数组,其中t表示数组e的长度,一般取21为宜。f表示对e进行最小二乘法拟合后的一维数组。
{f[i]=argmin[i+t∑j=i(f[j,j+t]−e[j,j+t])2],i∈(0,n−t)Efit(i)=f(i+t2) (3) 式中:Efit(i)为拟合后的边缘响应函数(ERF)。
对ERF进行拟合求导,并归一化计算(PSF)曲线。
P(x,y)=dEfit(i)max(dEfit(i)) (4) 1.2 缺陷定量评价模型
在特定检测系统和特定工艺下,得到实际CT图像中缺陷双边缘(PSF)曲线后,将其用于金属材料内部缺陷的定量检测方法中。在CT图像上设置一条通过缺陷中心的直线,该直线包含缺陷两端边缘,获取该直线所在位置对应的灰度值l[n]。采用1.1中分段最小二乘法拟合并求导的方法计算缺陷双边缘退化后的PSF曲线,左边缘曲线和右边缘曲线分别设为Vl和Vr。根据CT图像实际PSF计算缺陷边缘特征函数,计算公式为
R=Vl(r)P(x,y) (5) 此时缺陷边缘特征函数R可以根据边缘类型分为三类:
(1) 垂直平面边界特征函数Rvertical,其垂直于CT扫描平面,图像中灰度分布表现为有且仅有一个最大值,并且最大值附近数据急剧下降,呈现尖峰状,判断准则为num{maxRvertical}=1,如图 2所示。
(2) 斜边缘特征函数Rslope;其倾斜于X射线扫描平面,图像中灰度分布表现为方波数据,有n个最大值并且最大值附近数据急剧下降,判断准则为num{maxRslope}=n,如图 3所示。
(3) 任意曲面(不规则面)的边界特征函数Rsurface,如图 4所示。
根据上述缺陷边缘类型进行金属材料内部缺陷的定量反演:首先对过缺陷中心的直线灰度值l[n]进行求导,根据极值所在点确定缺陷边缘两端初始位置,设为a, b。取a, b中间点a+b2,将l[n]一分为二。分开进行定量反演,设其中一段为L(w),其中w表示一维数组长度,假定一长度估计值x=x0,x0值由该边导数极值而定,其退化前的灰度幅值为H=max[L(w)]∑P。设边界类型数据为Xold(w),计算方法如下:
(1) 当边界类型为垂直平面时
Xold (w)=H×ε(t),ε(t)={0,t<a,t>a+x01,a⩽t⩽a+x0 (6) (2) 当边界类型判断为斜边缘时
Xold(w)=H×ε(t)+H×[1−t−a−x0num{maxRslope }]ε(t)={0,t<a,t>a+x01,a⩽t⩽a+x0 (7) (3) 当边界类型判断为任意曲面(不规则面)的边界时
求max(Rsurface),对应的索引为y,取x≥y的所有数据设为R′ surface[x]
Xold (w)=Hε(t1)+R′surface [y+t2−a−x0]ε(t2)ε(t1)={0,t1<a,t1>a+x01,a⩽t1⩽a+x0ε(t2)={0,t2⩽a+x01,t2>a+x0 (8) 将Xold(w)与标准PSF曲线进行卷积,计算与L(w)的标准偏差S(w),x递增循环迭代。获得标准偏差S(w)的曲线,求最小偏差Smin(w)时的x,即为最佳长度值,代入Xold(w)中,获得最佳重构曲线Xnew(w),处理流程如图 5所示。
根据像素对应实际距离进行尺寸测量。设CT实际成像范围u×u,单位(mm),图像尺寸v×v,单位(像素数)。CT图像像素物理距离为u/v。设测量曲线起点与终点分别为(x1, y1)和(x2, y2)。根据最佳重构曲线Xnew(w),其实际长度应该为中间值的长度k(单位:mm)
{N=num(Xnew(w)>p),p=max[Xnew(w)]+min[Xnew(w)]2k=uv×N×√(x2−x1)2+(y2−y1)2max(|x2−x1|,|y2−y1|)+s (9) 式中:s为系统误差。
2. 结果与分析
2.1 实验方案及结果
图 6所示为本文使用的6 MeV高能工业CT系统,主要性能参数包括6 MeV能量的直线加速器,采用线阵探测器作为信号采集装置,具有608个通道数,垂直准直器开口为0.3 mm,相邻通道间隔为1.3 mm,水平准直器可调范围为0.25~5 mm,采用三代扫描方式,剂量率800 cGy/min·m,焦点直径≤2 mm,空间分辨率约为2 lp/mm,密度分辨率0.3%~1%,重建图像为4096×4096像素,重建视场半径为300 mm,切片厚度为1 mm,可微动1~10次。
选取图 7所示对比试块作为具体实验对象。试块为圆形薄片,直径为50 mm、厚度为1.2 mm,材质为高碳高铬钢,密度为7.75 g/cm3。采用电火花手段加工一系列不同直径的通孔来模拟金属材料内部小缺陷,通孔直径d分别为0.1,0.2,0.3,0.4,0.5 mm。
利用高能工业CT系统获取断层扫描图像,计算边缘响应函数(ERF),得到系统近似点扩散函数(PSF),如图 8所示。
2.2 实验分析与验证
为验证以上所讨论的金属材料内部缺陷定量检测方法的有效性,将本文方法与文献[21]和[22]使用的传统半高宽法的计算结果进行比较,对比分析两方法曲线随缺陷尺寸和切片厚度的变化规律。图 9为两种方法在切片厚度为0.5 mm时不同尺寸缺陷的真实值/测量值比较曲线。从图 9中可以看出,在缺陷较大的情况下,PSF近似计算法和半高宽法的测量结果相对吻合;但当缺陷尺寸小于0.4 mm以后,半高宽法测量值较真值迅速增大,而本文方法的测量结果更接近真值。图 10所示为直径0.4 mm人工缺陷,在不同切片厚度下,采用两种方法获得的缺陷真实值/测量值比较曲线。可以看出,缺陷一定时,切片厚度较大的情况下,两种测量方法结果相对吻合;但当切片厚度小于0.5 mm以后,半高宽法测量值较真值偏差迅速增大,而本文方法的测量结果更接近真值。上述两种现象是由于当缺陷尺寸接近或小于有效射束宽度(BW)时,CT图像上缺陷呈现出的灰度变化规律与缺陷尺寸呈非线性趋势,故采用半高宽法测量值显然已与实际不相符。对比结果显示:本方法的测量精度优于传统半高宽法。
接下来通过观察缺陷尺寸为0.3 mm的灰度分布来分析PSF近似法的实际性能。如图 11所示,提取经过缺陷中心的线性灰度分布曲线图,本研究方法的缺陷近似计算结果与CT图像中缺陷大小相吻合,利用近似边缘退化缺陷灰度分析反求出缺陷的实际尺寸,能够较好地控制检测误差。
3. 结论
提出了基于工业CT测量模型的金属材料内部缺陷定量检测方法。该方法利用工业CT图像金属材料边缘近似计算图像PSF曲线,绘制不同缺陷边缘特征函数曲线,利用近似PSF曲线和缺陷边缘特征函数曲线反求出缺陷的实际尺寸。分析研究和人工缺陷的实验检测结果均验证了该方法对金属材料缺陷定量检测的有效性。相比较于半高宽缺陷测量方法,本研究方法对缺陷检测精度有较大的提高,尤其是在缺陷尺寸接近或小于有效射束宽度(BW)时,能够保证较高的缺陷测量精度。
-
表 1 国内多轴多框架光电转台研制情况
Table 1. Development of domestic multi-axis and multi-frame photoelectric turntables
structure type manufacturer weight/kg volume/mm×mm load function stabilization precision/µrad two-axis four-frame CIOMP 80 490×650 visible light detection
focal distance: 20~500 mm
infrared detection
FOV(8~12 µm)
14.40×10.80、2.40×1.80
laser ranging
wavelength: 1154 µm
Distance: 20 km25 two-axis four-frame 014 35 358×508 visible light detection
focal distance: 15~300 mm
infrared detection
FOV(8~12µm)
24×180、30×2.2040 two-axis two-frame 618 25 280×455 visible light detection
focal distance: 10~150 mm
three modes: color, high resolution, low light100 two-axis four-frame CIOMP 200 600×850 visible light detection
focal distance: 240~1000 mm
infrared detection
FOV(3~5 µm)
30×2.25020 two-axis four-frame CIOMP 10 260×420 visible light detection
focal distance: 5.4~72 mm
two modes: color, high resolution100 表 2 国外多轴多框架光电转台研制情况
Table 2. Development of foreign multi-axis and multi-frame photoelectric turntables
structure type manufacturer weight/kg volume/mm×mm load function stabilization precision/µrad two-axis four-frame USA
Westing house32 384×596 visible light detection
focal distance: 20~280 mm
infrared detection FOV:
7.50×9.70、2.250×2.90
laser ranging distance: 10 km25 two-axis four-frame Israel
TOPLITE53 406×662 visible light detection
focal distance: 20~240 mm
infrared detection
FOV (8~120 µm, 3~5 µm)
18×240、3.90×5.10
1.30×1.70
laser Ranging
wave length: 1.54 µm
distance: 20 km25 two-axis four-frame Italy
Astro30 380×596 infrared detection
FOV(8~12 µm)
40×2.70、16×10.70
visible light detection:
10 times zoom25 three-axis stabilization Canada
mescam30 356×548 visible light detection
focal distance: 16~160 mm
infrared detection FOV:
230×170、2.30×1.70
laser ranging distance: 10 km35 two-axis four-frame France
Scarnoff35 356×548 visible light detection FOV:
24×16、50×3.30
infrared detection FOV:
24×180、30×2.20
laser ranging distance: 10 km35 two-axis four-frame Russia
GS-257.5 340×552 visible light detection
infrared detection
laser ranging
laser irradiation50 OIRU-100 OIRU-500 OIRU-1000 jitter performance (1−1000 Hz) 100 nrad 500 nrad 2 μrad IKA performance 1.2 mrad at 20 min 2.2 mrad at 20 min 6 mrad at 20 min mechanism size 9.5″ dia. × 5.5″ high 7″ dia. × 5″ high 3.1″ dia. × 3.5″ high mechanism mass 18 lbs 10 lbs 3 lbs electronics mass 2 lbs 2 lbs 2 lbs power 45 W 30 W 18 W gyro (3) Emcore 1.3k FOG (3) KVH DSP-1750D FOG (1) 2-Axis NG G2000 DTG ARS (3) ARS-16 (3) ARS-16 (3) ARS-15 acceleration 2 rad/s2 2 rad/s2 10 rad/s2 bandwidth (open-loop crossover) >100 Hz >100 Hz >150 Hz position resolution <1 μrad <1 μrad <1 μrad controller type digital (FPGA) digital (FPGA) digital (FPGA) -
[1] 胡浩军. 运动平台捕获、跟踪与瞄准系统视轴稳定技术研究[D]. 长沙: 国防科学技术大学, 2005: 12-15Hu Haojun. Line-of-sight stabilization of acquisition, tracking and pointing system on moving bed[D]. Changsha: National University of Defense Technology, 2005: 12-15 [2] Satyarthi S. Optical line-of-sight steering using gimbaled mirrors[C]//Proceedings of the SPIE 9076, Airborne Intelligence, Surveillance, Reconnaissance (ISR) Systems and Applications XI. 2014: 90760E. [3] 吕宏宇, 金刚石, 高旭辉. 两轴四框架机载光电平台稳定原理分析[J]. 激光与红外, 2015, 45(2):194-198. (Lü Hongyu, Jin Gangshi, Gao Xuhui. Stabilization analysis of airborne electro-optical platform with two-axis and four-gimbal[J]. Laser & Infrared, 2015, 45(2): 194-198 doi: 10.3969/j.issn.1001-5078.2015.02.017Lü Hongyu, Jin Gangshi, Gao Xuhui. Stabilization analysis of airborne electro-optical platform with two-axis and four-gimbal[J]. Laser & Infrared, 2015, 45(2): 194-198 doi: 10.3969/j.issn.1001-5078.2015.02.017 [4] 唐涛, 马佳光, 陈洪斌, 等. 光电跟踪系统中精密控制技术研究进展[J]. 光电工程, 2020, 47:200315. (Tang Tao, Ma Jiaguang, Chen Hongbin, et al. A review on precision control methodologies for optical-electric tracking control system[J]. Opto-Electronic Engineering, 2020, 47: 200315Tang Tao, Ma Jiaguang, Chen Hongbin, et al. A review on precision control methodologies for optical-electric tracking control system[J]. Opto-Electronic Engineering, 2020, 47: 200315 [5] 侯瑞博, 魏涛, 宋景. 航空光电侦察平台关键技术及其发展[J]. 电子元器件与信息技术, 2019(3):51-53,57. (Hou Ruibo, Wei Tao, Song Jing. Key technologies and its development of aeronautical photoelectric reconnaissance platform[J]. Electronic Component and Information Technology, 2019(3): 51-53,57 doi: 10.19772/j.cnki.2096-4455.2019.3.014Hou Ruibo, Wei Tao, Song Jing. Key Technologies and its development of aeronautical photoelectric reconnaissance platform[J]. Electronic Component and Information Technology, 2019(3): 51-53, 57 doi: 10.19772/j.cnki.2096-4455.2019.3.014 [6] 洪华杰, 王学武, 翁干飞. 光电侦察装备中的反射镜稳定技术[J]. 应用光学, 2011, 32(4):591-597. (Hong Huajie, Wang Xuewu, Weng Ganfei. Mirror stabilization in electro-optical reconnaissance system[J]. Journal of Applied Optics, 2011, 32(4): 591-597 doi: 10.3969/j.issn.1002-2082.2011.04.001Hong Huajie, Wang Xuewu, Weng Ganfei. Mirror stabilization in electro-optical reconnaissance system[J]. Journal of Applied Optics, 2011, 32(4): 591-597 doi: 10.3969/j.issn.1002-2082.2011.04.001 [7] Hilkert J M, Cohen S. Development of mirror stabilization line-of-sight rate equations for an unconventional sensor-to-gimbal orientation[C]//Proceedings of the SPIE 7338, Acquisition, Tracking, Pointing, and Laser Systems Technologies XXIII. 2009: 733803. [8] 王琦, 孙广利, 黎纯宁, 等. 基于半捷联方式的反射镜视轴稳定技术[J]. 红外与激光工程, 2015, 44(10):3070-3075. (Wang Qi, Sun Guangli, Li Chunning, et al. Inertial line-of-sight stabilization technique of semi-strapdown control using mirrors[J]. Infrared and Laser Engineering, 2015, 44(10): 3070-3075 doi: 10.3969/j.issn.1007-2276.2015.10.035Wang Qi, Sun Guangli, Li Chunning, et al. Inertial line-of-sight stabilization technique of semi-strapdown control using mirrors[J]. Infrared and Laser Engineering, 2015, 44(10): 3070-3075 doi: 10.3969/j.issn.1007-2276.2015.10.035 [9] 宋江鹏, 孙广利, 周荻, 等. 反射镜光电平台视轴稳定技术研究[J]. 红外与激光工程, 2015, 44(6):1904-1911. (Song Jiangpeng, Sun Guangli, Zhou Di, et al. Line-of-sight stabilization techniques for mirror electro-optical platform[J]. Infrared and Laser Engineering, 2015, 44(6): 1904-1911 doi: 10.3969/j.issn.1007-2276.2015.06.041Song Jiangpeng, Sun Guangli, Zhou Di, et al. Line-of-sight stabilization techniques for mirror electro-optical platform[J]. Infrared and Laser Engineering, 2015, 44(6): 1904-1911 doi: 10.3969/j.issn.1007-2276.2015.06.041 [10] Xia Yunxia, Bao Qiliang, Liu Zidong. A new disturbance feedforward control method for electro-optical tracking system line-of-sight stabilization on moving platform[J]. Sensors, 2018, 18: 4350. doi: 10.3390/s18124350 [11] Schneeberger T J, Barker K W. High-altitude balloon experiment: a testbed for acquisition, tracking, and pointing technologies[C]//Proceedings of the SPIE 1950 Acquisition, Tracking, and Pointing VII. 1993: 2-15. [12] Luniewicz M F, Murphy J, O'Neil E, et al. Testing the inertial pseudo-star reference unit[C]//Proceedings of the SPIE 2221, Acquisition, Tracking, and Pointing VIII. 1994: 638-649. [13] Eckelkamp-Baker D, Sebesta H R. Optical inertial reference unit for kilohertz bandwidth submicroradian optical pointing and jitter control: U. S. Patent 7227, 111[P]. 2007-06-05. [14] Gilmore J P, Luniewicz M F, Sargent D. Enhanced precision pointing jitter suppression system[C]//Proceedings of the SPIE 4632, Laser and Beam Control Technologies. 2002: 38-49. [15] Walter R E, Danny H, Donaldson J. Stabilized inertial measurement system (SIMS)[C]//Proceedings of the SPIE 4724, Laser Weapons Technology III. 2002: 57-68. [16] ATA(Applied Technology Associate). Optical inertial reference unit[EB/OL]. (2019-01-01)[2021-10-14].https://bluehalo.com/product/optical-inertial-reference-unit-oiru/. [17] Yue Ronggang, Wang Humei, Jin Ting, et al. Image motion measurement and image restoration system based on an inertial reference laser[J]. Sensors, 2021, 21: 3309. doi: 10.3390/s21103309 期刊类型引用(9)
1. 郭成龙,倪培君,齐子诚,付康. 基于斜孔散射校正板的锥束X射线CT散射校正方法. 强激光与粒子束. 2024(07): 80-88 . 本站查看
2. 徐凤丽,倪培君,齐子诚,张维国,付康. 活塞波浪形内冷油道位置的工业CT测量. 无损检测. 2022(08): 31-35+64 . 百度学术
3. 齐子诚,倪培君,张维国,姜伟,裘信国. 工业CT检测中小缺陷定量方法. 科学技术与工程. 2021(03): 958-964 . 百度学术
4. 贾晓亮,秦怡. 煤矿坑道钻探钻杆断裂智能CT检测研究. 煤炭科学技术. 2021(03): 135-140 . 百度学术
5. 刘丽东,钱承,倪培君,潘锋,郭淼. 无损检测新技术能力实验室认可现状与展望. 无损检测. 2021(09): 39-44 . 百度学术
6. 金献忠,陈建国,林力,张建波,付冉冉,陈颖娜,刘永丰. 激光剥蚀电感耦合等离子体质谱在金属材料缺陷分析中的应用. 电镀与涂饰. 2020(10): 665-669 . 百度学术
7. 齐子诚,张维国,唐盛明,郑颖. 基于CT图像变异系数的小缺陷定量不确定度评定. 计算机测量与控制. 2020(12): 262-266 . 百度学术
8. 朱凯然,王波,武晓朦,高建申. 原位岩土CT及其关键技术探讨. CT理论与应用研究. 2019(06): 641-652 . 百度学术
9. 葛浙东,戚玉涵,罗瑞,陈龙现,王艳伟,周玉成. 木材CT断层成像系统旋转中心校正方法. 林业科学. 2018(11): 164-171 . 百度学术
其他类型引用(8)
-