%腦部MRI圖像的三維重建及動畫顯示
%brainavi．m?
clear%清除內存?
clc%清除屏幕?

Figwin=figure（‘position‘?[50?50?450?450]‘Name‘??‘腦部三維重建及旋轉動畫演示?‘‘Numbertitle‘??‘off‘?‘Menubar‘‘none‘?）;?
%產生標題為”腦部?維重建及旋轉動畫演示”的圖形窗?
%口?
%%%1．讀入腦部MRl圖像%%%
load?mri??%調入腦部MRI圖像數據?
D=squeeze（D）;?%將D從4維轉換為3維?
Ds=smooth3（D）;%采用高斯低通濾波器對D進行平?
%滑，得Ds?
%%%2．腦部三維重建與旋轉動畫顯示%%%?
fv=isosurface（Ds20）;?
%腦部等值面抽取，閾值k=20，見公式（1）。fv是一個結構?
%數組，其中?
%fv．veaices為圖形的頂點信息；fv．faces為圖形的表面信?
%息?
fv2=isocaps（D5）;%腦部上蓋的等值面抽取，閩值k=5?
yuan=fv.vertices;?%令yuan為原腦部圖形的頂點信息?
yuan2=fv2.vertices;%令yuan2為原腦部上蓋圖形的頂點?
%信息?
N=length（yuan）;%N和N2分別為yuan和yuan2的?
%像素個數?
N2=length（yuan2）;?
xg=sum（yuan（:?1））/N;yg=sum（yuan（:?2））/N;?
zg=sum（yuan（:?3））/N;
xg2=sum（yuan2（:?1））/N2;
yg2=sum（yuan2（:?2））/N2;?
zg2=sum（yuan2（:?3））/N2;?
%求yuan和yuan2的質

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

?????文件?????213940??2012-03-06?21:14??基于Matlab的腦部MRI圖像三維重建與動畫顯示.pdf

?????文件???????3216??2012-03-07?14:37??brainavi.m

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

???????????????217156????????????????????2