Tomografie

Tomografie ( jině řecky τομή  - řez a γράφω - píšu) - získání vrstveného obrazu vnitřní struktury předmětu.

Typy tomografie

Anatomická tomografie

Anatomická neboli destruktivní tomografie (biotomie) je založena na fyzickém provedení zkoumaných úseků organismu s jejich následnou fixací pomocí chemikálií. Klasickými příklady anatomické tomografie jsou Pirogovovy řezy a obrazy histologických preparátů . Pro zachování tvaru karoserie při řezech se korpus fixuje např. zmrazením .

Rekonstrukční tomografie

Rekonstrukční nebo nedestruktivní tomografie - získávání tím či oním způsobem informace o rozložení zájmového parametru v objektu vyšší dimenze podél jeho projekcí nižší dimenze bez zničení objektu; antonymum anatomické tomografie. Rozsah konceptu zahrnuje analogovou rekonstrukční tomografii a počítačovou (počítačovou) tomografii .

Analogová rekonstrukční tomografie je rekonstrukční tomografie, která nepoužívá digitální, ale analogová výpočetní zařízení (například optická) k obnovení distribuce parametru objektu.

Metoda byla navržena pro rentgenové vyšetření francouzským lékařem Bocaillem a implementována jako přístroj (nazývaný „tomograf“) italským inženýrem Vallebona (a přibližně ve stejnou dobu i inženýry z jiných zemí) ve dvacátých letech a na počátku 1930, a byl založen na pohybu dvou ze tří složek radiografie ( rentgenová trubice , rentgenový film , předmět studia ). Tomograf umožnil získat jeden snímek – snímek vrstvy ležící ve zvolené hloubce studovaného objektu. Nejpoužívanější způsob střelby, kdy zkoumaný objekt zůstal nehybný a rentgenka a filmová kazeta se koordinovaně pohybovaly v opačných směrech. Při synchronním pohybu tubusu a kazety je na filmu zřetelná pouze potřebná vrstva, protože pouze její příspěvek k celkovému stínu zůstává vůči filmu nehybný, vše ostatní je rozmazané, téměř bez zásahu do rozboru výsledného obraz. Metoda se nazývá klasická nebo lineární tomografie . V současné době se podíl posledně jmenované metody ve výzkumu ve světě snižuje kvůli relativně nízkému informačnímu obsahu a vysoké radiační zátěži.

Panoramatická tomografie je široce používána v medicíně pro diagnostiku onemocnění dentoalveolárního systému . Díky pohybu zářiče a kazety s rentgenovým filmem po speciálních trajektoriích je vybrán obraz ve formě válcové plochy. To umožňuje získat snímek zobrazující všechny zuby pacienta.

Počítačová tomografie je obor matematiky , který vyvíjí matematické metody a algoritmy pro rekonstrukci vnitřní struktury objektu z projekčních dat – digitálních snímků objektu pořízených opakovaným prosvětlováním tohoto objektu v různých protínajících se směrech. Vnitřní struktura je obvykle reprezentována ve formě voxelu . Získání řady voxelů z řady promítacích obrazů se nazývá přímý tomografický problém . Do oboru výpočetní tomografie patří i řešení inverzního tomografického problému  - vytvoření libovolného projekčního pohledu na základě známé vnitřní struktury.

Počítačová tomografie je teoretickým základem počítačové tomografie, metody pro získávání vrstvených obrazů předmětu ve třech rovinách s možností jejich trojrozměrné rekonstrukce. Nejčastěji se počítačová tomografie týká rentgenové počítačové tomografie (CT).

Na rozdíl od rentgenového CT využívá magnetická rezonance (MRI) nízkoenergetické elektromagnetické vlny a při častém používání nepředstavuje pro pacienta nebezpečí. MRI a CT mají rozdíly a používají se v různých případech, nejsou zaměnitelné [1] .

Historie tomografie

Klasifikace typů tomografie

Vzájemná poloha zdroje sondovacího záření, objektu a detektoru

Z hlediska vzájemné polohy zdroje sondovacího záření, objektu a detektoru lze tomografické metody rozdělit do následujících skupin:

Rozměry studovaných objektů

Rozsah

Podle rozsahu aplikace existují:

Sondovací záření

Tomografické algoritmy

Je známo několik tisíc algoritmů, které se používají pro problémy výpočetní (rekonstrukční) tomografie. Lze je kombinovat do několika velkých hlavních skupin.

Od dob Abela, Radona, Weinsteina se používají analytické algoritmy inverzní transformace. Matematickým rysem těchto problémů je, že patří do třídy problémů položených špatně podle Hadamarda , které zpravidla souvisejí s Fredholmovými integrálními rovnicemi. Účinným prostředkem k jejich řešení s konečným počtem projekcí je metoda regularizace akademika A. N. Tichonova , kterou následně vyvinuli Phillips, Arsenin, Yaglom, Tanana a mnoho dalších.

Pro osově symetrické systémy se přímo používá inverzní Abelova transformace. Jeho diskrétní verzi poprvé použil Van Cittert na problém řešení za Rayleighovým limitem.

Pro dvourozměrné systémy popsané dvěma oddělitelnými proměnnými se používá elementární transformace Agravala a Sodha. Pro systémy se známou grupou symetrie udává Weinsteinova věta nejmenší počet průmětů, které postačují k přesné rekonstrukci systému.

Od 40. let 20. století (Tikhonov et al.) lze tomografické problémy pro 2- a 3-rozměrné objekty řešit numerickými metodami. Numerický diskrétní model soustavy integrálních rovnic je v konečném důsledku redukován zpravidla na speciální (nedostatečně nebo naopak přeurčenou a nekonzistentní) soustavu lineárních rovnic velkých rozměrů, navíc s rozměry od 3- a 4. - (pro dvourozměrnou tomografii) na 5- a 6-rozměrnou (pro trojrozměrnou tomografii). Čtyřrozměrná tomografie je známá v experimentální jaderné fyzice a fyzice svazku nabitých částic (Sandia National Laboratories, Brookhaven National Laboratory, CERN , M. V. Keldysh Research Center, Moskevský institut fyziky a technologie atd.).

Řešení takových soustav klasickými "exaktními" metodami (Gauss-Jordan apod.) je tedy nereálné vzhledem ke kubickému počtu prvků objektu =N M , kde N je charakteristická lineární velikost objektu, M je dimenze, velké výpočetní náklady (což dokazuje Klyuev-Kokovkin-Shcherbak teorém ). Například pro dvourozměrné úlohy v řádu 100×100 bude potřeba asi 1 bilion operací s nahromaděním zaokrouhlovacích chyb a pro trojrozměrné úlohy 100×100×100 asi 10 18 operací, což odpovídá doba cca 1 hodiny počítání na superpočítačích s výkonem desítek petaflopů.

Třída 1 je tedy výpočetně neuspokojivá. K jejich řešení se používají tři další třídy algoritmů:

První technické a biologické výpočetní introskopy-tomografy v SSSR (1940-50. léta) a první lékařské výpočetní tomografy v USA (70. léta) skutečně využívaly řadu verzí metody polského matematika Kaczmarze (1937), včetně sovětské matematik I. A. Bochek (1953, Moskevský fyzikální a technologický institut). Nositel Nobelovy ceny Cormack a Hounsfield tedy nazvali Kaczmarzův algoritmus, který použili (zajišťující dosažení bodu nejmenších čtverců) ART (1973); algoritmus sovětského matematika Tarasca (zajišťující dosažení bodu maximální pravděpodobnosti, 60. léta, IPPE, Obninsk) nazvali MART; použili také algoritmus japonského matematika Kuino Tanabe (1972), který je relaxační a superrelaxační verzí Kaczmarzova algoritmu. Často se používá Friedenův algoritmus (zajišťující dosažení maximálního bodu entropie). Stochastické metody pro výčet rovnic v projekcích (první z nich byla stochastická verze algoritmu I. A. Bochka, publikovaná v roce 1971) umožňují vyhnout se pravidelným artefaktům a výrazně zlepšit kvalitu obrazu.

Pokud je pro skenovací schémata s „tenkými paprsky“ systém rovnic relativně dobře podmíněn (výsledek rekonstrukce tedy není příliš citlivý na nevyhnutelné chyby v projekčním měření), pak pro skenování „tlustými paprsky“ (což je typické pro problematika NMR tomografie, ultrazvuk, PET, mikrovlnná introskopie Oshchepkov, elektrická proudová tomografie se systém rovnic ukazuje jako velmi špatně podmíněný, což vede k prudkému zpomalení přístupu iterací výše uvedených projekčních metod k řešení. K řešení takových systémů se používají metody A. V. Gorškova (MIPT) a S. Elsakova (SUSU), které se liší necitlivostí na špatnou podmíněnost soustav řešených rovnic a také nutným stochastickým výčtem rovnic. v nich absence pravidelných artefaktů a konečně míra konvergence (v praktických problémech) je o 2–3 řády vyšší, než bylo uvedeno dříve.

Pro nelineární rovnice a tomografii velkorozměrných objektů (trojrozměrné v medicíně, vědě a technice, 4-, 5-, 6-rozměrné v jaderné fyzice a fyzice plazmatu a svazku částic, v technologii urychlovačů), varianty Monte Carlo jsou efektivní metodou řešení v metrických prostorech velkých rozměrů.

Algoritmus sovětského a ruského matematika A. A. Abramova současná komprimace iterací na řešení a iterací na ortogonalizaci poskytuje záruku stabilní konvergence k řešení a zároveň velmi přesný odhad chyby a rychlosti rekonstrukce. Upozorňujeme, že u špatně podmíněných systémů jsou jako jeho elementární iterace doporučovány iterace druhého řádu (Gorshkov, Elsakova atd.), nikoli iterace prvního řádu (Kachmarz-Bocek, Tarasco, Frieden atd.) , nebo dokonce (je-li to nutné, v praktických problémech se zatím nesetkáme) iteraci 3. nebo vyšších řádů.

Všimněte si, že bychom neměli zbytečně používat iterace příliš vysokých řádů, protože výpočetní náklady pro ně s neomezeným nárůstem řádu iterací mají tendenci být kubické (v N**M) (jako v přímé Gauss-Jordanově inverzi).

K řešení výpočetních problémů soufázové ultrazvukové, mikrovlnné, SBMM a elektropotenciální tomografie je použit algoritmus akademika Lavrentieva.

Poznámky

  1. Davydov, D. 8 mýtů o zobrazování magnetickou rezonancí  / D. Davydov, A. Krjuchkov // Tinkoff Journal. - 2022. - 12. dubna.
  2. Jurij Erin. Vytvořena čtyřrozměrná elektronová tomografie . Živly (2. srpna 2010). Získáno 3. srpna 2010. Archivováno z originálu dne 24. srpna 2011.
  3. Kwon, O.-H. 4D elektronová tomografie ]  / O.-H. Kwon, A. H. Zewail // Věda . - 2010. - Sv. 328, č.p. 5986. - S. 1668-1673. - doi : 10.1126/science.1190470 .
  4. Tomografie. // Malá lékařská encyklopedie .. - M .  : Lékařská encyklopedie, 1991–96.
  5. Tomografie. // První pomoc .. - M .  : Velká ruská encyklopedie, 1994.
  6. Lékařská tomografie. // Encyklopedický slovník lékařských termínů .. - M .  : Sovětská encyklopedie, 1982–1984.

Literatura

Odkazy