from os.path import abspath
from PARRECread import PARRECread
from skimage.restoration import unwrap_phase
import pyqtgraph as pg
alpha = -10.3e-9
gamma = 0.267513e9
B0 = 3
echoTime = 19.5e-3
#
sonalleveNativeFile = abspath(
r"C:\Users\plettlk\DCGAN_Image_Completion\19-PLEK-RLES-05\30W-sonication-wo-motion-artefact_2019-11-12_13-55-06_Native.PAR",
)
sonalleveNativeData,Params,Dims = PARRECread(sonalleveNativeFile)
phaseMapsOriginal = sonalleveNativeData[:,:,0,0,0,1,:].swapaxes(0,2)
MMaps = sonalleveNativeData[:,:,0,0,0,0,:].swapaxes(0,2)
pg.image(MMaps, title='mag')
##
phaseMaps = np.zeros_like(phaseMapsOriginal)
for i in range(9,27):
phaseMaps[i,:,:] = unwrap_phase(phaseMapsOriginal[i,:,:], seed=100)
phaseMapsCropped = phaseMaps[9:27,49:113, 50:114]
pg.image(phaseMapsCropped)
#
phaseMapsCropped_recentered_QA_heated = np.zeros_like(phaseMapsCropped)
#Tofu phantom data
sonalleveNativeFile1 = abspath(
r"C:\Users\plettlk\DCGAN_Image_Completion\19-PLEK-RLES-03\baseline-scan-1_2019-10-11_11-14-17_Native.PAR",
sonalleveNativeFile2 = abspath(
r"C:\Users\plettlk\DCGAN_Image_Completion\19-PLEK-RLES-03\baseline-scan-2_2019-10-11_11-46-22_Native.PAR",
)
sonalleveNativeFile3 = abspath(
r"C:\Users\plettlk\DCGAN_Image_Completion\19-PLEK-RLES-03\baseline-scan-3_2019-10-11_12-20-01_Native.PAR",
)
# read temperature data from PARREC files
sonalleveNativeData1,Params,Dims = PARRECread(sonalleveNativeFile1)
phaseMaps1 = sonalleveNativeData1[:,:,0,0,0,1,:].swapaxes(0,2)
MMaps1 = sonalleveNativeData1[:,:,0,0,0,0,:].swapaxes(0,2)
phaseMaps1 = phaseMaps1[0:203,:,:]
MMaps1 = MMaps1[0:203,:,:]
sonalleveNativeData2,Params,Dims = PARRECread(sonalleveNativeFile2)
phaseMaps2= sonalleveNativeData2[:,:,0,0,0,1,:].swapaxes(0,2)
MMaps2 = sonalleveNativeData2[:,:,0,0,0,0,:].swapaxes(0,2)
phaseMaps2 = phaseMaps2[0:709,:,:]