只用生成一个全白的图就OK
在这里插入代码片
clear all;
cm = 0.01;
um = 1e-6;
nm = 1e-9;
mm = 1e-3;
lambda = 638.8*nm;
Pitch = 20*um;
M_a = 1200;
N_a = 950;
xa =[-N_a/2:N_a/2-1]*Pitch;
ya = [-M_a/2:M_a/2-1]'*Pitch;
[theta,rou]=cart2pol(xa,ya);
l1=-5 ;
l=0;
theta_offset = pi;
vox1 = exp(l1*j*(theta+theta_offset));
vox1 = (vox1);
%vox2 = exp(l*j*(theta+theta_offset));
vox2 = exp(l*j*(theta+theta_offset));
vox2 = (vox2);
Image_1=mat2gray(angle(vox1));
Image_2=mat2gray(abs(vox2));
% figure()
% imshow(Image_1);
% imshow(Image_2);
imwrite(Image_1,'fig_1.bmp');
imwrite(Image_1,'fig_2.bmp');
%水平拼接两幅图
C=[Image_1,Image_2]
[M_a,N_a]=size(C)
% figure()
% imshow(C)
imwrite(C,'055_-1_5.bmp')
%%%------------
%注释:这里生成的全白图时在螺旋相位图片上修改的,
%需要满足条件时l=0,Image_2=mat2gray(abs(vox2));
%这里用的是幅值abs而不是相角。然后生成就行了。
%我这里没有imshow fig_2,没显示,但是已经保存到我的bin
%文件里了。