NCL绘制泰勒图

youngman

泰勒图多用于评估模式结果相对于参考值（控制实验，OBS，再分析资料）的区别（相关性，方差，bias），当然，这个评估指标也可自己定义，可以评估多个不同case和不同变量。
NCL绘制最基本的泰勒图需要计算两个量：ratio 和 cc
ratio是方差之比：Case_Variance/Reference_Variance
cc是相关系数：cross correlation coef of Case to Reference
由这两个值定义泰勒图坐标系一个点的位置，然后绘图。
NCL绘制泰勒图另外需要load一个图像框架函数：taylor_diagram.ncl（见官网），参考官网一个简单绘图示例：

1. load "./source_taylor/taylor_diagram.ncl"
   
2. ; Assume the following have already been computed:
   
3. ;   _ratio are the ratio:  Case_Variance/Reference_Variance
   
4. ;   _cc    are the cross correlation coef of Case to Reference
   
5. ; All ratio values must be between 0 and 1.65 (change in taylor_diagram.ncl)
   
6. ; All cross correlation values are 0.0 to 1.0 [inclusive]
   
7. ;-------------------------------------------------------------
   
8. ; Cases [Model]
   
9. case      = (/"Case A", "Case B"/)
   
10. nCase     = dimsizes(case)                  ; # of Cases [Cases]
    
11. ; variables compared
    
12. var       = (/"SLP","Tsfc","Prc","Prc 30S-30N","LW","SW", "U300", "Guess"/)
    
13. nVar      = dimsizes(var)                   ; # of Variables
    
14. ; "Case A"
    
15. CA_ratio   = (/1.230, 0.988, 1.092, 1.172, 1.064, 0.966, 1.079, 0.781/)
    
16. CA_cc      = (/0.958, 0.973, 0.740, 0.743, 0.922, 0.982, 0.952, 0.433/)
    
17. ; "Case B"
    
18. CB_ratio   = (/1.129, 0.996, 1.016, 1.134, 1.023, 0.962, 1.048, 0.852 /)
    
19. CB_cc      = (/0.963, 0.975, 0.801, 0.814, 0.946, 0.984, 0.968, 0.647 /)
    
20. ; Put ratios and pattern correlations into arrays for plotting
    
21. ratio      = new((/nCase, nVar/),typeof(CA_cc))
    
22. cc         = new((/nCase, nVar/),typeof(CA_cc))
    
23. ratio(0,:) = CA_ratio
    
24. ratio(1,:) = CB_ratio
    
25. cc(0,:)    = CA_cc
    
26. cc(1,:)    = CB_cc
    
27. 
    
28. wks = gsn_open_wks("png","taylor1")
    
29. res                     = True
    
30. res@tiMainString        = "Example"
    
31. res@caseLabels          = case
    
32. res@varLabels           = var
    
33. res@varLabelsYloc       = 1.5                ; Move location of variable labels [default 0.45]
    
34. res@Colors              = (/"red","blue"/)
    
35. res@Markers             = (/16, 16/)         ; make all solid fill
    
36. res@markerTxYOffset     = 0.04               ; offset btwn marker & label
    
37. res@gsMarkerSizeF       = 0.01               ; marker size
    
38. res@txFontHeightF       = 0.015              ; text size
    
39. res@stnRad              = (/0.5, 1.5/)       ; additional standard radii
    
40. res@ccRays              = (/0.6, 0.9/)       ; correlation rays
    
41. res@centerDiffRMS       = True               ; RMS 'circles'
    
42. plot = taylor_diagram(wks, ratio, cc, res)
    

复制代码

taylor_diagram.ncl (18.52 KB, 下载次数: 88)

2020-5-11 17:59 上传

点击文件名下载附件

youngman

Datian 发表于 2020-7-29 09:42
楼主您好，请问在计算ratio和cc之前需要对原始数据进行标准化处理吗？

一般不需要吧，看你需求了，标准化可能是为了在不同量级的数据之间比较，模拟结果一般和观测量级不会差多少。另外NCL在计算相关系数时默认取距平了

jzcmwwbb

求助！想请问一下泰勒图如果想画半圆，需要改脚本里的哪些参数呢，我看到官网里有半圆的程序，但是那个是画三个变量的，我只有两个变量，但是变量有负值，需要画半圆。谢谢！

小鸣REN

楼主，你好。请问如果想把图上的1-8数字换成不同的标记表示（比方说，圆形、三角形、矩形等等），该怎么设置啊？在官网也没找到例子，请指教！

youngman

小鸣REN 发表于 2020-6-10 23:03
楼主，你好。请问如果想把图上的1-8数字换成不同的标记表示（比方说，圆形、三角形、矩形等等），该怎么设 ...

修改res@Markers =

小鸣REN

youngman 发表于 2020-6-11 08:57
修改res@Markers =

我试了，不行。您可能没明白我的意思，我是想画类似下图这样的，用不同的marker代表不同模式，然后不同颜色代表不同变量。但是我按照ncl官网的教程画出了只能是数字编号代表不同模式，那样看着比较乱。第二张图是我画的，都叠在一起了，看不清。

simple.

绘制泰勒图时，可以生成CAM Metrics的一个表格，请问楼主知道这代表什么意思，要怎么分析呢？谢谢

youngman

小鸣REN 发表于 2020-6-12 09:19
我试了，不行。您可能没明白我的意思，我是想画类似下图这样的，用不同的marker代表不同模式，然后不同颜 ...

在taylor_diagram.ncl脚本366行：ptRes@txFontHeightF 调为0，下面3行注释掉，可以不显示数字，在你的脚本里把 varLabels 相关行注释掉，可以不显示关于var的label注释。然后把数组模式维和变量维交换一下，应该就可以实现了

youngman

simple. 发表于 2020-6-13 00:02
绘制泰勒图时，可以生成CAM Metrics的一个表格，请问楼主知道这代表什么意思，要怎么分析呢？谢谢

没用过这个功能，不太清楚~ 什么样的表格？贴一下看看？

simple.

youngman 发表于 2020-6-13 10:14
没用过这个功能，不太清楚~ 什么样的表格？贴一下看看？

C:\Users\L\Desktop\微信图片_20200603214603.png

simple.

就是这个表格，刚才传错了。。。

小鸣REN

youngman 发表于 2020-6-13 10:11
在taylor_diagram.ncl脚本366行：ptRes@txFontHeightF 调为0，下面3行注释掉，可以不显示数字，在你的脚 ...

谢谢楼主，我试一下