matlab 角域重采样代码

以下是 MATLAB 中的角域重采样代码示例：

% 定义原始图像和重采样后的图像的角度范围
theta1 = linspace(0, 2*pi, 128);
theta2 = linspace(0, 2*pi, 256);

% 定义原始图像
f = zeros(128, 128);
f(64, 64) = 1;

% 计算原始图像的傅里叶变换
F = fftshift(fft2(f));

% 计算重采样后的傅里叶变换
F2 = zeros(256, 256);
for i=1:256
    for j=1:256
        x = (i-128)*cos(theta2(j)) + (j-128)*sin(theta2(j));
        y = -(i-128)*sin(theta2(j)) + (j-128)*cos(theta2(j));
        x = round(x) + 128;
        y = round(y) + 128;
        if (x>=1 && x<=128 && y>=1 && y<=128)
            F2(i, j) = F(x, y);
        end
    end
end

% 反傅里叶变换得到重采样后的图像
f2 = abs(ifft2(ifftshift(F2)));

% 显示原始图像和重采样后的图像
subplot(1,2,1);
imshow(f, []);
title('原始图像');
subplot(1,2,2);
imshow(f2, []);
title('重采样后的图像');

在此示例中，我们首先定义了原始图像和重采样后的图像的角度范围。然后，我们定义了原始图像并计算了其傅里叶变换。接下来，我们计算了重采样后的傅里叶变换，并使用反傅里叶变换得到了重采样后的图像。最后，我们显示了原始图像和重采样后的图像。