要画一个光路图,本来可以用proe,但是鼠标不好用,有些操作也忘了,用MATLAB画了个。下面是用到的图片。
但是三维坐标中显示彩色图片的目标没有搞定,做了个灰度图,然后用仿射程序将彩色图片贴到了二维灰度图中。
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]);%调节坐标轴比例
效果如下
之后又用了个仿射程序把衍射图样彩色化了