Matlab編程實現(xiàn)圖像放大

這些代碼是根據(jù)岡薩雷斯的《數(shù)字圖像處理》（第三版） 里的第二章的2.6.5闡述的原理編寫的，空間變換用的是反向映射，插值方法用的是雙線性插值。
PS：這是我?第一次寫博客，排版不美觀還請見諒！
初學(xué)數(shù)字圖像處理，覺得還挺有趣的，所以我才想寫下這篇博客，如果能幫到大家真是太好了！！！

注意： 這里處理的圖像是灰度圖像，若處理彩色圖像，輸出圖像會變成灰度圖像！！！
由于是原理性編程，不同于一般的直接使用matlab自帶函數(shù)處理圖像，所以運行耗費的時間比較長，但對于幫助理解圖像放大的原理可以起到很好的作用?

流程圖如下：

Created with Rapha?l 2.2.0開始讀入圖像并歸一化定義放大倍數(shù)定義輸出圖像大小掃描圖像，反向映射得(x,y)坐標(biāo)在原圖像內(nèi)？坐標(biāo)x或y是小數(shù)？進行雙線性插值掃描完畢？顯示圖像結(jié)束直接賦予該坐標(biāo)像素原圖像的灰度值yesnoyesnoyesno

代碼如下：

% 首先讀入需要放大的灰度圖像 I = imread('F:\數(shù)字圖像處理\textT.tif'); I_in = im2double(I); % 轉(zhuǎn)化為double類型，并歸一化 m = size(I_in,1); % 獲取原圖像大小 n = size(I_in,2);% 放大tn倍，放大幾倍修改tn值就好 tn = 2; T = [tn,0,0;0,tn,0;0,0,1]; % 仿射矩陣T% 定義輸出圖像 I_out I_out = ones(m*tn,n*tn); for i = 1:m*tn % 從圖像的左上角開始遍歷輸出圖像for j = 1:n*tn% 輸出圖像像素坐標(biāo)反向映射得到輸入圖像像素坐標(biāo)x和yz = [i,j,1]/T;x = z(1);y = z(2);% 下面判斷反向映射得到的坐標(biāo)(x, y)是否超出輸入圖像的大小，% if條件成立則沒有超出，% 條件不成立則不進行任何操作，定義I_out時已經(jīng)默認(rèn)該點灰度為1if (x>=1 & y>=1) & (x<=m & y<=n) % 然后再判斷x或y是否為小數(shù)，if條件成立則是小數(shù)if (rem(x,1))||(rem(y,1))% 坐標(biāo)為小數(shù)則需要插值% 先獲取該小數(shù)坐標(biāo)(x, y)鄰近的四個像素坐標(biāo) x0 = floor(x);x1 = ceil(x);y0 = floor(y);y1 = ceil(y);% 獲取鄰近坐標(biāo)的灰度值a = I_in(x0,y0);b = I_in(x0,y1);c = I_in(x1,y0);d = I_in(x1,y1);% 雙線性內(nèi)插g1 = a + (x-x0) * (c - a);g2 = b + (x-x0) * (d - b);I_out(i,j) = g1 + (y - y0) * (g2 - g1); else % x和y都為整數(shù)則直接賦予原圖像的灰度值I_out(i,j) = I_in(x,y);endendend end% 顯示原圖像和放大后的輸出圖像 figure;imshow(I); figure;imshow(I_out);

運行結(jié)果如下：

• 輸入圖像如下：
  
• 輸出圖像如下：
  
  這里圖像被放大了2倍，可以看到輸出圖像的“T”形邊緣被模糊了，這是插值的結(jié)果。
  下次有時間的話我會把旋轉(zhuǎn)變換和偏移變換的代碼也寫出來，歡迎大家捧場?！！！

總結(jié)

以上是生活随笔為你收集整理的Matlab编程实现图像放大的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。