©1994-2010ChinaAcademicJournalElectronicPublishingHouse.Allrightsreserved.()Vol.26No.4200812JournalofHubeiUniversityforNationalities(NaturalScienceEdition)Dec.2008:2008-09-11.:(1980-),(),.MATLAB(,445000):,MATLAB,,M.100.:;;MATLAB;:O242.1;TP391.41;R322.1+2:A:1008-8423(2008)04-0398-033DReconstructionofBloodVesselBasedonMATLABLIUBo(DepartmentofMathematics,HubeiUniversityforNationalities,Enshi445000,China)Abstract:Accordingtothebinaryimagefeaturesofbloodvesselslices,thethree-dimensionalrecon2structionprocessofbloodvesselwasdecomposedtoradiussearch,intersectionlocationandaxis-fittingthroughusingmatrixoperationandimageryprocessingofMATLABinstructions.AgeneralMprogrampackagewascompiledtoachievethegoalofautocalculationofwholecoursefromdataacquisitiontomod2elrendering.Atlast,3Dreconstructionfor100slicesofbloodvesseluseoftheMprogrampackagewasrealizedbymaking.Keywords:bloodvesselslices;binaryimage;MATLAB;three-dimensionalreconstruction,,,[1].,,,,.[2][3,4][5,6],,;,;.,,.[2,4,5][3][6],,;,;.,,.1,r,XOY©1994-2010ChinaAcademicJournalElectronicPublishingHouse.Allrightsreserved.=l/2,.XOZYOZ,.r.2,,,,:2.1(2001),,imread,.(1),,lh,00,=0,lh,(2).1,01,A,DA1,2=0Ak,:A(i,j)=1,A(i,j)D0,A(i,j)D(1)l=sum(a(k,))=min{sumA(m,),,sumA(m,))},mknsumA(m,)Ak,,.60,r=30.,area(3).2.2[7]:,;,(4).,(5).345Fig.3TheminimumaccumulatedvaluelinebylineFig.4ThemaximuminscribedcircleofsliceFig.5Coordinatesofpointofintersection,r,.A,DA19934:MATLAB©1994-2010ChinaAcademicJournalElectronicPublishingHouse.Allrightsreserved.(i,j)D,A(i,j)r312D,:(A(i,j)D)&&(A(i+r33sin,j+r33cos)D)=1=linspace(0,23pi,13)(2)A(i,j),.r3,.,r3r,,(2),.2.3,polyfitXOZYOZ(6,7),,.:x=1.254910-10z7-3.560410-8z6+3.473210-6z5-0.0001388z4+0.0027936z3-0.024197z2+0.0010758z+96.21y=3.527810-13z8-2.936910-10z7+8.08610-8z6-1.015710-5z5+0.0062806z4-0.01941z3+0.33812z2-2.1325z+261.2577(3):x(0,512),y(0,512),z(0,99).surface(3)r=30,,l8.6XOZ7YOZ8Fig.6ThecurvefittinginXOZplaneFig.7ThecurvefittinginYOZplaneFig.8Thesimulationgraphicof3Dreconstructionofbloodvessel3,,MATLAB,M,.100,MATLAB.,.,..:[1],.[J].,2002,19(5):54-58.[2],,.[J].,2002,19(5):22-28.[3],,.[J].,2002,19(5):29-34.[4],,.[J].,2002,19(5):35-40.[5],,.[J].,2002,19(5):41-46.[6],,.[J].,2002,19(5):47-53.[7],.[J].(A),1998,13(S1):86-90.004()26