專利名稱：高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取方法及系統(tǒng)的制作方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取方法及系統(tǒng)。
背景技術(shù)：
衛(wèi)星遙感具有覆蓋面大、持續(xù)時間長、實(shí)時性強(qiáng)、不受國界、地域限制等獨(dú)特優(yōu)勢，廣泛地應(yīng)用于資源開發(fā)、環(huán)境監(jiān)測、災(zāi)害研究、全球變化分析等領(lǐng)域，深受各國的高度重視。衛(wèi)星圖像的空間分辨率是衡量衛(wèi)星遙感能力的ー項主要指標(biāo)，也是衡量ー個國家航天遙感水平的重要標(biāo)志。提高衛(wèi)星觀測空間分辨率已成為衛(wèi)星工程技術(shù)研究熱點(diǎn)。衛(wèi)星獲取圖像過種中，有很多因素會導(dǎo)致圖像質(zhì)量下降，大氣擾動、運(yùn)動、散焦、傳輸和噪聲都會直接影響到圖像分辨率下降，特別是衛(wèi)星截荷所限使光學(xué)系統(tǒng)截止頻率有限高，及其CCD芯片像元尺寸有限小，限止了衛(wèi)星圖像空間高頻分量，使得圖像分辨率不夠高。根據(jù)光學(xué)傅里葉頻譜理論，光學(xué)系統(tǒng)存在截止頻率Cf=(D-1) (2fX)，其中D為等效透鏡直徑，I為C⑶芯片尺寸，f為透鏡焦距和入為光波長。若CXD芯片的像元尺寸為w，則按采樣定理，也有截止頻率為uw=l (2w)。被攝物中，只有同時小于Uw和Cf的空間頻率分量才能獲取并成像，若Cfデuw,則導(dǎo)致采樣資源或光學(xué)成像資源浪費(fèi)。設(shè)衛(wèi)星與被拍攝物的距離為R，則衛(wèi)星圖像的可分辨距離Ax=wR/f=AR(D_l)。若減少CXD芯片的像元尺寸提高IV并且光學(xué)截止頻率也隨之提高cf=uw，則可以提高圖像的空間分辨率(目前最小值為50 y m2)，但是CXD芯片像元尺寸w太小，信噪比太低，以至無法正常使用。因此，衛(wèi)星圖像的高頻分量缺失是ー個繞不開的科學(xué)難題。傳統(tǒng)的高分辨率(參考文獻(xiàn)1:J.L.Harris, Diffraction and resolving power, J.0pt.Soc.Amer., 54 (7): 931-133, 1964和又獻(xiàn)2:W.Lukosz, Optical systems with resolving power exceeding the classicallimit.J.0pt.Soc.Amer.，56 (11): 1463-1471，1966)是指對超出光學(xué)系統(tǒng)截止頻率 Cf 而被丟失的圖像高頻信息進(jìn)行恢復(fù)，這種方法稱為高分辨率復(fù)原技木。多數(shù)人認(rèn)為要準(zhǔn)確恢復(fù)截止頻率之外的頻譜信息是不可能的，并稱此為高分辨率神話。

當(dāng)可以獲得多幅同一場景的圖像序列時，可以建立數(shù)學(xué)模型:Si=Hsi+!^ i=l, 2，...，k，其中g(shù)i，Si, Hi分別表第i幀低分辨率圖像、高分辨率圖像和噪聲圖像，H表示各種導(dǎo)致圖像低分辨率的各種因素。通過多幀插值方法I# 見又獻(xiàn) 3:L.Zhang, X.ffu, An Edge-Guidea Image Interpolation Algorithmvia Directional Filtering and Data Fusion, IEEE Transaction on imageprocessing,15(8):2226-2238, 2006,文獻(xiàn) 4:D.Rajan D,S.Chaudhuri, Generalizedinterpolation and its application in super-resolution imaging, Imageand Vision Computing, 19(13):957-969, 2001,文獻(xiàn) 5:A.Sdnchez-Beato andG.Pajares, Non-1terative interpolation-based super-resolution minimizingaliasing in the reconstructed image, IEEE Trans.1mage Process., 17(10), pp.1817 - 1826,2008,文獻(xiàn) 6:S.Lertrattanapanich, N.K.BOSE, High resolution imageformation from low resolution frames using Delaunay triangulation, IEEETransaction on Image Processing，11(12):1427-1441，2002 和文獻(xiàn) 7F.Zhou，W.Yang,and Q.Liao, Interpolation-Based Image Super-Resolution Using MultisurfaceFitting, IEEE Transaction on Image Processing, 21 (7): 3312-28，2012)、利用先驗知識進(jìn)行優(yōu)化求解(參見文獻(xiàn) 8:X.Gaoj K.Zhang, D.Tao and X.Li，Image Super-ResolutionWith Sparse Neighbor Embedding,I IEEE Trans.0n Image Processing, Vol.21, No 7，pp.3194-3205，2012，文獻(xiàn) 9:Z.M.Wang and W.W.Wang, Fast and Adaptive Methodfor SAR Superresolution Imaging Based on Point Scattering Model and OptimalBasis Selection, IEEE Tran, on Image Processing，18(7):1477-1486，2009，文獻(xiàn) 10:A.Marquina and S.j.0sher，Image super-resolution by TV regularization andBregman iteration,Journal of Scientific Computing, vol.37，n0.3，pp.367 - 382，2008 和文獻(xiàn) 11:J.Yang，J.Wright，T.S.Huang and Y.Ma，“Image super-resolutionvia sparse representation, ”IEEE Trans.1mage Process.，vol.19，n0.11，pp.2861-2873，2010)、基于學(xué)習(xí)方法(參見文獻(xiàn) 11，文獻(xiàn) 12:T.Goto，Y.Kawamoto，Y.Sakuta，A.Tsutsui and M.^akuraij Learning-based Super-resolution Image Reconstruction onMult1-core Processor, IEEE Transactions on Consumer Electronics, 58(3):941 946，2012，文獻(xiàn) 13:P.Purkait and B.Chanda, Super Resolution Image ReconstructionThrough Bregman Iteration Using Morphologic Regularization,IEEE Trans.0nImage Processing, 21 (9):4029 4040，2012，文獻(xiàn) 14:P.P.Gajjar and M.V.Joshi，Newlearning based super-resolution:Use of DWT and IGMR-F prior,IEEE Trans,on Image Processing, Vol.19，N0.5，pp.1201-1213，2010，文獻(xiàn) 15:M.S.Crouse, R.D.Nowak, R.G Baraniuk, Wavelet-based statistical signal processing using niddenMarkov models,IEEE Transactions on Signal Processing, 46(4):886-902,1998 和文獻(xiàn) 16:M N Do, M.Vetterli, The contourlet transform: An efficient directionalmult1-resolution image representation,IEEE Transactions on Image Processing，14(12):2091-2106, 2005)等等(參見文獻(xiàn) 17:D.D.-Y Po and DO M.N.Do，Directionalmult1-scale modeling of images using the contourlet transform, IEEE Transactionson Image Processing, 15(6): 1610-1620，2006 和文獻(xiàn) 18:W.Dong, L Zhang, G.Shi andX.Wuj Image deblurring and superresolution by adaptive sparse domain selectionand adaptive regularization, IEEE Trans.1mage Process.，vol.20，n0.7，pp.533 -549，Jul.2011)獲取高分辨率圖像，改善欠釆樣而導(dǎo)致的圖像質(zhì)量退化。但在衛(wèi)星觀測攝像中，同一視場的多幀圖像釆集是極浪費(fèi)資源且難做到的，單幀圖像超分辨復(fù)原技術(shù)才是遙感圖像超分辨復(fù)原的關(guān)鍵技術(shù)，但至今沒有實(shí)質(zhì)性突破。為研究單幀圖像的高分辨率復(fù)原，我們考慮以下低分辨率圖像形成機(jī)制的數(shù)學(xué)模型.
若不考慮干擾因素，對于ー個成像系統(tǒng)其成像過程可用下式加以描述為:g (X，y) =p (X，y) *s (x, y)這里g(x，y)，p(x，y)，s(x，y)分別表示遙感圖像，視場和成像系統(tǒng)點(diǎn)擴(kuò)散函數(shù)，*表示卷積運(yùn)算。對圖像頻譜函數(shù)為:G(u，v)=P(u，v)S(u，V)，點(diǎn)擴(kuò)散函數(shù)的頻譜函數(shù)P(u，v)=l，|u|〈cf&|v|〈cf是有限帶寬的矩形窗。當(dāng)實(shí)際視場截止頻率Cs大于光學(xué)成像系統(tǒng)截止頻率Cf時，視場的Cf以外的高頻分量丟失，成為低分辨率圖像。傳統(tǒng)習(xí)慣認(rèn)為在光學(xué)成像系統(tǒng)截止頻率之外Cf頻譜是不可以復(fù)原的。但根據(jù)解析延拓定理:若解析在某一有限區(qū)間上為已知，就可以唯一地延拓到全部區(qū)域。也就說〃如果兩個解析函數(shù)在任一給定的區(qū)域上完全一致〃則它們一定在整體上完全一致〃即為同一函數(shù)(參見文19:E.B.Saff andA.D.Snider,Fundamentals of Complex Analysis with Applications to Engineeringand Science，2003，Pearson Education)。視場可以看作是ー個有界定義域上的函數(shù)，其譜函數(shù)是ー個解析函數(shù)。因此，按照解析延拓定理，可以由圖像頻譜數(shù)據(jù)G(u，v)=P(u，v)S(u，v), |u|〈cf&|v|〈cf,延拓到整個頻譜空間，截止頻率Cf= 早期研究從單巾貞圖像進(jìn)行高分辨率復(fù)原主要方法是頻譜外推(參見文獻(xiàn)20:H.Greenspan, C.H.Anderson, S.Akber, Image enhancement by nonlinear extrapolation in frequency space,IEEETrans.1m age Processing，vol.9，no6，pp.1035—1048,2000)，長捕球波函數(shù)法(參見文獻(xiàn)21:H.A.Brown, Effect of Truncation on Image Enhancement by Prolate SpheroidalFunctions,Journal of the Optical Society of America, Vol.59，no2，pp.228-229，1969)，迭加正弦模板法(參見文獻(xiàn) 22:S.Wadaks，T.Sato, Superresolution in IncoherentImaging System, Journal of the Optical Society of America，65(3):354-355，1975)、插值方法(參見文獻(xiàn)3)等超分辨復(fù)原技木。但是這些方法充分利用低分辨率圖像中隱含著圖像高分辨率信息，沒有理解和運(yùn)用解析延拓定理數(shù)學(xué)原理，從而無法或很少研究從低分辨圖像獲取高頻信息方法，因而效果十分有限(參見文獻(xiàn)23:S.C.Park, M.K.Park, M.Kj.Kang, buper-resolution image reconstruction: a technical overview, IEEE bigna丄Processing Magazine, Vol.20，n0.3，pp.21-36，May2003)o

發(fā)明內(nèi)容
本發(fā)明的目的在于提供一種高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取方法及系統(tǒng)，能夠在高頻頻譜數(shù)據(jù)缺失的情況下，快速高效地獲取高分辨率圖像復(fù)原的坐標(biāo)參數(shù)以供高分辨率圖像的復(fù)原。為解決上述問題，本發(fā)明提供一種高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取方法，包括:獲取待復(fù)原的高分辨率圖像的橫向或縱向像素點(diǎn)個數(shù)和一幅低分辨率圖像，根據(jù)所述像素點(diǎn)個數(shù)和低分辨率圖像獲取所述高分辨率圖像的低頻頻譜數(shù)據(jù)；根據(jù)所述高分辨率圖像的低頻頻譜數(shù)據(jù)獲取所述高分辨率圖像的補(bǔ)零法頻譜數(shù)據(jù)；對所述補(bǔ)零法頻譜數(shù)據(jù)作傅里葉變換以獲取高分辨率圖像的低頻頻譜數(shù)據(jù)補(bǔ)零法圖像；根據(jù)所述低頻頻譜數(shù)據(jù)補(bǔ)零法圖像獲取最佳奇異化算子，根據(jù)所述最佳奇異化算子獲取奇異函數(shù)，根據(jù)所述奇異函數(shù)獲取奇異譜函數(shù)；根據(jù)所述最佳奇異化算子并運(yùn)用點(diǎn)擴(kuò)散函數(shù)層析法獲取高分辨率圖像復(fù)原的坐標(biāo)參數(shù)。進(jìn)ー步的，在上述方法中，獲取待復(fù)原的高分辨率圖像的橫向或縱向像素點(diǎn)個數(shù)和一幅低分辨率圖像，根據(jù)所述像素點(diǎn)個數(shù)和低分辨率圖像獲取所述高分辨率圖像的低頻頻譜數(shù)據(jù)的步驟中，
一幅低分辨率圖像gl (i, j)，i, j=0, 1，...，I要復(fù)原到高分辨率圖像g(i，j)，i，j=0，l，...，N，N>>l，g(i，j)圖像的頻譜數(shù)據(jù)表示為6(1^，1^)，1^1^￡ Q, Q 為所述高分辨率圖像的頻譜空間，I表示低分辨率圖像的橫向或縱向像素點(diǎn)個數(shù)，N表示待復(fù)原的高分辨率圖像的橫向或縱向像素點(diǎn)個數(shù)，低分辨率圖像的頻譜數(shù)據(jù)表示為G1G^ ky)，
g(i，j)圖像的低頻頻譜數(shù)據(jù)表示為
權(quán)利要求
1.一種高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取方法，其特征在于，包括: 獲取待復(fù)原的高分辨率圖像的橫向或縱向像素點(diǎn)個數(shù)和一幅低分辨率圖像，根據(jù)所述像素點(diǎn)個數(shù)和低分辨率圖像獲取所述高分辨率圖像的低頻頻譜數(shù)據(jù)； 根據(jù)所述高分辨率圖像的低頻頻譜數(shù)據(jù)獲取所述高分辨率圖像的補(bǔ)零法頻譜數(shù)據(jù)； 對所述補(bǔ)零法頻譜數(shù)據(jù)作傅里葉變換以獲取高分辨率圖像的低頻頻譜數(shù)據(jù)補(bǔ)零法圖像； 根據(jù)所述低頻頻譜數(shù)據(jù)補(bǔ)零法圖像獲取最佳奇異化算子，根據(jù)所述最佳奇異化算子獲取奇異函數(shù)，根據(jù)所述奇異函數(shù)獲取奇異譜函數(shù)； 根據(jù)所述最佳奇異化算子并運(yùn)用點(diǎn)擴(kuò)散函數(shù)層析法獲取高分辨率圖像復(fù)原的坐標(biāo)參數(shù)。
2.如權(quán)利要求1所述的高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取方法，其特征在于，獲取待復(fù)原的高分辨率圖像的橫向或縱向像素點(diǎn)個數(shù)和一幅低分辨率圖像，根據(jù)所述像素點(diǎn)個數(shù)和低分辨率圖像獲取所述高分辨率圖像的低頻頻譜數(shù)據(jù)的步驟中， 一幅低分辨率圖像gl(i，j),i, j=0, 1，...，I要復(fù)原到高分辨率圖像g(i，j)，i，j=0，l，...，N，N>>l，g(i，j)圖像的頻譜數(shù)據(jù)表示為6(1^，1^)，1^1^￡ Ω，Ω 為所述高分辨率圖像的頻譜空間，I表示低分辨率圖像的橫向或縱向像素點(diǎn)個數(shù)，N表示待復(fù)原的高分辨率圖像的橫向或縱向像素點(diǎn)個數(shù)，低分辨率圖像的頻譜數(shù)據(jù)表示為G1G^ ky)，
3.如權(quán)利要求2所述的高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取方法，其特征在于，根據(jù)所述高分辨率圖像的低頻頻譜數(shù)據(jù)獲取所述高分辨率圖像的補(bǔ)零法頻譜數(shù)據(jù)的步驟中， 所述高分辨率圖像的補(bǔ)零法頻譜數(shù)據(jù)表示為G(kx，ky)P(kx, ky)，其中，
4.如權(quán)利要求3所述的高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取方法，其特征在于，對所述補(bǔ)零法頻譜數(shù)據(jù)作傅里葉變換以獲取高分辨率圖像的低頻頻譜數(shù)據(jù)補(bǔ)零法圖像的步驟中， 所述高分辨率圖像的低頻頻譜數(shù)據(jù)補(bǔ)零法圖像表示為
5.如權(quán)利要求4所述的高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取方法，其特征在于，根據(jù)所述低頻頻譜數(shù)據(jù)補(bǔ)零法圖像獲取最佳奇異化算子的步驟包括: 初始化:
6.如權(quán)利要求5所述的高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取方法，其特征在于，根據(jù)所述最佳奇異化算子獲取奇異函數(shù)的步驟包括: 根據(jù)差分方程Φ (i，j)*h(i，j)=δ (i，j)的零狀態(tài)的解獲取奇異函數(shù)h(i，j)。
7.如權(quán)利要求6所述的高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取方法，其特征在于，根據(jù)所述奇異函數(shù)獲取奇異譜函數(shù)的步驟中，奇異譜函數(shù)為
8.如權(quán)利要求7所述的高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取方法，其特征在于，根據(jù)所述最佳奇異化算子并運(yùn)用點(diǎn)擴(kuò)散函數(shù)層析法獲取高分辨率圖像復(fù)原的坐標(biāo)參數(shù)的步驟包括:初始化
9.一種高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取系統(tǒng)，其特征在于，包括: 低頻頻譜數(shù)據(jù)模塊，用于獲取待復(fù)原的高分辨率圖像的橫向或縱向像素點(diǎn)個數(shù)和一幅低分辨率圖像，根據(jù)所述像素點(diǎn)個數(shù)和低分辨率圖像獲取所述高分辨率圖像的低頻頻譜數(shù)據(jù)； 補(bǔ)零法頻譜數(shù)據(jù)模塊，用于根據(jù)所述高分辨率圖像的低頻頻譜數(shù)據(jù)獲取所述高分辨率圖像的補(bǔ)零法頻譜數(shù)據(jù)； 補(bǔ)零法圖像模塊，用于對所述補(bǔ)零法頻譜數(shù)據(jù)作傅里葉變換以獲取高分辨率圖像的低頻頻譜數(shù)據(jù)補(bǔ)零法圖像； 奇異譜函數(shù)模塊，用于根據(jù)所述低頻頻譜數(shù)據(jù)補(bǔ)零法圖像獲取最佳奇異化算子，根據(jù)所述最佳奇異化算子獲取奇異函數(shù)，根據(jù)所述奇異函數(shù)獲取奇異譜函數(shù)； 坐標(biāo)參數(shù)模塊，用于根據(jù)所述最佳奇異化算子并運(yùn)用點(diǎn)擴(kuò)散函數(shù)層析法獲取高分辨率圖像復(fù)原的坐標(biāo)參數(shù)。
10.如權(quán)利要求9所述的高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取系統(tǒng)，其特征在于，所述低頻頻譜數(shù)據(jù)模塊用于將一幅低分辨率圖像表示為gl(i，j),i, j=0, I,..., 1，將要復(fù)原到高分辨率圖像表示為g(i, j), I, J-O, I,..., N, N > > I, g(i, j)圖像的頻譜數(shù)據(jù)表示為G(kx, ky)，kx，ky e Ω，Ω為所述高分辨率圖像的頻譜空間，I表示低分辨率圖像的橫向或縱向像素點(diǎn)個數(shù)，N表示待復(fù)原的高分辨率圖像的橫向或縱向像素點(diǎn)個數(shù)，低分辨率圖像的頻譜數(shù)據(jù)表示為
11.如權(quán)利要求10所述的高分辨率圖像復(fù)原系統(tǒng)，其特征在于，所述補(bǔ)零法頻譜數(shù)據(jù)模塊將所述高分辨率圖像的補(bǔ)零法頻譜數(shù)據(jù)表示為G(kx，ky)P(kx, ky)，其中，
12.如權(quán)利要求11所述的高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取系統(tǒng)，其特征在于，所述補(bǔ)零法圖像模塊將所述高分辨率圖像的低頻頻譜數(shù)據(jù)補(bǔ)零法圖像表示為識U) = 1K(Uv)P(H)]。
13.如權(quán)利要求12所述的高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取系統(tǒng)，其特征在于，所述奇異譜函數(shù)模塊，用于 初始化:
14.如權(quán)利要求13所述的高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取系統(tǒng)，其特征在于，所述奇異譜函數(shù)模塊根據(jù)差分方程Φα，j)*h(i，j)=S (i，j)的零狀態(tài)的解獲取奇異函數(shù)h(i, j)。
15.如權(quán)利要求14所述的高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取系統(tǒng)，其特征在于，所述奇異譜函數(shù)模塊根據(jù)"(O,..)=刃沖,./)]茯取奇異譜函數(shù)。
16.如權(quán)利要求15所述的高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取系統(tǒng)，其特征在于，所述坐標(biāo)參數(shù)模塊，用于 初始化:
全文摘要
本發(fā)明提供了一種高分辨率圖像復(fù)原的坐標(biāo)參數(shù)獲取方法及系統(tǒng)。所述方法包括獲取待復(fù)原的高分辨率圖像的橫向或縱向像素點(diǎn)個數(shù)和一幅低分辨率圖像，根據(jù)所述像素點(diǎn)個數(shù)和低分辨率圖像獲取所述高分辨率圖像的低頻頻譜數(shù)據(jù)；根據(jù)所述高分辨率圖像的低頻頻譜數(shù)據(jù)獲取所述高分辨率圖像的補(bǔ)零法頻譜數(shù)據(jù)；對所述補(bǔ)零法頻譜數(shù)據(jù)作傅里葉變換以獲取高分辨率圖像的低頻頻譜數(shù)據(jù)補(bǔ)零法圖像；根據(jù)所述低頻頻譜數(shù)據(jù)補(bǔ)零法圖像獲取最佳奇異化算子，根據(jù)所述最佳奇異化算子獲取奇異函數(shù)，根據(jù)所述奇異函數(shù)獲取奇異譜函數(shù)；根據(jù)所述最佳奇異化算子并運(yùn)用點(diǎn)擴(kuò)散函數(shù)層析法獲取高分辨率圖像復(fù)原的坐標(biāo)參數(shù)，能夠在高頻頻譜數(shù)據(jù)缺失的情況下，快速高效地獲取高分辨率圖像復(fù)原的坐標(biāo)參數(shù)以供高分辨率圖像的復(fù)原。
文檔編號G06T5/50GK103116882SQ201310074160
公開日2013年5月22日 申請日期2013年3月7日 優(yōu)先權(quán)日2013年3月7日
發(fā)明者駱建華, 敬忠良 申請人:上海交通大學(xué)