圆形状、特殊的频谱图如何生成

% 读取输入图像
image = imread('input.png');

% 进行傅里叶变换
fft_image = fftshift(fft2(image));

% 取得幅度谱图
amplitude_spectrum = log(1 + abs(fft_image));
amplitude_spectrum = amplitude_spectrum / max(amplitude_spectrum(:));
amplitude_spectrum = im2uint8(amplitude_spectrum);

% 取得相位谱图
phase_spectrum = angle(fft_image);
phase_spectrum = im2uint8((phase_spectrum + pi) / (2 * pi));

% 显示幅度谱图和相位谱图
subplot(1,2,1);
imshow(amplitude_spectrum);
title('Amplitude Spectrum');
subplot(1,2,2);
imshow(phase_spectrum);
title('Phase Spectrum');