【正文】 因此還有： 式中，是像素點(diǎn)的特性函數(shù)，為沖擊響應(yīng)。MATLAB提供了離散傅里葉變換（DFT）、快速傅里葉變換（FFT）、離散余弦變換（DCT），連續(xù)小波變換（CWT）、離散小波變換（DWT）及其反變換等變換。MATLAB包括擁有數(shù)百個(gè)內(nèi)部函數(shù)的主包和三十幾種工具包（Toolbox）。圖像復(fù)原過(guò)程往往是對(duì)提高圖像質(zhì)量起著重要的作用的數(shù)字圖像處理方法。MATLAB的開(kāi)放性廣受用戶歡迎，除去內(nèi)部函數(shù)，MATLAB的所有主包文件和各種工具包都是可讀并可以修改的文件，通過(guò)對(duì)源程序的修改或添加，用戶可以構(gòu)造自己的專用工具包。（5）圖像復(fù)原與重建。由此可看出，如果已知退化圖像的傅里葉變換和“濾波”傳遞函數(shù)，則可以求得原始圖像的傅里葉變換，經(jīng)反傅里葉變換就可以求得原始圖像f(x,y)，這就是逆濾波法的基本原理。在實(shí)際系統(tǒng)中，維納濾波經(jīng)常用下式近似：，K為特殊常數(shù)。 F=fftshift(fft2(I))。imshow(uint8(I2))。LR算法是一種迭代非線性復(fù)原算法，它是從最大似然公式印出來(lái)的，圖像用泊松分布加以模型化的。)。subplot(2,2,3)。因此，盲去卷積算法就是那些不以PSF知識(shí)為基礎(chǔ)的圖像復(fù)原的方法。,0,V)。)。3圖像重建 圖像重建的概述圖像處理中的一個(gè)重要研究分支是物體圖像的重建，它被廣泛的應(yīng)用于檢測(cè)和觀察中，而這種重建方法一般是根據(jù)物體的一些橫截面部分的投影而進(jìn)行的。為得到圖像函數(shù)，我們必須進(jìn)行反變換運(yùn)算，即：這些結(jié)論很容易推廣到三維情形中。按照二維傅里葉反變換標(biāo)準(zhǔn)定義，有 作代換：，寫成極坐標(biāo)（R,θ）的形式： ① 利用傅里葉變換共軛對(duì)稱性，有：②令： ③則③可以表示為：， ④在②中，當(dāng)用FFT計(jì)算投影數(shù)據(jù)的傅里葉變換F(R,θ)時(shí)，投影數(shù)據(jù)g(ρ,θ)總被有限截?cái)?。imagesc(theta3,xp,R3)。subplot(2,2,2)。代數(shù)重建技術(shù)就是事先對(duì)未知圖像的各像素給予一個(gè)初始估值，然后利用這些假設(shè)數(shù)據(jù)去計(jì)算各射線穿過(guò)對(duì)象時(shí)可能得到的投影值（射影和），再用它們和實(shí)測(cè)投影值進(jìn)行比較，根據(jù)差異獲得一個(gè)修正值，利用這些修正值，修正各對(duì)應(yīng)射線穿過(guò)的諸像素值。申老師嚴(yán)謹(jǐn)?shù)闹螌W(xué)風(fēng)范、高度的責(zé)任心給我留下了極其深刻的印象。參考文獻(xiàn)[1] 龔聲蓉，劉純平，[M].北京:清華大學(xué)出版社,123144[2] 張德豐，：清華大學(xué)出版社，2012.[3] （MATLAB版）[M].北京:電子工業(yè)出版社,:123133[4] 張錚,倪紅霞,苑春苗，[M].北京:人民郵電出版社,:194196[5] 陳傳波,[M].北京：機(jī)械工業(yè)出版社,161176[6] 林旭梅，[M].山東：中國(guó)石油大學(xué)出版社，:4975[7] 張德豐,[J].中山大學(xué)學(xué)報(bào),(6)：4447[8] [M].西安:西安電子科技大學(xué)出版社,[9] ：控制工程篇[M].北京:電子工業(yè)出版社,[10] 羅軍輝，[M].北京:機(jī)械工業(yè)出版社,[11] 趙怡紅,:化學(xué)工業(yè)出版社,:38 125[12]：控制工程篇[M].北京：電子工業(yè)出版社，[13]劉維一，[J].，2002,13（4）：389400[14]::CDMD:[15][D].西安：西安電子科技大學(xué)出版社，2010,2125致謝 在論文完成之際，我的心情十分激動(dòng)，無(wú)法用言語(yǔ)表達(dá)。title(39。I1=iradon(R1,10)。[R2,xp]=radon(P,theta2)。由此，我們得出一個(gè)結(jié)論，即：三維圖像不能在得到部分投影數(shù)據(jù)的過(guò)程中局部的重建，而必須延遲到所有投影數(shù)據(jù)都獲得之后才能重建。即： 現(xiàn)在假設(shè)將函數(shù)投影到一條經(jīng)過(guò)旋轉(zhuǎn)的直線上，該直線的旋轉(zhuǎn)角度為θ。imshow(P,[])。subplot(221)。,7,10)。另外此法存在一些較嚴(yán)重的缺陷，一是噪聲放大問(wèn)題，二是對(duì)于恢復(fù)圖像中的不同部分，分別執(zhí)行多少迭代才合適的問(wèn)題。imshow (Blurred)。Blurred=imfilter(I,PSF,39。title(維納濾波復(fù)原圖)。 F0=fftshift(fft2(I1))。title(39。 退化圖像、原始圖像、噪聲都是均勻隨機(jī)場(chǎng)，噪聲的均值為零，且與圖像不相關(guān)。 一般來(lái)說(shuō)，圖像的生成可以簡(jiǎn)單地被描述為如下數(shù)學(xué)模型：f(x,y)是成像景物，H是綜合退化因子，g(x,y)是退化圖像。（3）圖像變換。矩陣是MATLAB的基本數(shù)據(jù)單位，它的指令表達(dá)式與數(shù)學(xué)、工程中常用的形式十分相似，故用MATLAB來(lái)解決問(wèn)題事件比用CFORTRAN等語(yǔ)言簡(jiǎn)捷方便得多。圖像處理中的一個(gè)重要的研究分支是圖像重建，其意義在于要檢測(cè)到獲得物體的內(nèi)部結(jié)構(gòu)圖像，而不會(huì)其造成任何物體上的損傷。 MATLAB的主要功能 MATLAB是一種用高級(jí)技術(shù)計(jì)算語(yǔ)言和交互式環(huán)境，它集算法開(kāi)發(fā)、數(shù)值計(jì)算、數(shù)據(jù)分析以及數(shù)據(jù)可視化為一體?？捎媚鏋V波復(fù)原、維納濾波復(fù)原等方法實(shí)現(xiàn)圖像復(fù)原。但在實(shí)際中用逆濾波法存在病態(tài)的情況：當(dāng)H(u,v)=0時(shí)，或非常小的數(shù)值點(diǎn)上，F(xiàn)(u,v)將變成無(wú)窮大或非常大的數(shù)。 ②①的維納濾波要求未退化圖像和噪聲的功率必須是已知的。 k=。title(39。當(dāng)下面這個(gè)迭代收斂時(shí)模型的最大似然函數(shù)就可以得到一個(gè)令人滿意的方程： *代表卷積，代表未退化圖像的估計(jì)，g和h和以前定義一樣。 V=.003。imshow (BN)。 在過(guò)去的20年里，一種盲去卷積的方法