====== Correction of B1 inhomogeneity artifacts in CEST images ====== {{ :b1_corr.png?direct&200|}} Here you find a B1-correction algorithm as published in **Windschuh et al. (2015), NMR Biomed., 28(5):529-37. [[http://onlinelibrary.wiley.com/doi/10.1002/nbm.3283/abstract|doi: 10.1002/nbm.3283]]**. Dowload zipped Matlab implementations [[https://github.com/cest-sources/B1_correction/archive/master.zip|here]] or find the package on [[https://github.com/cest-sources/B1_correction/]] ====== Tutorial ====== There are two functions, **Z_B1_correction** and **contrast_B1_correction**, which both require data of several CEST acquisitions with different B1. Both reconstruct either Z-spectrum stacks or contrast stacks at the B1 values given by **B1_output**. function [Z_stack_corr] = Z_B1_correction(Z_stack,rel_B1map,B1_input,B1_output,SEGMENT,fit_type,B1_input_index) % [Z_stack_corr] = Z_B1_correction(Z_stack,rel_B1map,SEGMENT_2D,B1_input,fit_type,B1_input_index,B1_output) % output: B1 corrected Z-stack (5D-stack) % input: Z_stack = 5D-stack of Z-spectra (y,x,z,offset,B1) % rel_B1map = relative B1map % B1_output = output B1 (scalar or vector) % SEGMENT = voxel/pixel mask for evaluation % fit_type = fit or interpolation type % B1_input_index = index vector of B1 samples that should be used % for correction % function [corr_img]=contrast_B1_correction(img,rel_B1map,B1_input,B1_output,SEGMENT_2D,fit_type,B1_input_index) % [corr_img]=contrast_B1_correction(img,rel_B1map,B1_input,B1_output,SEGMENT_2D,fit_type,B1_input_index) % output: B1 corrected images (2D or 3D-stack) (x,y,B1) % input: img = 3D-stack of images (x,y,B1) which are used for B1-correciton % rel_B1map = relative B1map % B1_output = output B1 (scalar or vector) % SEGMENT_2D = pixel mask for evaluation % fit_type = fit or interpolation type % B1_input_index = index vector of B1 samples that should be used % for correction {{ :logo_dkfz.png?direct&180|powered by DKFZ}}