NCL画图求助：画xy折线图时由于缺测值失败

里斯斯里

1. begin
   
2. _getvar = asciiread("Chengde_T_1", -1, "float")
   
3. y = new(dimsizes(_getvar), typeof(_getvar), -999.9)
   
4. y = _getvar
   
5. y@_FillValue = -999.9
   
6. time = ispan(1, dimsizes(y), 1)
   
7. 
   
8. res = True
   
9. res@xyMarkLineMode = "Lines"
   
10. res@xyDashPattern = 5
    
11. res@xyMarker = 0
    
12. res@gsnCenterString = "xyMarkLineMode = 'Lines'"
    
13. 
    
14. wks_1 = gsn_open_wks("png", "plot1")
    
15. plot_1 = gsn_csm_xy(wks_1, time, y, res)
    
16. 
    
17. res@xyMarkLineMode = "MarkLines"
    
18. res@gsnCenterString = "xyMarkLineMode = 'MarkLines'"
    
19. wks_2 = gsn_open_wks("png", "plot2")
    
20. plot_2 = gsn_csm_xy(wks_2, time, y, res)
    
21. end

复制代码

由于画图的数据存在缺测值（-999.9），在xy二维坐标图上能画散点图，但是连起来的折线图就是画不出来，请大佬帮忙看看，感激不尽。
数据格式和上图代码画的两张图如下所示：

灭火器

你这缺测值也太多了，画不出来才是正常现象……
非要画的话把有效值提取出来，或者插值吧。

里斯斯里

灭火器 发表于 2019-7-1 09:37
你这缺测值也太多了，画不出来才是正常现象……
非要画的话把有效值提取出来，或者插值吧。

是一个站点的时间序列，三分之二是缺测值，插值应该不合适，
另外想和模拟值作对比，单独提取有效值出来画图也不行
我想的是有没有什么设置能使数据里的有效值连成一条线（忽略缺测值）

灭火器

既然是连成一条线，那么把缺测值线性插值得到的图像是一样的。

里斯斯里

灭火器 发表于 2019-7-1 14:12
既然是连成一条线，那么把缺测值线性插值得到的图像是一样的。

懂了，谢谢了