decimating data points for image maps, otherwise redrawing takes forever

1 files changed, 19 insertions, 6 deletions

diff --git a/MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m b/MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m
index 579925d..e107cb2 100644
--- a/MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m
+++ b/MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m
@@ -14,27 +14,40 @@ Epos = 1/sqrt(2)*(EpR + 1i*EpL);
 Eneg = 1/sqrt(2)*(EpR - 1i*EpL);
 Imor = abs(Epos).^2 - abs(Eneg).^2;
 
+%% Maps are not so useful with too many pixels in the image
+% Here we limit number of image points
+Npx = 800;
+Npy=600;
+
+Ndx = length(z_1);
+Ndy = length(t_1);
+
+skip_xp = max(1, floor(Ndx/Npx) );
+skip_yp = max(1, floor(Ndy/Npy) );
+ind_x = 1:skip_xp:Ndx;
+ind_y = 1:skip_yp:Ndy;
+
 
 %% field propagation
 z_1=z_1*100; % z in cm
 t_1=t_1*1e6; % time now measured in uS
 figure(1)
-subplot(2,2,1); imagesc(z_1, (t_1), IdL_out_1.'); colorbar
+subplot(2,2,1); imagesc(z_1(ind_x), t_1(ind_y), IdL_out_1(ind_x, ind_y).'); colorbar
 xlabel('z (cm)')
 ylabel('t (uS)')
 zlabel('I_{dL}')
 title('I_{dL}')
-subplot(2,2,3); imagesc(z_1, (t_1), IpL_out_1.'); colorbar
+subplot(2,2,3); imagesc(z_1(ind_x), t_1(ind_y), IpL_out_1(ind_x, ind_y).'); colorbar
 xlabel('z (cm)')
 ylabel('t (uS)')
 zlabel('I_{pL}')
 title('I_{pL}')
-subplot(2,2,2); imagesc(z_1, (t_1), IdR_out_1.'); colorbar
+subplot(2,2,2); imagesc(z_1(ind_x), t_1(ind_y), IdR_out_1(ind_x, ind_y).'); colorbar
 xlabel('z (cm)')
 ylabel('t (uS)')
 zlabel('I_{dR}')
 title('I_{dR}')
-subplot(2,2,4); imagesc(z_1, (t_1), IpR_out_1.'); colorbar
+subplot(2,2,4); imagesc(z_1(ind_x), t_1(ind_y), IpR_out_1(ind_x, ind_y).'); colorbar
 xlabel('z (cm)')
 ylabel('t (uS)')
 zlabel('I_{pR}')
@@ -114,14 +127,14 @@ legend('before', 'after')
 
 %%
 subplot(2,2,2);
-imagesc(z_1, (t_1), Ieit.'); colorbar
+imagesc(z_1(ind_x), t_1(ind_y), Ieit(ind_x, ind_y).'); colorbar
 xlabel('z (cm)')
 ylabel('t (uS)')
 zlabel('I_{eit}')
 title('I_{eit}')
 
 subplot(2,2,4);
-imagesc(z_1, (t_1), Imor.'); colorbar
+imagesc(z_1(ind_x), t_1(ind_y), Imor(ind_x, ind_y).'); colorbar
 xlabel('z (cm)')
 ylabel('t (uS)')
 zlabel('I_{mor}')