?%1.重建數據的采集
?image1=imread（‘data1.jpg‘）;%使用imread（）函數讀入現有的n幅圖像
?image2=imread（‘data2.jpg‘）;
?image3=imread（‘data3.jpg‘）;
?image4=imread（‘data4.jpg‘）;
?image5=imread（‘data5.jpg‘）;
?image6=imread（‘data6.jpg‘）;
?D=cat（3image1image2image3image4image5image6）;%使用cat（）函數創建三維矩陣D
?%2.重建數據項預處理
?
?[x?y?z?D]=reducevolume（D[a?b?c]）;%使用reducevolume（）函數減少數據量，其中a，b，c為
???????????????????????????%x，y，z軸數據抽取的比例，根據數據情況自行定義
?D=smooth3（D）;%使用smooth（）函數對數據驚醒平滑處理
?%3.計算數據集在顯示平面累計投影
?
?fv=isosurface（xyzDisovalue）;%使用isosurface（）函數計算數據集在顯示平面累計投影，
????????????????????????????????%isovalue根據實際情況自行定義
?%4.構造三維體重建碎片
?p=patch（fv‘FaceColor‘‘yellow‘‘EdgeColor‘‘none‘）;%使用patch（）函數對碎片進行構造
????????????????????????????????????????????%并對圖像的顏色，光線進行定義，其中fv是
???????????????????????????????????????????%第三步中得到的
?%5.設置圖像的顏色，陰影及顯示效果
?colormap（gray）;%利用colormap（）函數為圖像定義顏色集
?view（3）;%利用view（）函數定義觀察者視角
?lighting?gouraud;%利用lighting定義顯示圖像的光線陰影
?axis?equal;%利用axis定義顯示圖像的軸
?daspect（[x?y?z]）;%使用daspect（）定義x，y，z軸的顯示比例

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1278??2011-04-25?11:07??頭部CT圖像三維體重建的過程\頭部CT圖像三維體重建的過程.m

?????目錄??????????0??2011-04-26?10:00??頭部CT圖像三維體重建的過程

-----------?---------??----------?-----??----

?????????????????1278????????????????????2