Pozitronová emisní tomografie je zobrazovací metoda nukleární medicíny. Tomografický obraz je rekonstruován díky současné detekci dvou fotonů, které byly vyzářeny při anihilaci pozitronu vzniklého při rozpadu radiofarmaka v pacientově těle. Názorný popis principu PET najdete na webu Principy PET
Úkolem následujícího cvičení je demonstrovat tvorbu signálu v pozitronové emisní tomografii. Potřebná data a kódy si stáhněte zde.
phantom(p_def);
s použitím definičních matic pro konstrukci fantomu, které jsou uložené v datovém souboru ph_def.mat. Dalši parametry funkce budou následující (ke každé veličině jsou v závorkách uvedeny hodnoty, které použijete, podrobnější vysvětlení je v sekciTvorba signálu a jeho měření). [3b.]
PETreconstruction(detPair,r,Nd);
(nachází se také v úvodním archivu). Spočtete a vizualizujte relativní chybu mezi rekonstruovaným fantomem z vašeho měření a originálním fantomem. [1b.]
Funkce které jste si stáhli na začátku cvičení neodevzdávejte.
Naprogramujte funkci, která vypočítá výstup měřícího prstence PET. Vstupem funkce nechť je čas podání radiofarmaka <latex>T_i</latex> [s], čas začátku měření <latex>T_b</latex> [s], čas konce měření <latex>T_e</latex> [s], počet detektorů Nd, poloměr detektorové kružnice r[px], rozptyl pozitronů <latex>\sigma^2</latex> [px], množství n[mol], poločas rozpadu radiofarmaka <latex>\tau_r</latex> [s] a fantom aktivit (relativní distribuce radiofarmaka v těle pacienta, pro absolutní roložení celý fantom znormujeme, aby součet byl 1 a vynásobíme celkovám počtem rozpadlých částic <latex>\Delta N</latex>). Výstupem pak bude matice současně aktivovaných detektorů P která má na souřadnici [i,j] počet současných aktivací (detekcí) elementů číslo i a j a celkový počet měřených rozpadů <latex>\Delta N</latex> . Získaná data si můžete rekonstruovat pomocí funkce activity=PETreconstruction(P,r,Nd) ze staženého balíčku.
Předpokládejme, že v čase <latex>T=0</latex> bylo vyrobeno radiofarmakum, které pak bylo v čase <latex>T_i</latex> podáno pacientovi.
Počet vyrobených molekul je <latex> N_0 = n \cdot N_A</latex>,kde <latex>N_A</latex> je Avogadrova konstanta <latex>N_A=6.0221415\cdot 10^{23} \textrm{mol}^{-1}</latex>. Pro určení počtu molekul dostupných v čase <latex>T_i</latex> musíme určit, jaká část radioizotopu se do této doby již rozpadla.
Uvažujme, že za čas <latex>dt</latex> se v látce rozpadne <latex>dN</latex> částic <latex> dN = -\lambda_r \cdot N dt \label{difeq}</latex>, kde <latex>\lambda_r</latex> je rozpadová konstanta <latex>\lambda_r=\frac{\ln{2}}{\tau_r}</latex>. Řešením rovnice je známá exponenciální závislost počtu izotopů na čase <latex> N_1=N_0\textrm{e}^{-\lambda_r t} </latex>.
Před začátkem měření je potřeba vyčkat, než se radioaktivní látka dostane do místa určení (často je podávána systémově, například intravenózně), proto zahájíme měření až v <latex>T_b</latex>. Naše úloha bude simulovat měření pomocí fludeoxyglukózy (FDG), která je analogem glukózy s inkorporovaným atomem fluóru 18. Protože je aktivní látka od času <latex>T_i</latex> v lidském těle, začne se, kromě rozpadu, také vylučovat. Přibližně 75% FDG je zachyceno ve tkáni a rozpadá se s <latex>\tau_r=110</latex> min, 25% látky se vyloučí ledvinami s <latex>\tau_v=16</latex> min. To je samozřejmě možné modelovat, ale vzhledem ke krátkému <latex>\tau_v</latex> a známému poměru zachyceného a vyloučeného radiofarmaka si můžeme práci zjednodušit (při zachování dostatečné přesnosti) a uvažovat, že je pro rozpad k dispozici jen 75% molekul radiofarmaka.
Každý z rozpadajících se atomů vyzáří pozitron, který se pohybuje náhodným směrem tkání a poté co ztratí svou kinetickou energii, anihiluje s elektronem hmoty. Při tom vzniknou dva fotony s energií přibližně 511 keV, vzdalující se od místa srážky na opačné strany po náhodně orientované přímce. Pokud dvojice fotonů poletí podél osy pacientova těla, nemá pro náš zobrazovací systém význam. Z geometrického uspořádání je jasné, že našeho myšleného prstencového detektoru dosáhne jen malá část záření.
Uvažujme pacienta vysokého 180 cm, z něhož pořizujeme přibližně 5 mm řez. Pokud by bylo radiofarmakum zachycováno rovnoměrně podél osy pacientova těla, pak v našem řezu zůstane přibližně 5/1800 částic. Toto zjednodušeni nám výrazně usnadní práci, při zachování řádové přesnosti.
Náš detektor je tenký prstenec kolmý k ose pacienta a z toho je jasně patrné, že nemůže zachytit fotony, jejichž dráha není téměř kolmá k pacientovu tělu. Takové dvojice fotonů minou detektor a jsou z našeho pohledu nepodstatné, protože nepřispívají k tvorbě obrazu. Definujme že rozdílový úhel (mezi rovinou detektoru a dráhou fotonů) může být maximálně <latex>1^\circ</latex> , detektor tedy zasáhnou jen 2/180 fotonů vzniklých v řezu.
Z první části dostáváme celkový počet detekcí v jedné vrstvě. Pokračovat budeme konfigurací detektoru a výpočtem (simulací) detekovaných rozpadů.
fspecial
a aplikujete ho na data pomocí imfilter
),r= \sqrt{ m^2 + n^2 } \cdot \frac{1}{2} </latex> kde m,n jsou velikosti obrázku/fantomu).
K určení průsečíků lze použít funkci lineAndCircle2intersectionAngles(a, b, c, rad, mid_point);
z uvodního balíčku. První tři parametry jsou vektory koeficientů jednotlivých přímek (v obecném tvaru <latex>a \cdot x + b \cdot y + c = 0</latex>). rad
je poloměr detektorového kruhu, mp
souřadnice středu fantomu. Funkce vrací průsečíky (ang1
, ang2
) s detektorovým kruhem jako úhel od počátku. Výše zmíněná funkce také může vykreslit jednotlivé čáry, pokud jí navíc dáme jako parametr i vstupní matici fantomu. Možný výstup pro bod se čtyřmi rozpady je na následujícím obrázku