Фазово-контрастная рентгеновская визуализация или фазово-чувствительная рентгеновская визуализация — это общий термин для различных технических методов, которые используют информацию об изменениях фазы рентгеновского луча , проходящего через объект , для создания его изображений. Стандартные методы рентгеновской визуализации, такие как рентгенография или компьютерная томография (КТ), основаны на уменьшении интенсивности ( затухания ) рентгеновского луча при прохождении образца , которое можно измерить непосредственно с помощью детектора рентгеновского излучения . Однако при фазово-контрастной рентгеновской визуализации фазовый сдвиг луча , вызванный образцом, не измеряется напрямую, а преобразуется в изменения интенсивности, которые затем могут быть зарегистрированы детектором. [1]
Помимо создания проекционных изображений , фазово-контрастная рентгеновская визуализация, как и обычная передача, может сочетаться с томографическими методами для получения трехмерного распределения реальной части показателя преломления образца. При применении к образцам, состоящим из атомов с низким атомным номером Z , фазово-контрастная рентгеновская визуализация более чувствительна к изменениям плотности в образце, чем обычная рентгеновская визуализация на основе пропускания . Это приводит к получению изображений с улучшенным контрастом мягких тканей . [2]
За последние несколько лет были разработаны различные методы фазово-контрастной рентгеновской визуализации, все из которых основаны на наблюдении интерференционных картин между дифрагированными и недифрагированными волнами. [3] Наиболее распространенными методами являются кристаллическая интерферометрия, визуализация на основе распространения, визуализация на основе анализатора, краевая подсветка и визуализация на основе решетки (см. ниже).
Первым, кто открыл рентгеновские лучи, был Вильгельм Конрад Рентген в 1895 году, когда он обнаружил, что они способны проникать через непрозрачные материалы. Он записал первый рентгеновский снимок, на котором запечатлена рука его жены. [4] Он был удостоен первой Нобелевской премии по физике в 1901 году «в знак признания выдающихся заслуг, которые он оказал, открыв замечательные лучи, впоследствии названные в его честь». [5] С тех пор рентгеновские лучи стали использовать в качестве инструмента для безопасного определения внутренней структуры различных объектов, хотя долгое время информация была получена путем измерения только передаваемой интенсивности волн, а информация о фазе не была получена. доступный.
Принцип фазово-контрастного изображения был впервые разработан Фрицем Цернике во время его работы с дифракционными решетками и видимым светом. [6] [7] Применение его знаний в микроскопии принесло ему Нобелевскую премию по физике в 1953 году. С тех пор фазово-контрастная микроскопия стала важной областью оптической микроскопии .
Перевод фазово-контрастной визуализации из видимого света в рентгеновские лучи занял много времени из-за медленного прогресса в улучшении качества рентгеновских лучей и недоступности рентгеновских линз. В 1970-х годах стало понятно, что синхротронное излучение , испускаемое заряженными частицами, циркулирующими в накопительных кольцах, построенных для экспериментов по ядерной физике высоких энергий, могло быть более интенсивным и универсальным источником рентгеновских лучей, чем рентгеновские трубки ; [8] это, в сочетании с прогрессом в развитии рентгеновской оптики, имело основополагающее значение для дальнейшего развития рентгеновской физики.
Пионерская работа по внедрению метода фазового контраста в рентгеновскую физику была представлена в 1965 году Ульрихом Бонсом и Майклом Хартом с факультета материаловедения и инженерии Корнелльского университета, Нью-Йорк. Они представили кристаллический интерферометр , изготовленный из большого и совершенного монокристалла . [9] Не менее 30 лет спустя японские учёные Ацуси Момосе , Тохору Такеда и коллеги приняли эту идею и усовершенствовали её для применения в биологической визуализации, например, путём увеличения поля зрения с помощью новых конфигураций установки и фазового методы поиска . [10] [11] Интерферометр Бонса-Харта обеспечивает на несколько порядков более высокую чувствительность в биологических образцах, чем другие методы фазового контраста, но он не может использовать обычные рентгеновские трубки, поскольку кристаллы принимают только очень узкую энергетическую полосу рентгеновского излучения. лучи (ΔE / E ~ 10 −4 ). В 2012 году Хан Вэнь и его коллеги сделали шаг вперед, заменив кристаллы нанометровыми фазовыми решетками. [12] Решетки расщепляют и направляют рентгеновские лучи в широком спектре, тем самым снимая ограничения на полосу пропускания источника рентгеновского излучения. Они обнаружили субнанорадианное преломление рентгеновских лучей в биологических образцах с помощью решеточного интерферометра Бонса-Харта. [12]
В то же время появились два дополнительных подхода к фазово-контрастному изображению с целью преодоления проблем кристаллической интерферометрии. Метод визуализации на основе распространения был впервые представлен группой Анатолия Снигирева голографии в 1948 году. [14]
в ESRF (Европейском центре синхротронного излучения) в Гренобле, Франция, [13] и был основан на обнаружении «полос Френеля», возникающих при определенных условиях. обстоятельства распространения в свободном пространстве. Экспериментальная установка состояла из линейного источника рентгеновского излучения, образца и детектора и не требовала каких-либо оптических элементов. Концептуально это было идентично замыслу революционной работы Денниса Габора поАльтернативный подход, называемый визуализацией на основе анализатора, был впервые исследован в 1995 году Виктором Ингалом и Еленой Беляевской в рентгеновской лаборатории в Санкт-Петербурге, Россия, [15] и Тимом Дэвисом и его коллегами из CSIRO (Организация научных и промышленных исследований Содружества). ) Отделение материаловедения и технологий в Клейтоне, Австралия. [16] Этот метод использует кристалл Брэгга в качестве углового фильтра для отражения на детектор только небольшой части луча, удовлетворяющего условию Брэгга . Важный вклад в развитие этого метода был внесен американскими исследовательскими группами Дина Чепмена, Чжун Чжуна и Уильяма Томлинсона, например, выделением дополнительного сигнала, вызванного сверхмалым угловым рассеянием [17] и первым КТ-изображение, полученное с помощью анализатора. [18] Альтернатива визуализации на основе анализатора, которая обеспечивает эквивалентные результаты без использования кристалла, была разработана Алессандро Оливо и его коллегами на синхротроне Elettra в Триесте, Италия. [19] Этот метод, называемый «краевым освещением», осуществляет точный отбор в направлении рентгеновских лучей, используя физический край самих пикселей детектора, отсюда и название. Позже Оливо в сотрудничестве с Робертом Спеллером из Университетского колледжа Лондона адаптировал метод для использования с обычными источниками рентгеновского излучения [20] , открыв путь к его использованию в клинических и других приложениях. Питер Манро (также из UCL) внес существенный вклад в развитие лабораторного подхода, продемонстрировав, что он практически не предъявляет требований к согласованности [21] и что, несмотря на это, он по-прежнему является полностью количественным. [22]
Последний подход, обсуждаемый здесь, - это так называемая визуализация на основе решеток, в которой используется эффект Тальбота , открытый Генри Фоксом Тэлботом в 1836 году . [23] Этот эффект самоизображения создает интерференционную картину после дифракционной решетки . На определенном расстоянии эта картина в точности напоминает структуру решетки и регистрируется детектором. Положение интерференционной картины можно изменить, поместив в луч объект, вызывающий фазовый сдвиг. Это смещение интерференционной картины измеряется с помощью второй решетки и с помощью определенных методов восстановления получается информация о действительной части показателя преломления. Так называемый интерферометр Талбота-Лау первоначально использовался в атомной интерферометрии , например, Джоном Ф. Клаузером и Шифаном Ли в 1994 году. [24] Первые интерферометры с рентгеновской решеткой, использующие синхротронные источники, были разработаны Кристианом Дэвидом и его коллегами из Институт Пауля Шеррера (PSI) в Виллингене, Швейцария [25] и группа Ацуши Момосе из Токийского университета. [26] В 2005 году, независимо друг от друга, группа Дэвида и Момосе включила компьютерную томографию в решеточную интерферометрию, что можно рассматривать как следующую веху в развитии визуализации на основе решеток. [27] [28] В 2006 году еще одним большим достижением стал перенос метода на основе решеток на обычные лабораторные рентгеновские трубки Францем Пфайффером и его коллегами, [29] что значительно расширило потенциал этого метода для клинического использования. Примерно два года спустя группе Франца Пфайффера также удалось извлечь дополнительный сигнал из своих экспериментов; так называемый «сигнал темного поля» был вызван рассеянием из-за пористой микроструктуры образца и давал «дополнительную и иным образом недоступную структурную информацию об образце в масштабе длины микрометра и субмикрометра». [30] В то же время Хан Вэнь и его коллеги из Национального института здравоохранения США разработали значительно упрощенную технику получения изображения рассеяния («темного поля»). Они использовали единственную проекцию сетки и новый подход к выделению сигналов, названный «однократный анализ Фурье». [31] Недавно было проведено много исследований по улучшению метода, основанного на решетках: Хан Вэнь и его команда проанализировали кости животных и обнаружили, что интенсивность сигнала темного поля зависит от ориентации сетки, и это связано с анизотропии костной структуры. [32]Они добились значительного прогресса в биомедицинских приложениях, заменив механическое сканирование решеток электронным сканированием источника рентгеновского излучения. [33] Поле фазоконтрастной КТ на основе решетки было расширено за счет томографических изображений темнопольного сигнала [34] и фазоконтрастной КТ с временным разрешением. [35] Кроме того, были опубликованы первые доклинические исследования с использованием фазово-контрастной рентгеновской визуализации на основе решеток. Марко Стампанони и его группа исследовали нативную ткань молочной железы с помощью «дифференциальной фазово-контрастной маммографии» [36] , а группа под руководством Дэна Статмана исследовала, как использовать визуализацию на основе решеток для мелких суставов руки. [37]
Совсем недавно значительный прогресс в визуализации на основе решеток произошел благодаря открытию Веном и его коллегами эффекта фазового муара [38] [39] . Это привело к интерферометрии, выходящей за пределы диапазона самоизображения Тэлбота, с использованием только фазовых решеток и обычных источников и детекторов. Фазовые решетки рентгеновского излучения могут быть изготовлены с очень малыми периодами, что позволяет получать изображения при низких дозах радиации и достигать высокой чувствительности.
Обычная рентгеновская визуализация использует падение интенсивности за счет ослабления, вызванного объектом в рентгеновском луче, а излучение рассматривается как лучи, как в геометрической оптике . Но когда рентгеновские лучи проходят через объект, изменяются не только их амплитуда, но и фаза. Вместо простых лучей рентгеновские лучи можно рассматривать как электромагнитные волны . Тогда объект можно описать его комплексным показателем преломления (см. [8] ):
Член δ представляет собой декремент действительной части показателя преломления, а мнимая часть β описывает показатель поглощения или коэффициент экстинкции. Заметим, что в отличие от оптического света действительная часть показателя преломления меньше, но близка к единице, это «связано с тем, что спектр рентгеновских лучей обычно лежит в высокочастотной стороне различных резонансов, связанных с связывание электронов». [8] Фазовая скорость внутри объекта больше скорости света c . Это приводит к иному поведению рентгеновских лучей в среде по сравнению с видимым светом (например, углы преломления имеют отрицательные значения), но не противоречит закону относительности , «который требует, чтобы только сигналы, несущие информацию, не перемещались быстрее, чем с» . сигналы движутся с групповой скоростью , а не с фазовой скоростью, и можно показать, что групповая скорость на самом деле меньше с » . [8]
Влияние показателя преломления на поведение волны можно продемонстрировать на примере волны, распространяющейся в произвольной среде с фиксированным показателем преломления n . Для простоты здесь предполагается монохроматическая плоская волна без поляризации . Волна распространяется в направлении, нормальном к поверхности среды, обозначенном в этом примере z (см. рисунок справа). Скалярная волновая функция в вакууме равна
Внутри среды угловое волновое число изменяется от k до nk . Теперь волну можно описать так:
где δkz — фазовый сдвиг, а e −β kz — коэффициент экспоненциального затухания, уменьшающий амплитуду волны E 0 . [8]
В более общих терминах общий фазовый сдвиг луча, распространяющегося на расстояние z , можно рассчитать с помощью интеграла
где λ — длина волны падающего рентгеновского луча. Эта формула означает, что фазовый сдвиг представляет собой проекцию декремента действительной части показателя преломления в направлении изображения. Это выполняет требование томографического принципа , который гласит, что «входные данные для алгоритма реконструкции должны быть проекцией величины f , которая передает структурную информацию внутри образца. Затем можно получить томограмму, которая отображает значение f ». [40] Другими словами, при фазово-контрастной визуализации карта реальной части показателя преломления δ(x,y,z) может быть восстановлена с помощью стандартных методов, таких как обратная проекция с фильтром , которая аналогична обычной рентгеновской компьютерной томографии. где можно получить карту мнимой части показателя преломления.
Чтобы получить информацию о составе образца, а именно о распределении плотности образца, необходимо связать измеренные значения показателя преломления с собственными параметрами образца, такое соотношение задается следующими формулами:
где ρ a — атомная плотность, σ a — сечение поглощения , k — длина волнового вектора и
где p — сечение фазового сдвига.
Вдали от краев поглощения (пиков в сечении поглощения из-за повышенной вероятности поглощения фотона, имеющего частоту, близкую к резонансной частоте среды), эффектами дисперсии можно пренебречь; это касается легких элементов ( атомный номер Z <40), которые являются компонентами тканей человека, и рентгеновских лучей с энергией выше 20 кэВ, которые обычно используются в медицинской визуализации. В предположении этих условий сечение поглощения приближенно определяется выражением
где 0,02 — константа, выраженная в barn , типичная единица площади поперечного сечения взаимодействия частиц, k — длина волнового вектора , k 0 — длина волнового вектора с длиной волны 1 Ангстрем и Z — атомный номер . [41] Действительная формула в этих условиях для сечения фазового сдвига:
где Z — атомный номер , k — длина волнового вектора , а r 0 — классический радиус электрона .
В результате получаются следующие выражения для двух частей комплексного показателя преломления:
Подстановка типичных значений тканей человека в приведенные выше формулы показывает, что δ обычно на три порядка превышает β в пределах диагностического рентгеновского диапазона. Это означает, что фазовый сдвиг рентгеновского луча, распространяющегося через ткань, может быть намного больше, чем потеря интенсивности, что делает фазово-контрастную рентгеновскую визуализацию более чувствительной к изменениям плотности в ткани, чем визуализацию поглощения. [42]
Благодаря пропорциональности
преимущество фазового контраста перед обычным контрастом поглощения даже возрастает с увеличением энергии. Кроме того, поскольку формирование фазово-контрастного изображения не связано с поглощением рентгеновских лучей в образце, поглощенную дозу потенциально можно уменьшить за счет использования более высоких энергий рентгеновских лучей. [29] [42]
Как уже говорилось выше, применительно к видимому свету действительная часть показателя преломления n может сильно отклоняться от единицы (n стекла в видимом свете колеблется от 1,5 до 1,8), тогда как для рентгеновских лучей в различных средах отклонение от единицы обычно имеет порядка 10 -5 . Таким образом, углы преломления, возникающие на границе двух изотропных сред, рассчитанные по формуле Снелла, также очень малы. Следствием этого является то, что углы преломления рентгеновских лучей, проходящих через образец ткани, не могут быть обнаружены напрямую и обычно определяются косвенно путем «наблюдения интерференционной картины между дифрагированными и недифрагированными волнами, создаваемой пространственными изменениями действительной части показателя преломления». ." [3]
Кристаллическая интерферометрия , иногда также называемая рентгеновской интерферометрией , является старейшим, но и самым сложным методом, используемым для экспериментальной реализации. Он состоит из трех светоделителей геометрии Лауэ , расположенных параллельно друг другу. (См. рисунок справа) Падающий луч, который обычно предварительно коллимируется и фильтруется монохроматором (кристаллом Брэгга), разделяется на первом кристалле (S) дифракцией Лауэ на два когерентных луча: опорный луч, который остается невозмущенным, и луч, проходящий через образец. Второй кристалл (Т) действует как передающее зеркало и заставляет лучи сходиться друг к другу. Два луча встречаются в плоскости третьего кристалла (А), который иногда называют кристаллом-анализатором, и создают интерференционную картину, форма которой зависит от разности оптических путей между двумя лучами, вызванной образцом. Эта интерференционная картина обнаруживается детектором рентгеновского излучения, расположенным за кристаллом анализатора. [9] [43]
Поместив образец на вращающийся столик и записав проекции под разными углами, можно получить трехмерное распределение показателя преломления и, следовательно, томографические изображения образца. [40] В отличие от методов, описанных ниже, с помощью кристаллического интерферометра измеряется сама фаза, а не какое-либо ее пространственное изменение. Для извлечения фазового сдвига из интерференционных картин; используется метод, называемый фазовым шагом или сканированием полос: в опорный луч вводится фазовращатель (имеющий форму клина). Фазовращатель создает прямые интерференционные полосы с регулярными интервалами; так называемые полосы несущей. При помещении образца в другой пучок несущие полосы смещаются. Фазовый сдвиг, вызванный образцом, соответствует смещению полос несущих. Несколько интерференционных картин регистрируются для разных смещений опорного луча, и путем их анализа можно извлечь информацию о фазе по модулю 2 π . [40] [43] Эта неоднозначность фазы называется эффектом фазовой развертки и может быть устранена с помощью так называемых «методов фазовой развертки». [44] Эти методы можно использовать, когда отношение сигнал/шум изображения достаточно высокое и изменение фазы не слишком резкое. [28]
В качестве альтернативы методу сканирования полос можно использовать метод преобразования Фурье для извлечения информации о фазовом сдвиге только с помощью одной интерферограммы, тем самым сокращая время экспозиции, но это имеет недостаток, заключающийся в ограничении пространственного разрешения расстоянием между несущими. бахрома. [45]
Рентгеновская интерферометрия считается наиболее чувствительной к фазовому сдвигу из четырех методов и, следовательно, обеспечивает самое высокое разрешение плотности в диапазоне мг/см 3 . [28] Но из-за его высокой чувствительности полосы, создаваемые образцом с сильным фазовым сдвигом, могут стать неразрешимыми; Чтобы преодолеть эту проблему, недавно был разработан новый подход, названный «когерентно-контрастная рентгеновская визуализация», где вместо фазового сдвига для контраста изображения имеет значение изменение степени когерентности, вызванное образцом. [46]
Общее ограничение пространственного разрешения этого метода дается размытием в кристалле-анализаторе, возникающим из-за динамической рефракции, т.е. угловое отклонение луча из-за преломления в образце усиливается в кристалле примерно в десять тысяч раз, поскольку путь луча внутри кристалла сильно зависит от угла его падения. Этот эффект можно уменьшить за счет утончения кристалла анализатора, например, при толщине анализатора 40 мкм было рассчитано разрешение около 6 мкм . В качестве альтернативы кристаллы Лауэ можно заменить кристаллами Брэгга , чтобы луч не проходил через кристалл, а отражался от поверхности. [47]
Другим ограничением метода является требование очень высокой стабильности установки; выравнивание кристаллов должно быть очень точным, а разница в длине пути между лучами должна быть меньше длины волны рентгеновских лучей; Для этого интерферометр обычно изготавливают из цельного блока кремния высокой степени совершенства путем вырезания двух канавок. При монолитном производстве очень важная пространственная когерентность решетки между всеми тремя кристаллами может поддерживаться относительно хорошо, но это ограничивает поле зрения небольшими размерами (например, 5 см х 5 см для 6-дюймового слитка), а также потому, что образец обычно размещаемый на одном из путей луча, размер самого образца также ограничен размером кремниевого блока. [9] [48] Недавно разработанные конфигурации, использующие два кристалла вместо одного, значительно увеличивают поле зрения, но еще более чувствительны к механическим нестабильностям. [49] [50]
Другая дополнительная трудность кристаллического интерферометра заключается в том, что кристаллы Лауэ фильтруют большую часть входящего излучения, что требует высокой интенсивности луча или очень длительного времени экспозиции. [51] Это ограничивает использование метода очень яркими источниками рентгеновского излучения, такими как синхротроны.
В соответствии с ограничениями на установку кристаллический интерферометр лучше всего подходит для получения изображений с высоким разрешением небольших образцов, которые вызывают небольшие или плавные градиенты фазы .
Чтобы обеспечить превосходную чувствительность кристаллической интерферометрии Бонса-Харта без некоторых основных ограничений, монолитные кристаллы были заменены нанометрическими рентгеновскими фазовыми решетками. [52] Первые такие решетки имеют периоды от 200 до 400 нанометров. Они могут расщеплять рентгеновские лучи по широкому энергетическому спектру обычных рентгеновских трубок. Основное преимущество этого метода заключается в том, что он использует большую часть входящих рентгеновских лучей, которые были бы отфильтрованы кристаллами. Поскольку используются только фазовые решетки, изготовление решеток менее сложное, чем методы, в которых используются абсорбционные решетки. Первый решеточный интерферометр Бонса-Харта (gBH) работал при энергии фотонов 22,5 кэВ и спектральной полосе пропускания 1,5%.
Входящий луч имеет форму щелей размером в несколько десятков микрометров, так что длина поперечной когерентности превышает период решетки. Интерферометр состоит из трех параллельных и равноотстоящих друг от друга фазовых решеток и рентгеновской камеры. Падающий луч дифрагируется на первой решетке периода 2P на два луча. Далее они дифрагируются на второй решетке периода P на четыре луча. Два из четырех сливаются на третьей решетке периода 2P. Каждый из них далее дифрагирует на третьей решетке. Множественным дифрагированным лучам позволяют распространяться на достаточное расстояние, так что разные порядки дифракции разделяются в камере. Существует пара дифрагированных лучей, которые совместно распространяются от третьей решетки к камере. Они мешают друг другу, создавая полосы интенсивности, если решетки слегка не совмещены друг с другом. Центральная пара дифракционных путей всегда одинакова по длине независимо от энергии рентгеновских лучей или угла падающего луча. Интерференционные картины от фотонов различной энергии и углов падения синхронизированы по фазе.
Изображаемый объект размещается вблизи центральной решетки. Абсолютные фазовые изображения получаются, если объект пересекает один из пары когерентных путей. Если оба пути проходят через объект в двух местах, разделенных поперечным расстоянием d, то обнаруживается изображение разности фаз Φ(r) - Φ(rd). Для получения фазовых изображений выполняется фазовый шаг одной из решеток. Изображение разности фаз Φ(r) - Φ(rd) можно интегрировать для получения изображения фазового сдвига объекта.
Этот метод достиг существенно более высокой чувствительности, чем другие методы, за исключением кристаллического интерферометра. [12] [53] Основным ограничением метода является хроматическая дисперсия дифракционной решетки, которая ограничивает его пространственное разрешение. Настольная система с рентгеновской трубкой с вольфрамовой мишенью, работающей при напряжении 60 кВ, будет иметь предельное разрешение 60 мкм. [12] Еще одним ограничением является то, что ширина рентгеновского луча составляет всего лишь десятки микрометров. Потенциальное решение было предложено в виде параллельной визуализации с несколькими щелями. [12]
Визуализация на основе анализатора (ABI) также известна как визуализация с дифракционным усилением , интроскопия с фазовой дисперсией и рентгенография с несколькими изображениями [54]. Его установка состоит из монохроматора (обычно одиночного или двойного кристалла, который также коллимирует луч) перед образец и кристалл-анализатор расположены в геометрии Брэгга между образцом и детектором. (См. рисунок справа)
Этот кристалл-анализатор действует как угловой фильтр для излучения, исходящего от образца. При попадании этих рентгеновских лучей на кристалл-анализатор условие брэгговской дифракции выполняется лишь для очень узкого диапазона углов падения. Когда рассеянные или преломленные рентгеновские лучи имеют углы падения вне этого диапазона, они вообще не отражаются и не вносят вклад в сигнал. Преломленные рентгеновские лучи в этом диапазоне будут отражаться в зависимости от угла падения. Зависимость отраженной интенсивности от угла падения называется кривой качания и является внутренним свойством системы формирования изображений, т. е. она представляет собой интенсивность, измеренную в каждом пикселе детектора, когда кристалл анализатора «раскачивается» (слегка повернут на угол θ) без присутствия объекта, поэтому его можно легко измерить. [54] Типичный угловой аксептанс составляет от нескольких микрорадиан до десятков микрорадиан и связан с полной шириной на полувысоте (FWHM) кривой качания кристалла.
Когда анализатор идеально совмещен с монохроматором и, таким образом, расположен на пике кривой качания, получается стандартная рентгеновская фотография с повышенным контрастом, поскольку нет размытия рассеянными фотонами. Иногда это называют «контрастом угасания».
Если же в противном случае анализатор ориентирован под малым углом (углом расстройки) по отношению к монохроматору, то рентгеновские лучи, преломленные в образце на меньший угол, будут отражаться меньше, а рентгеновские лучи, преломленные на больший угол, будут отражаться. более. Таким образом, контраст изображения основан на различных углах преломления образца. Для малых градиентов фазы угол преломления можно выразить как
где k — длина волнового вектора падающего излучения, а второй член в правой части — первая производная фазы в направлении дифракции. Поскольку измеряется не сама фаза, а первая производная фазового фронта, визуализация на основе анализатора менее чувствительна к низким пространственным частотам, чем кристаллическая интерферометрия, но более чувствительна, чем PBI.
В отличие от предыдущих методов визуализация на основе анализатора обычно предоставляет информацию о фазе только в направлении дифракции, но не чувствительна к угловым отклонениям в плоскости, перпендикулярной плоскости дифракции. Эта чувствительность только к одному компоненту фазового градиента может привести к неоднозначности в оценке фазы. [55]
Путем записи нескольких изображений под разными углами расстройки, то есть в разных положениях на кривой качания, получается набор данных, который позволяет получить количественную информацию о дифференциальной фазе. Существует несколько алгоритмов восстановления информации по кривым качания, некоторые из них предоставляют дополнительный сигнал. Этот сигнал возникает в результате сверхмалого угла рассеяния на субпиксельных структурах образца и вызывает угловое уширение луча и, следовательно, уширение формы кривой качания. На основе этого контраста рассеяния можно создать новый вид изображения, называемый изображением темного поля. [17] [54] [56]
Томографическую визуализацию с помощью визуализации на основе анализатора можно выполнить, зафиксировав анализатор под определенным углом и повернув образец на 360 ° во время сбора проекционных данных. Из одного и того же образца получают несколько наборов проекций с разными углами расстройки, после чего можно восстановить томографическое изображение. Если предположить, что кристаллы обычно выровнены так, что производная показателя преломления измеряется в направлении, параллельном томографической оси, полученное «изображение рефракционной КТ» показывает чистое изображение градиента вне плоскости.
Для визуализации на основе анализатора требования к стабильности кристаллов менее строгие, чем для кристаллической интерферометрии, но установка по-прежнему требует идеального кристалла-анализатора, для которого необходимо очень точно контролировать угол и размер кристалла-анализатора, а также ограничения, налагаемые на луч. должен быть параллельным, что также ограничивает поле зрения. Кроме того, как и в кристаллической интерферометрии, общее ограничение пространственного разрешения этого метода определяется размытием кристалла анализатора из-за эффектов динамической дифракции , но его можно улучшить, используя дифракцию скользящего падения для кристалла. [55]
Хотя этот метод в принципе требует монохроматического, высококоллимированного излучения и, следовательно, ограничен источником синхротронного излучения, недавно было показано, что метод остается осуществимым с использованием лабораторного источника с полихроматическим спектром, когда кривая качания адаптирована к спектральной линии K α . излучение материала мишени. [57]
Благодаря своей высокой чувствительности к небольшим изменениям показателя преломления этот метод хорошо подходит для визуализации образцов мягких тканей и уже применяется в медицинской визуализации, особенно в маммографии для лучшего обнаружения микрокальцификаций [1] и при исследованиях костного хряща. [58]
Визуализация на основе распространения (PBI) — наиболее распространенное название этого метода, но его также называют линейной голографией , визуализацией с усилением рефракции [59] или фазово-контрастной рентгенографией . Последнее название связано с тем, что экспериментальная установка этого метода в основном такая же, как и в обычной рентгенографии. Он состоит из линейного расположения источника рентгеновского излучения, образца и детектора рентгеновского излучения, и никаких других оптических элементов не требуется. Разница лишь в том, что детектор размещается не сразу за образцом, а на некотором расстоянии, поэтому преломленное образцом излучение может интерферировать с неизменившимся лучом. [13] Эта простая установка и низкие требования к стабильности обеспечивают большое преимущество этого метода перед другими методами, обсуждаемыми здесь.
При пространственно-когерентном освещении и промежуточном расстоянии между образцом и детектором создается интерференционная картина с «полосами Френеля»; т.е. полосы возникают при распространении в свободном пространстве в режиме Френеля , что означает, что для расстояния между детектором и образцом справедливо приближение формулы дифракции Кирхгофа для ближнего поля, уравнение дифракции Френеля . В отличие от кристаллической интерферометрии регистрируемые интерференционные полосы в PBI пропорциональны не самой фазе, а второй производной (лапласиану ) фазы волнового фронта. Поэтому метод наиболее чувствителен к резким изменениям декремента показателя преломления. Это приводит к более сильному контрасту, очерчивающему поверхности и структурные границы образца ( улучшение краев ) по сравнению с обычной радиограммой. [60] [61]
PBI можно использовать для повышения контрастности изображения поглощения, в этом случае информация о фазе в плоскости изображения теряется, но вносит вклад в интенсивность изображения ( улучшение краев изображения затухания). Однако возможно также разделить фазу и контраст затухания, т.е. восстановить распределение действительной и мнимой частей показателя преломления отдельно. Однозначное определение фазы волнового фронта ( фазовое восстановление ) можно реализовать, записывая несколько изображений на разных расстояниях детектор-образец и используя алгоритмы, основанные на линеаризации дифракционного интеграла Френеля, для восстановления распределения фазы, но этот подход страдает усиленный шум для низких пространственных частот и, следовательно, медленно меняющиеся компоненты не могут быть точно восстановлены. Существует еще несколько подходов к восстановлению фазы, хороший обзор их дан в [62] [63] .
Томографические реконструкции трехмерного распределения показателя преломления или «голотомография» реализуются путем вращения образца и записи для каждого угла проекции серии изображений на разных расстояниях. [64]
Для разрешения интерференционных полос требуется детектор с высоким разрешением, что практически ограничивает поле зрения этого метода или требует больших расстояний распространения. Достигаемое пространственное разрешение относительно высокое по сравнению с другими методами и, поскольку в пучке отсутствуют оптические элементы, в основном ограничивается степенью пространственной когерентности пучка. Как упоминалось ранее, для формирования полос Френеля ограничение на пространственную когерентность используемого излучения очень строгое, что ограничивает метод небольшими или очень удаленными источниками, но в отличие от кристаллической интерферометрии и визуализации на основе анализатора ограничение на временную когерентность , т.е. полихроматичность достаточно ослаблена. [55] Следовательно, метод можно использовать не только с синхротронными источниками, но и с полихроматическими лабораторными источниками рентгеновского излучения, обеспечивающими достаточную пространственную когерентность, такими как микрофокусные рентгеновские трубки . [60]
Вообще говоря, контрастность изображения, обеспечиваемая этим методом, ниже, чем у других обсуждаемых здесь методов, особенно если изменения плотности в образце невелики. Благодаря своей способности усиливать контраст на границах, он хорошо подходит для визуализации образцов волокон или пенопласта. [65] Очень важным применением PBI является исследование окаменелостей с помощью синхротронного излучения, которое позволяет выявить детали палеонтологических образцов, которые в противном случае были бы недоступны без разрушения образца. [66]
Визуализация на основе решеток (GBI) включает интерферометрию сдвига или рентгеновскую интерферометрию Тальбота (XTI) и полихроматическую интерферометрию в дальнем поле (PFI) . [38] С момента создания первого рентгеновского решеточного интерферометра, состоящего из двух фазовых решеток и кристалла-анализатора [25] , были разработаны различные несколько отличающиеся установки для этого метода; в дальнейшем основное внимание уделяется стандартному в настоящее время методу, состоящему из фазовой решетки и решетки анализатора. [26] (См. рисунок справа).
Метод XTI основан на эффекте Тальбота или «феномене самоизображения», который представляет собой эффект дифракции Френеля и приводит к повторению периодического волнового фронта после определенного расстояния распространения, называемого « длиной Тальбота ». Этот периодический волновой фронт может быть создан пространственно когерентным освещением периодической структуры, такой как дифракционная решетка , и если это так, то распределение интенсивности волнового поля на длине Тальбота точно напоминает структуру решетки и называется собственным изображением. [23] Также было показано, что шаблоны интенсивности будут создаваться при определенных дробных длинах Тальбота. На половине расстояния появляется то же распределение интенсивности, за исключением бокового смещения на половину периода решетки, в то время как на некоторых меньших дробных расстояниях Тальбота собственные изображения имеют дробные периоды и дробные размеры максимумов и минимумов интенсивности, которые становятся видимыми в распределении интенсивности. за решеткой так называемый ковер Талбота. Длина Тальбота и дробные длины могут быть рассчитаны, зная параметры освещающего излучения и освещаемой решетки, и, таким образом, дают точное положение максимумов интенсивности, которое необходимо измерить в GBI. [67] Хотя эффект Тальбота и интерферометр Тальбота были открыты и широко изучены с использованием видимого света, несколько лет назад они были продемонстрированы также и в режиме жесткого рентгеновского излучения. [68]
В GBI образец помещается до или за фазовой решеткой (линии решетки показывают незначительное поглощение, но значительный фазовый сдвиг), и, таким образом, интерференционная картина эффекта Тальбота модифицируется за счет поглощения, преломления и рассеяния в образце. Для фазового объекта с малым фазовым градиентом рентгеновский луч отклоняется на
где k — длина волнового вектора падающего излучения, а второй множитель в правой части — первая производная фазы в направлении, перпендикулярном направлению распространения и параллельном ориентации решетки. Поскольку поперечный сдвиг интерференционных полос линейно пропорционален углу отклонения, дифференциальная фаза волнового фронта измеряется в GBI так же, как и в ABI. Другими словами, угловые отклонения преобразуются в изменения локально передаваемой интенсивности. Выполняя измерения с образцом и без него, можно выявить изменение положения интерференционной картины, вызванное образцом. Период интерференционной картины обычно находится в диапазоне нескольких микрометров , что может быть удобно разрешено только детектором очень высокого разрешения в сочетании с очень интенсивным освещением (источником, обеспечивающим очень высокий поток) и, следовательно, ограничивает поле интерференции. просмотр значительно . [69] Именно по этой причине вторая решетка, обычно поглощающая, размещается на дробной длине Тальбота для анализа интерференционной картины. [26]
Решетка анализатора обычно имеет тот же период, что и интерференционные полосы, и, таким образом, преобразует локальное положение полос в изменение интенсивности сигнала на детекторе, который расположен сразу за решеткой. Чтобы отделить информацию о фазе от других составляющих сигнала, используется метод, называемый «ступенчатым изменением фазы». [27] Одну из решеток сканируют вдоль члена x g в поперечном направлении ; за один период решетки и для разных положений решетки снимается изображение. Сигнал интенсивности в каждом пикселе в плоскости детектора колеблется в зависимости от x g . Зарегистрированные колебания интенсивности могут быть представлены рядом Фурье , и путем записи и сравнения этих колебаний интенсивности с образцом или без него можно выделить отдельный дифференциальный фазовый сдвиг и сигнал поглощения относительно эталонного изображения. [27] Как и при визуализации на основе анализатора, также можно восстановить дополнительный сигнал, исходящий от сверхмалого угла рассеяния субпиксельными микроструктурами образца, называемый контрастом темного поля. [30] Этот метод обеспечивает высокое пространственное разрешение, но также требует длительного времени экспозиции.
Альтернативный подход — восстановление дифференциальной фазы с помощью муаровых полос . Они создаются как суперпозиция собственного изображения G1 и структуры G2 за счет использования решеток с одинаковой периодичностью и наклона G2 относительно G1 относительно оптической оси на очень небольшой угол (<<1). Эти муаровые полосы действуют как полосы несущей, поскольку они имеют гораздо больший интервал/период (меньшую пространственную частоту), чем полосы Тальбота, и, таким образом, фазовый градиент, вносимый образцом, может быть обнаружен как смещение муаровых полос. [26] С помощью Фурье-анализа муара также можно выделить сигнал поглощения и темнопольного сигнала. [70] При использовании этого подхода пространственное разрешение ниже, чем при использовании метода фазового шага, но общее время экспозиции может быть намного короче, поскольку дифференциально-фазовое изображение может быть получено только с одним муаровым узором. [71] Метод однократного анализа Фурье использовался в ранних визуализациях рассеяния на основе сетки [31] аналогично датчику волнового фронта Шака-Гартмана в оптике, что позволило провести первые исследования на живых животных. [72]
Методом устранения механического сканирования решетки и сохранения максимального пространственного разрешения является электронный фазовый шаг. [33] Он сканирует пятно источника рентгеновской трубки с помощью электромагнитного поля. Это заставляет проекцию объекта двигаться в противоположном направлении, а также вызывает относительное движение между проекцией и муаровыми полосами. Изображения смещаются в цифровом виде для выравнивания проекций. В результате проекция объекта остается неподвижной, а муаровые полосы перемещаются по ней. Этот метод эффективно синтезирует процесс фазового шага, но без затрат и задержек, связанных с механическими движениями.
С обоими этими методами фазовой экстракции томография применима путем вращения образца вокруг томографической оси, записи серии изображений с разными углами проекции и использования алгоритмов обратной проекции для восстановления трехмерных распределений действительной и мнимой частей рефракционной части. индекс. [27] [71] Количественная томографическая реконструкция темнопольного сигнала также была продемонстрирована для метода фазового шага [34] , а совсем недавно и для подхода с использованием муара. [70]
Также было продемонстрировано, что визуализация в темном поле с помощью решетчатого интерферометра может использоваться для извлечения информации об ориентации структурных деталей в субмикрометровом режиме за пределами пространственного разрешения системы обнаружения. Если рассеяние рентгеновских лучей в направлении, перпендикулярном линиям решетки, обеспечивает темнопольный контраст, то рассеяние в направлении, параллельном линиям решетки, приводит лишь к размытию изображения, которое не видно при низком разрешении. детектор. [31] Это внутреннее физическое свойство установки используется для извлечения ориентационной информации об угловом изменении локальной рассеивающей способности образца путем вращения образца вокруг оптической оси установки и сбора набора из нескольких темнопольных изображений. изображения, каждое из которых измеряет компонент рассеяния, перпендикулярный линиям решетки для этой конкретной ориентации. Это может быть использовано для определения местного угла и степени ориентации кости и может дать ценную информацию для улучшения исследований и диагностики заболеваний костей , таких как остеопороз или остеоартрит . [73] [74]
Стандартная конфигурация, показанная на рисунке справа, требует пространственной когерентности источника и, следовательно, ограничена источниками синхротронного излучения с высокой яркостью. Эту проблему можно решить, добавив третью решетку рядом с источником рентгеновского излучения, известную как интерферометр Тальбота-Лау . Эта решетка источника, которая обычно представляет собой решетку поглощения с пропускающими щелями, создает «массив индивидуально когерентных, но взаимно некогерентных источников». Поскольку решетка источника может содержать большое количество отдельных отверстий, каждое из которых создает достаточно когерентный виртуальный линейный источник, можно эффективно использовать стандартные генераторы рентгеновского излучения с размерами источника в несколько квадратных миллиметров и значительно увеличить поле зрения. [29]
Поскольку положение интерференционных полос, образующихся за светоделительной решеткой, не зависит от длины волны в широком диапазоне энергий падающего излучения, интерферометр в фазовой конфигурации все еще может эффективно использоваться с полихроматическим излучением. [27] Для конфигурации муарового узора ограничение на энергию излучения немного более строгое, поскольку конечная полоса пропускания энергии вместо монохроматического излучения приводит к уменьшению видимости муаровых полос и, следовательно, к снижению качества изображения, но умеренная полихроматичность нежелательна. все еще разрешено. [75] Большим преимуществом использования полихроматического излучения является сокращение времени экспозиции, и это недавно было использовано при использовании белого синхротронного излучения для реализации первой динамической (разрешенной по времени) фазово-контрастной томографии. [35]
Техническим барьером, который необходимо преодолеть, является изготовление решеток с большим удлинением и малыми периодами. Производство этих решеток из кремниевой пластины включает такие методы микрообработки, как фотолитография , анизотропное влажное травление , гальваника и формование . [76] Очень распространенным процессом изготовления рентгеновских решеток является LIGA , который основан на глубокой рентгеновской литографии и гальванике. Он был разработан в 1980-х годах для изготовления микроструктур с чрезвычайно высоким удлинением учеными из Технологического института Карлсруэ (KIT) . [77] Еще одним техническим требованием является стабильность и точное выравнивание и перемещение решеток (обычно в диапазоне нескольких нм), но по сравнению с другими методами, например кристаллическим интерферометром, это ограничение легко выполнить.
Проблема изготовления решеток была облегчена открытием эффекта фазового муара [38] , который позволил создать интерферометр с полнофазной решеткой, работающий с компактными источниками, называемый полихроматическим интерферометром дальнего поля (см. рисунок справа). Фазовые решетки изготовить проще по сравнению с упомянутыми выше решетками источника и анализатора, поскольку глубина решетки, необходимая для создания фазового сдвига, намного меньше, чем необходимая для поглощения рентгеновских лучей. Фазовые решетки с периодом 200–400 нанометров использовались для улучшения фазовой чувствительности в настольных сканерах PFI. [39] В PFI фазовая решетка используется для преобразования тонких интерференционных полос в широкую картину интенсивности в дистальной плоскости на основе эффекта фазового муара . Помимо более высокой чувствительности, еще одним стимулом для меньших периодов решетки является то, что латеральная когерентность источника должна составлять хотя бы один период решетки.
Недостатком стандартной установки GBI является чувствительность только к одной составляющей фазового градиента, то есть направлению, параллельному одномерным решеткам. Эта проблема была решена либо путем регистрации дифференциально-фазовых контрастных изображений образца в обоих направлениях x и y путем поворота образца (или решеток) на 90° [78] , либо путем использования двумерных решеток. [79]
Будучи дифференциально-фазовым методом, GBI не так чувствителен, как кристаллическая интерферометрия, к низким пространственным частотам, но из-за высокой устойчивости метода к механическим нестабильностям, возможности использования детекторов с большими пикселями и большим полем зрения и, что особенно важно, Важность, применимость к обычным лабораторным рентгеновским трубкам, визуализация на основе решеток является очень перспективным методом медицинской диагностики и визуализации мягких тканей. Первые медицинские применения, такие как доклиническое маммографическое исследование, демонстрируют большой потенциал будущего этого метода. [36] Помимо этого, GBI имеет применение в широкой области материаловедения, например, его можно использовать для улучшения проверки безопасности. [30] [80]
Краевое освещение (EI) было разработано на итальянском синхротроне (Elettra) в конце 90-х годов [19] как альтернатива ABI. Он основан на наблюдении, что при освещении только края пикселей детектора достигается высокая чувствительность к фазовым эффектам (см. рисунок).
Также в этом случае используется связь между углом преломления рентгеновских лучей и первой производной фазового сдвига, вызванного объектом:
Если рентгеновский луч вертикально тонкий и падает на край детектора, рефракция рентгеновских лучей может изменить статус отдельного рентгеновского луча с «обнаруженного» на «необнаруженный» и наоборот, эффективно играя ту же роль, что и кривая качания кристалла в ABI. Эта аналогия с ABI, уже наблюдавшаяся при первоначальной разработке метода, была формально продемонстрирована в работе [19] . [81] Фактически достигается тот же эффект – точный угловой отбор направления фотона; однако, хотя при визуализации на основе анализатора луч должен быть строго коллимированным и монохроматическим, отсутствие кристалла означает, что краевое освещение может быть реализовано с помощью расходящихся и полихроматических лучей, подобных тем, которые генерируются обычной рентгеновской трубкой с вращающимся анодом. . Это делается путем введения двух специально разработанных масок (иногда называемых масками с «кодированной апертурой» [20] ), одной непосредственно перед образцом и одной в контакте с детектором (см. рисунок).
Целью последней маски является просто создание нечувствительных областей между соседними пикселями, и ее использования можно избежать, если использовать специальную технологию детектора. Таким образом, конфигурация краевой подсветки реализуется одновременно для всех рядов пикселей зонального детектора. Такое множество отдельных бимлетов означает, что, в отличие от описанной выше реализации синхротрона, сканирование образца не требуется – образец помещается после маски образца и отображается за один кадр (два, если выполняется фазовый поиск [ 22] ). Хотя устройство, возможно, внешне напоминает решетчатый интерферометр, лежащий в его основе физический механизм иной. В отличие от других методов фазово-контрастной рентгеновской визуализации, краевое освещение является некогерентным методом, и фактически было доказано, что он работает как с пространственно, так и с некогерентными источниками без какой-либо дополнительной апертуры или коллимации источника. [22] [82] Например, обычно используются фокальные пятна размером 100 мкм, которые совместимы, например, с диагностическими маммографическими системами. Количественное восстановление фазы также было продемонстрировано с использованием (неколлимированных) некогерентных источников, показав, что в некоторых случаях могут быть получены результаты, аналогичные синхротронному золотому стандарту. [22] Относительно простая установка краевого освещения приводит к фазовой чувствительности, по крайней мере, сравнимой с другими методами фазово-контрастной рентгеновской визуализации, [83] приводит к ряду преимуществ, которые включают сокращение времени экспозиции при той же мощности источника, снижение дозы радиации, устойчивость к вибрациям окружающей среды и более легкий доступ к высокой энергии рентгеновского излучения. [83] [84] [85] [86] Кроме того, поскольку их соотношение сторон не особенно требовательно, маски дешевы, просты в изготовлении (например, не требуют рентгеновской литографии) и уже могут масштабироваться до больших площадей. Метод легко расширяется на фазовую чувствительность в двух направлениях, например, за счет реализации Г-образных апертур для одновременного освещения двух ортогональных краев в каждом пикселе детектора. [87] В более общем плане, хотя в простейшей реализации бимлеты сопоставляют отдельные строки пикселей (или пиксели), метод очень гибок, и, например, можно использовать разреженные детекторы и асимметричные маски [88] , а также компактные [89] и микроскопические методы. [90] системы могут быть построены. На данный момент метод успешно продемонстрирован в таких областях, как защитное сканирование, [91] биологическая визуализация, [83] [89] материаловедение, [92] палеонтология [93] [94] и других; Также была продемонстрирована адаптация к 3D (компьютерная томография). [93] [95]Помимо простого перевода для использования с обычными источниками рентгеновского излучения, реализация краевого освещения когерентным синхротронным излучением дает существенные преимущества, среди которых высокая производительность при очень высоких энергиях рентгеновского излучения [94] и высокое угловое разрешение. [96]
В контексте медицинской визуализации были выявлены четыре потенциальных преимущества фазового контраста: [42]
Количественное сравнение маммографии с фазовым и абсорбционным контрастом, в котором принимались во внимание реалистичные ограничения (доза, геометрия и экономия фотонов), пришло к выводу, что фазово-контрастная визуализация на основе решеток (интерферометрия Тальбота) не демонстрирует общей разницы между сигналом и улучшение шума по сравнению с контрастом поглощения, но производительность во многом зависит от задачи. [97] [98] Такое сравнение еще предстоит провести для всех методов фазового контраста, однако следующие соображения имеют решающее значение для такого сравнения:
Некоторые компромиссы проиллюстрированы на рисунке справа, где показано преимущество фазового контраста перед контрастом поглощения для обнаружения различных целей, важных для маммографии, в зависимости от размера цели. [97] Обратите внимание, что эти результаты не включают потенциальные выгоды от сигнала темного поля.
После предварительных лабораторных исследований, например, компьютерной томографии [100] и маммографии, [101] фазово-контрастная визуализация начинает применяться в реальных медицинских приложениях, таких как визуализация легких, [102] визуализация конечностей, [103] интра- визуализация оперативного препарата. [104] Начало применению фазово-контрастной визуализации in vivo положило новаторское исследование маммографии с синхротронным излучением, проведенное в Триесте, Италия. [105]
{{cite book}}
: CS1 maint: несколько имен: список авторов ( ссылка ){{cite book}}
: CS1 maint: несколько имен: список авторов ( ссылка ){{cite book}}
: CS1 maint: отсутствует местоположение издателя ( ссылка )