Moyamoya Ua Duab: 8 Kauj Ruam

2024 Tus sau: John Day | [email protected]. Kawg hloov kho: 2024-01-30 09:29

Moyamoya, "pa luam yeeb," yog kab mob uas tsis tshua muaj tshwm sim los ntawm kev thaiv cov hlab ntsha ntawm lub hauv paus ganglia, uas yog thaj tsam ntawm lub hauv paus ntawm lub hlwb. Tus kab mob no yog kab mob cerebrovascular zuj zus uas feem ntau cuam tshuam rau menyuam yaus. Cov tsos mob ntawm moyamoya suav nrog qhov pib mob stroke, tas li mini-strokes, nqaij leeg tsis muaj zog, tuag tes tuag taw, lossis qaug dab peg vim yog cov hlab ntsha loj zuj zus. Yog tsis kho, moyamoya yuav ua rau muaj teeb meem nrog kev hais lus, kev nkag siab tsis zoo, thiab kev nco qab tsis zoo. Hauv peb txoj haujlwm, peb yuav siv MATLAB txhawm rau txhawm rau ua tiav cov duab MRI lossis MRA siv ntau lub ntxaij lim dej kom txo suab nrov hauv cov duab los nrhiav thaj chaw cuam tshuam. Ib qho ntxiv, peb yuav siv qhov kev txhim kho kom pom thaj chaw cuam tshuam ntau dua. Ntxiv mus, peb yuav tom qab ntawd ntsuas kev ywj pheej t-kuaj txhawm rau txiav txim siab seb puas muaj qhov sib txawv ntawm cov hlab ntshav hauv lub hlwb ib txwm piv rau lub hlwb cuam tshuam los ntawm moyamoya.

Kauj Ruam 1: Nrhiav MRI thiab MRA Txheeb Xyuas Lub Hlwb Li Cas thiab Lub Hlwb cuam tshuam nrog Moyamoya

Cov duab no yog qhov ntsuas peb siv rau qhov haujlwm uas peb pom hauv online. Ob daim duab nrog cov hlab ntshav nyob hauv nruab nrab yog MRA scan, thaum ob daim duab yog MRI ntsuas.

Cov kev txuas hauv qab no yog qhov pom cov duab no:

static.cambridge.org/resource/id/urn:cambr…

static.cambridge.org/resource/id/urn:cambr…

Kauj Ruam 2: Thauj Cov Duab Los Rau MATLAB thiab Muab Cov Duab rau Hloov pauv los tso saib cov duab

Txhawm rau pib txheej txheem, pib los ntawm kev tshem lub qhov rai hais kom ua, kaw txhua qhov ua tau thiab cov duab uas tuaj yeem tau qhib thiab pom qhov hloov pauv uas twb tau muab rau hauv koj qhov chaw ua haujlwm.

Tom qab, tsim lub voj los ntawm 1 txog 2 siv cov lus txib i = [1: 2]

Tom qab, thauj cov MRA cov duab siv cov lus txib imread (sprintf ('filename%.filetype', i)) los nyeem cov duab los ntawm cov ntaub ntawv teev tseg los ntawm cov npe ntaub ntawv ua raws tus lej siv rau lub voj siv sprintf rau kev thauj khoom thiab muab nws mus rau qhov txawv txav.

Tom qab ntawd los tso saib cov duab hauv daim duab, siv qhov hais kom ua imshow (I).

Txhawm rau muab cov xim greormap, siv cov lus txib colormap (grey).

Txhawm rau tshem tawm cov xim tag nrho thiab hloov pauv 3D matrix rau cov duab rau hauv 2D, siv cov lus txib rgb2gray (I) thiab muab nws rau qhov sib txawv sib txawv.

Tom qab ntawd thauj cov duab MRI los ntawm kev siv cov lus txib yav dhau los tau teev tseg los yog nyeem (sprintf ('filename%.filetype', i)) thiab muab nws rau qhov hloov pauv tshiab

Rov hais dua rgb2gray hais kom ua nrog qhov hloov pauv tshiab siv rau MRI cov duab.

Yog tias tsim nyog, koj tuaj yeem hloov kho cov duab siv cov lus txib imresize (A, nplai) thiab muab faib rau qhov sib txawv sib txawv.

Kauj Ruam 3: Txhim Kho Cov Qauv Ntev hauv Kev Siv Duab Nrog Ntau Lub Lim Piam

Siv qhov hloov pauv tshiab, siv cov lus txib fibermetric (A) txhawm rau txhim kho cov tubular cov qauv hauv cov duab

Nrog rau qhov sib txawv dhau los, siv qhov hais kom ua histeq (B) txhawm rau txhim kho qhov sib npaug histogram los ntawm kev hloov pauv qhov siv ntawm cov duab thiab muab nws rau qhov hloov pauv tshiab.

Tso cov histogram siv cov lus txib imhist (B)

Tsim qhov hloov pauv tshiab los tsim qhov pib rau lub lim. Hauv qhov no muab qhov sib txawv yav dhau los> 0.875, lim tawm pixel siv nyob rau hauv tus nqi ntawm 0.875

Tom qab, tsim daim duab tshiab thiab siv cov lus txib imshow (A) los tso saib cov duab lim tshiab.

Kauj Ruam 4: Khiav 2D Median Lim

Siv kev hais kom ua medfilt2 (A, [m n]), khiav 2D lim dej nruab nrab, qhov twg txhua qhov pixel tso tawm muaj tus lej nruab nrab hauv mxn ciam teb nyob ib puag ncig cov pixel hauv cov duab nkag.

Tsim tus lej tshiab thiab siv imshow (A) los tso saib cov duab nruab nrab lim.

Kauj Ruam 5: Npog Daim Duab

Siv cov duab nruab nrab lim, siv qhov hais kom ua [labeledImage, numberOfBlots] = bwlabel (A) kom suav tus naj npawb ntawm cov kab dawb hauv daim duab

Tom qab ntawd, siv thaj tsam cov khoom siv ua haujlwm hauv xeev = thaj tsam cov khoom siv (daim ntawv loImage, 'Cheeb Tsam') txhawm rau xam thaj tsam ntawm txhua qhov blot lossis cov hlab ntsha

Muab txhua qhov chaw rau hauv ib qho sib txawv

Tom qab ntawd siv lwm qhov sib txawv, suav cov blots uas ntau dua 50 pixels

Tom qab, txheeb tawm ib qho blots uas qis dua 50 pixels hauv kev nqis nqis los siv kev hais kom ua [sortedAreas, sortedIndicies] = txheeb (thaj tsam, 'nqis')

Tom qab ntawd, siv lwm qhov sib txawv, siv cov lus txib ismember (labeledImage, sortedIndicies (1: numberToExtract)) kom rov qab cov kab nrog cov ntsiab lus ntawm labeledImage tau pom nyob hauv sortedIndicies los ntawm tus lej 1 mus rau tus naj npawb ntawm cov hlab ntshav kom rov qab muaj qhov laj thawj 1 (tseeb) lossis logical 0 (tsis tseeb).

Nrog rau qhov sib txawv hauv cov kauj ruam dhau los, nrhiav cov ntsiab lus uas muaj tseeb (qhov muaj txiaj ntsig> 0) thiab tsim cov lus sib cav los ua cov duab binary thiab muab nws hloov pauv tshiab.

Tsim tus lej tshiab thiab siv imshow (A) cov duab binary tshiab.

Tom qab ntawd, hloov cov duab siv cov lus txib ua tiav (A) thiab muab nws rau qhov sib txawv sib txawv.

Txhawm rau tsim lub ntsej muag duab, siv qhov hloov pauv tshiab nrog hais kom ua resizedimage.*Uint8 (invertedimage)

Tsim tus lej tshiab thiab siv imshow (A) los tso saib daim duab ntsej muag.

Txhawm rau tag nrho cov cai, nco ntsoov siv cov lus txib 'kawg' kom xaus tag nrho rau lub voj

Kauj Ruam 6: Xaiv MRA Scans rau Kev Xeem Txog Kev Txheeb Ze

Txhawm rau npaj rau kev ntsuas tus lej, xaiv MRA qhov kev tshuaj xyuas uas yuav siv rau qhov kev kuaj tus kheej t-test. Vim tias peb ob qhov kev kuaj yuav yog lub hlwb cuam tshuam rau Moyamoya, thiab lub hlwb zoo ib yam, xaiv qhov tsim nyog los ntsuas MRA ntawm txhua pab pawg.

Kauj Ruam 7: Xam Thaj Chaw Ntawm Cov Roj Ntshav Hauv Kev Npaj rau Kev Txheeb Zej Zog

Qhov kev ntsuas ntsuas yuav tsom mus rau qhov ntev lossis ntau npaum li cas ntawm cov hlab ntshav uas qhia ntawm MRA kev tshuaj xyuas. Yog li, peb yuav tsum suav thaj tsam ntawm cov hlab ntshav ua ntej sib piv.

Pib los ntawm kev lim cov MRAs ntawm lub hlwb ib txwm thiab suav tus nqi ntawm cov hlab ntshav. Txhawm rau ua qhov no, khiav ib lub voj. Vim tias muaj peb daim duab, qhov xwm txheej yuav yog kuv = [1: 3].

Qhib cov duab nrog cov lus txib nyeem thiab muab nws rau qhov sib txawv.

Tom ntej no, tsim ib qho if/else nqe lus nrog yog, lwm qhov hais kom ua. Rau qhov hais tawm, siv qhov hais kom ua me me (A, 3) == 3, qhov twg A yog qhov sib txawv siv los qhib rau cov duab, los tsim cov lus yog tias thaum twg qhov thib peb qhov loj ntawm cov array yog 3. Tom qab ntawd hloov cov duab rau 2D thiab tshem tawm cov xim siv qhov hais kom ua rgb2gray (A) thiab muab nws rau qhov hloov pauv tshiab. Siv cov lus txib imresize (A, [m n]) kom hloov kho cov duab. Hauv qhov no, peb hloov kho cov duab rau 1024 x 1024 matrix. Txhawm rau txhim kho cov qauv tubular ntawm daim duab, siv cov lus txib fibermetric dua, thiab muab nws tso rau qhov hloov pauv tshiab.

Cov hauv qab no yog rau lwm nqe lus. Yog tias daim duab tsis yog 3D matrix, peb xav hla kev hloov pauv. Ua tib yam li cov lus hais, tab sis tsis muaj rgb2gray (A) hais kom ua.

Tsim qhov hloov pauv tshiab, teeb tsa nws sib npaug ntawm qhov sib txawv los ntawm qib fibermetric ntau dua 0.15. Qhov no pib qhov duab rau qhov hnyav ntau dua 0.15.

Peb yuav rov ua cov kab ntawm cov lej los ntawm cov kauj ruam 4 thiab 5 ntawm cov lus qhia los ntawm kab lim nruab nrab kom txog thaum kab imshow (I). Tom qab, siv cov lus txib suav (Kuv (:)) txhawm rau ntxiv tag nrho cov pixels uas ua rau cov hlab ntshav, thiab muab nws rau qhov sib txawv sib txawv. Npe lub npe hloov pauv tshiab NormalBloodVessels (i) thiab teeb nws sib npaug rau qhov sib txawv los ntawm cov lej (Kuv (:)) hais kom ua. Qhov no ntxiv cov ntaub ntawv rau hauv matrix.

Xaus lub voj thiab rov ua dua tab sis rau MRAs ntawm lub hlwb cuam tshuam los ntawm Moyamoya. Lub npe qhov sib txawv ntawm qhov kawg MoyaMoyaBloodVessels (i) kom tsis txhob cuam tshuam nrog lub hlwb qub li MRAs.

Kauj Ruam 8: Khiav Ib Qhov Kev Xeem T-xeem

Txij li muaj ob qho kev kuaj ywj pheej, thiab cov neeg tsawg, ua qhov kev kuaj ywj siab t-kuaj.

Tsim thiab lub npe muaj nuj nqi uas ua qhov kev kuaj ywj siab t-kuaj los txiav txim seb cov ntshav cov ntshav hauv MRAs ntawm lub paj hlwb li qub puas sib npaug lossis tsis yog rau MRAs ntawm lub paj hlwb cuam tshuam los ntawm Moyamoya.

Zaub qhov kev xav uas tau teeb tsa rau kev sim los ntawm kev siv qhov hais kom ua disp ('X'). Ntawm kab thawj, tso saib, "Cov kev xav rau ob tus qauv t kuaj." Ntawm kab thib ob, tso saib, "H0 = Tus nqi ntawm cov hlab ntshav ntawm lub hlwb ib txwm sib npaug ntawm cov hlab ntshav ntawm lub hlwb nrog tus kab mob Moyamoya," los hais qhia qhov kev xav tsis meej. Ntawm kab thib peb, tso saib, "HA = Tus nqi ntawm cov hlab ntshav ntawm lub hlwb ib txwm tsis sib npaug ntawm cov hlab ntshav ntawm lub hlwb nrog tus kab mob Moyamoya." los hais txog lwm txoj kev xav.

Siv 95% kev ntseeg siab ib ntus thiab qhov piv txwv ntawm 3, suav cov qhab nia t los ntawm kev siv cov lus txib tinv ([0.025 0.975], 2) thiab muab rau qhov sib txawv t. Siv cov lus txib var (NormalBloodVessels) thiab var (MoyaMoyaBloodVessels) thiab muab lawv rau cov lej sib txawv los xam qhov sib txawv ntawm ob qho ntaub ntawv teev tseg.

Txheeb xyuas tias qhov sib txawv nyob ze rau qhov sib npaug lossis tsis. Ua qhov no los ntawm kev tsim ib qho if/else nqe lus nrog yog tias, lwm qhov hais kom ua. Txog qhov xwm txheej hauv daim ntawv yog, sau A / B == [0.25: 4], qhov twg A yog qhov sib txawv uas suav nrog qhov sib txawv ntawm cov hlab ntsha ib txwm muaj thiab B yog qhov sib txawv uas suav nrog qhov sib txawv ntawm Moyamoya cov hlab ntsha. 0.25 thiab 4 los ntawm kev kwv yees dav rau kev txiav txim siab seb qhov sib txawv sib npaug los tsis yog. Tom qab ntawd ua ob qhov qauv t kuaj nrog [h, p] = ttest2 (A, B, 0.05, 'ob qho tib si', 'sib npaug'), nrog A thiab B yog tib qho kev hloov pauv raws li tau hais ua ntej. Rau lwm nqe lus, siv [h, p] = ttest2 (A, B, 0.05, 'ob qho tib si', 'tsis sib xws') los ua ob qho qauv t kuaj hauv rooj plaub uas qhov sib txawv tsis sib npaug. Xaus qhov if/else nqe lus. Qhov no yuav suav p.

Tsim ib qho/lwm nqe lus uas yuav nthuav tawm qhov kev txiav txim siab raws li tus nqi p. Qhov xwm txheej rau yog tias nqe lus yuav yog p> 0.05. Txij li thaum peb ib txwm tsis ua raws li qhov kev xav tsis zoo thaum tus nqi p ntau dua 0.05, siv qhov tsis txaus siab ('X') hais kom ua "Vim tias p tus nqi ntau dua 0.05 peb tsis kam lees qhov kev xav tsis raug," thiab " Yog li peb tsis lees paub tias cov hlab ntshav ntawm lub hlwb ib txwm sib npaug ntawm lub hlwb nrog Moyamoya kab mob. " Hauv lwm nqe lus, txij li peb ib txwm tsis lees paub qhov kev xav tsis zoo thaum tus nqi p qis dua 0.05, siv qhov tsis txaus siab ('X') hais kom ua "Vim tias p tus nqi tsawg dua 0.05 peb tsis lees txais qhov kev xav tsis raug," thiab " Yog li peb tsis lees paub tias cov hlab ntshav ntawm lub hlwb ib txwm tsis sib npaug nrog lub hlwb nrog Moyamoya kab mob. " Xaus qhov if/else nqe lus.