clear all;clc;close all;

im1 = rgb2gray(imread(\'F:\sup.jpg\'));% 掩模，第一张图，白色的部分在后面的surf中透明化了
im1 = imresize(im1,[64 64]);

im2 = rgb2gray(imread(\'F:\im3.jpg\'));% 样品，第二张图
im2(:,end) = [];
im2 = imresize(im2,[64,64]);
im2 = flipud(im2);


im3 = rgb2gray(imread(\'F:\diff.jpg\'));% 衍射图样，第三张图
im3(513:end,:) = [];
im3(:,513:end) = [];
im3 = flipud(im3);

% figure;
% subplot(121);imshow(im1);
% subplot(122);imshow(im4);
mask = PIEmask(256,200,128,128);
mask = imresize(mask,[64 64]);
%%
figure(\'color\',\'white\');

N = size(im3);
[x z] = meshgrid(-N/2:N/2-1);
y = 0*ones(N,N);
a = surf(x,y,z,  double(im3) );   % 画衍射图
set(a,\'linestyle\',\'none\');  %隐藏网格
colormap(gray);

hold on;
N = size(im2);
[x z] = meshgrid(-N/2:N/2-1);
y = -200*ones(N,N);
b = surf(x,y,z,  double(im2) );  % 画样品图
set(b,\'linestyle\',\'none\');  %隐藏网格
colormap(gray);

hold on;
N = size(im1);
[x z] = meshgrid(-N/2:N/2-1);
y = -210*ones(N,N);
% c = surf(x,y,z,  double(im1) );
c = surf(x,y,z,  double(im1),\'FaceAlpha\',\'flat\',\'AlphaDataMapping\',\'scaled\',\'AlphaData\',1-mask);  % 画掩模图，并将中间部分透明化
% alpha(0.5)
set(c,\'linestyle\',\'none\');  %隐藏网格
colormap(gray);
%%  后面的部分就是画红色的‘光’了
[yy zz] = meshgrid(-500:-210,-32:32);
xx = 32*ones(size(yy));

p = surf(xx,yy,zz,double(255*ones(size(zz))),\'FaceAlpha\',\'flat\',...
\'AlphaDataMapping\',\'scaled\',...
\'AlphaData\',0.2*ones(size(zz)),...
\'FaceColor\',\'red\');
set(p,\'linestyle\',\'none\');  %隐藏网格

p = surf(-xx,yy,zz,double(255*ones(size(zz))),\'FaceAlpha\',\'flat\',...
\'AlphaDataMapping\',\'scaled\',...
\'AlphaData\',0.2*ones(size(zz)),...
\'FaceColor\',\'red\');
set(p,\'linestyle\',\'none\');  %隐藏网格
%==========================================
[xx yy] = meshgrid(-32:32,-500:-210);
zz = 32.*ones(size(xx));

p = surf(xx,yy,zz,double(255*ones(size(zz))),\'FaceAlpha\',\'flat\',...
\'AlphaDataMapping\',\'scaled\',...
\'AlphaData\',0.2*ones(size(zz)),...
\'FaceColor\',\'red\');
set(p,\'linestyle\',\'none\');  %隐藏网格

p = surf(xx,yy,-zz,double(255*ones(size(zz))),\'FaceAlpha\',\'flat\',...
\'AlphaDataMapping\',\'scaled\',...
\'AlphaData\',0.2*ones(size(zz)),...
\'FaceColor\',\'red\');
set(p,\'linestyle\',\'none\');  %隐藏网格
%%
t = 0:0.01:2*pi;
x = 25*cos(t\');x = repmat(x,1,300);
z = 25*sin(t\');z = repmat(z,1,300);

y = linspace(-210,-200,300);
y = repmat(y,629,1);

p = surf(x,y,z,x*0,\'FaceAlpha\',\'flat\',...
\'AlphaDataMapping\',\'scaled\',...
\'AlphaData\',0.5*ones(size(z)),...
\'FaceColor\',\'red\');
set(p,\'linestyle\',\'none\');  %隐藏网格
%%
s1 = -25*sqrt(2)/2;
s2 = -256;
t1 = linspace(s1,s2,100);
x = [];
t2 = linspace(-s1,-s2,100);
for k = 1:100;
    tmp = (linspace(t1(k),t2(k),200))\';
    x = [x tmp];
end

y = linspace(-200,0,100);
y = repmat(y,200,1);
z = linspace(25*sqrt(2)/2,255,100);
z = repmat(z,200,1);

p = surf(x,y,z,x*0,\'FaceAlpha\',\'flat\',...
\'AlphaDataMapping\',\'scaled\',...
\'AlphaData\',0.2*ones(size(z)),...
\'FaceColor\',\'red\');
set(p,\'linestyle\',\'none\');  %隐藏网格

p = surf(x,y,-z,x*0,\'FaceAlpha\',\'flat\',...
\'AlphaDataMapping\',\'scaled\',...
\'AlphaData\',0.2*ones(size(z)),...
\'FaceColor\',\'red\');
set(p,\'linestyle\',\'none\');  %隐藏网格

p = surf(z,y,x,x*0,\'FaceAlpha\',\'flat\',...
\'AlphaDataMapping\',\'scaled\',...
\'AlphaData\',0.15*ones(size(z)),...
\'FaceColor\',\'red\');
set(p,\'linestyle\',\'none\');  %隐藏网格

p = surf(-z,y,x,x*0,\'FaceAlpha\',\'flat\',...
\'AlphaDataMapping\',\'scaled\',...
\'AlphaData\',0.15*ones(size(z)),...
\'FaceColor\',\'red\');
set(p,\'linestyle\',\'none\');  %隐藏网格
%%

grid off;
axis off;

view([150 -256 50]);
% view([150 -100 0]);
daspect([10,3,10]);%调节坐标轴比例