Как многодиапазонные Rf-пульсы создаются на основе принципов физики мри
Table of Contents
Введение
Неустанное стремление к более быстрой и эффективной магнитно-резонансной визуализации (МРТ) привело к значительным достижениям в разработке импульсной последовательности. Среди наиболее эффективных инноваций является разработка многодиапазонных радиочастотных (РЧ) импульсов, которые позволяют одновременно возбуждать и считывать несколько срезов или областей. Этот метод, широко известный как одновременная многослойная (СМС) визуализация, особенно преобразует такие приложения, как функциональная МРТ, диффузионная визуализация и сердечный киносин. Однако проектирование этих импульсов требует глубокого понимания основных принципов физики МРТ, включая уравнение Лармора, уравнения Блоха, выбор градиентного среза и связь Фурье между формой радиочастотного импульса и возбужденным частотным профилем. Эта статья обеспечивает подробное исследование того, как многодиапазонные радиочастотные импульсы разработаны, физические принципы, которые управляют их поведением, инженерные проблемы, связанные с и передовые стратегии, используемые для достижения высококачественных, без артефактов изображений с резко сокращенным временем сканирования.
Основы РЧ возбуждения в МРТ
Уравнение Лармора и ядерная магнетизация
В основе МРТ лежит феномен ядерного магнитного резонанса. Ядра водорода (протоны) обладают магнитным моментом и угловым моментом (спином). При размещении в статическом магнитном поле B0 эти моменты предшествует направлению поля на частоте Лармора ω0 = γ B0, где γ — гиромагнитное отношение (приблизительно 42,58 МГц/Т для протонов). Чистая макроскопическая намагниченность M0 выравнивается с B0 по продольной оси (z-оси). RF-импульсы наносятся на частоту Лармора, чтобы отвести эту намагниченность от равновесия, генерируя поперечную намагниченность, которую можно обнаружить как сигнал МР. Количество угла наклона определяется амплитудой и продолжительностью RF-импульса (угол скольжения). В обычной визуализации одиночный RF-импульс возбуждает один срез, одновременно применяя линейное градиентное поле во время импульса, которое кодирует пространственное положение вдоль направления среза-выбора с уникальной
Выбор срезов и роль градиентных полей
Выбор среза опирается на использование градиента магнитного поля (например, Gz), приложенного вдоль оси ортогональной к желаемой плоскости среза. Этот градиент заставляет частоту Лармора изменяться линейно с положением: ω(z) = γ(B0 + Gz·z). RF-импульс с конечной полосой Δf, центрированной на некоторой частоте f0, будет возбуждать только те спины, резонансные частоты которых попадают в эту полосу пропускания. Толщина среза определяется отношением полосы импульса к амплитуде градиента (Δz = Δf/ (γ Gz)). Обычный RF-импульс — часто синхронный импульс — возбуждает одну прилегающую плиту ткани. Для многодиапазонной визуализации цель состоит в том, чтобы возбуждать несколько дискретных срезов одновременно с использованием одного RF-импульса. Это достигается путем проектирования RF-формы волны, частотный спектр которой содержит несколько различных полос, каждая из которых соответствует различному расположению с
Уравнение Блоха и проблема проектирования импульса РФ
Эволюция намагниченности под комбинированным влиянием полей RF и полей градиента регулируется уравнениями Блоха. Проектирование импульса RF по существу является обратной проблемой — при наличии желаемого пространственного профиля поперечного намагничивания (например, однородного возбуждения над несколькими тонкими срезами с нульами между ними), мы должны определить форму волны RF и модуляцию градиента, которые его производят. Для небольших углов сгиба отношения между оболочкой RF импульса B1(t) и полученным профилем возбуждения M xy(z) хорошо аппроксимированы линейным преобразованием Фурье: M xy(z) ≈ iγ M0 ⁇ B1(t) exp(-iγ ⁇ G z(t') dt'·z) dtt, где интегралом является траектория k-пространства в измерении срез-выбора. Для больших углов сгиба нелинейность уравнений Блоха требует более сложных методов, таких как преобразование Шиннар-Ле Ру
Принципы многодиапазонного RF-пульсового дизайна
Одновременное возбуждение нескольких срезов
Простейший подход к многодиапазонному возбуждению заключается в суммировании RF-импульсных оболочек, которые индивидуально возбуждали бы каждый желаемый срез. Например, если мы хотим возбуждать два среза в положениях z1 и z2 с толщиной δz, мы бы спроектировали два срез-селективных импульса (например, синхронизированных) с несущими частотами, смещенными на γ Gz·(z1 - z2), а затем суммировали их волновые формы. В режиме малого углового возбуждения общее возбуждение является линейным, поэтому профиль намагничивания представляет собой суперпозицию двух отдельных профилей. Эта прямая суперпозиция, однако, вводит рябь артефактов и нежелательных боковых долей в профиле, особенно если срезы находятся близко друг к другу. Более общая конструкция рассматривает многодиапазонный импульс как единую композитную форму волны, предназначенную для получения частотной реакции, содержащей несколько полос пропускания. Это достигается путем умножения базового срез-селективного импульса на сумму сложных экс
Частота и фазовая модуляция
Описанный подход модуляции представляет собой форму комбинированной амплитуды и фазовой модуляции. На практике форма импульсной волны RF является сложной (компоненты I и Q). При назначении отдельного смещения фазы каждой полосе мы можем уменьшить конструктивную интерференцию, которая может привести к увеличению пиковой мощности (и, следовательно, более высокой удельной скорости поглощения, SAR). Простейший многодиапазонный импульс, используя равномерное смещение фазы по полосам, часто дает высокую амплитуду пика. Общая стратегия заключается в использовании циклической фазы — назначении случайных или линейно увеличивающихся фаз для каждого среза возбуждения. Например, если у нас есть восемь одновременных срезов, мы можем использовать фазовое увеличение 180° между последовательными полосами для уменьшения амплитуды пика путем уничтожения когерентного сложения. Эта фазовая модуляция может быть выражена как B1−MB(t) = B1−base(t) · Σn exp(i[ωn t + φn], где φn являются тщательно выбранными фазами (часто оптимизированными
Crosstalk и Band разделились
Одной из основных проблем в конструкции многодиапазонного импульса является минимизация перекрестного намагничивания между близко расположенными срезами. Кросстолк проявляется как нежелательная намагниченность в смежных срезах, возникающая из-за конечной ширины профиля среза и спектрального перекрытия между полосами. Это может привести к загрязнению сигнала и уменьшению временного разрешения для исследований временных рядов. Три основные стратегии направлены на перекрестный нарез: (1) оптимизация крутизны профиля среза с использованием более длинных импульсов RF с более высокими продуктами полосы времени (например, использование большего количества боковых долей для импульсов синхронизации); (2) применение градиентных дробилок после импульса RF, который дефазирует остаточное поперечное намагничивание от соседних срезов; и (3) использование алгоритмов параллельной реконструкции изображений, которые используют различные модели чувствительности катушки для отдельных перекрывающихся сигналов. В современных реализациях многодиапазонные импульсы часто сочетаются с градиентной схемой с перевернутым CAIPIRINHA (контролируемое удаление в параллельных результатах визуализации при
Стратегии проектирования многоканальных пульсов РФ
Мультиполосные импульсы на основе Фурье
Для небольших углов скольжения (обычно ниже 30°) верно приближение Фурье уравнения Блоха, и импульсы могут быть сконструированы непосредственно в частотной области. Профиль целевого возбуждения представляет собой сумму прямоугольных функций, центрированных в нужных положениях среза. Обратное преобразование Фурье этого профиля дает комплексно-значный RF-импульс. Поскольку профиль содержит множество различных полос, полученная в результате форма волны домена времени представляет собой модулированную версию однополосного синхронного импульса. Этот метод является быстрым и интуитивно понятным, но он производит импульсы с высокой пиковой амплитудой (особенно когда многие срезы возбуждаются одновременно) и потенциально плохие срезные профили из-за звона Гиббса. Для смягчения звона прямоугольный профиль заменяется более гладкой функцией (например, увеличенной косинусной или Хамминговой окне), которая уменьшает боковые доли в импульсе. Компромисс немного увеличенной толщины среза или уменьшенного разделения полос. Несмотря на эти ограничения, мультиполосные импульсы на основе Фурье
Shinnar-Le Roux (SLR) - Трансформация для многополосных импульсов
Для больших углов скольжения (например, 90° или 180°) приближение Фурье не удается, поскольку уравнения Блоха становятся нелинейными. Преобразование СЛР является точным методом проектирования импульсов для любого угла скольжения путем сопоставления задачи с конструкцией цифровых фильтров. Многополосный импульс на основе СЛР проектируется путем сначала построения многополосного пространственного профиля — то есть шаблона намагниченности против z, который чередуется между нулем (нулями) и желаемым углом скольжения (пассаблями). Этот профиль разлагается на набор полиномов Чебышева, которые в свою очередь определяют RF-импульс. СЛР-импульсы могут достигать чрезвычайно резких переходов между полосами и очень плоскими полосами пропускания, но они требуют тщательной оптимизации длины фильтра и ряби. Для многополосного алгоритма СЛР можно расширить путем проектирования однополосного импульса с очень высоким продуктом с пропускной способностью во времени (для достижения тонких срезов) и затем модулирования его, как и раньше. Альтернативно, сам профиль может быть спроектирован с несколькими полос
Изменяемая селективная скорость возбуждения (VERSE)
VERSE — это методика, которая модифицирует амплитуду и силу градиента RF одновременно для уменьшения пиковой амплитуды RF или для сокращения длительности импульса при сохранении одного и того же профиля возбуждения. В многодиапазонной визуализации VERSE особенно ценна, поскольку многодиапазонные импульсы часто имеют высокие пиковые амплитуды, которые приближаются к аппаратным пределам и увеличивают SAR. Принцип VERSE заключается в применении изменяющейся во времени градиентной формы волны, которая медленнее пересекает k-пространство при высоких амплитудах RF и быстрее при низких амплитудах RF. Это эффективно распределяет энергию RF, понижая пик, но потенциально удлиняя импульс. Для многодиапазона VERSE может применяться к композитной форме волны или к каждой полосе отдельно. Более продвинутые версии используют VERSE для уменьшения перекрытия между оболочкой RF и градиентной формой волны, таким образом также уменьшая чувствительность импульса к несовершенствам градиента. При сочетании с параллельной
Оптимизация и итеративные методы
Современная многодиапазонная конструкция импульсов все больше опирается на численную оптимизацию для решения конкурирующих требований длительности импульса, SAR, точности угла поворота и подавления артефактов. Проблема часто ставится как ограниченная оптимизация: минимизация пиковой амплитуды RF (или общей энергии), подлежащая уравнению Блоха и ограничениям на профиль возбуждения при каждом положении среза. Для случаев с малым углом наклона может использоваться выпуклое квадратичное программирование; для случаев с большим углом наклона применяются нелинейные методы, такие как теория оптимального управления (например, алгоритм градиентного спуска или алгоритм Бройдена-Флетчера-Голдфарба-Шанно). Эти методы позволяют включать аппаратные ограничения (максимальный B1, максимальная скорость убоя), ограничения SAR (глобальные и локальные) и даже параллельные поля передачи (B1+ карты). Полученные импульсы часто имеют неправильные формы, которые не просто модулируются синхронностями, но они достигают значительно лучшей производительности с точки зрения качества профиля и SAR.
Практические ограничения и интеграция с параллельными изображениями
Управление уровнем поглощения (SAR)
SAR является основной проблемой для многодиапазонных импульсов, особенно при сверхвысоких напряжённостях поля (7 Т и выше). Мультидиапазонный импульс откладывается энергией на несколько срезов одновременно, поэтому и пиковый, и средний SAR увеличиваются примерно линейно с количеством полос, если не применяется оптимизация. Специфическая скорость поглощения регулируется интегралом времени квадрата RF-оболочки, а любое снижение амплитуды пика или длительности импульса напрямую уменьшает SAR. Для сохранения SAR в пределах FDA/IEC необходимы такие методы, как фазовый цикл, VERSE и оптимизация. Кроме того, использование параллельных передающих (pTx) систем, где приводятся в действие несколько независимых радиочастотных каналов с различными волновыми формами, позволяет пространственное распределение энергии RF для уменьшения локальных горячих точек SAR. Мультидиапазонная конструкция импульса pTx остается активной областью исследований, направленных на объединение ускорения одновременного мультирезультата с безопасностью параллельной передачи. Для клинических сканеров, работающих на 3 Т, многодиапазонные импульсы с
Комбинация с SENSE и GRAPPA
Одного многополосного возбуждения недостаточно для получения отдельных изображений из каждого среза, поскольку все срезы способствуют одному и тому же считывающему сигналу. Разделение срезов сигналов осуществляется в реконструкции изображения с использованием известных профилей чувствительности катушки приемника массива. Это аналогично параллельному ускорению изображения в направлении среза-кодирования. Такие методы, как SENSE (SENSitivity Encoding) или GRAPPA (GeneRalized Autocalibrating Partially Parallel Acquisition) расширяются для обработки одновременного срезания. В SENSE раскрытие псевдосрезанных изображений использует карты чувствительности катушки; в GRAPPA недостающие k-пространственные линии заполняются с использованием ядра, оцененного по данным автокалибровки. Комбинация многополосного возбуждения и параллельного изображения называется одновременным многослойным (SMS) изображением и широко реализуется в коммерческих последовательностях (например, Siemens "SMS", GE "Hyperband", Philips "MultiBand"). Коэффициент ускорения от SMS умножается на ускорение в
Применение многодиапазонных пульсов РФ
Функциональная МРТ (fMRI) и сети государств отдыха
Одно из самых ранних и наиболее эффективных применений многополосных импульсов - в фМРТ, особенно для визуализации, зависящей от уровня кислорода в крови (BOLD). Способность приобретать покрытие всего мозга с высоким временным разрешением (например, TR < 1s) произвела революцию в изучении динамики мозга. Многополосные коэффициенты ускорения 4-8 позволяют отбирать пробы всего мозга каждые 400-800 мс, захватывая быстрые колебания сетей состояния покоя и ответов, вызванных задачами. Конструкция многополосных импульсов для фМРТ часто подчеркивает низкие SAR и короткие TE, требуя высокоэффективных импульсов, которые работают под малыми углами скольжения (обычно 10-20 °). Стандартный подход - это многополосный импульс на основе Фурье с повышенным козиновым профилем, за которым следует перевернутый градиент CAIPIRINHA для улучшения разделения срезов. Последние разработки включают использование обобщенных срезов ускоренных двухэхо-приобретений, которые непосредственно извлекают выгоду из улучшенной статистической мощности, обеспечиваемой более высоким временным разрешением, что приводит к лучшему обнаружению функциональн
Диффузионная МРТ и трактография
Диффузионно-взвешенная визуализация (DWI) требует длительного времени сканирования для получения нескольких направлений градиента и b-значений. Многодиапазонные импульсы сыграли важную роль в создании диффузионной визуализации с высоким угловым разрешением (HARDI) и диффузионно-тензорной визуализации (DTI) клинически осуществимыми. Одновременно захватывающие 2-4 срезов, эффективное покрытие объема на TR увеличивается, уменьшая общее время получения на один и тот же фактор. Однако диффузионные градиенты сильны и могут вызывать вихревые токи и артефакты движения, которые усугубляются многодиапазонными возбуждениями. Поэтому конструкция RF-импульса для диффузии часто включает тщательное рассмотрение градиентного момента нулирования и использование двухдифференцированных спиновых эхо-эхо для смягчения вихревых токов. Кроме того, многодиапазонные диффузионные последовательности часто используют специальную «пирамидную» форму RF-импульса для минимизации распада T2* в течение длительного эхо-в
Кардиологическая визуализация и динамика в реальном времени
Для МРТ сердца требуется как высокое временное разрешение, так и широкий пространственный охват для изображения бьющегося сердца. Многополосные радиоимпульсы позволяют многослойное получение в пределах одной сердечной фазы, уменьшая количество удержания дыхания или позволяя свободно дышать. Для сбалансированных последовательностей с постоянным состоянием свободной прецессии (bSSFP), которые являются основой сердечного синего, конструкция многополосных импульсов должна обеспечивать устойчивое состояние во всех одновременно возбужденных срезах. Это достигается с помощью поезда идентичных многополосных импульсов с постоянным TR и тщательным испорчением остаточной поперечной намагниченности. Расширенные конструкции, такие как радиальная выборка «стоп-звезд» в сочетании с многополосным возбуждением, производят очень надежные изображения с пониженной чувствительностью к движению. Управление SAR особенно важно в сердечных приложениях из-за большой массы ткани и потенциала для радиочастотного нагрева вблизи грудной клетки, необходимы для соблюдения пределов безопасности при сохранении высоких углов переворачивания (40°-60°), типичных для bSSFP.
Будущие направления и новые технологии
Эволюция многодиапазонной конструкции RF-импульса продолжает раздвигать границы ускорения МРТ. Методы машинного обучения, в частности глубокие нейронные сети, применяются к обратной задаче проектирования импульсов. Эти модели могут изучать сложное отображение из желаемых многодиапазонных профилей в оптимальные RF-формы волн, учитывая ограничения, такие как SAR и аппаратная нелинейность. Ранние результаты показывают, что нейронные сети могут генерировать импульсы, которые превышают производительность обычных алгоритмов оптимизации с точки зрения точности профиля и длительности импульсов. Другим перспективным направлением является использование сверхнизкополевой МРТ (например, 0,55 Т), где более низкая частота Лармора уменьшает SAR, позволяя более высокие мультидиапазонные факторы без превышения нормативных ограничений. При сверхвысоких полях (10,5 Т и выше) комбинация многодиапазонных и pTx становится незаменимой, и исследователи активно работают над совместной оптимизацией как передающих, так и принимающих массивов для одновременного ускорения среза. Наконец, концепция «g-фактор-свободных» SMS с использованием адаптированных
Заключение
Многополосные радиочастотные импульсы представляют собой краеугольный камень современной ускоренной МРТ. Их конструкция глубоко укоренилась в фундаментальной физике магнитного резонанса, включая уравнение Лармора, уравнения Блоха и отношения преобразования Фурье между радиочастотными импульсами и пространственными профилями. Ключевые стратегии проектирования - от простой суперпозиции Фурье до продвинутой оптимизации и SLR-трансформаций - позволили перевести эту технику в рутинную клиническую практику. Интеграция с параллельной визуализацией и тщательным управлением SAR сделала одновременное многослойное получение безопасным и эффективным для широкого спектра применений, от МРТ мозга до сердечной визуализации. По мере того, как спрос на более быстрое, более высокое разрешение и более всеобъемлющую МРТ продолжает расти, физико-управляемая конструкция многополосных радиочастотных импульсов останется важной областью исследований и инноваций.