Картирование B1-поля в магнитно-резонансной томографии

Картирование пространственного распределения радиочастотного B1-поля представляет собой ключевую технологию современной магнитно-резонансной томографии, обеспечивающую количественную оценку пространственной неоднородности радиочастотного возбуждения. Значение метода возрастает по мере распространения высокопольных МР-систем (3 Тл и выше), внедрения параллельной передачи радиочастотных импульсов (pTx) и количественных методов визуализации, в которых неоднородность B1-поля напрямую определяет качество изображений и точность измеряемых тканевых параметров. За последние два года в области произошли заметные изменения: появились сверхбыстрые трёхмерные реализации классических методик, нейросетевые алгоритмы, предсказывающие карты B1+ по стандартным локализационным изображениям, а динамическая параллельная передача и универсальные импульсы перешли из стадии методических разработок в стадию клинического внедрения на системах 7 Тл.

Введение

Магнитно-резонансная томография как метод неинвазивной диагностики продолжает развиваться в направлении повышения пространственного и временного разрешения, а также улучшения количественных характеристик получаемых изображений. Одним из ключевых факторов, влияющих на качество МР-изображений и точность количественных измерений, является пространственная однородность радиочастотного B1-поля, создаваемого передающими катушками.

Неоднородность B1-поля проявляется в виде пространственных вариаций угла поворота намагниченности, что приводит к искажениям контрастности изображений, артефактам затенения и систематическим ошибкам в количественных измерениях релаксационных параметров. Проблема становится особенно значимой при использовании высокопольных МР-систем (3 Тл и выше): по мере роста частоты Лармора длина радиочастотной волны в тканях уменьшается и становится сопоставимой с размерами исследуемой области, вследствие чего эффекты диэлектрического резонанса и интерференции волн формируют выраженные зоны усиления и ослабления поля. На системах 7 Тл эти зоны затрагивают практически всё исследуемое пространство, а в абдоминальной области возможно формирование участков, где амплитуда B1+ приближается к нулю («зоны выпадения» сигнала).

Развитие технологий параллельной передачи радиочастотных импульсов (pTx) открывает новые возможности для коррекции неоднородностей B1-поля, однако требует знания пространственного распределения поля каждого передающего элемента. Это обстоятельство делает методы B1-картирования не просто желательными, а необходимыми для реализации полного потенциала современных МР-систем. Показательно, что принципиальная идея параллельной передачи — использование нескольких независимо управляемых передающих каналов для формирования заданного профиля возбуждения — была сформулирована как аналог параллельного приёма ещё в середине 2000-х годов, а её клиническая реализация потребовала почти двух десятилетий развития методов калибровки, среди которых B1-картирование занимает центральное место.

Количественные методы МРТ, такие как картирование времён релаксации T1 и T2, диффузионно-тензорная визуализация, магнитно-резонансная спектроскопия и методы, основанные на переносе насыщения (CEST и магнитизационный перенос), особенно чувствительны к неоднородностям B1-поля. В этих приложениях даже небольшие отклонения угла поворота от номинального значения могут приводить к существенным ошибкам в определении тканевых параметров, что критично для клинической диагностики и мониторинга терапевтического ответа.

Терминология: B1+ и B1−, относительные и абсолютные карты

В МР-физике принято различать две компоненты радиочастотного поля. Передающее поле B1+ — вращающаяся в направлении прецессии компонента поля, создаваемого передающей катушкой, — определяет фактический угол поворота намагниченности α, который при заданных форме и длительности импульса пропорционален амплитуде B1+. Приёмная чувствительность B1− характеризует способность катушки регистрировать сигнал из данной точки пространства и определяет неоднородность яркости изображения, но не угол поворота. В задачах, рассматриваемых в настоящей статье, под «картой B1» понимают, как правило, карту B1+, то есть пространственное распределение отношения фактического угла поворота к номинальному. Поправка на неоднородность приёмной чувствительности B1− обычно выполняется отдельно, в рамках стандартных алгоритмов коррекции поля смещения интенсивности (bias field correction), широко используемых в программных пакетах обработки нейровизуализационных данных.

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

Основные механизмы влияния анатомии на B1-поле (на примере головного мозга)

Геометрические параметры головы и распределение поля

Размеры черепной коробки и объём мозга влияют на пространственное распределение B1-поля. Увеличение общих размеров головы, как правило, сопровождается ростом неоднородности B1+: расстояние между поверхностью катушки и различными областями мозга изменяется, а вместе с ним — степень затухания и интерференции электромагнитного поля. Центральные структуры (таламус, базальные ганглии) демонстрируют более низкую вариабельность амплитуды поля, тогда как периферические зоны — в первую очередь лобная кора и нижние отделы височных долей — подвержены наибольшим отклонениям от номинального значения. Индивидуальные различия геометрии головы и положения пациента в катушке столь существенны, что фиксированный набор настроек радиочастотного возбуждения не может быть оптимален для всей популяции, и именно это обстоятельство обосновывает необходимость индивидуального B1-картирования либо использования специально сконструированных «универсальных» импульсов (см. ниже).

Диэлектрические свойства тканей

Диэлектрическая проницаемость и проводимость различных тканей мозга существенно различаются. Белое вещество содержит меньшую объёмную долю свободной воды по сравнению с серым, что определяет различия в затухании радиочастотного поля на высоких частотах. Ликвор, обладающий высокой проводимостью и диэлектрической проницаемостью, формирует границы раздела с выраженным контрастом диэлектрических свойств, на которых происходит перераспределение поля. В результате при ультравысоких полях центральные структуры нередко сохраняют интенсивность сигнала благодаря более однородным локальным диэлектрическим свойствам, тогда как кортикальные области демонстрируют артефакты затухания. Эти различия требуют адаптивных методов коррекции B1-карт в зависимости от типа ткани и анатомической области.

Возрастные изменения и объём ЦСЖ

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

Зависимость от напряжённости основного поля и анатомической области

Неоднородность B1+ не является исключительной проблемой исследований головного мозга и систем ультравысокого поля. Уже при 3 Тл существенные неоднородности наблюдаются в молочной железе: показано, что они негативно влияют на точность определения T1 методом переменного угла поворота (variable flip angle, VFA), а следовательно — на параметры динамического контрастного усиления, и что картирование B1 методом сдвига Блоха — Зигерта повышает точность и воспроизводимость измерений T1 в этой области. В области живота и таза при 7 Тл динамический диапазон амплитуды B1+ особенно велик, что делает задачу абсолютного картирования в теле человека существенно более сложной, чем в голове, и стимулирует разработку методик с расширенным динамическим диапазоном.

Методы B1-картирования на основе анализа амплитуды сигнала

Метод двойного угла (Double Angle Method, DAM) — классический подход к B1-картированию, основанный на сравнении интенсивностей сигнала, полученных при использовании радиочастотных импульсов с углами поворота α и 2α. Теоретическая основа метода — различная зависимость интенсивности сигнала от угла поворота: отношение интенсивностей зависит от косинуса угла α и, при полном восстановлении намагниченности, не зависит от T1.

Для устранения зависимости от времени релаксации T1 в методе двойного угла необходимо использовать длительные времена повторения (TR порядка нескольких T1), обеспечивающие практически полное восстановление продольной намагниченности между возбуждениями. Это ограничение существенно увеличивает время сканирования — в классических реализациях в десятки минут на объём, — делая исходный вариант метода малопригодным для рутинного клинического применения. Точность DAM ограничивается несколькими факторами: неидеальностью профиля среза в двумерных реализациях, чувствительностью к неоднородностям основного магнитного поля B0 и влиянием несовершенного подавления остаточной поперечной намагниченности. Эти ограничения особенно выражены в областях с выраженной неоднородностью магнитного поля.

В последние годы интерес к модификациям DAM сохраняется именно благодаря возможности реализации на стандартных, поставляемых производителем последовательностях: предложены варианты с эхо-планарным считыванием, а также двухугловые методики с коррекцией по T1, позволяющие получать карты B1 за клинически приемлемое время без разработки специальных импульсных последовательностей. Сравнение такого подхода с методами AFI и сдвига Блоха — Зигерта при 3 Тл в головном мозге показало его пригодность для введения поправки на неоднородность B1 при количественном картировании T1.

Метод фактического угла поворота (Actual Flip-angle Imaging, AFI) — усовершенствованный вариант двухугловых подходов, разработанный для преодоления основных ограничений классической схемы. Ключевой особенностью AFI является использование двух времён повторения, TR1 и TR2, где TR2 в несколько раз превышает TR1; два последовательных возбуждения в стационарном режиме позволяют определить фактический угол поворота по отношению сигналов независимо от T1 при условии эффективного подавления остаточной поперечной намагниченности (спойлинга).

Преимуществами метода AFI являются возможность трёхмерной реализации с хорошим анатомическим покрытием, повышенное отношение сигнал/шум по сравнению с двумерными мультисрезовыми методами и относительная нечувствительность к эффектам профиля среза. Недостатки метода — необходимость тщательной настройки спойлинга, требование точной регистрации изображений для устранения артефактов движения и ограничения по удельному поглощению энергии (SAR), связанные с применением интенсивных подавляющих градиентов. Время сканирования остаётся порядка минут, поэтому AFI нередко используется как «золотой» эталон для проверки и калибровки более быстрых методик: так, в работах 2025 года по картированию B1+ для массивов параллельной передачи абсолютные карты строятся путём калибровки относительно единственной AFI-карты, полученной при одновременной передаче по всем каналам.

Метод ускоренного 3D Look-Locker

Отдельное направление представляют подходы, основанные на модификации ускоренных трёхмерных последовательностей типа Look-Locker, первоначально разработанных для измерения времени релаксации T1. В рамках этого подхода предложены варианты без инвертирующих импульсов, в которых чередование радиочастотных импульсов с различной амплитудой позволяет получать карты угла поворота в течение времени порядка минуты при трёхмерном покрытии. Заявленным преимуществом является сочетание высокой скорости сканирования и приемлемой точности по сравнению с AFI; ограничения связаны со снижением точности на краях диапазона углов поворота и чувствительностью к несовершенному подавлению переходных процессов. Данные об этом подходе ограничены и нуждаются в независимой многоцентровой валидации.

Метод насыщающего TurboFLASH

Широкое распространение получили методики с предварительным насыщающим радиочастотным импульсом и быстрым считыванием по типу TurboFLASH (saturated TurboFLASH, satTFL). Номинальный угол поворота насыщающего импульса «кодируется» в отношении сигналов, полученных с насыщением и без него, а последующее быстрое считывание обеспечивает приемлемую продолжительность исследования. Метод применяется как для головы, так и для тела и обладает достаточно широким динамическим диапазоном, что делает его одним из референсных в абдоминальных исследованиях при 7 Тл.

Методы B1-картирования на основе фазы сигнала

Сдвиг Блоха — Зигерта

Метод, основанный на сдвиге Блоха — Зигерта (Bloch–Siegert shift, BSS), относится к фазовым методам: информация об амплитуде B1+ кодируется не в модуле, а в фазе сигнала. В последовательность вводят внерезонансный радиочастотный импульс, частота которого смещена относительно резонанса воды; такой импульс вызывает дополнительный фазовый сдвиг ядер, пропорциональный квадрату амплитуды B1+. Выполнение двух измерений с импульсом, смещённым на одинаковую частоту в противоположные стороны от резонанса, позволяет вычесть вклады неоднородности основного поля B0 и химического сдвига. В оригинальной работе показано, что получение карты с затратой порядка 25 секунд на срез количественно сопоставимо с двухугловой картой, требующей около 21 минуты, что и определило быстрое распространение метода.

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

DREAM и родственные методики

Метод DREAM (Dual Refocusing Echo Acquisition Mode) использует стимулированное эхо и сигнал спада свободной индукции, получаемые за один цикл повторения; отношение их амплитуд даёт фактический угол поворота, а получение данных за один цикл повторения обеспечивает очень короткое время измерения и малую чувствительность к движению. Метод широко применяется в исследованиях при 7 Тл, в том числе в вариантах с фазовым кодированием передающих каналов (PE-DREAM), позволяющих получать комплексные карты B1+ каждого канала массива и использовать их для расчёта kT-точечных импульсов, например при визуализации печени. В 2025 году предложен спиральный трёхмерный вариант (spiral 3DREAM), заменяющий декартово считывание ускоренным считыванием «стопкой спиралей» со схемой выборки CAIPIRINHA; по результатам фантомных и in vivo экспериментов он показал хорошее соответствие эталонным методам при уменьшении размытия изображений и карт угла поворота по сравнению с декартовым вариантом.

Родственным направлением является фазочувствительный метод, предложенный Моррелом (2008), а также ряд гибридных схем, сочетающих преимущества амплитудного и фазового кодирования. В литературе по картированию B1 принято подчёркивать, что ни одна из методик не является универсально оптимальной: выбор определяется анатомической областью, ожидаемым динамическим диапазоном поля, допустимыми SAR и временем сканирования.

Сравнительная характеристика основных методик B1-картирования

Картирование B1-поля в магнитно-резонансной томографии

Метод	Принцип кодирования	Основные преимущества	Основные ограничения
Двойной угол (DAM)	Амплитуда; отношение сигналов при α и 2α Амплитудный	Простота, возможность реализации на стандартных последовательностях	Длительные TR, продолжительное сканирование, чувствительность к профилю среза и B₀
Фактический угол поворота (AFI)	Амплитуда; два TR в стационарном режиме Амплитудный	3D-покрытие, хорошее отношение сигнал/шум, эталон для проверки	Время порядка минут, чувствительность к спойлингу, ограничения по SAR
Насыщающий TurboFLASH (satTFL)	Амплитуда; насыщение + быстрое считывание Амплитудный	Широкий динамический диапазон, быстрое получение	Требует аккуратной настройки насыщающего импульса
Сдвиг Блоха — Зигерта (BSS)	Фаза; внерезонансный импульс Фазовый	Высокая точность, устойчивость к B₀	Нагрев тканей (SAR), относительно длительное сканирование при большом объёме
DREAM / PE-DREAM / spiral 3DREAM	Амплитуда; стимулированное эхо и FID за один TR Амплитудный	Сверхбыстрое получение, малая чувствительность к движению	Ограниченный динамический диапазон, чувствительность к B₀ и T₂*
Фингерпринтинг (B1-MRF)	Подгонка к словарю сигналов Реконструктивный	Расширенный динамический диапазон в теле при 7 Тл	Сложность реконструкции, потребность в словаре и валидации

Современные инновационные подходы

Машинное обучение и вероятностные субпространственные модели

Одним из наиболее динамично развивающихся направлений стало применение машинного обучения и вероятностных субпространственных моделей для оценки B1-полей. Метод основан на предварительном изучении пространственных паттернов вариаций B1-поля на обучающих данных и последующем использовании этой информации для оценки B1-полей в новых экспериментальных данных.

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

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

Нейросетевое предсказание карт B1+ по локализационным изображениям

Следующим этапом стало прямое предсказание карт B1+ с помощью глубоких нейронных сетей. В работе 2025 года (Krueger и соавт.) для головного мозга при 7 Тл с восьмиканальной передающей катушкой применены комплекснозначные нейронные сети, оценивающие поканальные карты B1+ по многосрезовым локализационным изображениям, полученным в режиме циркулярной поляризации (CP+), при различных ориентациях срезов; целью являлось ускорение индивидуальной калибровки pTx. Среднее относительное расхождение по модулю составило около 2,7 ± 2,9% для трансверзальных, 1,8 ± 0,7% для сагиттальных и 1,3 ± 0,3% для корональных срезов в ткани мозга, а средняя абсолютная разность фаз — порядка 6,7°, 4,3° и 2,7° соответственно; на двух независимых добровольцах предсказанные карты были использованы для расчёта динамических kT-точечных импульсов с достижением однородного возбуждения во всём мозге.

Принципиальный вопрос для клинического внедрения — переносимость моделей между центрами. Недавнее исследование кросс-центровой обобщаемости показало, что свёрточные нейронные сети позволяют предсказывать карты B1+ на разных системах 7 Тл при использовании идентичных pTx-катушек, несмотря на различия моделей сканеров и аппаратной части; предсказанные карты оказались пригодны для динамической параллельной передачи, что рассматривается как основа быстрого, переносимого и масштабируемого решения. Ограничения подхода очевидны: сети обучаются на конкретной конфигурации катушки и в популяции здоровых добровольцев, поэтому перенос на другие катушки, анатомические области и патологические изменения требует отдельной проверки.

Фингерпринтинг, спиральные и реконструктивные подходы

Расширяется арсенал методик, использующих нестандартные принципы кодирования и реконструкции. Для абдоминальных исследований при 7 Тл предложена методика B1-MRF, основанная на последовательности MR-фингерпринтинга типа FLASH с высокой чувствительностью к B1+, которая обеспечивает большой динамический диапазон и сопоставлена с эталонными методами (PEX, satTFL, AFI, BSS). Для оценки массивов параллельной передачи предложен метод, использующий эквивалентность задачи алгоритму параллельной реконструкции ESPIRiT: он позволяет вычислять комплексные карты B1+ отдельных элементов с более высоким отношением сигнал/шум и расширенным динамическим диапазоном, а абсолютизация осуществляется калибровкой по одной AFI-карте. Предложены и трёхмерные последовательности «Sandwich», позволяющие получать карты B1+ всего мозга и сердца с ЭКГ-синхронизацией за одну задержку дыхания.

Прогнозирование локального SAR

Отдельное направление, непосредственно связанное с B1-картированием, — оценка локального SAR по картам B1+. В работе 2023 года (Gokyar и соавт.) для головы при 7 Тл с параллельной передачей трёхмерная свёрточная нейронная сеть использует карту B1+, карту с обращённой фазой и структурное МР-изображение для предсказания распределения локального SAR (усреднение по 10 г) и пикового пространственного SAR; такой подход потенциально позволяет уменьшить чрезмерные запасы безопасности, которые приходится закладывать при отсутствии индивидуальной оценки.

Клинические применения B1-картирования

Количественное T1-картирование

Одним из наиболее важных применений B1-картирования является коррекция неоднородности радиочастотного поля в методах количественного T1-картирования. Неточности определения угла поворота приводят к систематическим ошибкам оценки времени продольной релаксации, что критично для диагностических применений. Особенно чувствительны к ним методы, основанные на переменном угле поворота (VFA), и широко применяемая в исследованиях последовательность MP2RAGE; её чувствительность к B1+ при 7 Тл непосредственно влияет на межцентровую воспроизводимость T1 и на морфометрию гиппокампа, что показано в работе Haast и соавт. Для MP2RAGE разработаны схемы совместного использования карты B1+ и карт T1, снижающие B1-зависимость (Marques, Gruetter, 2013).

Валидационные исследования демонстрируют, что коррекция B1-неоднородности с использованием точных B1-карт приводит к значительному улучшению качества T1-карт как у здоровых добровольцев, так и у пациентов с патологией мозга. Важно повышение точности T1-измерений в областях с выраженной B1-неоднородностью — височных долях и основании мозга. Клиническая значимость точного T1-картирования проявляется в возможности более надёжной дифференциации патологических тканей, мониторинга терапевтического ответа и количественной оценки тканевых изменений; при высоких магнитных полях использование B1-коррекции становится обязательным условием воспроизводимости результатов.

Картирование плотности протонов

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

Параллельная передача радиочастотных импульсов

Развитие технологий параллельной передачи (pTx) делает B1-картирование необходимым для калибровки многоканальных передающих систем. Знание пространственного распределения B1-поля каждого передающего элемента требуется для расчёта оптимальных фаз и амплитуд радиочастотных импульсов. Специализированные применения pTx, такие как исследования спинного мозга при 7 Тл, требуют адаптации методов B1-картирования к нетрадиционным анатомическим областям; оптимизация интерферометрического кодирования обеспечивает достаточную точность карт для эффективной работы pTx-систем.

Принципиально новым этапом стало создание универсальных импульсов (universal pulses, UP) — калибровочно-независимых pTx-импульсов, рассчитанных на основе базы данных B1- и B0-карт популяции добровольцев. Валидация на двух коммерческих радиочастотных катушках показала, что при использовании UP-режима отклонение фактического угла поворота от номинального в среднем по популяции здоровых взрослых не превышало 11%, тогда как в режиме циркулярной поляризации оно превышало 25%. Для нейровизуализации универсальные импульсы признаются надёжной и простой в использовании альтернативой, тогда как индивидуально оптимизируемые «онлайн» импульсы пока остаются предметом дискуссии относительно их стабильного преимущества; при переходе 7 Тл к исследованиям тела, как отмечают авторы обзоров, субъектно-специфическое проектирование импульсов может стать неизбежным.

Следующим шагом являются динамические pTx-импульсы с индивидуальной оптимизацией. В 2025 году показано, что для трёхмерных последовательностей типа SPACE при 7 Тл индивидуально оптимизированные динамические pTx-импульсы улучшают отношение сигнал/шум и однородность во всех контрастах (T1, T2, FLAIR, DIR), а коэффициент вариации интенсивности в мозжечке — зоне, наиболее подверженной неоднородности, — снижается на 46% для FLAIR-изображений. В том же году опубликованы результаты клинического внедрения динамической параллельной передачи при 7 Тл с улучшением однородности и контраста на SPACE-последовательностях, а в 2026 году предложена схема быстрого онлайн-проектирования нескольких субъектно-специфических импульсов на основе универсальных импульсов субпопуляций для 3D SPACE и FLAIR: на 120 тестовых субъектах получено улучшение моделируемой амплитуды эхо-сигнала на 18% в областях с выраженной неоднородностью B0 и на 13% в мозге. Совокупный клинический контекст 7 Тл-МРТ в нейрорадиологии, включая роль динамической pTx и реконструкции на основе глубокого обучения, обобщён в обзоре, опубликованном в журнале RadioGraphics в 2025 году.

Методы переноса насыщения и спектроскопия

Эффекты магнитизационного переноса и химического обмена (CEST) зависят от среднеквадратичной амплитуды насыщающего поля B1rms, поэтому неоднородность B1 прямо искажает количественные CEST-параметры. Для 7 Тл разработан калибровочно-независимый подход PUSHUP, в котором универсальные импульсы, рассчитанные по базе B1-карт, гомогенизируют B1rms насыщения за счёт чередования нескольких насыщающих импульсов с различными B1-шиммами; это позволяет выполнять CEST-визуализацию всего мозга с 45 частотами насыщения при разрешении 1,6 мм менее чем за 4 минуты без индивидуальной калибровки. Аналогичные соображения относятся и к спектроскопической визуализации, в которой B1-карта используется для коррекции метаболических карт.

Технические ограничения

Влияние артефактов движения

Одним из основных технических ограничений методов B1-картирования является чувствительность к артефактам движения, особенно выраженная в методах, требующих совмещения нескольких изображений. Движение пациента между измерениями приводит к ошибкам вычисления отношений интенсивностей сигнала, что критично для точности B1-карт. Проблема особенно актуальна для AFI и двухугловых методов, где необходимо точное пространственное соответствие между изображениями, полученными с разными параметрами. Современные подходы к решению включают регистрацию изображений в реальном времени, использование сверхбыстрых однократных методик типа DREAM и разработку более устойчивых к движению последовательностей; для pTx отдельно исследуется влияние движения головы на точность ранее измеренных B1-карт и методы его компенсации.

Ограничения удельного поглощения энергии

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

Чувствительность к неоднородности основного поля

Многие методы B1-картирования чувствительны к неоднородностям основного магнитного поля B0, что может приводить к систематическим ошибкам в областях с выраженным градиентом поля. Проблема особенно актуальна вблизи границ воздух — ткань (околоносовые пазухи, основание черепа) и металлических имплантатов. Дифференциальные методы, такие как сдвиг Блоха — Зигерта, частично решают её вычитанием фазовых вкладов B0-неоднородности, однако остаточная чувствительность к B0 может ограничивать точность измерений в сложных анатомических областях. Для областей с металлическими конструкциями разработаны отдельные методики B1+-картирования на основе магнитудных изображений турбо-спин-эхо (TSE), менее подверженные влиянию искажений, вызванных металлом.

Ограничения пространственного разрешения

Достижение высокого пространственного разрешения в B1-картировании ограничивается компромиссом между разрешением, отношением сигнал/шум и временем сканирования. Большинство клинических протоколов используют относительно низкое пространственное разрешение (порядка 3–5 мм), что вполне допустимо ввиду плавного характера пространственных вариаций B1, но может быть недостаточным для коррекции в мелких анатомических структурах и вблизи границ тканей с резким контрастом диэлектрических свойств. Разработка ускоренных методов сбора данных, включая параллельную визуализацию, compressed sensing и спиральные траектории считывания, открывает возможности повышения разрешения B1-карт без пропорционального увеличения времени сканирования.

Обработка и анализ B1-карт

Методы сглаживания и фильтрации

Сырые B1-карты часто содержат шум, артефакты и нефизиологические вариации, не отражающие истинного распределения радиочастотного поля. Для минимизации распространения ошибок при использовании B1-карт в других количественных методах широко применяются сглаживание и фильтрация. Пространственное сглаживание снижает влияние шума и артефактов, однако чрезмерное сглаживание приводит к потере информации о локальных вариациях, поэтому необходим баланс между подавлением шума и сохранением пространственной точности. Физическая природа B1-поля — решения уравнений Максвелла в неоднородной среде — позволяет использовать его гладкость как регуляризующее условие: в некоторых реализациях (в частности, описанных в патентной литературе) для подготовки данных обучения нейросетей измеренную карту аппроксимируют полиномиальной или гармонической функцией по маске ткани, что подавляет разрывы и шум.

Современные подходы включают адаптивные методы фильтрации, учитывающие локальные характеристики изображения и анатомические границы. Они позволяют сохранять резкие переходы в B1-поле на границах тканей при одновременном подавлении шума в однородных областях.

Обеспечение качества B1-карт — критически важный аспект их клинического применения. Стандартные протоколы валидации включают сравнение с референсными методами, оценку воспроизводимости и анализ физиологической правдоподобности результатов. Для сопоставления методик часто используют фантом с известными свойствами в качестве эталона, а смещение каждого метода оценивают как мультипликативный поправочный коэффициент по гистограмме отношений карт; сравнительные исследования показывают, что точность и прецизионность методов существенно различаются и зависят от допустимого времени сканирования.

Фантомные исследования остаются основным способом валидации новых методов B1-картирования, обеспечивая контролируемые условия для оценки точности и воспроизводимости, однако результаты на фантомах не всегда полностью переносимы на клинические условия из-за различий диэлектрических свойств. Разработка стандартизированных протоколов контроля качества, включающих автоматизированное обнаружение артефактов и отклонений от ожидаемых паттернов распределения поля, остаётся актуальной задачей.

Клиническая стандартизация

Широкое внедрение B1-картирования в клиническую практику требует стандартизированных протоколов и рекомендаций для различных клинических применений: определения минимальных требований к качеству B1-карт для разных диагностических задач и руководств по выбору метода для конкретных приложений. Классификация методов на амплитудные и фазовые, сопоставление их параметров и ограничений, а также открытые образовательные ресурсы по количественной МРТ, включая курсы сообщества qMRLab, способствуют унификации терминологии и подходов.

Гармонизация методов B1-картирования между производителями МР-оборудования критически важна для воспроизводимости результатов многоцентровых исследований; разработка универсальных последовательностей и алгоритмов обработки способствует более широкому принятию этих технологий. Образовательные программы для медицинского персонала и интуитивно понятные интерфейсы интерпретации B1-карт будут способствовать интеграции методов в рутинную клиническую практику.

Заключение

B1-картирование представляет собой фундаментальную технологию, необходимую для реализации полного потенциала современных количественных методов МРТ и систем ультравысокого поля. Разнообразие доступных методов — от классических подходов на основе анализа амплитуды сигнала (DAM, AFI, satTFL) и фазовых методик (сдвиг Блоха — Зигерта, DREAM) до подходов с использованием фингерпринтинга, реконструктивных алгоритмов и машинного обучения — обеспечивает возможность выбора оптимального решения для конкретных клинических и исследовательских задач.

Современное состояние области характеризуется активным преодолением традиционных ограничений: длительного времени сканирования (сверхбыстрые трёхмерные и спиральные реализации, нейросетевое предсказание карт по локализаторам), чувствительности к артефактам движения и ограничений по SAR (методы предсказания локального SAR, универсальные и динамические pTx-импульсы). Особое значение приобретает переход параллельной передачи в клиническую практику на системах 7 Тл, где индивидуальная или популяционная калибровка B1 стала условием получения диагностически однородных изображений.

Клиническая значимость точного B1-картирования особенно очевидна в контексте количественной МРТ, где даже небольшие ошибки определения угла поворота могут приводить к существенным диагностическим ошибкам. Продолжающееся развитие методов коррекции неоднородности B1 и их интеграция в рутинные клинические протоколы будут способствовать повышению точности диагностики и мониторинга терапевтического ответа.

Дальнейшее развитие области, вероятно, будет характеризоваться автоматизацией, повышением скорости измерений, межцентровой переносимостью обученных моделей и интеграцией с адаптивными системами управления радиочастотным возбуждением. Эти достижения приближают реализацию персонализированной МРТ, в которой параметры сканирования автоматически оптимизируются для каждого пациента и клинической задачи.

♥ Если вам понравился эта публикация, поделитесь ей с друзьями и коллегами!

Ключевые цитирования

Brunheim S, Gratz M, Johst S, Bitz AK, Fiedler TM, Ladd ME, Quick HH, Orzada S. Fast and accurate multi-channel B1+ mapping based on the TIAMO technique for 7T UHF body MRI. Magn Reson Med. 2018;79(5):2652–2664. DOI: 10.1002/mrm.26925.
Yarnykh VL, Kisel AA, Khodanovich MY. Quantitative Multi-Parameter Mapping Optimized for the Clinical Routine. Front Neurosci. 2020;14:611194. DOI: 10.3389/fnins.2020.611194.
Boudreau M, Tardif CL, Stikov N, Sled JG, Lee W. B1 mapping for bias-correction in quantitative T1 imaging of the brain at 3T using standard pulse sequences. J Magn Reson Imaging. 2017;46(6):1673–1682. DOI: 10.1002/jmri.25692.
Kawashima M, Ohno Y, Kose K. Fast B1 Mapping Based on Double-Angle Method with T1 Correction Using Standard Pulse Sequence. J Med Phys. 2022;47(1):93–98. DOI: 10.4103/jmp.jmp_78_21.

Источники

Insko EK, Bolinger L. Mapping of the radiofrequency field. J Magn Reson A. 1993;103:82–85.
Yarnykh VL. Actual flip-angle imaging in the pulsed steady state: a method for rapid three-dimensional mapping of the transmitted radiofrequency field. Magn Reson Med. 2007;57(1):192–200.
Chung S, Kim D, Breton E, Axel L. Rapid B1+ mapping using a preconditioning RF pulse with TurboFLASH readout. Magn Reson Med. 2010;64(2):439–446.
Sacolick LI, Wiesinger F, Hancu I, Vogel MW. B1 mapping by Bloch–Siegert shift. Magn Reson Med. 2010;63(5):1315–1322.
Lutti A, Hutton C, Finsterbusch J, Helms G, Weiskopf N. Optimization and validation of methods for mapping of the radiofrequency transmit field at 3T. Magn Reson Med. 2010;64(1):229–238.
Nehrke K, Börnert P. DREAM — a novel approach for robust, ultrafast, multislice B1 mapping. Magn Reson Med. 2012;68(5):1517–1526.
Wiesinger F, Van de Moortele P-F, Adriany G, De Zanche N, Ugurbil K, Pruessmann KP. Potential and feasibility of parallel MRI at high field. NMR Biomed. 2006;19:368–378.
Padormo F, Beqiri A, Hajnal JV, Malik SJ. Parallel transmission for ultrahigh-field imaging. NMR Biomed. 2016;29:1145–1161.
Gras V, Vignaud A, Amadon A, Le Bihan D, Boulant N. Universal pulses: a new concept for calibration-free parallel transmission. Magn Reson Med. 2017;77:635–643.
Improved Bloch–Siegert based B1 mapping by reducing off-resonance shift. NMR Biomed. 2013 (PMID 23355474).
Bloch–Siegert B1-mapping improves accuracy and precision of longitudinal relaxation measurements in the breast at 3 T. PMC5201175.
Haast RAM, Lau JC, Ivanov D, Menon RS, Uludağ K, Khan AR. Effects of MP2RAGE B1+ sensitivity on inter-site T1 reproducibility and hippocampal morphometry at 7T. NeuroImage. 2021;224.
Lutz M, et al. B1-MRF: Large dynamic range MRF-based absolute B1+ mapping in the human body at 7T. Magn Reson Med. 2024. DOI: 10.1002/mrm.30242.
Niesen et al. Spiral 3DREAM sequence for fast whole-brain B1 mapping. Magn Reson Med. 2025. DOI: 10.1002/mrm.30282.
Felder et al. High dynamic range B1+ mapping for the evaluation of parallel transmit arrays. Magn Reson Med. 2025. DOI: 10.1002/mrm.30349.
Krueger F, Aigner CS, Lutz M, et al. Deep learning-based whole-brain B1+-mapping at 7T. Magn Reson Med. 2025;93(4):1700–1711. DOI: 10.1002/mrm.30359.
Cross-site generalization of CNN-based B1+ mapping in UHF MRI. PMC13051333.
Gokyar S, et al. Deep learning-based local SAR prediction using B1 maps and structural MRI of the head for parallel transmission at 7 T. Magn Reson Med. 2023. DOI: 10.1002/mrm.29797.
Piredda GF, Sleight E, Yu T, et al. Individually optimized dynamic parallel transmit pulses for 3D high-resolution SPACE imaging at 7T. Magn Reson Med. 2025;94:1616–1625.
Middlebrooks EH, Herrler J, Piredda GF, et al. Clinical implementation of dynamic parallel transmission in 7T brain MRI: improved homogeneity and contrast using SPACE sequences. AJNR Am J Neuroradiol. 2025.
Fast online 3D SPACE and FLAIR imaging at 7T using multiple subject-specific parallel transmission pulses based on subpopulation universal pulses. Magn Reson Med. 2026 (PMID 42003222).
Völzke Y, Akbey S, Löwen D, et al. Calibration-free whole-brain CEST imaging at 7T with parallel transmit pulse design for saturation homogeneity utilizing universal pulses (PUSHUP). Magn Reson Med. 2025;93:630–642.
Ultrahigh-field 7-T MRI in neuroradiology: a comprehensive clinical review. RadioGraphics. 2025 (DOI 10.1148/rg.250194).
Ready for Routine. Investigative Radiology. DOI: 10.1097/RLI.0000000000001252.
qMRLab. B1 mapping. Quantitative MRI mOOC (qmrlab.org/mooc/b1-mapping).
Marques JP, Gruetter R. New developments and applications of the MP2RAGE sequence — focusing the contrast and high spatial resolution R1 mapping. PLoS One. 2013;8(7):e69294.
Morrell GR. A phase-sensitive method of flip angle mapping. Magn Reson Med. 2008;60(4):889–894.