In this project, I show a couple of 3D vis techniques
[Ny,Nx,Nz] = size(data);
[x3d,y3d,z3d] = meshgrid(1:Nx,1:Ny,1:Nz);
s = data; % this if you want to normalize it -> 0.5*(data(:)-mean(data(:)))/std(data(:))+0.5;
s(s<0) = nan;
% s(s>1) = 1; % this is useful if you normalize it
hs = scatter3(x3d(:),y3d(:),z3d(:),0.5*s(:), ...