How to plot performance diagram (Roebber, 2009)?

Hi, I am trying to plot the dichotomous statictics of rainfall using the performance diagram based on the study of Roebber (2009). Here's the link of the paper https://www.atmos.albany.edu/daes/atmclasses/atm401/spring_2015/Roebber2009.pdf. I would appreciate if someone help me figure this one out. This is the sample output of the plot.

Thank you

[Dyuman Joshi]

What exactly is your question and how does it relate to MATLAB?

Answer 1

[Angelo Yeo]

bias = [0.3, 0.5, 0.8, 1.0, 1.3, 1.5, 2.0, 3.0, 5.0, 10.0];
[SR, POD] = meshgrid(linspace(0, 1, 100));
CSI = 1./(1./SR + 1./POD - 1);

figure;
temp_x = linspace(0, 1, 10);
for i_bias = 1:length(bias)
    plot(temp_x, temp_x * bias(i_bias), 'k--');
    hold on;
end
xlim([0, 1]); ylim([0, 1]);
contour(SR, POD, CSI, "ShowText", true, "LevelList", 0.1:0.1:0.9, "LabelSpacing", 500)
axis square
xlabel("Success Ratio (1-FAR)")
ylabel("Probability of Detection (POD)")

text(1, 0.3, "0.3")
text(1, 0.5, "0.5")
text(1, 0.8, "0.8")
text(1, 1, "1")
text(1/1.3, 1, "1.3", "VerticalAlignment", "bottom")
text(1/1.5, 1, "1.5", "VerticalAlignment", "bottom")
text(1/2.0, 1, "2.0", "VerticalAlignment", "bottom")
text(1/3.0, 1, "3.0", "VerticalAlignment", "bottom")
text(1/5.0, 1, "5.0", "VerticalAlignment", "bottom")
text(1/10.0, 1, "10.0", "VerticalAlignment", "bottom")
grid on;