| 查看: 1506 | 回复: 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
» 猜你喜欢
焦虑
已经有8人回复
-
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 最近卡在一步酰胺化合成很久，求大神们帮帮忙
已经有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
