From 569eba87754430827dcab097773d88919bd0639a Mon Sep 17 00:00:00 2001 From: AlmaBlomback <63227038+AlmaBlomback@users.noreply.github.com> Date: Fri, 4 Oct 2024 15:50:15 +0200 Subject: [PATCH] Update bsw_leakagecorr.py (#133) --- .../LeakageCorrection/bsw_leakagecorr.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/original/SR_TBG_BNIPhoenix_USA/LeakageCorrection/bsw_leakagecorr.py b/src/original/SR_TBG_BNIPhoenix_USA/LeakageCorrection/bsw_leakagecorr.py index aefe0360..12500371 100644 --- a/src/original/SR_TBG_BNIPhoenix_USA/LeakageCorrection/bsw_leakagecorr.py +++ b/src/original/SR_TBG_BNIPhoenix_USA/LeakageCorrection/bsw_leakagecorr.py @@ -53,7 +53,7 @@ def BSWleakagecorr(dR2s, nonenhance_map): for x in range(nX*nY*nZ): print("Voxel "+str(x)) if np.all(~np.isnan(dR2s_vec[x,:])) and np.all(~np.isinf(dR2s_vec[x,:])): - BSW_kvals,BSW_fitcov = curve_fit(BSWfunction,(dR2s_WBNE[x,:],dR2s_WBNE_integral),dR2s_vec[x,:],p0) + BSW_kvals,BSW_fitcov = curve_fit(BSWfunction,(dR2s_WBNE_avg,dR2s_WBNE_integral),dR2s_vec[x,:],p0) K1_dR2s[x] = BSW_kvals[0] K2_dR2s[x] = BSW_kvals[1] leakage = (-K2_dR2s[x]) * dR2s_WBNE_integral @@ -61,4 +61,4 @@ def BSWleakagecorr(dR2s, nonenhance_map): dR2s_leakagecorrected = np.zeros(dR2s.shape) dR2s_leakagecorrected = np.reshape(dR2s_BSW,dR2s.shape,order='F') - return dR2s_leakagecorrected \ No newline at end of file + return dR2s_leakagecorrected