Hisoblash anatomiyasining Bayes modeli - Bayesian model of computational anatomy

Hisoblash anatomiyasi (CA) ichidagi intizomdir tibbiy tasvir anatomik shakl va shaklni ko'rinadigan yoki ko'rinishda o'rganishga e'tibor qaratish yalpi anatomik ko'lami morfologiya. Ushbu soha keng tavsiflangan va tarkibida poydevor mavjud anatomiya, amaliy matematika va sof matematika, shu jumladan tibbiy tasvir, nevrologiya, fizika, ehtimollik va statistika. Bu tibbiy tasvirlash moslamalariga emas, balki tasvirlanadigan anatomik tuzilmalarga qaratilgan. Sub-maydonining markaziy yo'nalishi hisoblash anatomiyasi ichida tibbiy tasvir anatomik koordinata tizimlari bo'yicha ma'lumotlarni xaritalashdir, ko'pincha a ichida o'lchanadigan zich ma'lumot magnit-rezonansli tasvir (MRI). Suyuqlik dinamikasida ishlatiladigan harakat tenglamalariga o'xshash oqimlarni CA ga kiritish, tasvirni tahlil qilishda zich koordinatalar quyidagi tushunchadan foydalanadi. Lagrangian va Eylerian harakat tenglamalari. Dfefeomorfizmlarning Lagranjian va Eulerian oqimlariga asoslangan modellarda cheklov topologik xususiyatlar bilan bog'liq, masalan, ochiq to'plamlar saqlanib qolinadi, koordinatalar teskari xaritalashning o'ziga xosligi va mavjudligini anglatmaydi va bog'langan to'plamlar bog'lanib qoladi. Diffeomorfik usullardan foydalanish tezda o'sib bordi va Kristensendan keyin xaritalash usullari sohasida ustunlik qildi^[1]tez va nosimmetrik usullar mavjud bo'lgan asl qog'oz.^[2]^[3]

Asosiy statistik model

Shakllar manbasini deformatsiya qilinadigan shablonni ko'rsatadigan manba kanal modeli

{ displaystyle I doteq varphi cdot I _ { mathrm {temp}} in { mathcal {I}}}

va MRI sensori bilan bog'liq kanal chiqishi

{ displaystyle I ^ {D} in { mathcal {I}} ^ { mathcal {D}}}

Kontekstida hisoblash anatomiyasining markaziy statistik modeli tibbiy tasvir ning manba-kanal modeli bo'lgan Shannon nazariyasi; manba tasvirlarning deformatsiyalanadigan shablonidir ${ displaystyle I in { mathcal {I}}}$ , kanal chiqishi - bu kuzatiladigan narsalarga ega bo'lgan tasviriy sensorlar ${ displaystyle I ^ {D} in { mathcal {I}} ^ { mathcal {D}}}$ (rasmga qarang). Manba-kanal modelining ahamiyati shundaki, anatomik konfiguratsiyadagi o'zgarish tibbiy tasvirlarning sensorli o'zgarishlaridan ajratilgan holda modellashtirilgan. The Bayes nazariyasi model manba bo'yicha avvalgi bilan tavsiflanadi, ${ displaystyle pi _ { mathcal {I}} ( cdot)}$ kuni ${ displaystyle I in { mathcal {I}}}$ va kuzatiladigan narsaning shartli zichligi

{ displaystyle p ( cdot mid I) { text {on}} I ^ {D} in { mathcal {I}} ^ { mathcal {D}}}

shartli ${ displaystyle I in { mathcal {I}}}$ .

Deformatsiyalanadigan shablon nazariyasida rasmlar shablonga bog'langan bo'lib, shablonga ta'sir qiladigan guruh deformatsiyalari bilan; qarang hisoblash anatomiyasida guruh harakati Tasvir harakati uchun ${ displaystyle I (g) doteq g cdot I _ { mathrm {temp}}, g in { mathcal {G}}}$ , keyin guruhda oldingi ${ displaystyle pi _ { mathcal {G}} ( cdot)}$ tasvirlarda oldingi holatni keltirib chiqaradi ${ displaystyle pi _ { mathcal {I}} ( cdot)}$ , zichlik sifatida yozilgan log-posterior shaklni oladi

{ displaystyle log p (I (g) mid I ^ {D}) simeq log p (I ^ {D} mid I (g)) + log pi _ { mathcal {G}} (g).}

Quyidagi tasodifiy orbitaning modeli guruh elementlarini qanday yaratishni va shuning uchun oldingi taqsimotni tashkil etuvchi narsalarning tasodifiy purkashini belgilaydi.

Hisoblash anatomiyasining tasodifiy orbitali modeli

Shablonlar bo'yicha diffeomorfik guruh ta'siriga bog'liq bo'lgan miyaning orbitalari dastlabki teginish kosmik vektor maydonining tasodifiy hosil bo'lishiga bog'liq bo'lgan tasodifiy purkagich bilan geodezik oqimlarga bog'liq silliq oqim orqali tasvirlangan

{ displaystyle v_ {0} in V}

; yilda nashr etilgan.

The tasodifiy orbitaning modeli Hisoblash anatomiyasi birinchi marta paydo bo'lgan^[4]^[5]^[6] shakllar va shakllarning anatomik orbitasida tasvirlar manbasida tasodifiylikni keltirib chiqaradigan andozalar ustida ishlaydigan guruhning tasodifiyligi bilan bog'liq koordinatalar o'zgarishini modellashtirish va natijada tibbiy ko'rish vositalari orqali kuzatuvlar. Shunaqangi tasodifiy orbitaning modeli guruhdagi tasodifiylik tasvirdagi tasodifiylikni keltirib chiqaradigan maxsus evklid guruhi uchun guruh elementi bo'lgan ob'ektni aniqlash uchun tekshirildi ${ displaystyle g in { mathcal {G}}}$ maxsus evklid guruhi bo'lgan.^[7]

CA-da deformatsiyalanadigan shaklni o'rganish uchun hisoblash anatomiyasida ishlatiladigan yuqori o'lchovli diffeomorfizm guruhlari silliq oqimlar orqali hosil bo'ladi ${ displaystyle varphi _ {t}, t in [0,1]}$ oddiy diferensial tenglamani qondiradigan oqim maydonlarining Lagranj va Eylerian spetsifikatsiyasini qondiradigan:

Lagranj koordinatalari oqimini ko'rsatish

{ displaystyle x in X}

bog'liq vektor maydonlari bilan

{ displaystyle v_ {t}, t in [0,1]}

oddiy differentsial tenglamani qondirish

{ displaystyle { dot { varphi}} _ {t} = v_ {t} ( varphi _ {t}), varphi _ {0} = id}

.

{ displaystyle { frac {d} {dt}} varphi _ {t} = v_ {t} circ varphi _ {t}, varphi _ {0} = operator nomi {id} ;}

(Lagranj oqimi)

bilan ${ displaystyle v doteq (v_ {1}, v_ {2}, v_ {3})}$ vektor maydonlari ${ displaystyle { mathbb {R}} ^ {3}}$ deb nomlangan Evleriya zarrachalarning joylashish tezligi ${ displaystyle varphi}$ oqimning. Vektorli maydonlar funktsiya maydonidagi funktsiyalar bo'lib, ular silliq qilib modellashtirilgan Xilbert 1 uzluksiz hosilaga ega vektor maydonlari bilan bo'shliq. Uchun ${ displaystyle v_ {t} = { dot { varphi}} _ {t} circ varphi _ {t} ^ {- 1}, t in [0,1]}$ , oqimning teskarisi tomonidan berilgan

{ displaystyle { frac {d} {dt}} varphi _ {t} ^ {- 1} = - (D varphi _ {t} ^ {- 1}) v_ {t}, varphi _ { 0} ^ {- 1} = operator nomi {id},}

(Eulerianflow)

va ${ displaystyle 3 marta 3}$ Oqim uchun Jacobian matritsasi ${ displaystyle mathbb {R} ^ {3}}$ sifatida berilgan ${ displaystyle D varphi doteq chap ({ frac { kısmi varphi _ {i}} { qisman x_ {j}}} o'ng).}$

Vektor maydonlari teskari, diffeomorfizmlarning silliq oqishini ta'minlash uchun ${ displaystyle { mathbb {R}} ^ {3}}$ kosmosda kamida 1 marta doimiy farqlanadigan bo'lishi kerak^[8]^[9] ular Hilbert makonining elementlari sifatida modellashtirilgan ${ displaystyle (V, | cdot | _ {V})}$ yordamida Sobolev har bir element bo'lishi uchun teoremalarni joylashtirish ${ displaystyle v_ {i} H_ {0} ^ {3} da, i = 1,2,3,}$ 3 kvadrat bilan integrallanadigan lotinlarga ega. Shunday qilib ${ displaystyle (V, | cdot | _ {V})}$ bir martalik doimiy ravishda farqlanadigan funktsiyalarga muammosiz joylashtiring.^[8]^[9] Diffeomorfizm guruhi Sobolev normasida mutlaqo integrallanadigan vektor maydonlari bo'lgan oqimlar:

{ displaystyle operator nomi {Diff} _ {V} doteq chap { varphi = varphi _ {1}: { dot { varphi}} _ {t} = v_ {t} circ varphi _ {t}, varphi _ {0} = operator nomi {id}, int _ {0} ^ {1} | v_ {t} | _ {V} , dt < infty right }, }

(Diffeomorfizm guruhi)

qayerda ${ displaystyle | v_ {t} | _ {V} ^ {2} doteq int _ {X} Av_ {t} cdot v_ {t} dx}$ bilan ${ displaystyle A}$ chiziqli operator ${ displaystyle A: V mapsto V ^ {*}}$ RKHS normasini belgilash. Integral qachon bo'linmalar bo'yicha integratsiya bilan hisoblanadi ${ displaystyle Av}$ ikkilangan makondagi umumlashtirilgan funktsiya ${ displaystyle V ^ {*}}$ .

Riemann eksponent

In hisoblash anatomiyasining tasodifiy orbitali modeli, butun oqim diffeomorfizmni kodlovchi koordinatalarni hosil qiladigan dastlabki holatga keltiriladi. Dastlabki holatdan ${ displaystyle v_ {0}}$ keyin ga nisbatan geodezik joylashishni aniqlash Riemann metrikasi Euler-Lagranj tenglamasi oqimini hisoblash anatomiyasi hal qiladi. Geodeziyani dastlabki holatdan hal qilish ${ displaystyle v_ {0}}$ deb nomlanadi Riemann-eksponent, xaritalash ${ displaystyle operatorname {Exp} _ { mathrm {id}} ( cdot): V to operatorname {Diff} _ {V}}$ guruhga shaxsni aniqlashda.

Riemann eksponentligi qondiradi ${ displaystyle operatorname {Exp} _ { mathrm {id}} (v_ {0}) = varphi _ {1}}$ dastlabki holat uchun ${ displaystyle { dot { varphi}} _ {0} = v_ {0}}$ , vektor maydonining dinamikasi ${ displaystyle { dot { varphi}} _ {t} = v_ {t} circ varphi _ {t}, t in [0,1]}$ ,

klassik tenglama uchun diffeomorfik shakl impulsi ${ displaystyle int _ {X} Av_ {t} cdot w , dx}$ , ${ displaystyle Av in V}$ , keyin

{ displaystyle { frac {d} {dt}} Av_ {t} + (Dv_ {t}) ^ {T} Av_ {t} + (DAv_ {t}) v_ {t} + ( nabla cdot v ) Av_ {t} = 0 ;}

umumlashtirilgan tenglama uchun, keyin ${ displaystyle Av in V ^ {*}}$ , ${ displaystyle w in V}$

{ displaystyle int _ {X} { frac {d} {dt}} Av_ {t} cdot w , dx + int _ {X} Av_ {t} cdot ((Dv_ {t}) w- (Dw) v_ {t}) , dx = 0.}

U butun guruhga tarqaladi, ${ displaystyle varphi = operator nomi {Exp} _ { varphi} (v_ {0} circ varphi) doteq operator nomi {Exp} _ { mathrm {id}} (v_ {0}) circ varphi.}$ Qo'shimcha rasmda har bir namunadagi atrofdagi tasodifiy orbitalar tasvirlangan, ${ displaystyle m_ {0} in { mathcal {M}}}$ , identifikatsiyadagi dastlabki teginish kosmik vektor maydonini yaratish orqali oqimni tasodifiy ravishda hosil qilish ${ displaystyle v_ {0} in V}$ va keyin tasodifiy ob'ektni yaratish ${ displaystyle n doteq operator nomi {Exp} _ { mathrm {id}} (v_ {0}) cdot m_ {0} in { mathcal {M}}}$ .

Ikki o'lchovli katakchada joylashtirilgan sintez qilingan subkortikal tuzilmalarning tasodifiy purkagichini ko'rsatadigan rasm, bu sintez uchun momentum uchun ishlatiladigan o'ziga xos funktsiya farqini anglatadi.

Rasmda o'ngdagi multfilm orbitasida ko'rsatilgan, bu vektor maydonlarini tasodifiylashtirish natijasida hosil bo'lgan subkortikal manifoldlarning tasodifiy purkagichidir. ${ displaystyle v_ {0}}$ submanifoldlar tomonidan qo'llab-quvvatlangan. Tasodifiy orbitali model shakllar va tasvirlar oldingisini keltirib chiqaradi ${ displaystyle I in { mathcal {I}}}$ ma'lum bir atlas bilan shartlangan ${ displaystyle I_ {a} in { mathcal {I}}}$ . Buning uchun generativ model o'rtacha maydonni hosil qiladi ${ displaystyle I}$ shablon koordinatalarini tasodifiy o'zgartirish sifatida ${ displaystyle I doteq varphi cdot I_ {a}}$ , koordinatalardagi diffeomorfik o'zgarish geodezik oqimlar orqali tasodifiy hosil bo'ladi.

Ko'p atlasli orbitadagi modeldagi xaritani baholash

Tasodifiy orbitali model shakllar va tasvirlar oldingisini keltirib chiqaradi ${ displaystyle I in { mathcal {I}}}$ ma'lum bir atlas bilan shartlangan ${ displaystyle I_ {a} in { mathcal {I}}}$ . Buning uchun generativ model o'rtacha maydonni hosil qiladi ${ displaystyle I}$ shablon koordinatalarini tasodifiy o'zgartirish sifatida ${ displaystyle I doteq varphi cdot I_ {a}}$ , koordinatalardagi diffeomorfik o'zgarish geodezik oqimlar orqali tasodifiy hosil bo'ladi. Tasodifiy transformatsiyalar bo'yicha avvalgi ${ displaystyle pi _ { mathrm {Diff}} (d varphi)}$ kuni ${ displaystyle operatorname {Diff} _ {V}}$ oqim tomonidan indüklenir ${ displaystyle operatorname {Exp} _ { mathrm {id}} (v)}$ , bilan ${ displaystyle v in V}$ oldin Gauss tasodifiy maydoni sifatida qurilgan ${ displaystyle pi _ {V} (dv)}$ . Datchik chiqishidagi tasodifiy kuzatiladigan moslamalardagi zichlik ${ displaystyle I ^ {D} in { mathcal {I}} ^ {D}}$ tomonidan berilgan

{ displaystyle p (I ^ {D} mid I_ {a}) = int _ {V} p (I ^ {D} mid operatorname {Exp} _ { mathrm {id}} (v) cdot I_ {a}) pi _ {V} (dv) .}

Posteriori taxminiy maksimal (MAP) baholash zamonaviy uchun muhimdir statistik nazariya. Qiziqish parametrlari ${ displaystyle theta in Theta}$ kabi ko'plab shakllarni o'z ichiga oladi (i) kabi kasallik turlari neyrodejenerativ yoki neyro-rivojlanish kasalliklar, (ii) strukturaning turi, masalan, tasvirlarni segmentatsiyalash bilan bog'liq muammolarda kortikal yoki subkortikal tuzilmalar va (iii) populyatsiyalardan shablonni qayta tiklash. Kuzatilgan tasvirni hisobga olgan holda ${ displaystyle I ^ {D}}$ , MAP taxminiyligi eng yuqori darajani oshiradi:

{ displaystyle { hat { theta}} doteq arg max _ { theta in Theta} log p ( theta mid I ^ {D}).}

Bu shartli ehtimollarni hisoblashni talab qiladi ${ displaystyle p ( theta mid I ^ {D}) = { frac {p (I ^ {D}, theta)} {p (I ^ {D})}}}$ . Ko'p sonli atlas orbitasi modeli atlaslarning denumable to'plami ustida tasodifiylashadi ${ displaystyle {I_ {a}, a in { mathcal {A}} }}$ . Orbitadagi tasvirlardagi model aralashmaning ko'p modali taqsimoti shaklini oladi

{ displaystyle p (I ^ {D}, theta) = sum _ {a in { mathcal {A}}} p (I ^ {D}, theta mid I_ {a}) pi _ { mathcal {A}} (a) .}

Shartli Gauss modeli zich tasvirlarda aniq mos kelmaslik va diqqatga sazovor joylarni moslashtirish uchun juda tekshirildi.

Zich emage mosligi

Model ${ displaystyle I ^ {D} (x), x in X}$ shartli ravishda Gauss tasodifiy maydoni shartli, o'rtacha maydon sifatida, ${ displaystyle varphi _ {1} cdot I doteq I ( varphi _ {1} ^ {- 1}), varphi _ {1} Diff_ {V}} da$ . Bir xil xilma-xillik uchun so'nggi nuqta xato shartlari oxirgi nuqta atamasini beruvchi log-shartli (faqat o'rtacha maydonning funktsiyasi) rolini o'ynaydi:

{ displaystyle - log p (I ^ {D} mid I (g)) simeq E ( varphi _ {1}) doteq { frac {1} {2 sigma ^ {2}}} | I ^ {D} -I circ varphi _ {1} ^ {- 1} | ^ {2}.}

(Shartli-Gausscha)

Belgilangan joyga mos kelish

Model ${ displaystyle Y = {y_ {1}, y_ {2}, nuqta }}$ shartli ravishda o'rtacha maydon bilan Gauss ${ displaystyle varphi _ {1} (x_ {i}), i = 1,2, nuqtalar, varphi _ {1} in operator nomi {Diff} _ {V}}$ , belgilarga bog'liq bo'lmagan doimiy shovqin dispersiyasi. Log-shartli (faqat o'rtacha maydonning funktsiyasi) so'nggi nuqta atamasi sifatida qaralishi mumkin:

{ displaystyle - log p (I ^ {D} mid I (g)) simeq operator nomi {E} ( varphi _ {1}) doteq { frac {1} {2 sigma ^ {2 }}} sum _ {i} | y_ {i} - varphi _ {1} (x_ {i}) | ^ {2}.}

Bir nechta atlasga asoslangan xaritani segmentatsiya qilish

Ko'p atlas uchun tasodifiy orbitali model diffeomorfizmlarning guruh ta'siridan hosil bo'lgan ko'plab anatomik orbitalar bo'yicha birlashma sifatida shakllar orbitasini modellashtiradi, ${ displaystyle { mathcal {I}} = textstyle bigcup _ {a in { mathcal {A}}} displaystyle operatorname {Diff} _ {V} cdot I_ {a}}$ , har bir atlas shablonga va oldindan belgilangan segmentatsiya maydoniga ega ${ displaystyle (I_ {a}, W_ {a}), a = a_ {1}, a_ {2}, ldots}$ . parchalanishni MRG koordinatasining anatomik tuzilmalariga kiritish .. juftliklar voksel panjarasi ustida indekslanadi ${ displaystyle I_ {a} (x_ {i}), W_ {a} (x_ {i}), x_ {i} in X subset { mathbb {R}} ^ {3}}$ MRI tasviri va har bir voksel koordinatasining zich yorlig'i bilan. Parsellangan tuzilmalarning anatomik yorlig'i neyroanatomistlarning qo'lda belgilashidir.

Bayes segmentatsiyasi muammosi^[10] o'lchov berilgan ${ displaystyle I ^ {D}}$ o'rtacha maydon va parselatsiya bilan ${ displaystyle (I, W)}$ , anatomik yorliq ${ displaystyle theta doteq W}$ . mustg o'lchangan MRI tasviri uchun taxmin qilinadi. Kuzatiladigan narsalarning o'rtacha maydoni ${ displaystyle I ^ {D}}$ rasm shablonlardan birining tasodifiy deformatsiyasi sifatida modellashtirilgan ${ displaystyle I doteq varphi cdot I_ {a}}$ tasodifiy tanlangan, ${ displaystyle A = a}$ ,. Optimal diffeomorfizm ${ displaystyle varphi in { mathcal {G}}}$ yashiringan va tasodifiy tanlangan shablon tasvirining koordinatalari fon maydonida ishlaydi ${ displaystyle I_ {a}}$ . Bitta atlas berilgan ${ displaystyle a}$ , xulosa qilish ehtimoli modeli qo'shma ehtimol bilan belgilanadi ${ displaystyle p (I ^ {D}, W mid A = a)}$ ; bir nechta atlaslar bilan, ehtimollik funktsiyalarining birlashishi ko'p modali aralash modelini oldingi modellar bo'yicha o'rtacha o'rtacha qiymatini beradi.

Segmentatsiyani xaritasini baholash vositasi ${ displaystyle W_ {a}}$ maksimayzer hisoblanadi ${ displaystyle max _ {W} log p (W mid I ^ {D})}$ berilgan ${ displaystyle I ^ {D}}$ , bu aralashmani barcha atlaslarda o'z ichiga oladi.

{ displaystyle { hat {W}} doteq arg textstyle max _ {W} displaystyle log p (I ^ {D}, W) { text {with}} p (I ^ {D}) , W) = textstyle sum _ {a in { mathcal {A}}} displaystyle p (I ^ {D}, W mid A = a) pi _ {A} (a).}

Miqdor ${ displaystyle p (I ^ {D}, V)}$ bilan ko'p deformatsiyalanadigan atlaslardan kelib chiqish ehtimoli birlashtirilib hisoblab chiqiladi ${ displaystyle pi _ {A} (a)}$ kuzatilgan rasmning o'ziga xos shablon tasviridan rivojlanishining oldingi ehtimoli ${ displaystyle I_ {a}}$ .

Xarita segmentatsiyasini takroriy ravishda echish mumkin kutish-maksimallashtirish algoritm

{ displaystyle W ^ { text {new}} doteq arg max _ {W} int log p (W, I ^ {D}, A, varphi) , dp (A, varphi o'rtada W ^ { text {old}}, men ^ {D}).}

Populyatsiyalar va EM algoritmining hajm shablonlari xaritasini baholash

Aholidan empirik ravishda shablonlarni yaratish bu intizom uchun hamma joyda mavjud bo'lgan asosiy operatsiya bo'lib, Bayman statistikasiga asoslangan bir necha usul submanifoldlar va zich tasvir hajmlari uchun paydo bo'ldi. ${ displaystyle I ^ {D_ {1}}, I ^ {D_ {2}}, dots}$ muammo zich tasvirlar orbitasida shablonni taxmin qilishda ${ displaystyle I in { mathcal {I}}}$ . Ma protsedurasi boshlang'ich gipertemplatni oladi ${ displaystyle I_ {0} in { mathcal {I}}}$ boshlang'ich nuqtasi sifatida va noma'lum bo'lgan diffeomorfizm ostida shablonni orbitada modellashtiradi ${ displaystyle I doteq varphi _ {0} cdot I_ {0}}$ , log-koordinatalari taxmin qilinadigan parametrlari bilan ${ displaystyle theta doteq v_ {0}}$ giper-shablonning geodezik xaritasini aniqlash ${ displaystyle operatorname {Exp} _ { mathrm {id}} (v_ {0}) cdot I_ {0} = I in { mathcal {I}}}$ .

In Hisoblash anatomiyasining Bayesiya tasodifiy orbitasi modeli kuzatilgan MRI tasvirlari ${ displaystyle I ^ {D_ {i}}}$ o'rtacha maydonga ega shartli ravishda Gauss tasodifiy maydoni sifatida modellashtirilgan ${ displaystyle varphi _ {i} cdot I}$ , bilan ${ displaystyle varphi _ {i}}$ shablonning tasodifiy noma'lum o'zgarishi. Xaritani baholash muammosi noma'lum shablonni taxmin qilishdir ${ displaystyle I in { mathcal {I}}}$ kuzatilgan MRI tasvirlarini hisobga olgan holda.

Ma ning zich tasvir uchun protsedurasi dastlabki gipertemplatni oladi ${ displaystyle I_ {0} in { mathcal {I}}}$ boshlang'ich nuqtasi sifatida va noma'lum bo'lgan diffeomorfizm ostida shablonni orbitada modellashtiradi ${ displaystyle I doteq varphi _ {0} cdot I_ {0}}$ . Kuzatiladigan narsalar shartli tasodifiy maydonlar sifatida modellashtirilgan, ${ displaystyle I ^ {D_ {i}}}$ a shartli-Gausscha o'rtacha maydon bilan tasodifiy maydon ${ displaystyle varphi _ {i} cdot I doteq varphi _ {i} cdot varphi _ {0} cdot I_ {0}}$ . MAP tomonidan aniq taxmin qilinadigan noma'lum o'zgaruvchi giper-shablonni xaritalashdir ${ displaystyle varphi _ {0}}$ , Bayes protsedurasi orqali birlashtirilgan boshqa xaritalar noqulay yoki yashirin o'zgaruvchilar deb hisoblanadi. Bu yordamida amalga oshiriladi kutish-maksimallashtirish algoritm.

Orbita modeli noma'lum oqimlarni log-koordinatalariga bog'lash orqali foydalaniladi ${ displaystyle v_ {i}, i = 1, dots}$ Riemann geodezik jurnali va eksponent uchun hisoblash anatomiyasi identifikatsiyadagi teginish fazosidagi dastlabki vektor maydoni, shunday qilib ${ displaystyle operatorname {Exp} _ { mathrm {id}} (v_ {i}) doteq varphi _ {i}}$ , bilan ${ displaystyle operatorname {Exp} _ { mathrm {id}} (v_ {0})}$ giper-shablonni xaritalash.MAP-ni baholash muammosi paydo bo'ladi

{ displaystyle max _ {v_ {0}} p (I ^ {D}, theta = v_ {0}) = int p (I ^ {D}, theta = v_ {0} mid v_ {) 1}, v_ {2}, nuqta) pi (v_ {1}, v_ {2}, nuqta) , dv}

EM algoritmi xaritalashni parametrlashtiruvchi vektor-maydon koordinatalarini to'liq ma'lumot sifatida qabul qiladi, ${ displaystyle v_ {i}, i = 1, dots}$ va shartli-kutishni takroriy ravishda hisoblang

{ displaystyle { begin {case} Q ( theta = v_ {0}; theta ^ { text {old}} = v_ {0} ^ { text {old}}) & = - operator nomi {E } ( log p (I ^ {D}, theta = v_ {0} mid v_ {1}, v_ {2}, dots)) mid I ^ {D}, theta ^ { text {old }}) & = - | ({ bar {I}} ^ { text {old}} - I_ {0} circ operatorname {Exp} _ { mathrm {id}} (v_ {0) }) ^ {- 1}) { sqrt { beta ^ { text {old}}}} | ^ {2} - | v_ {0} | _ {V} ^ {2} end { holatlar}}}

Q-funktsiyani maksimal darajada oshiradigan yangi shablonni hisoblang, sozlash

{ displaystyle theta ^ { text {new}} doteq v_ {0} ^ { text {new}} = arg max _ { theta = v_ {0}} Q ( theta; theta ^) { text {old}} = v_ {0} ^ { text {old}}) = - left | ({ bar {I}} ^ { text {old}} - I_ {0} circ operatorname {Exp} _ { mathrm {id}} (v_ {0}) ^ {- 1}) { sqrt { beta ^ { text {old}}}} right | ^ {2} - | v_ {0} | _ {V} ^ {2}}

Rejim qiymatlari uchun kutilgan qiymatlarni yangilash uchun taxminiy rejimni hisoblang:

{ displaystyle v_ {i} ^ { text {new}} = arg max _ {v: { dot { varphi}} = v circ varphi} - int _ {0} ^ {1} | v_ {t} | _ {V} ^ {2} , dt- | I ^ {D_ {i}} - I_ {0} circ operatorname {Exp} _ { mathrm {id}} (v_ {0} ^ { text {old}}) ^ {- 1} circ operatorname {Exp} _ { mathrm {id}} (v) ^ {- 1} | ^ {2} .i = 1,2, nuqta}

{ displaystyle beta ^ { text {new}} (x) = sum _ {i = 1} ^ {n} | D operator nomi {Exp} _ { mathrm {id}} (v_ {i} ^ { text {new}}) (x) |, { text {with}} { bar {I}} ^ { text {new}} (x) = { frac { sum _ {i = 1 } ^ {n} I ^ {D_ {i}} circ operatorname {Exp} _ { mathrm {id}} (v_ {i} ^ { text {new}}) | D operatorname {Exp} _ { mathrm {id}} (v_ {i} ^ { text {new}}) (x) |} { beta ^ { text {old}} (x)}}}

Adabiyotlar

^ Kristensen, GE; Rabbitt, R.D .; Miller, M.I. (1996-02-01). "Katta deformatsiyaning kinematikasidan foydalangan holda deformatsiyalanadigan shablonlar". Rasmni qayta ishlash bo'yicha IEEE operatsiyalari. 5 (10): 1435–1447. Bibcode:1996ITIP .... 5.1435C. doi:10.1109/83.536892. PMID 18290061.
^ Ashburner, J. (2007 yil iyul). "Tasvirni tezkor diffeomorfik ro'yxatdan o'tkazish algoritmi". NeuroImage. 38 (1): 95–113. doi:10.1016 / j.neuroimage.2007.07.007. PMID 17761438.
^ Avants, B. B .; Epshteyn, K. L .; Grossman, M .; Gee, J. C. (2008-02-01). "O'zaro bog'liqlik bilan simmetrik diffeomorfik tasvirni ro'yxatdan o'tkazish: keksa va neyrodejenerativ miyaning avtomatlashtirilgan yorlig'ini baholash". Tibbiy tasvirni tahlil qilish. 12 (1): 26–41. doi:10.1016 / j.media.2007.06.004. ISSN 1361-8423. PMC 2276735. PMID 17659998.
^ Miller, Maykl; Banerji, Ayananshu; Kristensen, Gari; Joshi, Sarang; Xaneya, Navin; Grenander, Ulf; Mateych, Larissa (1997-06-01). "Hisoblash anatomiyasidagi statistik usullar". Tibbiy tadqiqotlarda statistik usullar. 6 (3): 267–299. doi:10.1177/096228029700600305. PMID 9339500.
^ U. Grenander va M. I. Miller (2007-02-08). Naqsh nazariyasi: vakillikdan xulosaga qadar. Oksford universiteti matbuoti. ISBN 9780199297061.
^ M. I. Miller va S. Mori va X. Tang va D. Tvard va Y. Jang (2015-02-14). Bayesian bir nechta Atlasning o'zgaruvchan shablonlari. Miya xaritasi: Entsiklopedik ma'lumotnoma. Akademik matbuot. ISBN 9780123973160.
^ Srivastava, S .; Miller, M. I .; Grenander, U. (1997-01-01). Byorns, Kristofer I.; Datta, Bisva N.; Martin, Klayd F.; Gilliam, Devid S. (tahrir). ATR uchun maxsus evklid guruhlari bo'yicha ergodik algoritmlar. Tizimlar va boshqarish: asoslar va ilovalar. Birkxauzer Boston. 327-350 betlar. CiteSeerX 10.1.1.44.4751. doi:10.1007/978-1-4612-4120-1_18. ISBN 978-1-4612-8662-2.
^ ^a ^b P. Dyupuy, U. Grenander, M.I. Miller, Diffeomorfizm oqimlari bo'yicha echimlarning mavjudligi, Amaliy matematikaning chorakligi, 1997 y.
^ ^a ^b Trouvé, A. (1995). "Action de groupe de dime infinie et reconnaissance de formes". Comptes Rendus de l'Académie des Sciences, Série I (frantsuz tilida). 321 (8): 1031–1034.
^ Tang, Syaoyin; Oishi, Kenichi; Fariya, Andreiya V.; Xillis, Argye E.; Albert, Merilin S.; Mori, Susumu; Miller, Maykl I. (2013-06-18). "Ko'p atlasli tasodifiy orbitadagi modelda Bayesian parametrlarini baholash va segmentatsiya qilish". PLOS ONE. 8 (6): e65591. Bibcode:2013PLoSO ... 865591T. doi:10.1371 / journal.pone.0065591. PMC 3688886. PMID 23824159.

[1] Kristensen, GE; Rabbitt, R.D .; Miller, M.I. (1996-02-01). "Katta deformatsiyaning kinematikasidan foydalangan holda deformatsiyalanadigan shablonlar". Rasmni qayta ishlash bo'yicha IEEE operatsiyalari. 5 (10): 1435–1447. Bibcode:1996ITIP .... 5.1435C. doi:10.1109/83.536892. PMID 18290061.

[2] Ashburner, J. (2007 yil iyul). "Tasvirni tezkor diffeomorfik ro'yxatdan o'tkazish algoritmi". NeuroImage. 38 (1): 95–113. doi:10.1016 / j.neuroimage.2007.07.007. PMID 17761438.

[3] Avants, B. B .; Epshteyn, K. L .; Grossman, M .; Gee, J. C. (2008-02-01). "O'zaro bog'liqlik bilan simmetrik diffeomorfik tasvirni ro'yxatdan o'tkazish: keksa va neyrodejenerativ miyaning avtomatlashtirilgan yorlig'ini baholash". Tibbiy tasvirni tahlil qilish. 12 (1): 26–41. doi:10.1016 / j.media.2007.06.004. ISSN 1361-8423. PMC 2276735. PMID 17659998.

[:1-4] Miller, Maykl; Banerji, Ayananshu; Kristensen, Gari; Joshi, Sarang; Xaneya, Navin; Grenander, Ulf; Mateych, Larissa (1997-06-01). "Hisoblash anatomiyasidagi statistik usullar". Tibbiy tadqiqotlarda statistik usullar. 6 (3): 267–299. doi:10.1177/096228029700600305. PMID 9339500.

[:9-5] U. Grenander va M. I. Miller (2007-02-08). Naqsh nazariyasi: vakillikdan xulosaga qadar. Oksford universiteti matbuoti. ISBN 9780199297061.

[:19-6] M. I. Miller va S. Mori va X. Tang va D. Tvard va Y. Jang (2015-02-14). Bayesian bir nechta Atlasning o'zgaruvchan shablonlari. Miya xaritasi: Entsiklopedik ma'lumotnoma. Akademik matbuot. ISBN 9780123973160.

[7] Srivastava, S .; Miller, M. I .; Grenander, U. (1997-01-01). Byorns, Kristofer I.; Datta, Bisva N.; Martin, Klayd F.; Gilliam, Devid S. (tahrir). ATR uchun maxsus evklid guruhlari bo'yicha ergodik algoritmlar. Tizimlar va boshqarish: asoslar va ilovalar. Birkxauzer Boston. 327-350 betlar. CiteSeerX 10.1.1.44.4751. doi:10.1007/978-1-4612-4120-1_18. ISBN 978-1-4612-8662-2.

[:22-8] P. Dyupuy, U. Grenander, M.I. Miller, Diffeomorfizm oqimlari bo'yicha echimlarning mavjudligi, Amaliy matematikaning chorakligi, 1997 y.

[:4-9] Trouvé, A. (1995). "Action de groupe de dime infinie et reconnaissance de formes". Comptes Rendus de l'Académie des Sciences, Série I (frantsuz tilida). 321 (8): 1031–1034.

[10] Tang, Syaoyin; Oishi, Kenichi; Fariya, Andreiya V.; Xillis, Argye E.; Albert, Merilin S.; Mori, Susumu; Miller, Maykl I. (2013-06-18). "Ko'p atlasli tasodifiy orbitadagi modelda Bayesian parametrlarini baholash va segmentatsiya qilish". PLOS ONE. 8 (6): e65591. Bibcode:2013PLoSO ... 865591T. doi:10.1371 / journal.pone.0065591. PMC 3688886. PMID 23824159.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]