% function J = myzoom2(I,k1,k2)
% myzoom1.m: 基于雙線性插值的圖像放大程序，其中，k1,k2分別 
%    是行、列的放大倍數。J為放大輸出圖像
%          
%  I = imread('D:\MATLAB7.0\images\FigProb9.20(left).jpg');
%  I = double(I);
clear all
I = [1 2 3;4 5 6;7 8 9];

%  if k1 < 1 || k2 < 1
%      error('input k1 or k2 is error');
%  end

k1 = 8.9; k2 = 3.6;
di = 1/k1; dj = 1/k2;
[r,c] = size(I);
si = round(k1*r);
sj = round(k2*c);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 將行、列放大分別放大K1,k2倍
for i = 1:r
    for j = 1:c
        if i == r && j ~= c
            J(round(k1*i),round(k2*(j-1)+1)) = I(i,j);           
        elseif j == c && i ~= r
            J(round(k1*(i-1)+1),round(k2*j)) = I(i,j);
        elseif i == r && j == c
            J(round(k1*i),round(k2*j)) = I(i,j);
        else
            J(round(k1*(i-1)+1),round(k2*(j-1)+1)) = I(i,j);
        end
    end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 對放大矩陣 J 進行雙線性插值

return