(MRI, Python) K-space에서 image로 바꾸는 방법

K-space data는 numpy라이브러리에서 제공하는 numpy array의 자료형으로 저장되어 있다.

 

이 K-space data인 숫자 행렬들로 부터 영상을 보기 위해 다음과 같은 작업을 한다.

 

import numpy as np
from numpy import fft
import matplotlib.pyplot as plt

path = './T2_7_C1M1_1.npy'

k_space = np.load(path)

unshifted_image = fft.fft2(k_space)

image = fft.fftshift(unshifted_image)

fig = plt.figure(figsize=(15, 6))

ax1 = plt.subplot(1,3,1);
ax1.set_title('K space')
plt.imshow(abs(k_space), cmap='gray')

ax2 = plt.subplot(1,3,2)
ax2.set_title('Fourier transform')
plt.imshow(abs(unshifted_image), cmap='gray')

ax3 = plt.subplot(1,3,3)
ax3.set_title('Fourier transform and Shifting')
plt.imshow(abs(image), cmap='gray')

이 결과는 다음과 같다.