图像的傅里叶变换的频谱特征 三（平移，旋转，相位的重要性）

最新推荐文章于 2025-03-18 12:21:34 发布
最新推荐文章于 2025-03-18 12:21:34 发布
阅读量1w 收藏 43
点赞数 9
分类专栏： 数字图像处理 文章标签： 图像处理
原文链接：https://blog.csdn.net/daduzimama/article/details/80394596
版权
本文探讨了图像傅里叶变换的频谱特征，重点分析了图像平移和旋转对频谱及相位的影响，通过Matlab代码演示了平移不会改变频谱，但会改变相位，而旋转则会使频谱随之旋转。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

图像傅里叶变换的频谱特征 三

 

6，平移和旋转

图像的平移并不会影响图像的频谱，同时，图像的相位会随着图像的旋转而旋转。

 

Part I 平移和旋转对频谱的影响

下面我用矩形的频谱图来说明图像中矩形的平移并不会对频谱有丝毫的影响。

 

Matlab代码：

 


 
 
  1. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     clear all
    
    
    
   
   
  2. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     close all
    
    
    
   
   
  3. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %% Author: J27
    
    
    
   
   
  4. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     % Jesus love you!
    
    
    
   
   
  5. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  6. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Isize = 
     
     512;
    
    
    
   
   
  7. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Rwidth = 
     
     50;
    
    
    
   
   
  8. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Rlength = 
     
     3*Rwidth;
    
    
    
   
   
  9. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  10. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect = zeros(Isize);
    
    
    
   
   
  11. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect(floor((Isize - Rlength)/
     
     2) + 
     
     1:floor((Isize - Rlength)/
     
     2) + Rlength,...
    
    
    
   
   
  12. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
          (floor(Isize - Rwidth)/
     
     2) + 
     
     1:floor((Isize - Rwidth)/
     
     2) + Rwidth) = 
     
     1;
    
    
    
   
   
  13. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  14. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     1,
     
     1)
    
    
    
   
   
  15. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshowpair(Irect,log(abs(fftshift(fft2(Irect)))+
     
     1),
     
     'montage')
    
    
    
   
   
  16. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  17. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect = zeros(Isize);
    
    
    
   
   
  18. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect(floor((Isize - Rlength)/
     
     2) + 
     
     150 + 
     
     1:floor((Isize - Rlength)/
     
     2) + 
     
     150 + Rlength,...
    
    
    
   
   
  19. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
          (floor(Isize - Rwidth)/
     
     2) + 
     
     1 + 
     
     200:floor((Isize - Rwidth)/
     
     2) + Rwidth + 
     
     200) = 
     
     1;
    
    
    
   
   
  20. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  21. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     1,
     
     2)
    
    
    
   
   
  22. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshowpair(Irect,log(abs(fftshift(fft2(Irect)))+
     
     1),
     
     'montage')
    
    
    
   
   
  23. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  24. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect = zeros(Isize);
    
    
    
   
   
  25. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect(floor((Isize - Rlength)/
     
     2) + 
     
     1:floor((Isize - Rlength)/
     
     2) + Rlength,...
    
    
    
   
   
  26. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
          (floor(Isize - Rwidth)/
     
     2) + 
     
     1 - 
     
     80:floor((Isize - Rwidth)/
     
     2) + Rwidth - 
     
     80) = 
     
     1;
    
    
    
   
   
  27. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  28. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     1,
     
     3)
    
    
    
   
   
  29. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshowpair(Irect,log(abs(fftshift(fft2(Irect)))+
     
     1),
     
     'montage')

 

再比如下面这个例子：

 

再来看看频谱随着矩形的旋转而旋转相同的角度。

Matlab代码：


 
 
  1. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Isize = 
     
     512;
    
    
    
   
   
  2. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Rwidth = 
     
     50;
    
    
    
   
   
  3. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Rlength = 
     
     3*Rwidth;
    
    
    
   
   
  4. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect = zeros(Isize);
    
    
    
   
   
  5. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  6. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect(floor((Isize - Rlength)/
     
     2) + 
     
     1:floor((Isize - Rlength)/
     
     2) + Rlength,...
    
    
    
   
   
  7. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
          (floor(Isize - Rwidth)/
     
     2) + 
     
     1:floor((Isize - Rwidth)/
     
     2) + Rwidth) = 
     
     1;
    
    
    
   
   
  8. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  9. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irot = imrotate(Irect, 
     
     15, 
     
     'crop', 
     
     'bilinear');
    
    
    
   
   
  10. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     1,
     
     1)
    
    
    
   
   
  11. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshowpair(Irot,log(abs(fftshift(fft2(Irot)))+
     
     1),
     
     'montage')
    
    
    
   
   
  12. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  13. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect = zeros(Isize);
    
    
    
   
   
  14. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect(floor((Isize - Rlength)/
     
     2) + 
     
     1:floor((Isize - Rlength)/
     
     2) + Rlength,...
    
    
    
   
   
  15. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
          (floor(Isize - Rwidth)/
     
     2) + 
     
     1:floor((Isize - Rwidth)/
     
     2) + Rwidth) = 
     
     1;
    
    
    
   
   
  16. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  17. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irot = imrotate(Irect, 
     
     45, 
     
     'crop', 
     
     'bilinear');
    
    
    
   
   
  18. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     1,
     
     2)
    
    
    
   
   
  19. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshowpair(Irot,log(abs(fftshift(fft2(Irot)))+
     
     1),
     
     'montage')
    
    
    
   
   
  20. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  21. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect = zeros(Isize);
    
    
    
   
   
  22. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect(floor((Isize - Rlength)/
     
     2) + 
     
     1:floor((Isize - Rlength)/
     
     2) + Rlength,...
    
    
    
   
   
  23. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
          (floor(Isize - Rwidth)/
     
     2) + 
     
     1:floor((Isize - Rwidth)/
     
     2) + Rwidth) = 
     
     1;
    
    
    
   
   
  24. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  25. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irot = imrotate(Irect, 
     
     90, 
     
     'crop', 
     
     'bilinear');
    
    
    
   
   
  26. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     1,
     
     3)
    
    
    
   
   
  27. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshowpair(Irot,log(abs(fftshift(fft2(Irot)))+
     
     1),
     
     'montage')

 

 

Part II 平移和旋转对相位的影响

     先用一个简单的例子来说明图像相位的作用（所用图像为cameraman），在图像的频域分析和滤波中，相位是常常被忽略的。虽然相位分量的贡献很不直观，但是它恰恰很重要。相位是频谱中各正弦分量关于原点的位移的度量。

上面的小实验充分说明了，看似无用的，且常常被忽略的相位，在DFT的频域中起到了多么重要的作用（注意区分实部和虚部（直角坐标系）VS 频谱和相位（极坐标系）！）。

Matlab代码：


 
 
  1. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     close all;
    
    
    
   
   
  2. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     clear all;
    
    
    
   
   
  3. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     I = im2double(imread(
     
     'cameraman.tif'));
    
    
    
   
   
  4. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     DFT = fft2(I);
    
    
    
   
   
  5. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     absF = abs(DFT);
    
    
    
   
   
  6. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Phi = atan2(imag(DFT),real(DFT));
    
    
    
   
   
  7. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     2,
     
     3,
     
     1);
    
    
    
   
   
  8. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(I,[]);
    
    
    
   
   
  9. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     title(
     
     'Image(Cameraman)');
    
    
    
   
   
  10. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     2,
     
     3,
     
     2);
    
    
    
   
   
  11. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(log(fftshift(absF) + 
     
     1),[]);
    
    
    
   
   
  12. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     title(
     
     'Spectrum');
    
    
    
   
   
  13. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     2,
     
     3,
     
     3);
    
    
    
   
   
  14. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(fftshift(Phi),[]);
    
    
    
   
   
  15. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     title(
     
     'Phase');
    
    
    
   
   
  16. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  17. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     F = absF.*exp(
     
     1j*Phi);
    
    
    
   
   
  18. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Rebuild = im2uint8(real(ifft2(F)));
    
    
    
   
   
  19. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     2,
     
     3,
     
     4);
    
    
    
   
   
  20. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(Rebuild,[]);
    
    
    
   
   
  21. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     title(
     
     'Rebuild with orginal specturm and phase');
    
    
    
   
   
  22. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     % rebuild with spectrum only
    
    
    
   
   
  23. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     F = absF.*exp(
     
     1j*
     
     1);
    
    
    
   
   
  24. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Rebuild = im2uint8(real(fftshift(ifft2(F))));
    
    
    
   
   
  25. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     2,
     
     3,
     
     5);
    
    
    
   
   
  26. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(Rebuild,[]);
    
    
    
   
   
  27. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     title(
     
     'Rebuild with spectrum only');
    
    
    
   
   
  28. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     % rebuild with Phase only
    
    
    
   
   
  29. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     g = 
     
     1.*exp(
     
     1j*Phi);
    
    
    
   
   
  30. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     g = im2uint8(real(ifft2(g)));
    
    
    
   
   
  31. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     2,
     
     3,
     
     6);
    
    
    
   
   
  32. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(g,[]);
    
    
    
   
   
  33. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     title(
     
     'Rebuild with Phase only');

     接下来我们再来看看图像在空间域中的移位和旋转对相位有什么影响。下图中，左边一列是图像，中间一列是频谱，右边一列是相位图。你必须意识到，通过肉眼，你很难从相位图中得到什么有用的信息。

(上图中最后一行打错了，不是“旋转改变了相位”而是“平移改变了相位”)

Matlab代码：


 
 
  1. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     clear all
    
    
    
   
   
  2. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     close all
    
    
    
   
   
  3. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %% Author: J27
    
    
    
   
   
  4. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     % Jesus love you!
    
    
    
   
   
  5. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  6. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     % shifting and rotation
    
    
    
   
   
  7. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Isize = 
     
     512;
    
    
    
   
   
  8. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Rwidth = 
     
     50;
    
    
    
   
   
  9. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Rlength = 
     
     3*Rwidth;
    
    
    
   
   
  10. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  11. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect = zeros(Isize);
    
    
    
   
   
  12. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect(floor((Isize - Rlength)/
     
     2) + 
     
     1:floor((Isize - Rlength)/
     
     2) + Rlength,...
    
    
    
   
   
  13. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
          (floor(Isize - Rwidth)/
     
     2) + 
     
     1:floor((Isize - Rwidth)/
     
     2) + Rwidth) = 
     
     1;
    
    
    
   
   
  14. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Idft = fft2(Irect);
    
    
    
   
   
  15. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     3,
     
     1);
    
    
    
   
   
  16. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(Irect);
    
    
    
   
   
  17. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     3,
     
     2);
    
    
    
   
   
  18. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(log(abs(fftshift(Idft))+
     
     1),[]);
    
    
    
   
   
  19. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     3,
     
     3);
    
    
    
   
   
  20. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow((atan2(imag(Idft),real(Idft))),[]);
    
    
    
   
   
  21. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  22. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irot = imrotate(Irect, 
     
     45, 
     
     'crop', 
     
     'bilinear');
    
    
    
   
   
  23. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Idft = fft2(Irot);
    
    
    
   
   
  24. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     3,
     
     4);
    
    
    
   
   
  25. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(Irot);
    
    
    
   
   
  26. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     3,
     
     5);
    
    
    
   
   
  27. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(log(abs(fftshift(Idft))+
     
     1),[]);
    
    
    
   
   
  28. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     3,
     
     6);
    
    
    
   
   
  29. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(atan2(imag(Idft),real(Idft)),[]);
    
    
    
   
   
  30. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  31. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect = zeros(Isize);
    
    
    
   
   
  32. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Irect(floor((Isize - Rlength)/
     
     2) + 
     
     150 + 
     
     1:floor((Isize - Rlength)/
     
     2) + 
     
     150 + Rlength,...
    
    
    
   
   
  33. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
          (floor(Isize - Rwidth)/
     
     2) + 
     
     1 + 
     
     200:floor((Isize - Rwidth)/
     
     2) + Rwidth + 
     
     200) = 
     
     1;
    
    
    
   
   
  34. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Idft = fft2(Irect);
    
    
    
   
   
  35. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     3,
     
     7);
    
    
    
   
   
  36. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(Irect);
    
    
    
   
   
  37. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     3,
     
     8);
    
    
    
   
   
  38. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(log(abs(fftshift(Idft))+
     
     1),[]);
    
    
    
   
   
  39. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     3,
     
     3,
     
     9);
    
    
    
   
   
  40. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(atan2(imag(Idft),real(Idft)),[]);

 

Matlab代码：


 
 
  1. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     I = im2double(imread(
     
     'cameraman.tif'));
    
    
    
   
   
  2. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Idft = fft2(I);
    
    
    
   
   
  3. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     figure;
    
    
    
   
   
  4. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     2,
     
     3,
     
     1);
    
    
    
   
   
  5. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(I,[]);
    
    
    
   
   
  6. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     2,
     
     3,
     
     2);
    
    
    
   
   
  7. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(log(abs(fftshift(Idft))+
     
     1),[]);
    
    
    
   
   
  8. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     2,
     
     3,
     
     3);
    
    
    
   
   
  9. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(atan2(imag(Idft),real(Idft)),[]);
    
    
    
   
   
  10. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  11. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Iswap = fftshift(I);
    
    
    
   
   
  12. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Idft = fft2(Iswap);
    
    
    
   
   
  13. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     2,
     
     3,
     
     4);
    
    
    
   
   
  14. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(Iswap,[]);
    
    
    
   
   
  15. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     2,
     
     3,
     
     5);
    
    
    
   
   
  16. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(log(abs(fftshift(Idft))+
     
     1),[]);
    
    
    
   
   
  17. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subplot(
     
     2,
     
     3,
     
     6);
    
    
    
   
   
  18. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(atan2(imag(Idft),real(Idft)),[]);

                                                                                （全文完）

                                                                                 谢谢收看！

数字图像处理 --- 图像的傅里叶变换的频谱特征 一（周期性，能量分布，fftshift，交错性）
松下J27录放机
05-02 6万+
图像傅里叶变换的频谱特征 傅里叶变换在一维信号处理中的地位是显著的，是不可撼动的，然后傅里叶变换在图像处理领域中的应用似乎稍逊一筹，黯然失色。究其原因，我想了很久，请允许我用非官方的，不正规的，但却通俗易懂的方式说一下。 一句话概括就是：要化繁为简（DSP），不要弄巧成拙（DIP）。 前半句说的是DF......
数字图像处理 --- 图像的傅里叶变换的频谱特征 三（平移，旋转，相位的重要性）
松下J27录放机
05-21 3万+
图像傅里叶变换的频谱特征 三 6，平移和旋转 图像的平移并不会影响图像的频谱，同时，图像的相位会随着图像的旋转而旋转。 Part I 平移和旋转对频谱的影响 下面我用矩形的频谱图来说明图像中矩形的平移并不会对频谱有丝毫的影响。 Matlab代码： clear all close all %%.........
对傅里叶变换FFT性质的理解 平移 旋转 缩放
Maxwave13的博客
08-31 1万+
网上很多关于FFT变换的应用，但是对其性质的探索并不多。 我做了一个简单的实验，来探讨平移，缩放，旋转的性质。 由于我浅薄的知识，如有不对的地方，请各位大佬指正。 平移：FFT具有平移不变性，在空域平移图像，频域的信号不发生变换。 旋转：旋转同一性，空域图像的旋转，也会带动频域图像的旋转。 缩放：在空域缩小图像，频域的信号会相应的缩小信息量，视觉上体现为放大频谱图。 代码实现 先写一些准备代码，这里要注意，需要将fft2变换的结果，进行后续处理，方便展示。 1.abs傅里叶变换的结果是复数，需要取模转为
【图像处理】基于图像特征检测——使用相位拉伸变换附Matlab代码
最新发布
Matlab_jiqi的博客
03-18 751
图像特征检测是计算机视觉领域一项至关重要的任务，它旨在提取图像中具有辨识性的显著性信息，用于后续图像分析、识别、匹配以及三维重建等应用。传统图像特征检测方法，如边缘检测、角点检测等，在处理某些复杂场景或噪声干扰下表现出局限性。近年来，基于相位拉伸变换 (Phase Stretch Transform, PST) 的特征检测技术由于其独特的优势，在图像特征提取方面展现出强大的潜力。本文将深入探讨基于图像特征检测的相位拉伸变换技术，分析其原理、优点以及应用，并对其未来发展方向进行展望。一、相位拉伸变换的原理。
图像的傅里叶变换频谱特性分析
01-31
关于傅立叶频谱分析，权威，对于深度了解傅立叶变换有重要意义
图像傅里叶变换
热门推荐
Rachel Zhang的专栏
06-01 19万+
冈萨雷斯版里面的解释非常形象：一个恰当的比喻是将傅里叶变换比作一个玻璃棱镜。棱镜是可以将光分解为不同颜色的物理仪器，每个成分的颜色由波长（或频率）来决定。傅里叶变换可以看作是数学上的棱镜，将函数基于频率分解为不同的成分。当我们考虑光时,讨论它的光谱或频率谱。同样, 傅立叶变换使我们能通过频率成分来分析一个函数。 Fourier theory讲的就是：任何信号（如图像信号）都可以表示成
图像处理离散傅里叶变换频谱相位谱幅度谱关系
01-12
例如，通过仅改变相位谱，可以保持幅度谱不变，实现图像的旋转或平移。另一方面，幅度谱的改变可能导致图像的亮度变化，但不会影响图像的基本结构。 四、MATLAB实现 MATLAB 提供了 `fft2` 函数来计算二维离散傅里叶...
数字图像处理实验二 图像傅里叶变换 实验报告 附代码 结果 数据
12-16
实验报告——数字图像处理：图像傅里叶变换 在图像处理领域，傅里叶变换是一种重要的工具，它能够将图像从空间域转换到频域，从而揭示图像的频率特性。本实验旨在深入理解图像频域处理的意义，掌握二维傅里叶变换的...
图像的二维傅里叶变换频谱图特点研究.rar
07-11
图像的二维傅里叶变换频谱图特点研究，以及基于halcon得到的一些图谱分析。
如何理解 图像傅里叶变换的频谱图
程序员，他们想的是什么?他们想的永远都是技术，他们崇尚的也永远都是技术。
03-23 7262
本文转载于松下J27，原文分为三章讲解，在此合并转载 文一：https://blog.csdn.net/daduzimama/article/details/80109139 文二：https://blog.csdn.net/daduzimama/article/details/80394570 文三：https://blog.csdn.net/daduzimama/article/details/80394596 如何理解 图像傅里叶变换的频谱图1、周期性2、高低频率的分布3、频谱图的能量分布4、纵横“交
图像的傅里叶变换
whaosoft143ai的博客
08-18 5129
比如南非世界杯时，南非人吹的呜呜主拉的声音太吵了，那么对现场的音频做傅立叶变化（当然是对声音的数据做），会得到一个展开式，然后找出呜呜主拉的特征频率，去掉展开式中的那个频率的sin函数，再还原数据，就得到了没有呜呜主拉的嗡嗡声的现场声音。他们是可逆的，想不到吧，乱七八糟的东西也有规律了。对于傅里叶变换的逆操作，使用OpenCV的函数就是返回结果=cv2.idft(原始数据)，然后计算幅度函数仍然是返回值=cv2.magnitude(参数1，参数2)，numpy.fft.ifftshift​​​​​​​。..
数字图像处理 --- 图像的傅里叶变换的频谱特征 二（方向性）
松下J27录放机
05-21 1万+
图像傅里叶变换的频谱特征 二 很多人都不了解图像（二维）频谱中的每一点究竟代表了什么，有什么意义? 简而言之：二维频谱中的每一个点都是一个与之一一对应的二维正弦/余弦波。 5，方向性（direction） 在二维频谱图中的任意“一对亮点”（注意：频谱的对称性），都在相应的空间域有一个与之相对应的二维正弦波。亮点在二维频谱中的位置决定了与之对应的正弦波的频率和方向。......
图像的二维傅里叶变换频谱图特点研究
weixin_30563917的博客
10-13 2117
一、先放一些相关的结论： 1、傅里叶变换的幅值称为傅里叶谱或频谱。 2、F(u)的零值位置与“盒状”函数的宽度W成反比。 3、卷积定理：空间域两个函数的卷积的傅里叶变换等于两个函数的傅里叶变换在频率域中的乘积。f(t)*h(t) <=> H(u)F(u) 4、采样定理：如果以超过函数最高频率的两倍的取样率来获得样本，连续的带限函数可以完全地从它的样本集来恢复。 5、严重的混淆...
图像处理之图像傅里叶变换
qq_44111805的博客
08-21 4万+
傅里叶变换是在以时间为自变量的“信号”与频率为自变量的“频谱”函数之间的某域研究中较复杂的问题在频域中变得简单起来，从而简化其分析过程；当自变量“时间”或“频率”为连续形式和离散形式的不同组合时，就可以形成各种不同的傅里叶变换对，即“信号”与“频谱”的对应关系。即傅里叶变换的物理意义是将图像的灰度分布函数变换为图像的频率分布函数，傅里叶逆变换是将图像的频率分布函数变换为灰度分布函数。函数fft()和fft2()的关系为fft2(X)=fft(fft(X).‘).’。，坐标原点是低频，向外是高频。
图像处理:傅里叶变换
forfuture3513的博客
03-19 4323
图像处理：傅里叶变换
[openCV]图像的傅里叶频谱
zhoufan的专栏
05-21 1万+
1.图像的傅里叶频谱的意义之前的博文其实已经归纳过这方面的内容了。我们常用的图像平滑处理，其实就是一个低通滤波，一定程度上去除高频信号，可以使得图像变得柔和（也就是平滑）。但是，在去除周期性噪声时候，空间域内的滤波（卷积）就不是那么好操作了。所以，这里时候，无论是理解起来方便，还是其他原因，都需要在频域内进行滤波。
二维傅里叶变换旋转matlab仿真
12-28
### 实现二维傅里叶变换旋转的Matlab仿真 为了在 MATLAB 中实现二维傅里叶变换并对其进行旋转仿真的操作，可以通过以下方法完成。首先，加载一张灰度图像，并对其应用 `fft2` 函数来进行二维傅立叶变换[^1]。 接着，通过调整频率域中的相位角来模拟旋转效果。这涉及到创建一个与原始图像大小相同的网格，并在此基础上构建相应的相位因子矩阵。最后再利用逆傅里叶变换 (`ifft2`) 将修改后的频谱转换回空间域以观察到旋转的效果。 以下是具体的实现过程： #### 加载图片并显示原图 ```matlab % 读取测试用的灰度图像 gray_image = imread('cameraman.tif'); figure; imshow(uint8(gray_image)); title('Original Image'); ``` #### 执行二维傅里叶变换 ```matlab % 对输入图像做双精度浮点数类型的强制转换后再求解其FFT fft_image = fftshift(fft2(double(gray_image))); figure; imagesc(log(abs(fft_image)+eps)); colormap gray; colorbar; title('Fourier Transform Magnitude Spectrum (log scale)'); axis equal tight; ``` #### 构建用于旋转的角度参数以及对应的相位因子 假设要顺时针方向旋转 θ 度，则有如下关系： \[ e^{-j \cdot 2\pi / N} \] 其中 \(N\) 是图像尺寸的一边长度（假定为正方形），\( j=\sqrt{-1}\) 表示虚单位。对于任意给定点 `(u,v)` 的位置，在新坐标系下的对应点可通过下述公式得到新的索引值 `(u',v')`: \[ u' = cos(\theta)\times u + sin(\theta)\times v \\ v' = -sin(\theta)\times u + cos(\theta)\times v \] 因此可以根据上述定义构造出整个图像范围内的相位偏移量数组。 ```matlab angle_degrees = 30; % 设置想要旋转的角度 [M, N] = size(gray_image); [U,V] = meshgrid(-1)/2), ... -floor(M/2):ceil((M-1)/2)); rotation_matrix = [cosd(angle_degrees) sind(angle_degrees);... -sind(angle_degrees) cosd(angle_degrees)]; UV_prime = rotation_matrix * [U(:)', V(:)']'; U_rotated = reshape(UV_prime(1,:), M, N); V_rotated = reshape(UV_prime(2,:), M, N); phase_factor = exp(-1i*2*pi*(U_rotated./N + V_rotated./M)); rotated_fft = double(fft_image .* phase_factor); ``` #### 进行反向傅里叶变换查看结果 ```matlab inverse_transformed = real(ifft2(ifftshift(rotated_fft))); figure; subplot(1,2,1); imshow(uint8(inverse_transformed)), title('Rotated Image via FFT'); subplot(1,2,2); fft_after_rotation = fftshift(fft2(double(inverse_transformed))); imagesc(log(abs(fft_after_rotation)+eps)); colormap gray; colorbar; title('Magnitude Spectrum After Rotation'); axis equal tight; ``` 这段程序展示了如何通过对频谱施加特定角度的相位变化来达到图像旋转的目的。值得注意的是这里仅考虑理想情况下的纯几何变换而不涉及任何物理意义层面的影响因素。
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值