1 UVOD Slikovne tehnike so postale temeljno orodje v sodobni medicini. Njihova uporaba pomaga v zdravniški praksi na različne načine, od diagnosticiranja zlomov kosti, odkritja rakastih tvorb do slikovno vodenih kirurških posegov. V uporabi so številne slikovne tehnike, kjer vsaka od njih prikazuje slikani objekt na drugačen način. Čeprav različne slikovne modalnosti, dobljene z različnimi tehnikami slikanja, v osnovi prikazujejo različne informacije, med njimi ostajajo tudi nekatere podobnosti. Na primer, medicinske slike so bolj ali manj pošumljene svetlostne slike na črnem ozadju, ki predstavlja okoliški zrak. Torej je na tako rekoč vseh medicinskih slikah opazen zunanji obris med slikanim objektom in okolico. Z uvedbo slikovnih postopkov [1], kot so: magnetna resonanca (MR), ultrazvok (US), računalniška tomografija (CT) in druge, je omogočen boljši vpogled v zgradbo anatomskih struktur. Pomemben korak postopkov razčlenjevanja medicinskih slik je definicija področja, na podlagi katere sliko razčlenjujemo. Razčlenjevanje slik lahko temelji na analizi prostora značilk slikovnih elementov, analizi homogenosti povezanih slikovnih elementov, iskanju robov med različnimi območji v prostoru slike ali na analizi fizikalnih lastnosti površin, ki jih predstavlja slika [2]. Magnetnoresonančno (MR) slikanje [3, 4, 5] spada med anatomske slikovne tehnike in je v nasprotju s CT za človeško telo neinvazivna slikovna tehnika. Nastane kot odziv vodikovih atomov v telesu na zunanje magnetno vzbujanje. Slikano telo se mora nahajati v močnem stalnem magnetnem polju, tako da gostota magnetnega polja doseže od 0.04 do 4 Tesla. Obstaja več različnih postopkov magnetnoresonančnega zajemanja slik, ki dajo različno informacijo o slikanem tkivu in tako predstavljajo različne modalnosti. Uveljavljene modalnosti MR slikanja so T1, T2 in PD. MR slikanje odlikuje dober kontrast med različnimi tipi tkiv in je zato primerno za podroben prikaz anatomije. Uporablja se predvsem za diagnosticiranje tumorjev in drugih nepravilnosti možganskega tkiva. Za Prejet 31. maj, 2012 Odobren 4. junij, 2012 130 FINKŠT pravilno vrednotenje rezultatov razčlenjevanja MR-slik je treba imeti tako ročno obrisana območja kot tudi območja, ki jih obrišejo različni postopki. Glede na primerjavo subjektivne presoje zdravnikov specialistov in samodejnih postopkov razčlenjevanja lahko vrednotimo samodejne postopke razčlenjevanja MR slik. Spoznavanje področja uporabe vzorcev objektov temelji na samodejnem razvrščanju (rojenju, samoorganizaciji) vzorcev iz dane množice in njihovem označevanju glede na pripadnost rojem (angl. clustering) [6]. S tem dano množico vzorcev preslikamo v učno množico vzorcev, ki jo pregleda in prilagodi še kompetentna oseba – učitelj oziroma strokovnjak za določeno področje uporabe. V našem primeru je namen omenjenih postopkov ugotavljanje rojenja vzorcev. S temi postopki preslikamo oziroma razbijemo dano množico vzorcev v podmnožice, kjer preslikava temelji na merjenju podobnosti vzorcev. Pri tem pričakujemo, da so v posameznem roju med seboj čim bolj podobni vzorci. Podobnost med dvema vzorcema merimo s preslikavo, ki primerjana vzorca preslika v realno število. Podobnost med vzorci, ki so predstavljeni s krajevnimi vektorji značilk, ponavadi merimo z razdaljami, med katerimi se najpogosteje uporablja Evklidova razdalja. Pri teh razdaljah manjše realno število pomeni večjo podobnost med primerjanima vzorcema. V obeh postopkih, tako k-tih povprečij (angl. k-means (kM)) in mehkih c-tih povprečij (angl. fuzzy c-means (FcM)) je treba na začetku nastaviti središče rojev [7]. Vhodni vektorji (podatkovne točke) so nato dodeljeni v enega od obstoječih rojev glede na kvadrat evklidske razdalje od roja, ki je najbližji. Povprečje (sredino) vsakega roja nato izračunamo tako, da dopolnimo središče roja. Ta posodobitev je posledica spremembe v pripadnosti za vsak roj. Združevanje vhodnih vektorjev in posodobitev središč rojev ponavljamo, dokler ni več sprememb vrednosti katerih koli središč rojev. V zadnjem času se postopek c-tih povprečij uporablja v nenadzorovanih postopkih rojenja na podlagi vidnih rezultatov pri samodejnem razčlenjevanju medicinskih slik [7,8,9,10]. Rojenje z mehkim postopkom c-tih povprečij [11,12,13] uspešno uporabljamo na številnih področjih, kot so astronomija, geologija, medicinsko slikanje in razčlenjevanje slik. Omenjena mehka razčlenjevalna metoda ima precejšnje prednosti, ker lahko ohrani veliko več informacij iz prvotne slike kot klasične razčlenjevalne metode [11]. Dva preprosta pristopa za začetno nastavitev središča roja sta naključna izbira začetne vrednosti ali izbira prvih k-vzorcev iz podatkovnih točk. Različne množice začetnih vrednosti so izbrane (zunaj podatkovnih točk) in izbrana je množica, ki je najbližja optimalni. Vendar pa je testiranje različnih začetnih vrednosti množic praktično neizvedljivo zaradi velikega števila rojev [9]. V literaturi [7,14,15,16] so bile predlagane različne metode. Rojenje je učinkovito orodje pri razčlenjevanju medicinskih slik in za nadaljnje načrtovanje zdravljenja. V MR slikah možganov je možgansko tkivo zapletena struktura in s tem je pravilna diagnoza številnih bolezni možganov v veliki meri odvisna od natančnega razčlenjevanja treh možganskih tkiv, in sicer belega tkiva (WM), sivega tkiva (GM) in cerebrospinalne tekočine (CSF). 2 RAZČLENJEVANJE SLIK Proces razčlenjevanja je razdelitev slike na neprekrivajoča se območja, tako da vsako območje ustreza kriteriju homogenosti [17], ki jo definiramo takole: naj bo ( , )f x y množica vrednosti značilk v različnih točkah koordinatnega sistema slike in ()P predikat homogenosti, definiran nad skupino povezanih slikovnih elementov. Razčlenjevanje je razdelitev slike ( , )f x y na n povezanih podmnožic oziroma območij 1 2 ( , , , ), n S S S tako da velja: 1 ( , ) ; 0, . n i i j i S f x y S S i j      (1) Predikat homogenosti ima za vsako območje i S vrednost ( ) 1 i P S  , medtem ko je vrednost predikata za povezani sosednji območji ( ) 0 i j P S S  . Ta definicija je splošna za vse tipe slik. Tako predikat homogenosti pomeni pri svetlostnih sivinskih slikah kriterij homogenosti svetlostne intenzitete, pri globinskih slikah pa kriterij homogenosti globinskega reliefa površin. Pri vektorskih slikah pa predikat definira homogenost barve ali teksture. V literaturi zasledimo veliko postopkov razčlenjevanja slik, vendar noben od njih ni optimalen za razčlenjevanje vseh tipov slik, pa tudi vse metode ne dajejo enakih rezultatov razčlenjevanja za isti tip slike. Podroben pregled metod razčlenjevanja sivinskih slik je podan v literaturi [17]. Postopki, razviti za določen tip slik, na primer svetlostne sivinske slike, niso vedno uporabni za magnetnoresonančne slike. Zato si v takšnem primeru pomagamo pri razčlenjevanju slike z modelom nastanka slike. Izbira primerne metode razčlenjevanja je zahtevna, ker ni splošne metode vrednotenja razčlenjene slike. Pavlidis [18] meni, da je razčlenjevanje vizualnih slik problem psihofizičnega zaznavanja. Vsak matematični postopek je ponavadi dopolnjen s hevristiko, ki vključuje znanje o vhodni sliki. Na splošno obstajata dva pristopa razčlenjevanja slik [17]: klasični in mehki matematični pristop. Opazno je pomanjkanje metod za objektivno vrednotenje rezultatov razčlenjevanja. Objektivna primerjava metod razčlenjevanja je bila podana za PRIMERJAVA RAZČLENJEVANJA MAGNETNORESONANČNIH SLIK Z UPORABO POSTOPKOV K-TIH IN MEHKIH …. 131 globinske slike [19]. V tem primeru je primerjalna referenca za razčlenjevanje globinskih slik dobro definirana. Bistvo razčlenjevalnih postopkov je združevanje delov slike z dovolj sorodnimi lastnostmi v en oziroma v isti objekt [6]. Obstaja množica razčlenjevalnih postopkov, katerih uporabnost je od primera do primera različna. Ni pa univerzalnega razčlenjevalnega postopka. V našem primeru smo uporabili za razčlenjevanje slik postopek rojenja [6]. Najbolj znana postopka rojenja sta postopek k-tih povprečij in mehki postopek c-tih povprečij. Obstaja še vrsta drugih razčlenjevalnih metod, ki pa že poznajo in upoštevajo določene lastnosti anatomije razčlenjenih slik. Postopke rojenja [6] na splošno razdelimo po načelu: objektivni, grafično-teoretični, hierarhični, ali po vrsti vzorcev: deterministični, statistični, mehki. V nadaljevanju bosta uporabljeni neskončni družini objektivnih funkcij postopkov rojenja, imenovani metoda rojenja k-tih povprečij in mehka metoda rojenja c-tih povprečij. Vse to kaže velik napredek tako teorije kot uporabe teh dveh postopkov v zadnjih letih [12]. 2.1 Postopek k-tih povprečij Postopek k-tih povprečij, ki ga uvrščamo med parametrične postopke rojenja, temelji na minimizaciji kriterijske funkcije J . V vsaki iteraciji izračuna srednjo intenziteto za vsak razred in potem razčleni sliko tako, da vsak slikovni element dodeli v razred, ki ima intenziteti slikovnega elementa najbližjo srednjo vrednost intenzitete. Ena najpogosteje uporabljanih kriterijskih funkcij je vsota kvadratov razdalj med vzorci v roju (enačba (2)) in aritmetična sredina roja (enačba (3)): 2 1 1 cN N i j j i J     x c (2) 1 1 N j ij N   c x , (3) kjer je x n-dimenzionalni vzorec, c N število rojev, j S j-ti roj vzorcev, j c aritmetična sredina vzorcev j-tega roja, j N pa je število vzorcev v roju j S . Postopek k-tih povprečij lahko strnemo v naslednje korake: (1) V prvi iteraciji izberemo , c N N N je število vseh vzorcev, središč rojev 1 2 (1), (1), , (1). cN c c c (2) V k-ti ponovitvi postopka razdelimo vzorce v c N rojev. Vzorec x uvrstimo v roj j, če: ( ) ( ) 1,2, , . j i c x m k x m k i N i j        (3) Izračunamo nova središča rojev ( 1), 1,2, , . j c k j N c (4) Če je ( 1) ( ) j j k k c c za 1,2, , , c j N je postopek končan. Uporaba tega postopka na splošno zahteva eksperimentiranje z različnimi začetnimi središči rojev in če vnaprejšnje število rojev ni znano, tudi z različnim številom rojev. Začetna središča rojev smo poiskali s pomočjo hierarhičnega postopka rojenja: najnižji nivo v hierarhiji rojenja je sistem, kjer je vsak vzorec sam zase v svojem roju, najvišji nivo pa rojenje vseh vzorcev v en sam roj [9]. Na vsakem nivoju rojenja, začenši z najnižjim, združimo dva najbolj podobna roja v en sam večji roj. Združena roja se na poznejših – višjih nivojih ne razdružita več. Postopek združevanja smo prekinili pri nivoju s tremi roji (tri vrste možganskih tkiv) in središča teh rojev smo uporabili kot začetna središča v postopku k-tih povprečij. 2.2 Mehki postopek c-tih povprečij Z mehkim rojenjem želimo podatke razdeliti na več mehkih rojev. Tako vsak vzorec pripada posameznemu mehkemu roju z določeno stopnjo pripadnosti, ki se izračuna glede na razdaljo vzorca do središč posameznih mehkih rojev, ki jih postopek določa v vsakem koraku sproti. Metoda mehkega rojenja c-tih povprečij temelji na minimizaciji kriterijske funkcije, ki jo podaja enačba (4) 2 1 1 ( ) cN N m MR ij i j j i J      x c (4) V enačbi (4) je m izbrano realno število, pri čemer mora veljati 1 .m  N je število vzorcev, c N je število rojev, ij  je vrednost pripadnostne funkcije j-tega roja za i-ti vzorec . i x j c označuje središče roja j, i j x c pa razdaljo, s katero je določena podobnost med središčem roja j c in vzorcem . i x V praksi najpogosteje uporabljamo Evklidsko razdaljo, ki jo podaja enačba (5). Tx x x  (5) Mehko rojenje se izvaja iterativno, pri čemer se v vsakem koraku sproti izračunavajo vrednosti pripadnostnih funkcij ij  , kot kaže enačba (6), in središča rojev j c , kot kaže enačba (7). 132 FINKŠT   2 1 1 1 , c ij mN i j k i k x c x c              (6) pri čemer za vsak {1, }, Ni velja 1 1 cN ijj        1 1 . N m ij i i j N m ik i x za i         c (7) V enačbah (6) in (7) je i indeks posameznega vzorca, j in k pa označujeta indeks posameznega mehkega roja oziroma njegovega središča. Vrednost parametra m določa mehkost oz. ostrost porazdelitve pripadnostnih funkcij v prostoru. Za skrajni primer 1m  velja, da so pripadnostne funkcije, ki so definirane nad prostorom, degenerirane v ostre stopnje pripadnosti. Pripadnostne funkcije lahko torej v vsaki točki prostora zavzamejo samo eno od vrednosti  0,1ij  . Iz enačbe (6) vidimo, da je 1 ij   , če je norma i j x c za i-ti vzorec in središče j-tega roja najmanjša glede na vsa središča rojev. Za preostale roje pa velja 0 ij   , kjer je {1, , } c k N . Za drugo skrajnost, m  , pa velja, da so pripadnostne funkcije degenerirano zmehčane, kar pomeni, da so vrednosti pripadnostnih funkcij po vsem prostoru enake: 1/ ij c N  , in sicer za vsak {1, , } c j N . V praksi sta najpogosteje uporabljeni vrednosti parametra 1,25m  ali 2m  [20]. Postopek mehkega rojenja c-tih povprečij je opisan v nadaljevanju. (1) Izberemo število rojev c N in parameter m ter določimo začetno matriko pripadnosti [(0) ] ij   . (2) V k-ti iteraciji določimo središča rojev j c za 1, , c j N glede na ( )k . (3) Določimo novo matriko pripadnosti ( 1)k  . (4) Če je ( 1) ( )k k     , končamo postopek, sicer se vrnemo na 2. korak. 3 REZULTATI Ker razčlenjevalni postopek začetne točke izbere naključno, tudi končni rezultati razčlenjevanja niso bili vedno povsem enaki, bili pa so si zelo podobni. Neponovljiva je bila izbira posameznih segmentov za posamezna tkiva. Na slikah možganov, ki so posnete z različnimi modalnostmi in v različnih okoliščinah, smo iskali tri območja različnih tkiv. Ta tkiva so belo možgansko tkivo (WM), sivo možgansko tkivo (GM) in možganska tekočina (CSF). Po primerjavi dobljenih rezultatov s sliko, ki je predstavljala normo razčlenjevanja (ang. ground truth), smo se odločili, katera razčlenjena površina na sliki pripada kateremu tkivu. Slika 1 prikazuje MR sliko in referenčno razčlenjeno sliko. V raziskave smo vključili 30 MRI slik možganov modalnosti T1 v axialni projekciji [21]. 3.1 Ocenjevanje delovanja razčlenjevalnih postopkov Ocenjevanje učinkovitosti in pravilnosti delovanja razčlenjevalnih postopkov je pomembno, saj lahko s pomočjo ocen izberemo najučinkovitejši postopek. Za primerjavo posameznih postopkov pa moramo dobiti posamezne postopke v posameznih primerih delovanja. Najbolje je, da so te ocene preproste skalarne vrednosti, ki so primerljive med seboj. Torej rešujemo problem transformacije ujemanja razčlenjene slike z objektom v realnosti v preprosto številko. Povedano z drugimi besedami: za primerjavo rezultatov potrebujemo definirano metriko [22], ki pa jo je zaradi zapletene mnogodimenzionalne narave razčlenjevalnih podatkov težko definirati. Poleg problema o metriki pa smo postavljeni še pred en problem - pogosto ni idealno razčlenjene slike, ki bi jo lahko uporabili kot referenco za popolnoma pravilno razčlenjeno sliko. Približki idealne razčlenjene slike, ki so ponavadi uporabljeni, so v večini primerov narejeni kot ročno razčlenjena slika. Sliko je ponavadi ročno razčlenil strokovnjak za območje, ki ga prikazuje slika. Slika 1: Vhodna slika za in referenčno razčlenjena slika PRIMERJAVA RAZČLENJEVANJA MAGNETNORESONANČNIH SLIK Z UPORABO POSTOPKOV K-TIH IN MEHKIH …. 133 To pomeni, da je za nastanek razčlenjevalnih norm potrebnega veliko napornega in monotonega dela. Zaradi lastnosti tega dela pa tudi pri ročno razčlenjenih slikah nastajajo napake. Ena od metod, kako te napake odpraviti, je uporaba več strokovnjakov [22], ki razčlenijo isto sliko, in uporaba statističnih metod, s katerimi dobljene rezultate združimo v boljšo sliko. Ena od metod za dobivanje normalnih razčlenitev za določeno sliko pa ne uporabi naravne slike, ampak umetno generirano [22]. V tem primeru generiramo umetno sliko, ki vsebuje objekte, ki jih bomo potem iskali. Iz tako generirane slike simuliramo nastanek slike, ki jo običajno dobimo iz diagnostičnih naprav (rentgen, MR, CT, ultrazvok,...), ker bolj ali manj dobro poznamo karakteristike prenosne poti original – naprava za zajem slik – slika. Tako lahko bolj ali manj natančno posnemamo popačenja, ki jih na sliko vneseta naprava in metoda zajema. Rezultat te metode je par slik – prva je vhodni podatek v naš postopek, in druga slika, ki je idealno razčlenjena, se uporabi za preverjanje pravilnosti delovanja postopka. Pri ocenjevanju smo primerjali rezultate, ki smo jih dobili z obema postopkoma, z rezultati, pridobljenimi iz norme razčlenjevanja. Ocenjevanje in merjenje učinkovitosti je bilo zelo poenostavljeno, ker smo imeli na voljo idealne poravnane slike. Tako smo morali preprosto primerjati ujemanje deleža slikovnih elementov na sliki, ki je bil rezultat našega razčlenjevanja, s slikovnimi elementi norme razčlenjevanja. Kot merilo za oceno je bila vzeta preslikava p, ki jo prikazuje enačba (8). ( , )z x y je število slikovnih elementov razčlenjene slike z omenjenima postopkoma, ( , ) n z x y pa pomeni število slikovnih elementov norme razčlenjevanja. ( , ) ( , ) ( , ) ( , ) n n z x y z x y p z x y z x y    (8) 3.2 Rezultati učinkovitosti obeh postopkov Sklep, do katerega lahko pridemo pri pregledu rezultatov ujemanja razčlenjevanja, podaja tabela 1. Omenimo, da smo imeli nabor 20 slik. Mehki postopek c-tih povprečij je učinkovitejši od postopka k-tih povprečij. Pri postopku c-tih povprečij smo si izbrali vrednost parametra m=2. Pri tej vrednosti parametra m smo dobili najboljše rezultate. Prav tako ugotavljamo, da smo dobili za belo možgansko tkivo pri razčlenjevanju najboljše ujemanje, najmanj pa se je ujemala možganska tekočina (SCF) glede na referenčno razčlenjeno sliko. Tabela 1: Učinkovitost razčlenjevalnih postopkov k-tih povprečij in mehkih c-tih povprečij Tkivo/Postopek Belo možgansko tkivo Sivo možgansko tkivo Možganska tekočina (CSF) k-tih povprečij 0,8446 0,6644 0,7045 mehkih c-tih povprečij 0,9402 0,8295 0,7815 a) b) c) d) e) f) Slika 2: Rezultat razčlenjevanja magnetnoresonančne slike modalnosti T1 s postopkom k-tih povprečij a) belo tkivo, c) sivo tkivo in e) možganska tekočina in s postopkom mehkih c- tih povprečij b) belo tkivo, d) sivo tkivo in f) možganska tekočina 134 FINKŠT 4 RAZPRAVA Prispevek opisuje izvedbo postopkov in primerjavo za razčlenjevanje magnetnoresonančnih slik (MR) možganov. Naš postopek razčlenjevanja je zasnovan tako, da po začetni nastavitvi rojev opravlja celotno razčlenjevanje popolnoma avtomatsko, brez interakcije z uporabnikom. Popolna avtomatizacija je v smislu razbremenitve rutinskega dela uporabnika-zdravnika zaželena, vendar pa hkrati tudi pomeni, da s tem popolnoma izgubimo nadzor nad potekom razčlenjevanja. Tako napake, ki jo naredimo pri razčlenjevanju v določenem koraku, ne moremo več odstraniti, ampak se skozi vse nadaljnje faze razpoznavanja le še povečuje. Zato tudi predlagamo, da medicinske razčlenjevalne postopke raje zasnujemo kot polavtomatske. S to zahtevo vnesemo veliko večji nadzor in povečamo pravilnost razčlenjevanja. Tako omogočimo, da strokovnjak v kritičnih trenutkih poseže v razčlenjevanje in dobljene rezultate zavrne ali pa jih popravi. Ta interakcija je še posebej pomembna, kadar imamo po razčlenjevanju združene objekte. Problem reševanja združenih objektov je izjemno zahteven. Če pa omogočimo poseg strokovnjakov, lahko le-ti razveljavijo razčlenjevanje tega dela slike ter zahtevajo ponovitev ali pa na tem delu slike opravijo razčlenjevanje ročno. V omenjenih postopkih omenimo še problem, da po končanem razčlenjevanju iste slike nismo dobili vsakič istih rezultatov. Ista vrsta tkiva je namreč v zaporednih zagonih metode razvrščena v poljubni razred, kar je odvisno od začetne izbire točk. 5 SKLEP Razčlenjevalne metode lahko ob predpostavki, da so uporabljene pravilno, bistveno olajšajo diagnosticiranje in druge medicinske raziskave. Med raziskavami smo spoznali, da se metode za razčlenjevanje medicinskih slik ne razlikujejo bistveno od metod, namenjenih za razčlenjevanje drugih slik. Kljub tej podobnosti med metodami pa obstajajo med njimi tudi bistvene razlike. Razčlenjevalne metode za medicinske slike morajo imeti zelo nizek odstotek napačno razpoznanih objektov (bolan ali deformiran organ ne sme biti razpoznan kot zdrav), saj je lahko vsaka takšna napaka usodna. Takšna razčlenjevalna metoda pa mora biti v vseh fazah načrtovanja plod dela strokovnjakov s področja digitalne obdelave slik in s področja medicine.