傅里叶变换全息图

最新推荐文章于 2024-07-03 22:21:00 发布

fpga和matlab

最新推荐文章于 2024-07-03 22:21:00 发布

阅读量4.3k

点赞数 10

分类专栏： MATLAB 板块2:图像-特征提取处理

本文链接：https://blog.csdn.net/ccsss22/article/details/113789787

版权

MATLAB 同时被 2 个专栏收录

1606 篇文章 1630 订阅

订阅专栏

板块2:图像-特征提取处理

317 篇文章 114 订阅

订阅专栏

该博客详细介绍了如何使用傅里叶变换处理全息图像。首先，从20120810191738_K2fjT.jpeg图片中获取256*256的灰度图像，并展示。接着，通过傅里叶变换得到物体的频谱，并展示了干涉记录。最后，经过复振幅处理和逆傅里叶变换，实现了全息图像的再现。

摘要由CSDN通过智能技术生成

function [varagout]=C(varagin)
%傅里叶变换全息图

clc;
close all;

I=zeros(256);
I=imresize(double(rgb2gray(imread('20120810191738_K2fjT.jpeg'))),[256,256]);;
subplot(2,2,1);imshow(I);
title('256*256灰度图像')
%定义接收屏幕
L=500;
M=256;
[x,y]=meshgrid(linspace(-L/2,L/2,M));
%先设参考波为平面波
r0=1;b=800000;a=800000;
R=r0.*exp(j.*2.*pi.*(b.*y+a.*x));
%对I进行傅里叶变换,并平移
G=fftshift(fft2(fftshift(I)));
subplot(2,2,2);imshow(abs(G),[]);title('物体的频谱')
%记录面上的光场 U=G+R;
%记录面上的光强
I1=G.*conj(R)+conj(G).*R;
subplot(2,2,3);
imshow(I1);title('CCD靶面的干涉记录')
%假如显影后，复振幅的透过率正比于光强
t1=2;
%负片
r1=-1;
t=t1+r1.*(G.*conj(R)+conj(G).*R);
%全息再现过程
U=ifft2(I1);
subplot(2,2,4);imshow(abs(U));title('再现后的图像')