Merge branch 'master' into 'main'

Centre k-space slab

See merge request sapme1/optimalgridsampling!7

ComputeMostHexagonalSampling.m

@@ -238,6 +238,8 @@ if nargin < 4 || isempty(ah)
 else
     axes(ah);
 end
+kz = kz-pi/sliceSpacing;
+
 scatter(ky, kz, 'MarkerEdgeColor', [0.5 0.5 0.5], 'MarkerFaceColor', [0.5 0.5 0.5], 'LineWidth', 5);
 
 hold on;
@@ -246,14 +248,14 @@ 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);
+yline(-pi/sliceSpacing,'r','LineWidth',3, 'Color', slabBorderColour); 
+yline(pi/sliceSpacing, 'r', 'LineWidth', 3, 'Color', slabBorderColour);
 
 % Plot zigzag line to indicate sampling chronology
 plot(ky(1:ceil(numel(ky)/3)), kz(1:ceil(numel(ky)/3)), 'LineWidth', 3, 'Color', [253 218 13]/255); % yellow
 
 axis equal;
-ylim(2*pi/sliceSpacing*[-0.2 1.2]);
+ylim(2*pi/sliceSpacing*[-0.7 0.7]);
 xlim([min(ky) max(ky)]);
 xlabel('k_y [rad/m]'); ylabel('k_z [rad/m]');
 hold off;