Фазово-контрастная рентгеновская съемка или фазово-чувствительная рентгеновская съемка — это общий термин для различных технических методов, которые используют информацию об изменениях фазы рентгеновского луча , проходящего через объект, для создания его изображений. Стандартные методы рентгеновской съемки, такие как рентгенография или компьютерная томография (КТ), основаны на уменьшении интенсивности рентгеновского луча ( затухании ) при прохождении через образец , которое можно измерить напрямую с помощью рентгеновского детектора . Однако при фазово-контрастной рентгеновской съемке сдвиг фазы луча, вызванный образцом, не измеряется напрямую, а преобразуется в изменения интенсивности, которые затем могут быть зарегистрированы детектором. [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 (European Synchrotron Radiation Facility) в Гренобле, Франция, [13] и был основан на обнаружении «френелевских полос», которые возникают при определенных обстоятельствах при распространении в свободном пространстве. Экспериментальная установка состояла из встроенной конфигурации источника рентгеновского излучения, образца и детектора и не требовала никаких оптических элементов. Она была концептуально идентична установке революционной работы Денниса Габора поАльтернативный подход, называемый визуализацией на основе анализатора, был впервые исследован в 1995 году Виктором Ингалом и Еленой Беляевской в рентгеновской лаборатории в Санкт-Петербурге, Россия, [15] и Тимом Дэвисом и коллегами из CSIRO (Commonwealth Scientific and Industrial Research Organisation) Division of Material Science and Technology в Клейтоне, Австралия. [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 . Это приводит к иному поведению рентгеновских лучей в среде по сравнению с видимым светом (например, углы преломления имеют отрицательные значения), но не противоречит закону относительности , «который требует, чтобы только сигналы, несущие информацию, не распространялись быстрее c . Такие сигналы движутся с групповой скоростью , а не с фазовой скоростью, и можно показать, что групповая скорость на самом деле меньше 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 — константа, заданная в барнах , типичной единице площади поперечного сечения взаимодействия частиц, k — длина волнового вектора , k 0 — длина волнового вектора с длиной волны 1 Ангстрем , а Z — атомный номер . [41] Действительная формула при этих условиях для поперечного сечения сдвига фаз имеет вид:
где Z — атомный номер , k — длина волнового вектора , а r 0 — классический радиус электрона .
Это приводит к следующим выражениям для двух частей комплексного показателя преломления:
Подстановка типичных значений человеческой ткани в приведенные выше формулы показывает, что δ обычно на три порядка больше, чем β в диагностическом рентгеновском диапазоне. Это подразумевает, что сдвиг фазы рентгеновского пучка, распространяющегося через ткань, может быть намного больше, чем потеря интенсивности, что делает фазово-контрастную рентгеновскую визуализацию более чувствительной к изменениям плотности в ткани, чем абсорбционная визуализация. [42]
Из-за пропорциональности
Преимущество фазового контраста над обычным абсорбционным контрастом даже растет с увеличением энергии. Более того, поскольку формирование фазового контрастного изображения не связано по сути с поглощением рентгеновских лучей в образце, поглощенная доза может быть потенциально снижена за счет использования более высоких энергий рентгеновских лучей. [29] [42]
Как упоминалось выше, относительно видимого света, действительная часть показателя преломления n может сильно отклоняться от единицы (n стекла в видимом свете колеблется от 1,5 до 1,8), в то время как отклонение от единицы для рентгеновских лучей в различных средах обычно составляет порядка 10−5 . Таким образом, углы преломления, возникающие на границе между двумя изотропными средами, рассчитанные по формуле Снеллиуса, также очень малы. Следствием этого является то, что углы преломления рентгеновских лучей, проходящих через образец ткани, не могут быть обнаружены напрямую и обычно определяются косвенно путем «наблюдения интерференционной картины между дифрагированными и недифрагированными волнами, создаваемыми пространственными изменениями действительной части показателя преломления». [3]
Кристаллическая интерферометрия , иногда также называемая рентгеновской интерферометрией , является старейшим, но и самым сложным методом, используемым для экспериментальной реализации. Он состоит из трех светоделителей в геометрии Лауэ , выровненных параллельно друг другу. (См. рисунок справа) Падающий луч, который обычно коллимируется и фильтруется монохроматором (кристаллом Брэгга) до этого, разделяется на первом кристалле (S) дифракцией Лауэ на два когерентных луча, опорный луч, который остается невозмущенным, и луч, проходящий через образец. Второй кристалл (T) действует как передающее зеркало и заставляет лучи сходиться друг к другу. Два луча встречаются в плоскости третьего кристалла (A), который иногда называют кристаллом-анализатором, и создают интерференционную картину, форма которой зависит от оптической разности хода между двумя лучами, вызванной образцом. Эта интерференционная картина обнаруживается рентгеновским детектором за кристаллом-анализатором. [9] [43]
Поместив образец на вращающийся столик и записав проекции под разными углами, можно получить трехмерное распределение показателя преломления и, таким образом, томографические изображения образца. [40] В отличие от методов, описанных ниже, с помощью кристаллического интерферометра измеряется сама фаза, а не ее пространственное изменение. Чтобы получить сдвиг фазы из интерференционных картин, используется метод, называемый фазовым шагом или сканированием полос: в опорный луч вводится фазовращатель (в форме клина). Фазовращатель создает прямые интерференционные полосы с регулярными интервалами; так называемые несущие полосы. Когда образец помещается в другой луч, несущие полосы смещаются. Сдвиг фазы, вызванный образцом, соответствует смещению несущих полос. Записывается несколько интерференционных картин для различных сдвигов опорного луча, и, анализируя их, можно извлечь информацию о фазе по модулю 2 π . [40] [43] Эта неоднозначность фазы называется эффектом обертывания фазы и может быть устранена с помощью так называемых «методов развертывания фазы». [44] Эти методы могут использоваться, когда отношение сигнал/шум изображения достаточно высокое, а изменение фазы не слишком резкое. [28]
В качестве альтернативы методу сканирования полос можно использовать метод преобразования Фурье для извлечения информации о сдвиге фаз с помощью только одной интерферограммы, что сокращает время экспозиции, но у этого метода есть недостаток, заключающийся в ограничении пространственного разрешения расстоянием между полосами носителя. [45]
Рентгеновская интерферометрия считается наиболее чувствительной к фазовому сдвигу из 4 методов, следовательно, обеспечивая наивысшее разрешение по плотности в диапазоне мг/см3 . [ 28] Но из-за ее высокой чувствительности полосы, созданные образцом с сильным фазовым сдвигом, могут стать неразрешимыми; для преодоления этой проблемы недавно был разработан новый подход, называемый «когерентно-контрастная рентгеновская визуализация», где вместо фазового сдвига для контрастности изображения имеет значение изменение степени когерентности, вызванное образцом. [46]
Общее ограничение пространственного разрешения этого метода обусловлено размытием в кристалле анализатора, которое возникает из-за динамической рефракции, т. е. угловое отклонение луча из-за рефракции в образце усиливается примерно в десять тысяч раз в кристалле, поскольку путь луча внутри кристалла сильно зависит от угла падения. Этот эффект можно уменьшить, утоньшив кристалл анализатора, например, при толщине анализатора 40 мкм было рассчитано разрешение около 6 мкм . В качестве альтернативы кристаллы Лауэ можно заменить кристаллами Брэгга , так что луч не проходит через кристалл, а отражается от поверхности. [47]
Другим ограничением метода является требование очень высокой стабильности установки; выравнивание кристаллов должно быть очень точным, а разница в длине пути между пучками должна быть меньше длины волны рентгеновских лучей; для достижения этого интерферометр обычно изготавливается из высокосовершенного цельного блока кремния путем вырезания двух канавок. Благодаря монолитному производству очень важная пространственная когерентность решетки между всеми тремя кристаллами может поддерживаться относительно хорошо, но это ограничивает поле зрения небольшим размером (например, 5 см x 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 является чувствительность только к одному компоненту градиента фазы, который является направлением, параллельным решеткам 1-D. Эта проблема была решена либо путем записи изображений контраста дифференциальной фазы образца в обоих направлениях x и y путем поворота образца (или решеток) на 90° [78], либо путем использования двумерных решеток. [79]
Будучи дифференциально-фазовой техникой, GBI не так чувствителен, как кристаллическая интерферометрия, к низким пространственным частотам, но из-за высокой устойчивости метода к механическим нестабильностям, возможности использования детекторов с большими пикселями и большим полем зрения и, что крайне важно, применимости к обычным лабораторным рентгеновским трубкам, визуализация на основе решеток является очень многообещающим методом для медицинской диагностики и визуализации мягких тканей. Первые медицинские применения, такие как доклиническое маммографическое исследование, показывают большой потенциал для будущего этой техники. [36] Помимо этого, GBI имеет применение в широкой области материаловедения, например, его можно использовать для улучшения досмотра. [30] [80]
Edge-illumination (EI) была разработана на итальянском синхротроне (Elettra) в конце 90-х годов [19] как альтернатива ABI. Она основана на наблюдении, что при освещении только края пикселей детектора достигается высокая чувствительность к фазовым эффектам (см. рисунок).
В этом случае также используется связь между углом преломления рентгеновских лучей и первой производной фазового сдвига, вызванного объектом:
Если рентгеновский луч вертикально тонкий и падает на край детектора, рефракция рентгеновских лучей может изменить статус отдельного рентгеновского луча с «обнаруженного» на «необнаруженный» и наоборот, фактически играя ту же роль, что и кривая качания кристалла в ABI. Эта аналогия с ABI, уже наблюдавшаяся при первоначальной разработке метода [19] , была недавно формально продемонстрирована. [81] По сути, достигается тот же эффект — точная угловая селекция направления фотона; однако, в то время как при визуализации на основе анализатора луч должен быть высококоллимированным и монохроматичным, отсутствие кристалла означает, что краевое освещение может быть реализовано с помощью расходящихся и полихроматических пучков, подобных тем, которые генерируются обычной рентгеновской трубкой с вращающимся анодом. Это делается путем введения двух своевременно разработанных масок (иногда называемых масками с «кодированной апертурой» [20] ), одна непосредственно перед образцом, а другая в контакте с детектором (см. рисунок).
Цель последней маски — просто создать нечувствительные области между соседними пикселями, и ее использования можно избежать, если использовать специализированную технологию детектора. Таким образом, конфигурация краевой подсветки одновременно реализуется для всех рядов пикселей площадного детектора. Это множество отдельных бимлетов означает, что, в отличие от реализации синхротрона, обсуждавшейся выше, не требуется сканирования образца — образец помещается ниже маски образца и отображается за один снимок (два, если выполняется фазовое восстановление [22] ). Хотя установка, возможно, поверхностно напоминает установку решетчатого интерферометра, лежащий в основе физический механизм отличается. В отличие от других методов фазово-контрастной рентгеновской визуализации краевая подсветка является некогерентной техникой, и на самом деле было доказано, что она работает как с пространственно, так и с временно некогерентными источниками, без какой-либо дополнительной апертуры или коллимации источника. [22] [82] Например, обычно используются фокусные пятна размером 100 мкм, которые совместимы, например, с диагностическими маммографическими системами. Количественное восстановление фазы также было продемонстрировано с (неколлимированными) некогерентными источниками, показав, что в некоторых случаях можно получить результаты, аналогичные синхротронному золотому стандарту. [22] Относительно простая установка краевой подсветки приводит к фазовой чувствительности, по крайней мере, сравнимой с другими методами фазово-контрастной рентгеновской визуализации, [83] приводит к ряду преимуществ, которые включают в себя сокращенное время экспозиции для той же мощности источника, сокращенную дозу облучения, устойчивость к вибрациям окружающей среды и более легкий доступ к высокой энергии рентгеновского излучения. [83] [84] [85] [86] Более того, поскольку их соотношение сторон не особенно требовательно, маски дешевы, просты в изготовлении (например, не требуют рентгеновской литографии) и уже могут быть масштабированы для больших площадей. Метод легко распространяется на фазовую чувствительность в двух направлениях, например, посредством реализации L-образных апертур для одновременного освещения двух ортогональных краев в каждом пикселе детектора. [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}}
: |journal=
проигнорировано ( помощь ){{cite book}}
: CS1 maint: отсутствует местоположение издателя ( ссылка )