基于傅里叶变换的图像三维旋转实验Matlab代码

以下是基于傅里叶变换的图像三维旋转实验Matlab代码：

% 读取图片
I = imread('lena.jpg');
% 转换为灰度图
I = rgb2gray(I);
% 图像大小
[N, M] = size(I);
% 三维旋转角度
theta_x = 0.1;
theta_y = 0.2;
theta_z = 0.3;

% 傅里叶变换
F = fft2(I);

% 生成傅里叶变换矩阵
u = 0:(M-1);
v = 0:(N-1);
idx = find(u>M/2);
u(idx) = u(idx)-M;
idy = find(v>N/2);
v(idy) = v(idy)-N;
[V, U] = meshgrid(v, u);
D = sqrt(U.^2 + V.^2);
F_angle = atan2(V, U);

% 三维旋转变换
Fx = F_angle + theta_x*D;
Fy = F_angle + theta_y*D;
Fz = F_angle + theta_z*D;

% 傅里叶反变换
F_rotate = interp2(V, U, F, Fx, Fy, 'cubic');
F_rotate = interp2(V, U, F_rotate, Fz, F_angle, 'cubic');
I_rotate = real(ifft2(F_rotate));

% 显示结果
figure;
subplot(1, 2, 1);imshow(I);title('原图');
subplot(1, 2, 2);imshow(I_rotate);title('旋转后图像');

其中，`lena.jpg` 是需要旋转的图像，可以自行更换。`theta_x`、`theta_y`、`theta_z` 分别表示绕 x 轴、y 轴、z 轴旋转的角度，可以自行调整。最后通过傅里叶反变换得到旋转后的图像，并显示出来。