modify call to B0 mapping scripts to use off-resonance estimation fro…

…m MRIFieldmaps.jl
BRAIN-TO · Jun 19, 2024 · 9e5db82 · 9e5db82
1 parent 34c9f0f
commit 9e5db82
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 2 deletions.
diff --git a/example/recon_config_joss_demo.jl b/example/recon_config_joss_demo.jl
@@ -18,7 +18,7 @@ params_general = Dict{Symbol,Any}();
 params_general[:gamma] = 42577478;
 
 # General options for reconstruction script:
-params_general[:do_load_maps] = true                 # if true, reloads B₀/SENSE maps instead of recalculating
+params_general[:do_load_maps] = false                 # if true, reloads B₀/SENSE maps instead of recalculating
 params_general[:do_save_recon] = true                # if true, saves reconstruction and all auxiliary image data (maps) as NIfTI files
 params_general[:do_plot_recon] = true                # if true, plots intermediate debugging and output recon figures (needs graphics, not recommended in multi-thread mode due to PyPlot)
 params_general[:do_process_map_scan] = true          # if true, compute sensitivity and B₀ maps from reconstructed Cartesian scan   

diff --git a/src/GIRFReco.jl b/src/GIRFReco.jl
@@ -24,6 +24,7 @@ using MRICoilSensitivities
 using RegularizedLeastSquares
 using Flux
 using MRIReco
+using MRIFieldmaps
 
 # using MosaicViews
 

diff --git a/src/utils/utils.jl b/src/utils/utils.jl
@@ -981,7 +981,7 @@ function run_cartesian_recon(params_general)
         b0_maps, times, out = b0map(
             reshape(cartesian_reco.data[:, :, :, :, 1, 1], (200, 200, num_slices,1, 2)),
             (TE1 / 1000.0, TE2 / 1000.0); 
-            smap = ones(size(b0_maps)),
+            smap = ones(ComplexF32,size(b0_maps,1),size(b0_maps,2),size(b0_maps,3),1),
             track =true,
             chat = true,
             chat_iter = 2