Make units consistent (metres)

ComputeMostHexagonalSampling.m

@@ -8,14 +8,12 @@ function [packingDensity, L_optimum, isDegenerate, L_CAIPI_A, packingDensityCAIP
 % IN
 %   R_inplane               in-plane undersampling factor
 %
-%   FOV_inplane             in-plane FOV
-%                           [same unit as sliceSpacing, sliceThickness & sliceGap]
+%   FOV_inplane             in-plane FOV [m]
 %
-%   sliceSpacing            distance between simultaneously excited slices
+%   sliceSpacing            distance between simultaneously excited slices [m]
 %                           (only used if either of the other input arguments
 %                           {nSlicesSimultaneously, nSlicesTotal, sliceThickness, sliceGap}
 %                           unavailable)
-%                           [same unit as FOV_inplane, sliceThickness & sliceGap]
 %
 %   nSlicesSimultaneously   number of simultaneously excited slices
 %                           (only needed if sliceSpacing empty or blipped-CAIPIRINHA options are required)
@@ -23,13 +21,11 @@ function [packingDensity, L_optimum, isDegenerate, L_CAIPI_A, packingDensityCAIP
 %   nSlicesTotal            total number of slices
 %                           (only needed if sliceSpacing empty)
 %
-%   sliceThickness          slice thickness
+%   sliceThickness          slice thickness [m]
 %                           (only needed if sliceSpacing empty)
-%                           [same unit as FOV_inplane, sliceSpacing & sliceGap]
 %
-%   sliceGap                gap between slices
+%   sliceGap                gap between slices [m]
 %                           (only needed if sliceSpacing empty)
-%                           [same unit as FOV_inplane, sliceSpacing & sliceThickness]
 %
 %   doPlot                  Plot grid and packing density curve {false}
 %
@@ -63,8 +59,8 @@ function [packingDensity, L_optimum, isDegenerate, L_CAIPI_A, packingDensityCAIP
 %                           figure (fh{2})
 %
 %
-% Example: [packingDensity, L_optimum] = ComputeMostHexagonalSampling(1, 192, 19.2, 3);
-%          [packingDensity, L_optimum] = ComputeMostHexagonalSampling(1, 192, [], 3, 3, 4, 15.2)
+% Example: [packingDensity, L_optimum] = ComputeMostHexagonalSampling(1, 192e-3, 19.2e-3, 3);
+%          [packingDensity, L_optimum] = ComputeMostHexagonalSampling(1, 192e-3, [], 3, 3, 4e-3, 15.2e-3)
 %
 % Author: Maria Engel
 % (c) Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, United Kingdom
@@ -97,8 +93,8 @@ if nargout>7
 end
 
 % in-plane spacing of lines in k-space
-deltak_inplane = 2*pi*R_inplane/FOV_inplane;
-kzSlabThickness = 2*pi/sliceSpacing;
+deltak_inplane = 2*pi*R_inplane/FOV_inplane;    % in rad/m
+kzSlabThickness = 2*pi/sliceSpacing;            % in rad/m
 
 if nargin < 4
     nSlicesSimultaneously = ceil(kzSlabThickness*R_inplane/deltak_inplane);

MakeExampleGridsPlot.m

@@ -2,7 +2,7 @@ clear all;
 
 baseFolder = pwd;
 
-doSave = true;
+doSave = false;
 saveFolder = pwd;
 if doSave
     dfws = get(0,'DefaultFigureWindowStyle');
@@ -11,7 +11,7 @@ end
 
 Nsim = 5;
 R = 4.4/Nsim;
-FOVinplane = 220;
+FOVinplane = 0.22;
 
 sliceSpacing = FOVinplane/2/Nsim;
 

MakeGridAndAlpGif.m

@@ -13,7 +13,7 @@ end
 
 Nsim = 6;
 R = 1;
-FOVinplane = 220;
+FOVinplane = 0.22;
 
 sliceSpacing = FOVinplane/2/Nsim;
 

MakeOptDvsRPlot.m

@@ -8,9 +8,9 @@ Rtot_start = 1;
 Rtot_end = 12;
 Rtot_increment = 0.001;
 
-FOVinplane = 220;
+FOVinplane = 0.22;
 
-doSave = true;
+doSave = false;
 saveFolder = pwd;
 if doSave
     dfws = get(0,'DefaultFigureWindowStyle');