Skip to content
Snippets Groups Projects
Commit 68cc33a6 authored by Tamas Borbath's avatar Tamas Borbath
Browse files

Minor visualization changes in DF_processing. lineBroadeningExchangingPeaks.m...

Minor visualization changes in DF_processing. lineBroadeningExchangingPeaks.m and plot_mean_std_of_spectra_DF.m
parent 953b56db
No related branches found
No related tags found
No related merge requests found
...@@ -89,10 +89,14 @@ plotPpm = true; ...@@ -89,10 +89,14 @@ plotPpm = true;
legend('DF_{5.75} k=6.76s^{-1}','DF_{6.80} k=2.346s^{-1}','DF_{8.20} k=9.32s^{-1}','DF_{8.30} k=13.8s^{-1}','DF_{8.50} k=3.31s^{-1}', 'NAA k=0.74s^{-1}'); legend('DF_{5.75} k=6.76s^{-1}','DF_{6.80} k=2.346s^{-1}','DF_{8.20} k=9.32s^{-1}','DF_{8.30} k=13.8s^{-1}','DF_{8.50} k=3.31s^{-1}', 'NAA k=0.74s^{-1}');
if plotPpm if plotPpm
xlim([5.7 8.6]) xlim([5.7 8.6])
xlim([4.6 8.6])
else else
xlim([400 1600]) xlim([400 1600])
end end
ylim([-0.001 0.1]) ylim([-0.001 0.25])
set(gca,'xdir','reverse')
xlabel('\delta (ppm)')
ylabel('Signal (arb. u.)')
% results of FWHM % results of FWHM
%DF58 = 2.17 Hz %DF58 = 2.17 Hz
......
...@@ -10,6 +10,7 @@ TEs = [24; 32; 40; 52; 60]; ...@@ -10,6 +10,7 @@ TEs = [24; 32; 40; 52; 60];
colors = {[0 0.7 0]; [1 0 0]; [0 0 1]; [0 0 0]; [0 0.5 0.5];}; colors = {[0 0.7 0]; [1 0 0]; [0 0 1]; [0 0 0]; [0 0.5 0.5];};
numberOfTEs = length(TEs); numberOfTEs = length(TEs);
offsetStep = 2.7e-5; offsetStep = 2.7e-5;
offsetStep = 1.e-5;
offSetBase = offsetStep * (numberOfTEs + 1); offSetBase = offsetStep * (numberOfTEs + 1);
for indexTE = 1:numberOfTEs for indexTE = 1:numberOfTEs
...@@ -99,14 +100,14 @@ end ...@@ -99,14 +100,14 @@ end
%% %%
figure(1) figure(1)
stdAlpha = 0.4; % stdAlpha = 0.4;
hold on % hold on
h = area(ppmScale,[(meanSpectrum - stdSpectrum) + offset (stdSpectrum * 2)]); % h = area(ppmScale,[(meanSpectrum - stdSpectrum) + offset (stdSpectrum * 2)]);
h(1).FaceColor = 'none'; % h(1).FaceColor = 'none';
h(2).FaceColor = colorSpectrum; % h(2).FaceColor = colorSpectrum;
h(2).FaceAlpha = stdAlpha; % h(2).FaceAlpha = stdAlpha;
h(1).EdgeColor = 'none'; % h(1).EdgeColor = 'none';
h(2).EdgeColor = 'none'; % h(2).EdgeColor = 'none';
hold on hold on
pSpectrum = plot(ppmScale,meanSpectrum + offset, 'LineWidth',2, 'Color', colorSpectrum); pSpectrum = plot(ppmScale,meanSpectrum + offset, 'LineWidth',2, 'Color', colorSpectrum);
...@@ -114,7 +115,8 @@ xlim([5.5 9.5]); ...@@ -114,7 +115,8 @@ xlim([5.5 9.5]);
set(gca,'xDir','reverse') set(gca,'xDir','reverse')
set(gca,'LineWidth',1); set(gca,'LineWidth',1);
set(gca,'ytick',[]) set(gca,'ytick',[])
xlabel('[ppm]'); xlabel('\delta (ppm)');
ylabel('Signal (arb.u.)');
title('Spectra mean and standard deviation') title('Spectra mean and standard deviation')
% title(['Spectra TE ' num2str(TE) ' ms']) % title(['Spectra TE ' num2str(TE) ' ms'])
end end
\ No newline at end of file
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment