У меня есть МР-изображение, и я выбрал данные вдоль спиральной траектории с помощью фазового картирования. Чтобы было понятно, алгоритм, который я использовал, был следующим:
(1) Получить местоположение (kx, ky) по спиральной траектории. Единицей измерения kx является выборка / метр или радиан / метр. (2) Рассчитайте фазовую карту (PM) с помощью exp (-j * 2 * pi * (kx x + ky y)), где (x, y) - местоположение изображения (в метрах). (3) Умножить изображение с PM (поэлементное умножение). Это умножение дает сдвинутое по фазе изображение (PMI). (4) Возьмите преобразование Фурье (fft2 в MATLAB) PMI. Это преобразует PMI в k-пространство сдвинутого по фазе изображения (KspPMI). (5) Соберите отсчеты по эталону (kx, ky), используемому для всех KspPMI.
Затем k-пространство проходит алгоритм построения сетки для реконструкции. Я получаю изображение обратно, но некоторые значения интенсивности и усиления усиливаются, а также усиление увеличивается с увеличением количества образцов, которые я беру. Может кто-нибудь сказать мне, почему это происходит?
Основной вопрос, который я здесь хотел задать, заключается в том, что мои спиральные данные в k-пространстве имеют гораздо большую амплитуду, чем предполагалось. Я пропускаю какие-либо необходимые нормализации?
Спасибо.