Image Processing Resim Döndürme

clc;clear;close all;
I = imread('cuval.jpg');
I = rgb2gray(I);
[H,W] = size(I);
R = ceil(sqrt((W^2)+(H^2)));
RI = uint8(zeros(R,R));
[RH,RW] = size(RI);
midx = round(W/2);
midy = round(H/2);
angle = (2*pi);

for x=1:H
for y=1:W
rotx = ceil(sin(angle) * (x - midx) + cos(angle) * (y - midy))+midx;
roty = ceil(cos(angle) * (x - midy) - sin(angle) * (y - midx))+midy;
if ((rotx < RW) && (roty < RH) && (rotx > 1) && (roty > 1))
RI(rotx,roty) = I(x,y);
end
end
end

imshow(RI);