boxplot疑问

lysbaba

作的一个boxplot，代码片段如下：
  res@tmXBMode = "Manual"
  res@trXMinF=-6
  res@tmXBTickStartF  = 0
  res@tmXBTickSpacingF=24
  plot=boxplot(wks,xtime,plotdata,opti,res,llres)
由于x轴用预报时效，不方便使用，想将其转换为日期，
代码片段设置为：
res@tmXBMode = "Explicit"
res@tmXBValues      = xtime(0:40:4)
res@tmXBLabels     =xtimelabel(0:40:4)
res@tmXBMinorValues=xtime

xtime(0:40:4)就是(/0,24,48,72,96,120,144,168,192,216,240/)
xtimelabel是xtime转换的日期数组。
运行时提示

boxplot: Fatal:Number of XB Labels does not match number of boxes, exiting
请画过同类图的大虾帮忙。

longlivehj

boxplot这个绘图函数是开放源码的，有兴趣可以看一下。
如果tmXBMode是Explicit，那么tmXBLabels的个数必须和box的个数相等。至于你贴出来的图中为何x轴标签个数少于box个数，是因为你用Manual覆盖了boxplot程序中默认的Explicit，使程序中辛辛苦苦生成的tmXBLabels失去了作用。

不过，仍然可以考虑用下面的方式来实现你要的功能：
xtimelabel2 = new(dimsizes(xtimelabel), "string")
xtimelabel2(:) = ""
xtimelabel2(0:40:4) = xtimelabel(0:40:4)

res@tmXBMode = "Explicit"
res@tmXBValues = xtime
res@tmXBLabels = xtimelabel2

plot=boxplot(wks,xtime,plotdata,opti,res,llres)

这样做的一个副作用就是，暂时无法体现出major和minor，还需要一些美化的操作！
你试试看上面的方式行不行先！

lysbaba

佩服！！！
xtimelabel2(:) = ""
xtimelabel2(0:40:4) = xtimelabel(0:40:4)

这个方法肯定行的通。

clannand

你好，请问你调用shea_util.ncl有这个问题吗fatal:The result of the conditional expression yields a missing value. NCL can not determine branch, see ismissing function
fatal:["Execute.c":8575]:Execute: Error occurred at or near line 2551 in file ~/plot/shea_util.ncl

fatal:["Execute.c":8575]:Execute: Error occurred at or near line 33 in file box.ncl

彦祖

clannand 发表于 2016-12-31 22:43
你好，请问你调用shea_util.ncl有这个问题吗fatal:The result of the conditional expression yields a mis ...

我也出现这个报错，将横坐标赋值后就好了，res@tmXBValues = （/.../）