尝试绘制一种类似泰勒图降水评分图

qxsw2016

def plot_ts(tsarray):
    for ts in tsarray:
        xs = ts/(1+ts)  #zero
        precision = arange(xs+0.001,1.01,0.01)
        pod = precision/(precision*(1+1/ts)-1)
        plot(precision,pod,'b--',lindwidth=1)
        xlim(0,1)
        ylim(0,1)

        text(2*xs-0.02,2*xs-0.02,'ts='+str(ts),fontsize=16)

def plot_bias(biasarray):
    precision =arange(0,1.01,0.1)
    n = len(biasarray)
    for i in range(n):
        bias = biasarray[i*precision
        plot(precision,bias,'--',lindwidth=1)
        xlim(0,1)
        ylim(0,1)

    bias = precision
    plot(precision,bias,'-',lindwidth=1)
    xlim(0,1)
    ylim(0,1)
   
    for i in range(n):
        if(i<4):
            text((1*0.775)/biasarray[i+0.1,0.89,'bias='+str(biasarray[i),coordinates='figure',fontsize=16)
        else:
            text(0.91,0.775*biasarray[i+0.1,'bias='+str(biasarray[i),coordinates='figure',fontsize=16)

#plot
tss = arange(0.1,1.0,0.1)
figure(figsize=[800,800)
ax=axes(position=[0.13,0.11,0.775,0.775,aspect='equal',clop=False)
plot_ts(tss)
biass = array([5,2.5,1.67,1.25,0.8,0.6,0.4,0.2)
plot_bias(biass)

precision01 = 0.3
pod01 = 0.5
s01 = scatter(precision01,pod01,'S',color = 'r')
precision10 = 0.4
pod10 = 0.6
s02 = scatter(precision10,pod10,'S',color = 'g')
legend([s01,s02)
xlabel(u'成功率',fontname=u'黑体',fontsize=16)
ylabel(u'命中率',fontname=u'黑体',fontsize=16)

一大碗年糕

MeteoInfo

很好的分享贴

任性小王子

不错不错，非常感谢楼主分享