clc,clear v=csvread('vdata.txt'); % 數(shù)據(jù) 這里的數(shù)據(jù)矩陣為v(5*27*6),可以看成一本27頁紙，每頁繪制了5*6的網(wǎng)格，% 然后27頁紙疊在一起。當(dāng)你理解本圖繪制后，數(shù)據(jù)可以隨意制作% 讀取文件，也可以讀取一部分 v=reshape(v,5,27,6); % 語法是 A = reshape（A，m，n）； 或者 A = reshape（A，[m,n]）;% 都是將A 的行列排列成m行n列。另外 reshape是按照列取數(shù)據(jù)的 [x,y,z]=meshgrid(1:27,1:5,1:6); % 三維網(wǎng)格化 h=figure(1); set(h,'name','取單切片') subplot(221) slice(x,y,z,v,[],[1],[]); % 切片函數(shù) y=1 shading interp % 平滑處理 set(gca,'zdir','reverse'); % ‘’都是屬性調(diào)整 axis equal % 將橫軸縱軸的定標(biāo)系數(shù)設(shè)成相同值 ,即單位長度相同 grid on subplot(222) slice(x,y,z,v,[],[2],[]); % 切片函數(shù) y=2 shading interp % 平滑處理 colormap('jet') % 顏色屬性 set(gca,'zdir','reverse'); % ‘’都是屬性調(diào)整 axis equal grid on subplot(223) slice(x,y,z,v,[],[3],[]); shading interp set(gca,'zdir','reverse'); axis equal grid on subplot(224) slice(x,y,z,v,[],[4],[]); shading interp set(gca,'zdir','reverse'); axis equal grid on% 2.全空間立體切片 h2=figure(2); set(h2,'name','全空間切片','MenuBar','none','ToolBar','none') slice(x,y,z,v,[1:2:27],[2 3 4],[2 3 4 5]) % x方向1-27，以2位單元，y方向 分成3分，z方向分成4分 shading interp % 圖像平滑 colorbar % 顏色條 colormap('jet') % 顏色屬性 set(gca,'zdir','reverse'); % 屬性調(diào)整 axis equal grid on box on % 顯示坐標(biāo)區(qū)輪廓 顯示圍繞當(dāng)前坐標(biāo)區(qū)的框輪廓% 3.立體包絡(luò)圖 h3=figure(3); subplot(221); set(h3,'name','定值包絡(luò)立體圖','MenuBar','none','ToolBar','none') set(gcf,'InvertHardcopy','off') fw=150; % fw是作主要可變參數(shù)，控制最外層包絡(luò)面的值 %%此值為最外層包絡(luò)面取值 fv=isosurface(x,y,z,v,fw); % 從三維體數(shù)據(jù)中提取等值面數(shù)據(jù) % isovalue 中指定的等值面值處的體數(shù)據(jù) V 計算等值面數(shù)據(jù)。等值面連接具有指定值的點，與等高線連接仰角相同的點的方式大致相同 p=patch(fv); % 創(chuàng)建一個或多個填充多邊形 set(p,'facecolor','b','edgecolor','none'); patch(isocaps(x,y,z,v, fw), 'FaceColor', 'interp', 'EdgeColor', 'none');% isocaps 計算等值面端帽幾何圖 % colorbar % colormap('jet') box on daspect([1,1,1]) % 坐標(biāo)軸尺寸變化為1:1:1 view(3) % 觀看角度 三維視角 set(gca,'zdir','reverse','color',[0.2,0.2,0.2]); camlight % 在照相機(jī)坐標(biāo)系中創(chuàng)建或移動光源對象 camproj perspective % 設(shè)置或查詢投影類型 lighting phong % 指定光照算法 axis equal grid on % title(['最外層表面的值為： ' , num2str(fw)]);% 三維投影 hold on subplot(222); set(h3,'name','定值包絡(luò)立體圖','MenuBar','none','ToolBar','none') set(gcf,'InvertHardcopy','off') fw=150; % fw是作主要可變參數(shù)，控制最外層包絡(luò)面的值 %%此值為最外層包絡(luò)面取值 fv=isosurface(x,y,z,v,fw); % 從三維體數(shù)據(jù)中提取等值面數(shù)據(jù) %isovalue 中指定的等值面值處的體數(shù)據(jù) V 計算等值面數(shù)據(jù)。等值面連接具有指定值的點，與等高線連接仰角相同的點的方式大致相同 p=patch(fv); % 創(chuàng)建一個或多個填充多邊形 set(p,'facecolor','b','edgecolor','none'); patch(isocaps(x,y,z,v, fw), 'FaceColor', 'interp', 'EdgeColor', 'none');% isocaps 計算等值面端帽幾何圖 % colorbar % colormap('jet') box on daspect([1,1,1]) % 坐標(biāo)軸尺寸變化為1:1:1 view(90,0); %在yoz面上的投影 set(gca,'zdir','reverse','color',[0.2,0.2,0.2]); camlight % 在照相機(jī)坐標(biāo)系中創(chuàng)建或移動光源對象 camproj perspective % 設(shè)置或查詢投影類型 lighting phong % 指定光照算法 axis equal grid on % title(['最外層表面的值為： ' , num2str(fw)]);hold on subplot(223); set(h3,'name','定值包絡(luò)立體圖','MenuBar','none','ToolBar','none') set(gcf,'InvertHardcopy','off') fw=150; % fw是作主要可變參數(shù)，控制最外層包絡(luò)面的值 %%此值為最外層包絡(luò)面取值 fv=isosurface(x,y,z,v,fw); % 從三維體數(shù)據(jù)中提取等值面數(shù)據(jù) %isovalue 中指定的等值面值處的體數(shù)據(jù) V 計算等值面數(shù)據(jù)。等值面連接具有指定值的點，與等高線連接仰角相同的點的方式大致相同 p=patch(fv); % 創(chuàng)建一個或多個填充多邊形 set(p,'facecolor','b','edgecolor','none'); patch(isocaps(x,y,z,v, fw), 'FaceColor', 'interp', 'EdgeColor', 'none');% isocaps 計算等值面端帽幾何圖 ....................................

?

D202

?

總結(jié)

以上是生活随笔為你收集整理的【三维重建】基于切片的三维重建MATLAB仿真的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。