Реконструкция самолетов МРТ ⇐ Python

[Гость]

In the following code i read a MRI scan with 3D images.The scans are axial slices. By executing the following code I try to obtain the sagittal and the coronal planes. My problem is that there is some much noise and I would like to know what I do wrong

from load_MRI_FLAIR import scan import matplotlib.pyplot as plt from scipy import ndimage import numpy as np print(scan.shape) #(22, 512, 512) z,y,x = scan.shape SLICE_NUM = np.random.randint(0,22) print("SLICE_NUM:",SLICE_NUM) coronal = ndimage.zoom(scan, (x/z,z/y,1), order=5) print("coronal.shape",coronal.shape) #coronal.shape (512, 22, 512) sagittal = ndimage.zoom(scan, (x/z,1,z/y), order=5) print("sagittal.shape",sagittal.shape) # sagittal.shape (512, 512, 22) scans =[scan[SLICE_NUM,:,:],coronal[:,SLICE_NUM,:],sagittal[:,:,SLICE_NUM]] plt.figure(figsize=(10,5)) for i in range(3): plt.subplot(1,3,i+1),plt.imshow(scans,cmap='gray') plt.show()




Источник: https://stackoverflow.com/questions/776 ... nstruction