КАРТИРОВАНИЕ ГРАДИЕНТА ВОСПРИИМЧИВОСТИ

1. Способ определения карты градиента восприимчивости; способ включает получение (300) массива магнитно-резонансных данных о k-пространстве, указанный массив данных содержит данные о магнитно-резонансном эхе, и трансформацию (306) полученного массива данных о k-пространстве в массив данных изображения, где для заданного направления k-пространства способ включает
a) осуществление одномерного дискретного преобразования Фурье массива данных изображения вдоль заданного направления k-пространства; указанное одномерное преобразование Фурье ведет к набору членов Фурье, каждый член Фурье из указанного набора ассоциирован с указателем k-пространства вдоль заданного направления k-пространства,
b) выбор (314) значения усечения, непрерывно меняющегося между -N/2 и N/2 вдоль указанного направления k-пространства, где N представляет собой размер массива данных о k-пространстве в указанном направлении k-пространства,
c) рекурсивную генерацию (310) обновленного массива данных изображения посредством определения члена Фурье из набора членов Фурье, с которым ассоциирован числовой указатель k-пространства, соответствующий указанному значению усечения, и вычитание определяемого члена Фурье из массива данных изображения, генерируемого на непосредственно предшествующей стадии генерации массива данных изображения; указанное вычитание ведет к новому обновленному массиву данных изображения,
d) вычисление (312) величины нового обновленного массива данных изображения,
e) повторение стадий с b) до d) с увеличением или уменьшением значения усечения,
f) вычисление (322) для заданного направления k-пространства эхо-сдвига по изменениям величин между новыми обновленными массивами данных изображения,
g) вычисление (326) градиента восприимчивости вдоль заданного направления k-пространства по указанному эхо-сдвигу.
2. Способ по п.1, в котором
- в первом цикле на стадии b) значение усечения непрерывно увеличивают от нуля до половины общего числа указателей k-пространства массива данных о k-пространстве вдоль указанного направления k-пространства, что ведет на стадии c) к первым новым обновленным массивам данных изображения и на стадии d) - к первым величинам,
- во втором цикле на стадии b) значения усечения непрерывно уменьшают от нуля до минус половины общего числа указателей k-пространства массива данных о k-пространстве вдоль указанного направления k-пространства, что ведет на стадии c) ко вторым новым обновленным массивам данных изображения и на стадии d) - ко вторым величинам,
- на стадии f) эхо-сдвиг вычисляют по разности между суммой всех вторых величин всех вторых новых обновленных массивов данных изображения и суммой всех первых величин всех первых новых обновленных массивов данных изображения.
3. Способ по п.2, в котором после каждой отдельной стадии вычисления первой или второй величины указанную, фактически вычисленную первую или вторую величину прибавляют (312) к уже вычисленной сумме всех первых или вторых величин, что ведет к новой, уже вычисленной сумме всех первых или вторых величин.
4. Способ по п.1, который дополнительно содержит повторение стадий с a) до g) для всех остальных направлений k-пространства.
5. Способ по п.1, в котором направление k-пространства представляет собой одно из направления градиента считывания, направления фазового кодирования или направления кодирования среза.
6. Способ по п.1, в котором трехмерные и одномерные дискретные преобразования Фурье представляют собой быстрые преобразования Фурье.
7. Способ по п.1, в котором в случае осуществления стадии c) в первый раз способ включает генерацию (302) начального массива данных изображения посредством трехмерного преобразования Фурье полученного массива данных о k-пространстве, где в этом случае на стадии c) предыдущую стадию генерации массива данных изображения задают посредством генерации указанного начального массива данных изображения.
8. Способ по п.1, в котором полученный массив данных содержит данные о магнитно-резонансном эхе для множества эхо с различными эхо-временами, где карту градиента восприимчивости определяют отдельно для каждого эхо-времени, что ведет к временной карте градиента восприимчивости, где способ дополнительно содержит вычисление общей карты градиента восприимчивости посредством комбинации всех определяемых временных карт градиентов восприимчивости.
9. Способ по п.8, в котором множество эхо представляет собой двойное эхо.
10. Способ по п.8, в котором регистрация массива магнитно-резонансных данных о k-пространстве для одного из множества эхо осуществляют с использованием ультракороткой эхо-временной последовательности.
11. Компьютерный программный продукт, который содержит исполняемые компьютером инструкции для осуществления способа по любому из предыдущих пп.1-10.
12. Устройство для определения карты градиента восприимчивости объекта посредством магнитно-резонансной визуализации, где устройство содержит систему регистрации данных для получения (300) массива магнитно-резонансных данных о k-пространстве, указанный массив данных содержит данные о магнитно-резонансном эхе, где устройство выполнено с возможностью преобразования (306) полученного массива данных о k-пространстве в массив данных изображения, где устройство дополнительно выполнено с возможностью
a) осуществления одномерного дискретного преобразования Фурье массива данных изображения вдоль заданного направления k-пространства; указанное одномерное преобразование Фурье ведет к набору членов Фурье, каждый член Фурье из указанного набора ассоциирован с указателем k-пространства вдоль заданного направления k-пространства,
b) выбора (314) значения усечения непрерывно между -N/2 и N/2 вдоль указанного направления k-пространства, где N представляет собой размер массива данных о k-пространстве в этом направлении,
c) рекурсивной генерации (310) обновленного массива данных изображения посредством определения члена Фурье из набора членов Фурье, с которым ассоциирован числовой указатель k-пространства, который соответствует указанному значению усечения, и вычитания определяемого члена Фурье из массива данных изображения, генерируемого на непосредственно предшествующей стадии генерации массива данных изображения; указанное вычитание ведет к новому обновленному массиву данных изображения,
d) вычисления (312) величины нового обновленного массива данных изображения,
e) повторения стадий с b) до d) с увеличением или уменьшением значения усечения,
f) вычисления (322) для заданного направления k-пространства эхо-сдвига по изменениям величин между новыми обновленными массивами данных изображения,
g) вычисления (326) градиента восприимчивости вдоль заданного направления k-пространства по указанному эхо-сдвигу.