T-test cutoff

Hi all,

I am currently working with a data set that will be displayed on a map. After I calculated the t-value I decided that It will be better to visualize only the values that are on the above my 90% or 2.32. How can I make this cutoff on my code without using if statement? Thanks in advance.

figure (3)
h= pcolor(lon_gridded,lat_gridded,tstat);
set (h, 'Edgecolor', 'none');
hold on
plot (coastlines(:,1), coastlines(:,2))
shading flat 
colormap(jet(9))
%caxis([bottom top]);
xlim([lon_start lon_end]);ylim([lat_start lat_end])

Answer 1

[Asvin Kumar]

Try adding the following lines before plotting the values:

tstatMod = tstat; 
tstatMod(tstat<=2.32) = NaN; 

And then modify your pcolor command to:

h = pcolor(lon_gridded,lat_gridded,tstatMod); 

The following articles might help you understand Logical Indexing better:

https://blogs.mathworks.com/steve/2008/01/28/logical-indexing/

https://www.mathworks.com/company/newsletters/articles/matrix-indexing-in-matlab.html

https://blogs.mathworks.com/loren/2013/02/20/logical-indexing-multiple-conditions/