%腦部MRI圖像的三維重建及動(dòng)畫(huà)顯示
%brainavi．m?
clear%清除內(nèi)存?
clc%清除屏幕?

Figwin=figure（‘position‘?[50?50?450?450]‘Name‘??‘腦部三維重建及旋轉(zhuǎn)動(dòng)畫(huà)演示?‘‘Numbertitle‘??‘off‘?‘Menubar‘‘none‘?）;?
%產(chǎn)生標(biāo)題為”腦部?維重建及旋轉(zhuǎn)動(dòng)畫(huà)演示”的圖形窗?
%口?
%%%1．讀入腦部MRl圖像%%%
load?mri??%調(diào)入腦部MRI圖像數(shù)據(jù)?
D=squeeze（D）;?%將D從4維轉(zhuǎn)換為3維?
Ds=smooth3（D）;%采用高斯低通濾波器對(duì)D進(jìn)行平?
%滑，得Ds?
%%%2．腦部三維重建與旋轉(zhuǎn)動(dòng)畫(huà)顯示%%%?
fv=isosurface（Ds20）;?
%腦部等值面抽取，閾值k=20，見(jiàn)公式（1）。fv是一個(gè)結(jié)構(gòu)?
%數(shù)組，其中?
%fv．veaices為圖形的頂點(diǎn)信息；fv．faces為圖形的表面信?
%息?
fv2=isocaps（D5）;%腦部上蓋的等值面抽取，閩值k=5?
yuan=fv.vertices;?%令yuan為原腦部圖形的頂點(diǎn)信息?
yuan2=fv2.vertices;%令yuan2為原腦部上蓋圖形的頂點(diǎn)?
%信息?
N=length（yuan）;%N和N2分別為yuan和yuan2的?
%像素個(gè)數(shù)?
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的質(zhì)

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

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

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

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

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