Master

ComputeMostHexagonalSampling.m

@@ -117,15 +117,12 @@ LArray = unique(sort([LArray L_CAIPI_B]));
 isBlippedCAIPI_A = ((ceil(LArray)==LArray) + (LArray<=nSlicesSimultaneously)) > 1;
 isBlippedCAIPI_B = ismember(LArray, L_CAIPI_B);
 
-colourCAIPI_A = [217, 83, 25]/255; % red
-colourCAIPI_B = [148, 0, 211]/255; % purple
-colourOptimum = [34, 139, 34]/255; % green
-
 howFar2Look = 10; % TODO: Find out what to put here!
 d_min = zeros(size(LArray));
 
 if doDebug
     figure('Name','Sampling grid and alps');
+    set(gcf, 'Position', [680   693   621   285])
     ah = subplot(1,2,1);
 end
 
@@ -144,19 +141,20 @@ for iL = 1:numel(LArray)
     if doDebug
         [kyArray, kzArray] = GetGridPoints(3*Lmax, deltak_inplane, sliceSpacing, LArray(iL));
         PlotGrid(kyArray, kzArray, sliceSpacing, ah);
-
+        title(sprintf('L = %.2f      %s = %.2f',LArray(iL), char(hex2dec('0394')), nSlicesSimultaneously/LArray(iL)));
+        
         subplot(1,2,2); hold on;
         plot(LArray(iL),packingDensityArray(iL),'.','Color','k');
 
         if isBlippedCAIPI_A(iL)
-            plot(LArray(iL),packingDensityArray(iL),'x','Color',colourCAIPI_A,'MarkerSize',10,'LineWidth',3);
+            plot(LArray(iL),packingDensityArray(iL),'x','Color',GetSamplingPatternColour(1),'MarkerSize',10,'LineWidth',3);
         elseif isBlippedCAIPI_B(iL)
-            plot(LArray(iL),packingDensityArray(iL),'x','Color',colourCAIPI_B,'MarkerSize',10,'LineWidth',3);
+            plot(LArray(iL),packingDensityArray(iL),'x','Color',GetSamplingPatternColour(2),'MarkerSize',10,'LineWidth',3);
         end
         
         xlabel('L'); ylabel('Packing density [%]');
-        xlim([Lmin Lmax]); ylim([min(packingDensityArray) 90.7]);
-        sgtitle(sprintf('L = %.2f',LArray(iL)));
+        xlim([Lmin Lmax]); ylim([0 90.7]);
+        axis square;
 
         %% Make GIF
         if ~isempty(gifName)
@@ -213,11 +211,11 @@ if doPlot
     plot(LArray, packingDensityArray, 'LineWidth', 3, 'Color', 'k', ...
         'DisplayName', 'All feasible sampling patterns');
     plot(L_CAIPI_A, packingDensityCAIPI_A, 'x', 'MarkerSize', 10, 'LineWidth', 3, ...
-        'Color', colourCAIPI_A, 'DisplayName', 'Blipped-CAIPIRINHA - A');
+        'Color', GetSamplingPatternColour(1), 'DisplayName', 'Blipped-CAIPIRINHA - A');
     plot(L_CAIPI_B, packingDensityCAIPI_B, 'x', 'MarkerSize', 10, 'LineWidth', 3, ...
-        'Color', colourCAIPI_B, 'DisplayName', 'Blipped-CAIPIRINHA - B');
+        'Color', GetSamplingPatternColour(2), 'DisplayName', 'Blipped-CAIPIRINHA - B');
     plot(L_optimum, packingDensity, 'x', 'MarkerSize', 10, 'LineWidth', 3, ...
-        'Color', colourOptimum, 'DisplayName', 'Optimum sampling');
+        'Color', GetSamplingPatternColour(0), 'DisplayName', 'Proposed method');
     legend;
     xlabel('L');
     ylabel('Packing density [%]');
@@ -244,6 +242,9 @@ scatter(ky, kz, 'MarkerEdgeColor', [0.5 0.5 0.5], 'MarkerFaceColor', [0.5 0.5 0.
 
 hold on;
 
+ind0 = find(ky==0);
+scatter(ky(ind0:ind0+1), kz(ind0:ind0+1), 'MarkerEdgeColor', [0 0 0], 'MarkerFaceColor', [0 0 0], 'LineWidth', 5);
+
 slabBorderColour = [100 149 237]/255; % light blue
 yline(0,'r','LineWidth',3, 'Color', slabBorderColour); 
 yline(2*pi/sliceSpacing, 'r', 'LineWidth', 3, 'Color', slabBorderColour);