Компенсация плотности дискретизации с использованием быстрого преобразования Фурье деконволюции

Основная информация

• ID статьи: 2510.14873
• Название: Sampling Density Compensation using Fast Fourier Deconvolution
• Авторы: Rui Luo, Peng Hu, Haikun Qi (ShanghaiTech University)
• Классификация: physics.med-ph
• Дата публикации: 16 октября 2025 г. (препринт arXiv)
• Ссылка на статью: https://arxiv.org/abs/2510.14873

Аннотация

Функции компенсации плотности (DCF) широко применяются при реконструкции некартезианской МРТ, как при прямой реконструкции с использованием неравномерного быстрого преобразования Фурье (NUFFT), так и при итеративной реконструкции недодискретизированных данных. Современные передовые методы требуют затратных десятков итераций, что является одним из основных препятствий для широкого применения эффективной некартезианской МРТ. В данной работе предложен эффективный неитеративный метод расчета DCF для произвольных некартезианских k-пространственных траекторий с использованием быстрого преобразования Фурье деконволюции. Результаты моделирования показывают, что метод может генерировать DCF для трёхмерной некартезианской реконструкции примерно за 20 секунд, достигая ускорения на порядок величины по сравнению с передовыми методами при сохранении аналогичного качества реконструкции.

Исследовательский контекст и мотивация

Определение проблемы

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

Значимость проблемы

1. Качество реконструкции: DCF критически важна как для прямой NUFFT-реконструкции, так и для итеративной реконструкции
2. Ускорение сходимости: DCF может ускорить сходимость итеративной реконструкции и реконструкции на основе глубокого обучения путём улучшения числа обусловленности задачи
3. Клиническое применение: Эффективный расчет DCF является ключевым узким местом для клинического применения некартезианской МРТ

Ограничения существующих методов

1. Метод диаграммы Вороного: Дорогостоящий в вычислительном отношении и численно нестабильный, особенно для трёхмерных траекторий
2. Итеративные методы: Классический итеративный метод, предложенный Pipe и Menon, требует десятков итераций, каждая из которых занимает от нескольких секунд до нескольких минут
3. Методы оптимизации: Хотя более точные, время выполнения увеличивается на два порядка величины по сравнению с итеративными методами

Исследовательская мотивация

Разработка быстрого неитеративного метода расчета DCF, особенно для трёхмерных некартезианских схем дискретизации, с целью преодоления вычислительного узкого места.

Основные вклады

1. Предложен неитеративный метод расчета DCF на основе быстрого преобразования Фурье деконволюции (FFD)
2. Достигнуто ускорение на 1-2 порядка величины, время расчета DCF для трёхмерных траекторий сокращено с примерно 10 минут до менее чем 20 секунд
3. Сохранено или незначительно улучшено качество реконструкции при одновременном создании более гладких DCF
4. Предоставлено универсальное решение, применимое к произвольным некартезианским k-пространственным траекториям
5. Реализована открытая версия, способствующая воспроизводимости метода и его широкому применению

Подробное описание метода

Определение задачи

Для заданной схемы дискретизации некартезианского k-пространства $K = \{k_i\}_{i=1}^{N_k}$ найти функцию компенсации плотности $D(k)$ таким образом, чтобы функция рассеяния точки (PSF) взвешенной схемы дискретизации приблизительно равнялась импульсной функции в пределах поля зрения: $P(x) \approx \delta(x)$ для $\|x\| < L$.

Основная теоретическая база

1. Моделирование задачи

Процесс дискретизации представляется как: $S_1(k) = \text{III}(k) \cdot S_0(k)$

где $\text{III}(k) = \sum_{i=1}^{N_k} \delta(k - k_i)$ — последовательность импульсов.

Взвешенное k-пространство: $S_2(k) = D(k) \cdot \text{III}(k) \cdot S_0(k)$

Взвешенная схема дискретизации (WSP): $E(k) = D(k) \cdot \text{III}(k)$

Соответствующая PSF: $P(x) = \mathcal{F}^{-1}\{E(k)\}$

2. Стратегия разложения PSF

Начальная оценка PSF $\hat{P}(x)$ разлагается на: $\hat{P}(x) = \hat{P}_{\text{in}}(x) + \hat{P}_{\text{out}}(x)$

где:

• $\hat{P}_{\text{in}}(x) = \hat{P}(x) \cdot W(x)$
• $\hat{P}_{\text{out}}(x) = \hat{P}(x) \cdot (1-W(x))$

$W(x)$ — оконная функция, равная нулю при $\|x\| \geq L$.

3. Оптимальная взвешенная схема дискретизации

Оптимальная WSP получается путём деконволюции: $E^*(k) = \hat{E}(k)/\hat{E}_{\text{in}}(k)$

Это обеспечивает:

• $P^*_{\text{in}}(x) = \mathcal{F}^{-1}\{1\} = \delta(x)$
• $P^*_{\text{out}}(x) = \mathcal{F}^{-1}\{\hat{E}_{\text{out}}(k)/\hat{E}_{\text{in}}(k)\}$

Технические инновации

1. Оптимизация оконной функции

Используется параметризованная оконная функция $W(x) = 1 - \|\bar{x}\|^p$, где $\bar{x} = x/L$, с определением оптимального параметра формы посредством поиска минимакса: $p^* = \arg\min_p \{\max_{i_{\text{test}}} \|P^*_{\text{out}}(x)\|/P^*_0\}$

Оптимальный параметр определён методом Монте-Карло: $p^* = 2.4$.

2. Обеспечение численной стабильности

Используется одномерная DCF в качестве начального приближения: $\hat{D}(k_i) = \|k_{i+1} - k_i\|_2 \cdot \|k_i\|_2^{N_d-1}$

где $N_d$ — размерность k-пространства.

3. Быстрое преобразование Фурье деконволюции

Основной расчет реализуется через FFD, избегая итеративного процесса и напрямую решая оптимальную DCF.

Экспериментальная установка

Набор данных

Используются комплексные двумерные/трёхмерные цифровые фантомы, включающие:

• Эллиптические оболочечные структуры
• Сердцевидные структуры
• Сферы различных размеров
• Фазовые карты, созданные белым шумом и пространственной низкочастотной фильтрацией
• Размер матрицы: 256×256×256
• Поле зрения: 500 мм

Траектории дискретизации

Тестируются четыре типа некартезианских траекторий:

1. Двумерные траектории: Variable Density Spiral (VdSpiral), Rosette
2. Трёхмерные траектории: Cones, Yarnball

Методы сравнения

В качестве эталона используется передовой метод компенсации плотности дискретизации для трёхмерной дискретизации Zwart и др., объединяющий:

• Базовую итеративную структуру Pipe
• Оптимальную функцию ядра Johnson
• Эффективный метод сеточной свёртки

Показатели оценки

1. Качество реконструкции:
   • Нормализованная среднеквадратичная ошибка (NRMSE)
   • Индекс структурного сходства (SSIM)
2. Вычислительная эффективность: время выполнения $T_{\text{exe}}$
3. Качество PSF: полная ширина на половине максимума (FWHM)

Детали реализации

• Язык программирования: Python 3.12.8
• Библиотека БПФ: FINUFFT
• Оборудование: 12-ядерный процессор 4.9 ГГц (Intel® Core™ i7-12700)
• Нормализация с нулевым средним и единичной дисперсией перед реконструкцией

Результаты экспериментов

Основные результаты

1. Сравнение скорости вычисления

2. Сравнение качества реконструкции

Ключевые находки

1. Анализ качества DCF

• DCF, созданная предложенным методом, более гладкая, в то время как эталонный метод демонстрирует значительные колебания
• Обе методики имеют одинаковую FWHM PSF (1.5× размер пикселя), что указывает на эквивалентное пространственное разрешение

2. Качество реконструированного изображения

• Реконструированные изображения не содержат явных искажений и размытия
• Присутствуют только незначительные артефакты типа колец Гиббса, вызванные усечением k-пространства
• Предложенный метод в большинстве случаев достигает лучших показателей NRMSE и SSIM

3. Прорыв в вычислительной эффективности

• Время расчета DCF для трёхмерных траекторий сокращено с примерно 10 минут до менее чем 20 секунд
• Достигнуто ускорение на 1-2 порядка величины
• Обеспечена возможность эффективной трёхмерной некартезианской реконструкции

Связанные работы

Эволюция традиционных методов

1. Метод диаграммы Вороного (Rasche и др., 1999): Интуитивный, но вычислительно дорогостоящий
2. Итеративные методы (Pipe & Menon, 1999): Заложили основу современного расчета DCF
3. Оптимизация функции ядра (Johnson & Pipe, 2009): Улучшена точность реконструкции, но скорость остаётся низкой
4. Сеточная свёртка (Zwart и др., 2012): Повышена эффективность итераций, но по-прежнему требуется множество итераций

Позиционирование вклада данной работы

По сравнению с существующими методами, данная работа впервые реализует:

• Неитеративный расчет DCF
• Ускорение на порядок величины
• Сохранение или улучшение качества реконструкции
• Применимость к произвольным некартезианским траекториям

Заключение и обсуждение

Основные выводы

1. Прорыв в эффективности: Достигнуто ускорение расчета DCF на 1-2 порядка величины
2. Гарантия качества: Сохранено или незначительно улучшено качество реконструкции, созданы более гладкие DCF
3. Универсальность: Применимо к произвольным двумерным/трёхмерным некартезианским k-пространственным траекториям
4. Практичность: Время расчета трёхмерной DCF сокращено до менее чем 20 секунд, что соответствует требованиям клинического применения

Ограничения

1. Форма оконной функции: Текущая параметризованная форма $W(x) = 1 - \|\bar{x}\|^p$ может быть не оптимальной
2. Оптимизация параметров: Предположение об оптимальности параметров независимо от размерности и траектории может требовать дальнейшей проверки
3. Численная стабильность: Численная стабильность операции деконволюции зависит от качества начальной оценки DCF
4. Валидация на реальных данных: Метод проверен только на смоделированных данных, требуется дальнейшее тестирование на реальных МРТ-данных

Направления будущих исследований

1. Оптимизация оконной функции: Исследование более оптимальных форм оконной функции и адаптивного выбора параметров
2. Валидация на реальных данных: Проверка эффективности метода на реальных МРТ-данных
3. Параллельная оптимизация: Дальнейшая оптимизация параллельных вычислений для повышения скорости
4. Интеграция в приложения: Интеграция с современными МРТ-конвейерами реконструкции

Углубленная оценка

Преимущества

1. Значительный технологический прорыв: Впервые реализован неитеративный расчет DCF, решена давно существующая проблема вычислительной эффективности
2. Прочная теоретическая база: Теоретическая база, основанная на разложении PSF и деконволюции, разумна и инновационна
3. Тщательный экспериментальный дизайн: Тестирование на множестве траекторий, количественные показатели оценки, сравнение с передовыми методами
4. Высокая практическая ценность: Значительное ускорение обеспечивает методу сильный потенциал клинического применения
5. Вклад в открытый исходный код: Обещание открытого кода способствует распространению метода и воспроизводимости

Недостатки

1. Недостаточно глубокий теоретический анализ: Отсутствуют теоретические гарантии сходимости и оптимальности
2. Упрощённый выбор параметров: Предположение об независимости оптимальных параметров оконной функции от размерности и траектории может быть чрезмерно упрощённым
3. Отсутствие реальных данных: Валидация только на смоделированных данных, отсутствие тестирования на реальных МРТ-данных
4. Робастность к шуму: Недостаточное обсуждение робастности метода к шуму
5. Ограниченный диапазон сравнения: Основное сравнение с одним эталонным методом, следует рассмотреть больше вариантов сравнения

Влияние

1. Академическая ценность: Предоставляет новую теоретическую базу и практический метод для расчета DCF
2. Клиническое значение: Значительное повышение вычислительной эффективности может способствовать клиническому применению некартезианской МРТ
3. Распространение технологии: Открытая реализация будет способствовать широкому применению метода в МРТ-сообществе
4. Последующие исследования: Предоставляет новые идеи для дальнейших исследований в смежных областях

Сценарии применения

1. Клиническая МРТ: Требуется быстрый расчет DCF для реконструкции МРТ в реальном или квазиреальном времени
2. Исследовательские приложения: Обработка больших объёмов данных некартезианской МРТ и разработка алгоритмов
3. Трёхмерная визуализация: Особенно подходит для вычислительно интенсивной трёхмерной некартезианской реконструкции
4. Множество типов траекторий: Применимо к расчету DCF для различных схем некартезианской дискретизации

Библиография

В данной работе цитируются ключевые работы в области расчета DCF, включая новаторские работы Pipe & Menon, оптимизацию функции ядра Johnson & Pipe, метод сеточной свёртки Zwart и др., обеспечивающие прочную теоретическую базу и эталоны сравнения для исследования.