| 查看: 1727 | 回复: 3 | ||
[求助]
大神们，求帮我看下我这个数字菲涅尔全息图的代码有什么事问题 已有1人参与
|
|
我这个主要思路是，用点云算法，myfresnel函数做菲涅尔变化。参考光为平面波，源图片的各个点模拟成点光源。算出来的结果不太对。原图片和结果图片在附件里 大家帮我看看，如果觉得太费眼，给我指点一些用matlab编写全息的书或者论文资料，也感激不尽。信息光学的书我看过了，不过那上面大多是大段原理公式，想找些有实例的书和资料。谢谢 % ------------------------------------------------------------------------- % 功能：读入物光光波振幅分布 [fname,pname] = uigetfile('*.*','读入物光光波振幅分布图象'); if pname name = strcat(pname,fname); Picture =im2bw( imread(name)); figure,imshow(Picture);title('物光光波振幅分布图像'); picture = Picture(:,  ;O = picture; end clear fname pname name picture Picture; % 物光光波传播至全息面上的复振幅分布 (菲涅耳衍射） wl = 632.8e-9; % 单色光波长(m) z = 0.5; % 全息记录距离(m) k = 2*pi/wl; %%波数 detaO = 2.65e-4; % 物光波场像素间隔(m) detaH = 1.325e-4; % 衍射场光波像素间隔 y=repmat((512:-1:1)',1,512)-256.5; R=10*exp(j*k*(z*cos(pi*10/180)+y*detaH*sin(pi*10/180)));%%平面参考光复振幅 tic; uo=myfresnel(O,detaO,detaH,z,wl);% 物光光波传播至全息面上的复振幅分布(菲涅耳衍射） toc u=uo+R; %%在全息图上，物光波和参考光的复振幅相加 I=abs(u).^2; figure,imshow(I,[]); % ------------------------------------------------------------------------- function [ U ] = myfresnel(Object,DETA_O,DETA_H,Z,WL ) %%菲涅尔全息 %坐标轴和矩阵索引的转换：U(256.5-y,x+256.5)等于求xy坐标系中（n-256.5,256.5-m） k = 2*pi/WL; U=zeros(512,512); Uo=0; for x=-255.5:255.5 for y=-255.5:255.5 for x0=-127.5:127.5 for y0=-127.5:127.5 r=((DETA_H*x-DETA_O*x0)^2+(DETA_H*y-DETA_O*y0)^2+Z^2)^(1/2); Uo=Uo+(Object(128.5-y0,128.5+x0)/r)*exp(j*k*r); end end U(256.5-y,x+256.5)=Uo; end end  ABCDEF.jpg |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 结果.fig
2014-05-22 10:15:50, 1.61 M
» 猜你喜欢
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有4人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
什么时候开奖？
已经有10人回复
-
2026国自然函评费到账
已经有16人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
让我中一个面上吧！
已经有14人回复
-
今天放榜吗？
已经有16人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 最近卡在一步酰胺化合成很久，求大神们帮帮忙
已经有27人回复
- 求大神们，帮我看看出啥问题了
已经有1人回复
- 求亲们帮我看看这个到底怎么反应，愁死我了。
已经有13人回复
- 急求答案，可是好穷没金币分给大家，大神们帮我看下这道题该怎么解？小女在此谢过啦！
已经有3人回复
- 求大神帮检查下一个简单的UDF
已经有13人回复
- 想组装一台电脑，求懂电脑的高手-大神帮帮忙
已经有11人回复
- 关于2H-MoS2模型建立，文献中不是很清楚，求大神进来帮帮我
已经有0人回复
- 酚醛树脂发泡的问题，求大神指导！
已经有1人回复
- 15考研，精细或药学未定，浙江上海未定，本科化学，求帮忙分析我该考哪所学校！！
已经有8人回复
- 我想问下，，我是一名的动物药学的学生，就要考研了。求有经验的大神指导。。。
已经有6人回复
- 过来看过来观，片片大餐接着上~请求各位植物大神大仙们帮我看看这些是什么植物呢~~
已经有20人回复
- 求各位大神帮我看下一个小小的程序，只需要你一点儿时间。就会帮我很大的忙
已经有15人回复
- 求大神帮我看看，我的计算思路有没有问题？
已经有5人回复
- 请大神们看下我这个反应，有可能得到不？
已经有8人回复
smy1982
木虫 (正式写手)
- 物理EPI: 1
- 应助: 21 (小学生)
- 金币: 3128.7
- 散金: 6748
- 红花: 12
- 帖子: 853
- 在线: 450.4小时
- 虫号: 305139
- 注册: 2006-12-09
- 性别: GG
- 专业: 光学信息获取与处理
2楼2014-05-23 10:16:37
smy1982
木虫 (正式写手)
- 物理EPI: 1
- 应助: 21 (小学生)
- 金币: 3128.7
- 散金: 6748
- 红花: 12
- 帖子: 853
- 在线: 450.4小时
- 虫号: 305139
- 注册: 2006-12-09
- 性别: GG
- 专业: 光学信息获取与处理
【答案】应助回帖
感谢参与，应助指数 +1
|
clear clc Obj=imread('Obj.bmp'); Obj=double(Obj); [M,N]=size(Obj); % M,N are the pixel number of the object. Lamda=0.6328; % um, Wavelength of the input light DDist=150.0e+3; % um, Diffraction distance to the output plane % M=512; N=512; % Number of the sample array dx=10.0; dy=10.0; % um, Sampling space figure(1); colormap(gray); imshow(Obj);%物体 % Set coordinate of the sample array [x,y]=meshgrid(-M*dx/2:dx  M-1)*dx/2,-N*dy/2:dyN-1)*dy/2); SphFunct=exp(i*pi*(x.^2+y.^2)/(Lamda*DDist)); % The spherical wave function % Simulate the Fresnel diffraction Obj_F=fftshift(fft2(fftshift(Obj))); %Fourier transform of the slit function. SphFunct_F=fftshift(fft2(fftshift(SphFunct))); %Fourier transform of the spherical function. FrDiffract=fftshift(ifft2(fftshift(Obj_F.*SphFunct_F))); % Inverse Fourier transform figure(2); FrDiffract_I=abs(FrDiffract).^2; FrDiffract_I=FrDiffract_I/max(max(FrDiffract_I)); imshow(FrDiffract_I,[0 0.5]);%物光波 % Display the diffraction pattern a=40*pi/180;b=90*pi/180; R=1000*exp(i*2*pi/Lamda*DDist*sqrt(1-cos(a).^2-cos(b).^2)).*exp(i*2*pi/Lamda*(x*cos(a)+y*cos(b)));%%平面参考光复振幅 Hol=abs(R+FrDiffract).^2; % FrDiffract_I=FrDiffract.*conj(FrDiffract); % Intensity distribution Hol=Hol./max(max(Hol)); figure(3); imshow(Hol,[0 0.5]); %全息图 %%  全息图.jpg  物体.jpg  物体150毫米的衍射图像.jpg |
3楼2014-05-23 11:14:13
smy1982
木虫 (正式写手)
- 物理EPI: 1
- 应助: 21 (小学生)
- 金币: 3128.7
- 散金: 6748
- 红花: 12
- 帖子: 853
- 在线: 450.4小时
- 虫号: 305139
- 注册: 2006-12-09
- 性别: GG
- 专业: 光学信息获取与处理
|
clear clc Obj=imread('Obj.jpg'); Obj=double(Obj(:,:,1)); [M,N]=size(Obj); % M,N are the pixel number of the object. Lamda=0.6328; % um, Wavelength of the input light DDist=100.0e+3; % um, Diffraction distance to the output plane % M=512; N=512; % Number of the sample array dx=10.0; dy=10.0; % um, Sampling space figure(1); colormap(gray); imshow(Obj);%物体 % Set coordinate of the sample array [x,y]=meshgrid(-M*dx/2:dx M-1)*dx/2,-N*dy/2:dyN-1)*dy/2); SphFunct=exp(i*pi*(x.^2+y.^2)/(Lamda*DDist)); % The spherical wave function % Simulate the Fresnel diffraction Obj_F=fftshift(fft2(fftshift(Obj))); %Fourier transform of the slit function. SphFunct_F=fftshift(fft2(fftshift(SphFunct))); %Fourier transform of the spherical function. FrDiffract=fftshift(ifft2(fftshift(Obj_F.*SphFunct_F))); % Inverse Fourier transform figure(2); FrDiffract_I=abs(FrDiffract).^2; FrDiffract_I=FrDiffract_I/max(max(FrDiffract_I)); imshow(FrDiffract_I,[0 0.5]);%物光波 % Display the diffraction pattern a=40*pi/180;b=90*pi/180; R=100000*exp(i*2*pi/Lamda*DDist*sqrt(1-cos(a).^2-cos(b).^2)).*exp(i*2*pi/Lamda*(x*cos(a)+y*cos(b)));%%平面参考光复振幅 Hol=abs(R+FrDiffract).^2; % FrDiffract_I=FrDiffract.*conj(FrDiffract); % Intensity distribution Hol=Hol./max(max(Hol)); figure(3); imshow(Hol,[0 0.5]); %全息图 %%% 你的物体最好保存为灰度的图像》  obj.jpg  全息图1.jpg |
4楼2014-05-23 11:25:39
40