菱形种子法实现相位解包裹

laughvn

1. clc; clear; close all;
   
2. % 创建一个模拟的包裹相位图像（这里只是一个简单的示例）  
   
3. G = 6;
   
4. N = 128;
   
5. phi0 = peaks(N)*G; %模拟初始相位
   
6. wrapped_phase = angle(exp(1i*phi0));         %包裹相位
   
7.   
   
8. % 选择一个种子点（例如图像中心）  
   
9. seed_x = 20;  
   
10. seed_y = 20;  
    
11. 
    
12. tic
    
13. % 调用菱形种子法函数进行相位解包裹  
    
14. unwrapped_phase = diamond_seed_unwrap(wrapped_phase, seed_x, seed_y);  
    
15. toc
    
16. 
    
17. % 显示结果（包裹相位和解包裹相位）  
    
18. figure;  
    
19. subplot(1, 2, 1); imagesc(phi0); title('Wrapped Phase');  
    
20. subplot(1, 2, 2); imagesc(unwrapped_phase); title('Unwrapped Phase');

复制代码

laughvn

512*512的包裹相位，运行时间 0.864690 秒