资源简介

Gerchberg-Saxton算法实现对图像相位的恢复,输入两张图片其中一张作为目标振幅。是对GS算法很好理解的一个程序。

资源截图

代码片段和文件信息

clear
clc
A=imread(‘C:\wo.jpg‘);  %从文件读取图象数据
B=imresize(A[512 512]);   %重新设定图像大小 
C=im2double(rgb2gray(B));  % 把灰度图像的数据类型转换成双精度浮点类型(有效数字可达15位,精确度更高)
a=imread(‘C:\gaosi.jpg‘);
b=imresize(a[512512]);
b=rgb2gray(b);
c=im2double(b); 
C0=abs(C);                 % 设定初始振幅
MB=abs(c);                 % 目标振幅
figuremesh(MB);title(‘已知目标振幅分布‘)
for p=1:600              % 迭代循环100次
    C=fftshift(C);   %将傅里叶变换的DC分量移到频谱中心区
    D=ifft2(C);     %二维傅里叶逆变换    
    D1=fftshift(D

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1132  2014-12-15 21:58  Gerchberg-  Saxton(GS)算法.m

----------- ---------  ---------- -----  ----

                 1132                    1


评论

共有 条评论

相关资源