振動模糊星圖的PSF函數(shù)估計及復(fù)原算法研究*

張　濤，王宏力，陸敬輝，姜　偉

（火箭軍工程大學(xué)控制工程系，西安710025）

振動模糊星圖的PSF函數(shù)估計及復(fù)原算法研究*

張濤，王宏力*，陸敬輝，姜偉

（火箭軍工程大學(xué)控制工程系，西安710025）

針對彈體角振動引起的星圖成像質(zhì)量退化的問題，提出了一種基于自適應(yīng)形態(tài)學(xué)PSF估計方法。首先通過研究角振動對星點像移的影響得到振動PSF模板，避免將模糊過程簡單近似為一維勻速直線運動和一維快速簡諧運動得到PSF；其次，為了消除或減弱彈體角振動引起的模糊，提出了一種改進約束條件的最小平方的星圖復(fù)原算法；最后進行仿真驗證，仿真結(jié)果表明本文方法能較好的估計隨機振動的PSF，提高了模糊星圖復(fù)原的精度和魯棒性，復(fù)原后的峰值信噪比和結(jié)構(gòu)相似比相比模糊星圖分別提高了23.3%和42.6%。

圖像處理；PSF估計；形態(tài)學(xué)；角振動；星圖復(fù)原

星敏感器作為一種具有高測量精度的空間姿態(tài)光學(xué)敏感器［1］，可為中遠程彈道導(dǎo)彈提高準確的空間方位和姿態(tài)基準［2］。對于捷聯(lián)安裝的星敏感器，由于其固連在導(dǎo)彈彈體上，在中制導(dǎo)實時調(diào)姿和變軌過程中會產(chǎn)生受迫振動，使得星圖模糊，星圖信噪比降低，進而影響星體質(zhì)心提取的精度［3］。因此如何復(fù)原振動模糊星圖就成為一個亟待解決的問題。

傳統(tǒng)的機械式、光學(xué)式、電子式像移補償法可以有效地抑制運動模糊，但存在體積大、成本高、結(jié)構(gòu)復(fù)雜、應(yīng)用不靈活等缺陷，而且難以適應(yīng)彈體高動態(tài)、隨機振動造成的模糊［4］。星圖復(fù)原技術(shù)通過對星圖進行后期處理消除振動模糊，且成本低，是解決振動模糊星圖問題的一個重要途徑。

星圖復(fù)原的關(guān)鍵技術(shù)主要有參數(shù)識別和復(fù)原濾波器設(shè)計兩項，前者利用不同先驗知識獲取適當(dāng)?shù)耐嘶瘏?shù)，即點擴展函數(shù)PSF，后者則根據(jù)這些參數(shù)估計真實星圖［5］。目前，對于PSF的估計，韓禮［4］、謝冰［6］等學(xué)者在研究中選用的模糊尺度和方向是通過計算機仿真得到的，實際應(yīng)用價值不高，且將模糊過程近似為一維勻速直線運動和一維快速簡諧運動來處理得到PSF，具有很大的局限性。對于復(fù)原濾波器的選擇，吳小娟［7］、段小慶［8］等分別引入維納濾波器和約束最小二乘濾波器，但維納濾波器工作時會增強低頻段中偏高的頻率，在視覺上表現(xiàn)為一些不必要的細節(jié)增強；而約束最小二乘濾波器工作時在高頻區(qū)Sn（u，v）/Sf（u，v）取較大的值，噪聲較強時復(fù)原效果較差。

針對上述分析中PSF估計和復(fù)原濾波器選擇存在的局限性，本文首先從振動模糊星圖本身得到PSF模板，提出了一種適用范圍更廣的自適應(yīng)形態(tài)學(xué)PSF估計方法，可以較好解決隨機振動過程中星圖退化的先驗知識無法獲得的難題；其次改進復(fù)原方法，提出一種改進約束條件的最小二乘星圖復(fù)原算法，可以達到更好的復(fù)原效果；最后，選取峰值信噪比和結(jié)構(gòu)相似性作為星圖復(fù)原程度的客觀評價參數(shù)，對復(fù)原效果進行了分析。

1　星圖振動PSF估計

1.1基本原理

自適應(yīng)形態(tài)學(xué)方法的基本原理是從星敏感器拍攝的星圖中估計振動PSF［9］。美國噴氣推進實驗室的C.C.Liebe提出了星敏感器針孔模型［10］，如圖1所示。從成像模型可得知，如果一顆恒星的成像記錄下了曝光時間內(nèi)的光軸的振動模糊信息，那么從星圖中恒星的成像就能精確獲得振動模糊PSF。

圖1中OPxPyP為像平面坐標系，OP為坐標原點，是星敏感器光學(xué)系統(tǒng)的光軸與圖像傳感器感光面的交點，xP軸與yP軸相互垂直，分別平行于感光面的行列方向。OSxSySzS為星敏感器坐標系，原點OS位于星敏感器光學(xué)系統(tǒng)的光心，zS軸為光軸指向，垂直于圖像傳感器的感光面，xS軸、yS軸分別于像平面坐標系的xP軸、yP軸平行，方向按右手定則確定。星敏感器捷聯(lián)安裝在導(dǎo)彈上，其安裝方式如圖2所示。

圖1　成像模型

圖2　安裝方式

圖2中Obxbybzb為彈體坐標系，星敏感器是安裝在彈體坐標系的yb軸上，且xS與zb軸平行、yS軸與xb平行、zS軸與yb軸重合。

由星敏感器的安裝方式以及導(dǎo)彈姿態(tài)角的定義可知［11］：當(dāng)導(dǎo)彈偏航振動Δψ時，相當(dāng)于星敏感器繞坐標系OSxSySzS的yS軸旋轉(zhuǎn)Δψ，這樣將使第i顆恒星在像平面上的成像Pi沿像平面坐標系的xP軸方向和yP軸方向移動至。假設(shè)像Pi在星敏感器坐標系OSxSySzS中的坐標為（xi，yi，-f），當(dāng)導(dǎo)彈偏航振動Δψ角后，像Pi和Pi′在星敏感器坐標系OSxSySzS中的坐標分別為（xiψ，yiψ，ziψ）和（x′iψ，y′iψ，z′iψ），則

因此，像偏移量為：

由式（3）可以看出，偏航角振動對星敏感器成像平面上星點的x軸、y軸方向上的坐標均產(chǎn)生影響。

當(dāng)導(dǎo)彈的俯仰角振動Δφ時，相當(dāng)于星敏感器繞坐標系OSxSySzS的zS旋轉(zhuǎn)Δφ角度。像Pi和在星敏感器坐標系 OSxSySzS中的坐標分別為（xiφ，yiφ，ziφ）和（x′iφ，y′iφ，z′iφ），則

因此，導(dǎo)彈的俯仰角Δφ引起的像Pi沿像平面坐標系的yP軸方向的像偏移量為：

由式（6）可以看出，俯仰角振動對星點的y軸方向上的坐標產(chǎn)生影響，而對x軸方向上的坐標不產(chǎn)生影響。

類似于俯仰角振動Δφ對像移的影響，可分析得滾動角振動Δγ只對星點的x軸方向上的坐標產(chǎn)生影響，而對y軸方向上的坐標不產(chǎn)生影響，則導(dǎo)彈的滾動角振動Δγ引起的像Pi沿像平面坐標系的xP軸方向的像偏移量為：

在曝光時間T0時刻，俯仰角抖動為φ（t），偏航角抖動為ψ（t），滾轉(zhuǎn)角抖動為γ（t）時，則偏移量為：

聯(lián)立式（3）、式（6）、式（7）得：

1.2提取算法

本文通過星敏感器拍攝的星圖，用自適應(yīng)形態(tài)學(xué)方法提取載體的PSF模板，較好的獲得復(fù)雜的星圖振動PSF。主要流程有星圖預(yù)處理、星圖二值化、形態(tài)學(xué)濾波（腐蝕和膨脹運算）、振動目標軌跡連通域標記和提取四大部分，提取算法的詳細步驟如下：

①星圖預(yù)處理是對原始星圖的數(shù)據(jù)線性化和灰度均衡化，以便為后續(xù)的星圖復(fù)原提供可靠量化和非灰度畸變的星圖信息。

②星圖二值化主要用來實現(xiàn)星點振動軌跡與背景噪聲分開，選用的方法是動態(tài)迭代閾值法［12］，二值化后用四連通判據(jù)（星點為凸形目標）去除星圖中大量單個或者少數(shù)像元構(gòu)成的噪點。

③利用形態(tài)學(xué)腐蝕和膨脹運算［13］消除二值化后星像內(nèi)部可能存在的隨機小洞和周邊可能存在的凹坑和凸刺，腐蝕運算在去除孤立噪聲的同時，星點的振動軌跡能夠更好保存下來，膨脹運算可以對腐蝕造成的軌跡失真和斷裂進行修補，對斷裂區(qū)域重新連接，確保振動軌跡所在星圖區(qū)域的完整。

④形態(tài)學(xué)濾波后將星圖中的連通區(qū)域初步確定為候選區(qū)并分別進行標記，由于成像噪聲、量化噪聲等會對PSF模板影響，因此通過比較各個標記窗口尺寸和能量信息，選取尺寸/能量最大的窗口作為振動PSF模板。

2　基于改進約束條件的最小平方復(fù)原方法

由文獻［6-7］可知，約束最小平方星圖復(fù)原方法中的二階微分對星圖中的孤立噪聲有較好的抑制作用，減小了H（u，v）對噪聲的敏感性，但二階微分對星圖本身的細節(jié)也有很強的響應(yīng)，在平滑星圖的同時也削弱了星圖的細節(jié)，基于上述原因，引入Sobel算子對振動模糊星圖進行一階微分處理。因為一階微分的約束條件將減少星圖細節(jié)的損失并抑制星圖模糊，彌補了單獨采用二階微分的約束最小平方復(fù)原算法的不足。算法表達式如式（10）所示：

3　復(fù)原星圖質(zhì)量評價方法

峰值信噪比PSNR（Peak Signal-to-Noise Ratio）和結(jié)構(gòu)相似比SSIM（Structure Similarity Ratio）是客觀評價星圖質(zhì)量最常用的評價標準之一［14］。PSNR 和SSIM的表達式如下：

其中 μ、σ和σxy分別為星圖的均值、標準差和協(xié)方差，C1、C2、C3為正實數(shù)常數(shù)。

4　實驗結(jié)果與對比分析

4.1實驗參數(shù)設(shè)定

星敏感器曝光時間為0.005 s，分辨率為1 024× 1 024，視場角為8°×8°，像元大小15 μm，視軸指向為（113°，30°），高斯半徑為5 pixel。振動模糊過程中加入均值為0，方差為0.001的高斯白噪聲，由于彈體上星敏感器的成像相機信噪比較高，Sn（u，v）/Sf是否精確對復(fù)原影響較小，因此k取值為0.001。

4.2仿真實驗

首先通過加入振幅為0.005 rad，頻率在20 Hz～2 000 Hz之間變化的隨機振動對標準星圖進行模糊，設(shè)定仿真星圖幀數(shù)為n=25得到星點振動軌跡，然后根據(jù)1、2節(jié)的PSF提取算法提取模糊星圖的PSF，實驗結(jié)果如圖3所示。

圖3　PSF提取過程

由圖3（f）標記的窗口可知，窗口3尺寸/能量最大，因此選擇其作為振動PSF模板。

利用圖4提取到的PSF模板，對經(jīng)過相同振動參數(shù)處理的模糊星圖進行復(fù)原，復(fù)原方法選擇改進約束條件的最小平方算法，當(dāng)?shù)螖?shù)k=87時得到圖5，從圖5可以看出振動模糊星圖得到了很好的復(fù)原。

圖4　PSF模板

圖5　復(fù)原星圖

此外還進行了對比實驗，對比組分別采用了常用的維納濾波法、最小平方法，由于三種復(fù)原方法均能使模糊星圖復(fù)原的很好，限于篇幅問題，就不將維納濾波法和最小平方法復(fù)原后的星圖在本文上體現(xiàn)。由于從主觀上無法評價這三種方法哪種方法復(fù)原效果更好，因此本文選用峰值信噪比PSNR和結(jié)構(gòu)相似比指數(shù)SSIM作為客觀評價指標，實驗數(shù)據(jù)統(tǒng)計如表1所示。

表1　實驗數(shù)據(jù)統(tǒng)計

從表1的實驗數(shù)據(jù)可以分析得到：振動會使得星圖的峰值信噪比和結(jié)構(gòu)相似比下降，而復(fù)原能較好的提高峰值信噪比和結(jié)構(gòu)相似比。維納濾波法、最小平方法以及本文方法這三種方法復(fù)原后得到的SSIM的值相差不大，但本文方法得到的值最高，比模糊星圖提高了42.6%；而復(fù)原得到的PSNR的值相差較大，比較結(jié)果為：維納濾波法＜最小平方法＜本文方法，本文方法相比模糊星圖提高了23.3%。綜上可知，本文選用的復(fù)原方法優(yōu)于維納濾波法和最小平方法。

5　結(jié)束語

本文通過研究彈體角振動對星點像移的影響得到載體振動PSF模板，并從模糊星圖中自適應(yīng)估計隨機振動的PSF，避免了將模糊過程簡單近似為一維勻速直線運動和一維快速簡諧運動得到PSF。改進約束條件的最小平方的星圖復(fù)原算法較好的減弱彈體角振動引起的模糊，這為下一步星圖識別以及姿態(tài)確定打下基礎(chǔ)，對提高星光制導(dǎo)的性能有著重要意義。

［1］陳雪芬，康國華.基于ITSI優(yōu)選算法的星敏感器導(dǎo)航星庫建立［J］.傳感技術(shù)學(xué)報，2014，28（3）：381-386.

［2］孫高飛，張國玉，劉石，等.高精度背景可控星圖模擬器設(shè)計［J］.紅外與激光工程，2015，44（7）：2196-2199.

［3］黃博.臨近空間高動態(tài)星圖仿真技術(shù)研究［D］.武漢：華中科技大學(xué)，2012.

［4］韓禮，張力軍，蔡洪，等.基于維納濾波的運動模糊星圖復(fù)原仿真分析［J］.光學(xué)技術(shù)，2014，40（4）：316-320.

［5］王輝，吳欽章.基于圖像質(zhì)量評價的自動圖像復(fù)原技術(shù)［J］.傳感技術(shù)學(xué)報，2012，25（7）：930-935.

［6］謝冰.基于TDICCD像移分析的PSF估計及圖像復(fù)原［D］.河北：燕山大學(xué)，2005.

［7］吳小娟，王新龍.星圖運動模糊及其復(fù)原方法［J］.北京航天航空大學(xué)學(xué)報，2011，37（11）：1338-1342.

［8］段小慶，王宏力.運動模糊圖像復(fù)原技術(shù)在星圖中的應(yīng)用［J］.微計算機信息，2008，24（11）：296-298.

［9］紀峰，李翠，常霞，等.基于改進閾值函數(shù)的自適應(yīng)圖像去噪方法［J］.傳感技術(shù)學(xué)報.2014，27（3）：351-354.

［10］Liebe C C.Accuracy Performance of Star Trackers-Atutorial［J］. IEEE Transactions on Aerospace Systems and Electronic System，2002，38（2）：587-589.

［11］王洪濤，羅長州，王渝，等.彈載星敏感器動態(tài)測角的分析［J］.現(xiàn)代防御技術(shù)，2010，29（1）：37-41.

［12］唐曉強.自適應(yīng)形態(tài)學(xué)圖像邊緣檢測算法的研究［D］.新疆：新疆大學(xué)，2010.

［13］顧茂松，楊小岡.自適應(yīng)邊緣檢測方法的研究［J］.計算機工程，2007，33（5）：185-186.

［14］鄧勤.衛(wèi)星顫振對遙感成像質(zhì)量影響研究及實驗驗證［D］.浙江：浙江大學(xué)，2013.

張濤（1992-），男，江西上饒人，碩士研究生，主要從事星敏感器測量精度研究；

王宏力（1965-），男，陜西鳳翔人，教授，博士，主要研究方向為導(dǎo)航制導(dǎo)與控制、復(fù)合制導(dǎo)、天文/慣性導(dǎo)航、故障診斷及可靠性評估，775128700@qq.com。

Research on PSF function Estimation and Restoration Algorithm for Vibration Fuzzy Star Image*

ZHANG Tao，WANG Hongli*，LU Jinghui，JIANG Wei
（Dept of Control Engineering，Rocket Force University of Engineering，Xi’an 710025，China）

This paper aim to solve the problem as to quality degradation of star mapping by adaptive morphology based PSF estimation method.Firstly，a vibration PSF template was obtained through studying the influence of angular vibration on star point displacement，the fuzzy process was avoided from being simply approximate to one-dimensional uniform linear motion and one-dimensional fast simple and harmonic motion，then PSF was obtained；secondly，a star image restoration algorithm of the least square based on improved constraint conditions was proposed to eliminate or weaken the fuzzy caused by angular vibration of the projectile body；finally，simulation verification was made.The result shows that the improved method could better estimate the randomly vibrating PSF and improve the accuracy and robustness of the restoration of fuzzy star image.Compared with the fuzzy star image，the peak signalto-noise ratio and the structure similarity ratio after restoration were improved by 23.3%and 42.6%，respectively.

image processing； PSF estimation； morphology； angular vibration； star map restoration

TP391.41

A

1004-1699（2016）08-1234-05

EEACC：6140；613510.3969/j.issn.1004-1699.2016.08.019

項目來源：陜西省科學(xué)技術(shù)研究發(fā)展基金項目（2014MJ2-6107）

2016-01-19修改日期：2016-03-29