Томография

Томогра́фия (др.-греч. τομή — сечение) — метод неразрушающего послойного исследования внутренней структуры объекта посредством его многократного просвечивания в различных пересекающихся направлениях.

Терминологические вопросы

См. также Интроскопия, Линейная томография, Томосинтез.

В классической трактовке под томографией понимается метод рентгенологического исследования, с помощью которого можно производить снимок слоя, лежащего на определённой глубине исследуемого объекта. Он был предложен Бокажем через несколько лет после открытия рентгеновских лучей и был основан на перемещении двух из трёх компонентов рентгенографии (рентгеновская трубка, рентгеновская плёнка, объект исследования). Наибольшее распространение получил метод съёмки, при котором исследуемый объект оставался неподвижным, а рентгеновская трубка и кассета с плёнкой согласованно перемещались в противоположных направлениях. При синхронном движении трубки и кассеты, четким на пленке получается только необходимый слой, потому что только его вклад в общую тень остаётся неподвижным относительно плёнки, всё остальное смазывается, почти не мешая проводить анализ полученного изображения. В настоящее время доля последнего метода в исследованиях в мире уменьшается в связи со своей относительно малой информативностью. В России в связи с дороговизной и недостаточной укомплектованностью медицинских учреждений современным диагностическим оборудованием и высокой заболеваемостью туберкулезом данный метод остается широко распространенным и актуальным. В настоящее время данный метод получил название классическая томография или линейная томография.

Вычислительная томография — область математики, занимающаяся разработкой математических методов и алгоритмов восстановления внутренней структуры объекта по проекционным данным (цифровым снимкам объекта, сделанным с разных точек). Внутренняя структура как правило представляется в воксельной форме. Получение массива вокселей по массиву проекционных снимков называется прямой томографической задачей. К области вычислительной томографии так же относится и решение обратной томографической задачи - формирование произвольного проекционного вида на основании известной внутренней структуры.

См. также преобразование Радона и экспоненциальное преобразование Радона.

Также вычислительной томографией может называться практическая область деятельности, занимающаяся томографией с использованием этих численных методов.

Компьютерная томография (КТ) — то же, что вычислительная томография. Чаще всего под КТ подразумевается томография рентгеновская.

Анатомическая томография — основана на получении срезов тканей человека с их последующей фиксацией с помощью химических веществ и регистрация их на фотоплёнку. Классическими примерами анатомической томографии являются пироговские срезы и изображения гистологических препаратов.

Реконструктивная томография — получение тем или иным способом распределения интересующего параметра в объекте большей размерности по его проекциям меньшей размерности без разрушения объекта; антоним анатомической томографии. В объём понятия входят вычислительная (она же компьютерная, она же цифровая) и аналоговая реконструктивная томографии.

Аналоговая реконструктивная томография — реконструктивная Т., использующая для восстановления распределения параметра объекта не цифровые, а аналоговые вычислительные устройства (например, оптические).

Неразрушающая томография — то же, что и реконструктивная томография. Антоним анатомической томографии.

История томографии

В этом разделе не хватает ссылок на источники информации. Информация должна быть проверяема, иначе она может быть поставлена под сомнение и удалена. Вы можете отредактировать эту статью, добавив ссылки на авторитетные источники. Эта отметка установлена 12 мая 2011.

• До ХХ в. Математики Фредгольм и Абель исследуют свойства семейства интегральных уравнений, позже ставших основой томографии.
• В 1895 г. В. К. Рентген открывает проникающие «Х-лучи», позже названные благодарными физиками и врачами его именем — «рентгеновские».
• В 1905 г. Я. Ван-Циттерт осуществил томографическое измерение распределение яркости далёкой звезды по радиусу как численное обратное преобразование Абеля.
• В 1917 году австрийский математик Иоганн Радон предложил способ обращения интегрального преобразования, впоследствии получившего его имя (см. преобразование Радона), благодаря которому стало возможно восстанавливать изначальную функцию, зная её преобразование. Однако в то время работа Радона не попала в поле зрения исследователей и вскоре была незаслуженно забыта современниками.
• В 20-х гг. ХХ в. французский врач Бокаж изобрёл и запатентовал томографический механический сканер, который должен был оставлять на рентгенограмме неразмытым только заданный слой тела пациента. Это называлось «рентгеновская планиграфия», а также «биотомия», а позже было названо «классическая томография».
• В 1930 г. А. Валлебона, как сообщают^[кто?], «разработал принцип послойного рентгенологического исследования (томографии)».
• В 1934 г. В. И. Феоктистов создал первый действующий рентгеновский томограф.
• В 1937 г. польский математик Качмаж опубликовал алгоритм, который впоследствии был использован Кормаком и Хаунсфилдом без ссылки на автора.
• В 1937 г. И. Раби открыл новое явление — ядерный магнитный резонанс (ЯМР) в изолированном ядре.
• В 1938 г. А. Ощепков изобрёл СВЧ-интроскопию.
• В 1941 г. А. Н. Тихонов изобрёл метод регуляризации, сделавший возможным реконструкцию при неточных проекциях.
• 1941—1945 гг. Рентгеновские, гамма- и нейтронные интроскопы с вычислительной томографической обработкой по Тихонову осуществлены в СССР для задач дефектоскопии в авиационной и пушечной отраслях промышленности, а к концу войны — и для контроля объёмного распределения процессов в ядерных реакторах.
• В 1944 г. Е. К. Завойский открыл новое явление — электронный парамагнитный резонанс (ЭПР).
• В 1946 г. Ф. Блох и Э. Парселл повторили открытие И. Раби в конденсированных средах.
• В 1953 г. И. А. Бочек изобрёл стохастическую версию алгоритма Качмажа, избавившую реконструкции от регулярных артефактов и значительно увеличивший качество изображений.
• В 1953 г. советский математик Вайнштейн доказал теорему о связи минимального достаточного количества проекций с группой симметрии объекта, резко упростившую томографию.
• В 1960 г. В. А. Иванов изобрёл ЯМР-томографию (внутривидение на основе ядерного магнитного резонанса).
• В 1963 году американский физик А. Кормак повторно (но отличным от Радона способом) решил задачу томографического восстановления, а в 1969 году английский инженер-физик Г. Хаунсфилд из фирмы EMI Ltd. сконструировал «ЭМИ-сканер» (EMI-scanner) — первый компьютерный рентгеновский томограф, чьи клинические испытания прошли в 1972 году. В 1979 году Кормак и Хаунсфилд «за разработку компьютерной томографии» были удостоены Нобелевской премии по физиологии и медицине.
• А в 2003 за изобретение метода магнитно-резонансной томографии, на основе открытия Реймонда Дамадьяна, Нобелевскую премию по физиологии и медицине получили Питер Мэнсфилд и Пол Лотербур.
• В 1991 г. в одном из крупнейших медицинских учреждений России НЦПЗ РАМН был установлен 1-ый отечественный МР-томограф "Образ-1", произведенный Научно Производственной Фирмой "Аз" организованной в 1988 году как инновационное предприятие в области разработки и производства отечественных магнитно-резонансных томографов, до сих пор успешно эксплуатируется.
• В 2010 году создана так называемая четырёхмерная электронная томография — техника визуализирования динамики трёхмерных объектов во времени. Эта техника позволяет наблюдать за пространственно-временными характеристиками микрообъектов^[1][2].

Классификация видов томографии

Взаиморасположение источника зондирующего излучения, объекта и детектора

С точки зрения взаиморасположения источника зондирующего излучения, объекта и детектора томографические методы могут быть разделены на следующие группы:

• трансмиссионные — регистрируется зондирующее внешнее излучение, прошедшее через пассивный (неизлучающий) объект, частично ослабляясь при этом (тень от объекта);
• эмиссионные — регистрируется излучение, выходящее из активного (излучающего) объекта с некоторым пространственным распределением источников излучения;
• комбинированные трансмиссионно-эмиссионные (люминесцентные, акустооптические и оптоакустические и др.) — регистрируется вторичное излучение от источников, распределённых по объёму объекта и возбуждённых внешним излучением;
• эхозондирование — регистрируется зондирующее внешнее излучение, отражённое от внутренних структур пассивного объекта.

Размеры исследуемых объектов

• Микроуровень (микротомография) — исследуются объекты размером с отдельную клетку.
• Уровень объектов, соизмеримых с человеческим телом (от отдельного органа или лабораторной мыши до самолёта).
• Макроуровень — атмосферные явления (облака, циклоны, торнадо), планеты и звёзды.

Сфера применения

По сфере применения выделяют:

• медицинскую томографию (как вид медицинской визуализации и медицинской диагностики);
• промышленную (техническую) томографию (как вид дефектоскопии);
• томографию макрообъектов.

Зондирующее излучение

• Томография с использованием звуковых волн (в том числе сейсмических):
  • ультразвуковая томография (УЗТ);
  • сейсмическая томография.
• Томография с использованием электромагнитного излучения:
  • радионуклидная эмиссионная томография (гамма-излучение);
    • однофотонная эмиссионная томография (ОФЭКТ);
    • двухфотонная эмиссионная или позитронно-эмиссионная томография (ПЭТ);
  • рентгеновская томография;
    • рентгеновская компьютерная томография (КТ, РКТ);
  • оптическая (лазерная) томография (ОТ);
  • томография в радиодиапазоне.
• Томография с использованием электромагнитных полей:
  • магнитно-резонансная томография (МРТ);
  • электро-импедансная томография.
• Томография с использованием элементарных частиц:
  • нейтронная томография;
  • электронная и позитронная томография;
  • протонная томография;
  • нейтринная томография.

Томографические алгоритмы

Этот раздел статьи следует викифицировать. Пожалуйста, оформите его согласно правилам оформления статей.

В этом разделе не хватает ссылок на источники информации. Информация должна быть проверяема, иначе она может быть поставлена под сомнение и удалена. Вы можете отредактировать эту статью, добавив ссылки на авторитетные источники. Эта отметка установлена 12 мая 2011.

Известны несколько тысяч алгоритмов, применяемых для задач вычислительной (реконструктивной) томографии. Их можно объединить в несколько больших основных групп.

Со времён Абеля, Радона, Вайнштейна применялись алгоритмы аналитического обратного преобразования. Математической особенностью этих задач является то, что они принадлежат классу некорректно по Адамару поставленных задач, как правило, родственных интегральным уравнениям Фредгольма. Эффективным средством их решения при конечном числе проекций является метод регуляризации академика А. Н. Тихонова, развитый впоследствии Филлипсом, Арсениным, Ягломом, Тананой и многими др.

Для осесимметричных систем применяют непосредственно обратное преобразование Абеля. Его дискретная версия впервые была применена Ван-Циттертом для задачи разрешения сверх предела Рэлея.

Для 2-мерных систем, описываемых 2-мя разделяющимися переменными, применяют элементарное преобразование Агравала и Содха. Для систем с известной группой симметрии теорема Вайнштейна указывает наименьшее число проекций, достаточных для точной реконструкции системы.

С 40-х годов (Тихонов и др.) томографические задачи для 2- и 3-мерных объектов поддаются решению численными методами. Численная дискретная модель системы интегральных уравнений сводится, в конечном итоге, как правило, к особенной (недоопределённой либо, напротив, переопределённой и несовместной) системе линейных уравнений большого размера, причём с размерностью от 3-х и 4-х (для 2-мерной томографии) до 5- и 6-мерной (для 3-мерной томографии). В экспериментальной ядерной физике и физике пучков заряженных частиц известна 4-мерная томография (Sandia Nat.Lab., Broockhaiwen Nat.Lab., CERN, Исследовательский центр им. М. В. Келдыша, МФТИ и др.).

Таким образом, решение таких систем классическими «точными» методами (Гаусса-Жордана и т. п.) нереально вследствие кубически по числу элементов объекта =N^M, где N — характерный линейный размер объекта, M — размерность, больших вычислительных затрат (что доказано теоремой Клюева — Коковкина — Щербака). Например, для 2-мерных задач порядка 100х100 потребуется порядка 1 трлн операций с накоплением погрешностей округления, а для 3-мерных 100х100х100 — порядка 10¹⁸ операций, что соответствует времени порядка 1 часа счёта на рекордных современных в мире многопетафлопных супер-ЭВМ. Итак, класс 1 вычислительно неудовлетворителен.

Для их решения применяют 3 иных класса алгоритмов.

Класс 2. Безытерационное обратное преобразование разложения проекций по ортогональным функциям (Фурье, Чебышёва, Котельникова, Хартли, Уолша, Радемахера и др.).

Класс 3. Регуляризация по Тихонову (или без неё до заранее оцененного предела сверхразрешения Косарева) в сочетании с итерационными методами многомерного поиска — спуска, Монте-Карло и др.

Класс 4. Регуляризация по Тихонову (или без неё до заранее оцененного предела сверхразрешения Косарева) в сочетании с итерационными проекционными алгоритмами. Все проекционные алгоритмы базируются на теореме математика Банаха (г. Львов) о сжимающих отображениях. Важным их достоинством является гарантированная и устойчивая сходимость итераций. Ещё более важным их достоинством для многомерной томографии является радикально более низкая вычислительная трудоёмкость — квадратичная по N**M. Для вышеуказанных параметров объекта в 2-мерном случае это пропорционально 100 млн операций и числу итераций, то есть порядка 1 часа счёта на рядовой современной ПЭВМ (для итераций первого порядка) и порядка 1 сек. для итераций второго порядка. В 3-мерном случае (100х100х100) это пропорционально 1 трлн операций и числу итераций, то есть порядка 1 сек. (если первого порядка) или порядка 1 миллисекунды (если второго порядка) на супер-ЭВМ.

Первые технические и биологические вычислительные интроскопы-томографы в СССР (40-е — 50-е гг.) и первые медицинские вычислительные томографы в США (70-е гг.) фактически использовали ряд версий метода польского математика Качмажа (1937 г.), в том числе советского математика И. А. Бочека (1953 г., МФТИ). Так, награждённые Нобелевской премией Кормак и Хаунсфилд использованный ими алгоритм Качмажа (обеспечивающий достижение точки наименьших квадратов) называли ART (1973 г.), алгоритм советского математика Тараско (обеспечивающий достижение точки максимума правдоподобия, 60-е гг., ФЭИ, г. Обнинск) они назвали MART, также они использовали алгоритм японского математика Куино Танабе (1972 г.), являющийся релаксационной и сверхрелаксационной версией алгоритма Качмажа. Часто используется алгоритм Фридена (обеспечивающий достижение точки максимума энтропии). Стохастические методы перебора уравнений в проекциях (первым из таких была стохастическая версия алгоритма И. А. Бочека, опубликованная в 1971 г.) позволяют избежать регулярных артефактов и значительно улучшить качество изображения.

Если для схем сканирования «тонкими лучами» система уравнений сравнительно хорошо обусловлена (следовательно, результат реконструкции мало чувствителен к неизбежным погрешностям измерений проекций), то для сканирования «толстыми лучами» (что характерно для задач ЯМР-томографии, УЗИ, ПЭТ, СВЧ-интроскопии Ощепкова, электротоковой томографии, система уравнений оказывается очень плохо обусловленной. Это приводит к резкому замедлению приближения итераций вышеупомянутых проекционных методов к решению. Для решения таких систем используют методы А. В. Горшкова (МФТИ) и С. Елсакова (ЮУрГУ), отличающиеся нечувствительностью к плохой обусловленности решаемых систем уравнений, а также, за счёт необходимого стохастического перебора уравнений в них, отсутствием регулярных артефактов, и, наконец, скоростью сходимости (в практических задачах) на 2-3 порядка большей, чем указанные ранее.

Для нелинейных уравнений и томографии объектов большой размерности (3-мерной в медицине, науке и технике, 4-, 5-, 6-мерной в ядерной физике и физике плазмы и пучков заряженных частиц, в ускорительной технике) эффективным методом решения являются варианты метода Монте-Карло в метрических пространствах большой размерности.

Алгоритм советского и российского математика А. А. Абрамова (МФТИ) одновременных сжимающих итерации к решению и итерации к ортогонализации обеспечивает гарантию устойчивой сходимости к решению и заодно весьма точную оценку погрешности и скорости реконструкции. Укажем, что в плохо обусловленных системах в качестве его элементарных итераций рекомендуются не итерации первого порядка (Качмажа-Бочека, Тараско, Фридена и т. п.), а второго порядка (Горшкова, Елсакова и др.), или даже (в случае необходимости, пока не встреченной в практических задачах) итерации 3-го или большего порядков.

Заметим, что не следует без необходимости использовать итерации слишком высоких порядков, так как вычислительные затраты на них при неограниченном увеличении порядка итерации стремятся к кубическим (по N**M) (как у прямого обращения Гаусса-Жордана).

Для решения вычислительных задач синфазных УЗ-, СВЧ-, СБММ- и электропотенциальной томографии используют алгоритм академика Лаврентьева.

Примечания

1. ↑ Юрий Ерин Создана четырехмерная электронная томография  (рус.). Элементы.ру (2 августа 2010). Архивировано из первоисточника 24 августа 2011. Проверено 3 августа 2010.
2. ↑ Oh-Hoon Kwon, Ahmed H. Zewail 4D Electron Tomography (англ.) // Science. — 2010. — Т. 328. — № 5986. — С. 1668—1673.

Литература

В Викисловаре есть статья «томография»

• A. C. Kak, M. Slaney Principles of Computerized Tomographic imaging. (IEEE Press, NY 1988)
• Хорнак Дж. П. Основы МРТ (1996—1999)
• Cormack A.M. Early two-dimensional reconstruction and recent topics stemming from it // Nobel Lectures in Physiology or Medicine 1971—1980. — World Scientific Publishing Co., 1992. — p. 551—563
• Hounsfield G.N. Computed Medical Imaging // Nobel Lectures in Physiology or Medicine 1971—1980. — World Scientific Publishing Co., 1992. — p. 568—586
• Lauterbur P.C. All science is interdisciplinary — from magnetic moments to molecules to men // Les Prix Nobel. The Nobel Prizes 2003. — Nobel Foundation, 2004. — p. 245—251
• Mansfield P. Snap-shot MRI // Les Prix Nobel. The Nobel Prizes 2003. — Nobel Foundation, 2004. — p. 266—283
• Мэнсфилд П. Быстрая магнитно-резонансная томография // Успехи физических наук, 2005, т. 175, № 10, с. 1044—1052 (перевод на русский)
• Дьячкова С. Я., Николаевский В. А. Рентгеноконтрастные средства. — Воронеж, 2006.
• Важенин А. В., Ваганов Н. В. Медицинско-физическое обеспечение лучевой терапии. — Челябинск, 2007.
• Левин Г. Г., Вишняков Г. Н. Оптическая томография. — М.: Радио и связь, 1989. — 224 с.
• Тихонов А. Н., Арсенин В. Я., Тимонов А. А. Математические задачи компьютерной томографии. — М.: Наука, Гл. ред. физ.-мат. лит., 1987. — 160 с.
• Тихонов А. Н., Гончарский А. В., Степанов В. В., Ягола А. Г. Численные методы решения некорректных задач. — М.: Наука, Гл. ред. физ.-мат. лит., 1990. — 232 с.
• Наттерер Ф. Математические аспекты компьютерной томографии. — М.: Мир, 1990. — 288 с.
• Васильев М. Н., Горшков А. В. Аппаратно-программный комплекс GEMMA и томографический метод измерения многомерных функций распределения в траекторном и фазовом пространствах при диагностике пучков заряженных частиц. // Приборы и техника эксперимента. — 1994. № 5. — С.79-94. // Перевод на англ.: Instruments and Experimental Techniques. — V.37. № 5. Part 1. 1994. -P.581-591.
• Горшков А. В. Пакет программ REIMAGE для существенного улучшения разрешения изображений при обработке данных физического эксперимента и метод нахождения неизвестной аппаратной функции. 26.01.94. // Приборы и техника эксперимента. — 1995. № 2. — С.68-78. // Перевод на англ.: Instruments and Experimental Techniques. — V.38. № 2. 1995. — P.185-191.
• Москалёв И. Н., Стефановский А. М. Диагностика плазмы с помощью открытых цилиндрических резонаторов. — М.: Энергоатомиздат, 1985.
• Хермен Г. Восстановление изображений по проекциям: Основы реконструктивной томографии. — М.: Мир, 1983. — 352 с.
• Горелова Л. Е. Нобелевская премия. Нобель, Мечников, Рентген. // Сайт ММА им. И. М. Сеченова.

Ссылки

• Группа Вычислительной томографии, Институт теоретической и прикладной механики Сибирского отделения РАН