【正文】 果。當時，此時維納濾波就變成了逆濾波；當時，表明維納濾波避免了逆濾波中出現的對噪聲過多的放大作用 ,也就是說圖像維納復原比逆濾波復原效果好。 %從正面觀察彩色立方體 rgbcube(10,0,10)。%從對角線觀察彩色立方體 %********** 索 引 圖 像 的 顯 示 和 轉 換*********f=imread(39。)。 %f 是 RGB真彩圖像 %rgb圖像轉換成 8色索引圖像，不采用抖動方式 [X1,map1]=rgb2ind(f,8,39。)。 %采用抖動方式轉換到 8色索引圖像 [X2,map2]=rgb2ind(f,8,39。)。%顯示效果要好一些 g=rgb2gray(f)。 %將灰色圖像經過抖動處理，轉換打二值圖像 figure,imshow(g)。 %顯示抖動處理后的二值圖像 程序運行結果： * 彩色立方體原圖 不采用抖動方式轉換到 8 色索引圖像采用抖動方式轉換到 8 色索引圖像 灰度圖像抖動處理后的二值圖像 (2)彩色空間轉換 f=imread(39。)。 %轉換到 NTSC彩色空間 %f 是 RGB真彩圖像 ntsc_image=rgb2ntsc(f)。% 顯 示 亮 度 信 息figure,imshow(ntsc_image(:,:,2))。%顯示色差信息 %轉換到 HIS彩色空間 hsi_image=rgb2hsi(f)。% 顯示色度信息figure,imshow(hsi_image(:,:,2))。%顯示亮度信息 程序運行結果： 原圖 轉換到 NTSC彩色空間 顯示亮度信息顯示色差信息顯示色差信息 轉換到 HIS彩色空間 顯示色差信息顯示飽和度信顯示亮度信息 (3)彩色變換 f=imread(39。)。image39。 %打開 ice窗口對圖像進行調整 %在窗口中執(zhí)行以下操作： %a）得到圖像的補色 %b）拖動映射曲線，對圖像顯示效果進行修改 %c）在顏色通道中選中某一顏色，然后對映射曲線進行修改 程序運行結果 (1)： 全彩色圖片 ICE窗口 它的補色 ICE窗口 拖動映射曲線，圖像的顯示效果 ICE 窗口 f2=imread(39。)。% 在CMYK 彩色空間內打開圖像選擇 Red ICE 窗口 g2=ice(39。,f2,39。,39。)。 程序運行結果 (2) 原圖 在 CMYK彩色空間內打開圖像 調整映射參數后 ICE 窗口 (4)彩色空間濾波 %********************彩色圖像平滑 ********************** fc=imread(39。)。h=rgb2hsi(fc)。 S=h(:,:,2)。 %轉換到 HIS彩色空間 %色度分量 %飽和度分量 %亮度分量 %25 25 的方形平滑算子 w=fspecial(39。,25)。replicate39。%對亮度分量進行平滑操作 H1=cat(3,H,S,I_filtered)。 f=min(f,1)。 %彩 色圖像平 滑處理 %******************* 彩色圖像 銳化********************** fb=imread(39。)。181。%拉普拉斯算子 g=imfilter(fb,lapmask,39。)。 %圖像銳化后與原圖像相減 fen=imsubtract(fb,g)。%顯示差值圖像 程序運行結 果： 原圖彩色圖像平滑處理 拉普拉斯算子處理 銳化處理 (5)彩色圖像分割 %*****************彩色邊緣檢測 ************ closeall。clc。39。%f 如圖 64所示 %直接計算彩色邊緣 VG，并得到對單獨彩色平面的二維梯度求和形成的梯度 PPG [VG,A,PPG]=colrad(f)。 %顯示彩色梯度圖像 figure,imshow(A,[])。 %顯示分量梯度圖像 g=VGPPG。 程序運行結果： %顯示差值圖像 彩色梯度圖像彩色梯度對應的角度圖像 分量梯度圖像兩種梯度差值圖像 %*****************彩色圖像分割 ************************* f=imread(39。)。 %取得選定區(qū)域圖像的紅色分量 %手工選定圖中紅色區(qū)域 R=immultiply(mask,f(:,:,1))。%取得選定區(qū)域圖像的藍色分量 B=immultiply(mask,f(:,:,3))。 %取得選定區(qū)域的 RGB圖像 figure,imshow(g)。 %選定區(qū)域的大小 I=reshape(g,M*N,3)。 % 找到 mask 中 的 像 素 在 重 排 后 的 圖 像 g 中 位 置I=double(I(idx,1:3))。 %計算協(xié)方差矩陣 C 主對角線上的元素，即方差 %sd=sqrt(d)%求得這些元素的平方根，即均方差 均方差 %用歐式距離對圖像進行彩色分割，閾值為 25 E25=colorseg(39。,f,25,m)。E251(E25==0)=0figure,imshow(E251,[])。mahalanobis39。figure,imshow(M25,[])。 程序代碼如下 %思路：把原圖想轉換為 128位索引圖 closeall。clc。D:\Picture\39。figure,imshow(f)。 % 將 f 轉 換 成 灰 度 圖 像figure,imshow(f1)。 figure,imshow(x)。map(:,1)=。map(1,1)=0。 程序運行結果： 原 RGB圖像 灰度圖像 索引圖修改索引表后的圖像 (2) 節(jié)的（ 4）是將圖像轉換到 HIS 空間后再對亮度分量進行平滑操作，試直接對原圖像的 RGB 三個分量分別進行平滑操作，并計算兩種操作方法的圖像差值。D:\Picture\39。figure,imshow(f)。 S=f(:,:,2)。 %色度分量 %飽和度分量 %亮度分量 %25 25 的方形平滑算子 w=fspecial(39。,25)。replicate39。 % 對色度分量進行平滑操作S_filtered=imfilter(S,w,39。)。replicate39。 %對亮度分量進行平滑操作 %將 H、 S、 I三個分量重新組合起來 Ha=cat(3,H_filtered,S_filtered,I_filtered)。 %轉化成 double 類型 figure,imshow(Ha)。 %轉換到 HIS彩色空間 H=h(:,:,1)。 I=h(:,:,3)。average39。 I_filtered=imfilter(I,w,39。)。fb=hsi2rgb(H1)。 figure,imshow(fb)。figure,imshow(fen)。 (1)通過這個實驗我了解了 RGB 圖像以及索引圖像在 MATLAB 中的存儲和詢問方式 (2)了解 RGB 空間到其他彩色空間的轉換方式和轉換方法 (3)知道了一些彩色圖像的空間濾波方法 (4)了解在 RGB 向量空間中檢測彩色邊緣的方法 (5)了解在 RGB向量空間中對圖像進行分割的方法 遇到的問題：在做思考題 (2)的時候，計算兩種平滑處理的得到圖像的差值的中程序出現了問題，經過看出錯信息，認真檢查程序后，發(fā)現是因為直接平滑處理的圖像是 unit8 類型的，而轉換到 HIS 空間后，得到的圖像是 double 類型 ，二者的類型不同，導致使用 imsubtract函數的時候出現了問題。 第五篇：遙感圖像處理實驗報告 遙感圖像處理實驗報告 班級 11資環(huán)姓名學號實驗專題實驗室 F樓機房成績評定教師簽字 專 題 一 ： DEM 圖 像 進 行 彩 色 制圖 2（敘述 制圖過程并把自己處理結果加載到本文檔里） 專 題 二 ： TM 與 SPOT 數 據 融合 3（敘述該過程并處理結果加載到本文檔里。 4（敘述該過程并處理結果加載到本文檔） 專 題 四 ： 切 取 某 研 究 區(qū) 域 的 操作 6（主要是快速制圖。 8（其過程將其結果加載到本檔里）。 給生成的彩色圖像添加圖名、格網、比例尺、灰度條、等高線及數值等信息。結果如下圖 2 實驗專題：專題二： TM與 SPOT數據融合 主圖象窗口選擇 TransformImageSharpeningHSV，從一個打開的彩色圖像中選擇三個波段進行變換。 將 HSV 圖像重新轉換為 RGB 圖像。 加載最終圖像，并保存結果。 在主菜單 Registration 里的 SelectGCP(GroundControlPoints)來選擇地面控制點，并調整誤差。 圖像鑲嵌。 5 、在 PixelBasedMosaic 對 話 框 ， 選 擇ImportImportFilesandEditProperties，調整羽化、背景等參數。結果如圖： 實驗專題：專題四：切取某研究區(qū)域的操作 加載 一多 波 段彩 色圖 像， 在 主圖 像窗 口中 ， 選擇OverlayRegionofInterest。 通過選擇 BasicToolsMaskingBuildMask，建立掩膜。結果如下： 實驗專題：專題五：地圖制圖的方法 打開顯示要輸出的圖像。 設置圖名、投影并添加 mapkey等信息，保存圖像。 選擇相應的數字高程模型（ DEM）輸入文件及所需要的DEMResolution（像元數）復選框。值越高，放大越多。 選取最佳視覺位置保存。 選擇 ClassificationSupervised需要的分類方法，以平行六面體 (Parallelepiped)為例。在“ MaxstdevfromMean”文本框里，鍵入一個數值。 執(zhí)行，并保