# -*- coding: utf-8 -*-
"""
Created on Mon Aug 19 10:07:26 2019

@author: txx
"""

import numpy as np
import matplotlib.pyplot as plt

ori_img1=np.zeros((218,170,24))
ori_img1=ori_img1.astype(np.complex64)

ori_img=np.load('D:/e14079s3_P09216.7_image_slice22.jpg.npy')

#ori_img1 =ori_img[:,:,0:12]+1j*ori_img[:,:,12:24] # 实部+虚部
for i in range(0,24,2):
    ori_img1[:,:,int(i/2)] =ori_img[:,:,i]+1j*ori_img[:,:,i+1] # 实部+虚部
    #fshift = np.fft.ifftshift(ori_img1)
    
    
plt.title('K space')
plt.imshow(np.abs(ori_img1[:,:,7]),cmap='gray')
plt.show()
    
image=np.fft.ifftn(ori_img1,axes=(-3,-2))
#image=np.fft.ifftn(ori_img1,axes=(-3,-2,-1))
plt.title('image') # 不能出现中文
plt.imshow(np.abs(image[:,:,7]),cmap = 'gray') # 显示复数图像
plt.show()

# 验证图像
v_image=np.fft.ifft2(ori_img1[:,:,7])
plt.title('fft2')
plt.imshow(np.abs(v_image),cmap='gray')
plt.show()