图像旋转 matlab程序

clear;clc
init=imread('原始路径程序');             %读入原始图像
[R, C,m] = size(init); % 获取图像大小
res = 255*ones( R,  C,m); % 构造结果矩阵。每个像素点默认初始化为255（白色）
alfa = -15 * 3.1415926 / 180.0; % 旋转角度
tras = [cos(alfa) -sin(alfa) 0; sin(alfa) cos(alfa) 0; 0 0 1]; % 旋转的变换矩阵
for i = 1 : R
    for j = 1 : C
        temp = [i; j; 1];
        temp = tras * temp;% 矩阵乘法
        x = uint16(temp(1, 1));
        y = uint16(temp(2, 1));
        % 变换后的位置判断是否越界
        if (x <= R) & (y <= C) & (x >= 1) & (y >= 1)
            res(i, j,:) = init(x, y,:);
        end
    end
end;
imshow(uint8(res));  % 显示图像