NCL绘制温度图出现异样

dashuhecanghais

各位同学好，我用NCL的gsn功能绘制EC数据温度图，为什么会出现这个形状的图的？
有人知道怎么解决吗？

SunJiaming

发一下代码看看

guoguohh

可能没有给数据赋坐标属性

dashuhecanghais

f = addfile("test.nc","r")
temperature = f->temperature
wks = gsn_open_wks("png",get_script_prefix_name())
gsn_define_colormap(wks,"whBlGrYeRe")
res = True
res@cnFillOn = True
res@cnInfoLabelOn = True
res@gsnSpreadColors = True

res@mpMinLatF = 5
res@mpMaxLatF = 60
res@mpMinLonF = 70
res@mpMaxLonF = 140

plot = gsn_csm_contour_map(wks,temperature(0,:,:),res)

dashuhecanghais

dashuhecanghais 发表于 2023-6-30 14:48
f = addfile("test.nc","r")
temperature = f->temperature
wks = gsn_open_wks("png",get_script_prefix ...

谢谢各位大神了

LeslieTse

貌似是数据的问题。用ncl_filedump命令看一下你的数据类型是不是short类型，如果是的话需要用函数short2flt转换一下。
https://www.ncl.ucar.edu/Documen ... ted/short2flt.shtml

dashuhecanghais

LeslieTse 发表于 2023-7-2 11:09
貌似是数据的问题。用ncl_filedump命令看一下你的数据类型是不是short类型，如果是的话需要用函数short2flt ...

看了一下，数据类型是type，不过数据不是全球的而是中国区的，这会不会有影响呢？

dashuhecanghais

guoguohh 发表于 2023-6-29 09:59
可能没有给数据赋坐标属性

数据有坐标属性的，会不会是gsn调用有隐藏的要求呢？

guoguohh

转一下数据格式，建议画图设置加上
  res@gsnDraw              = False        ; don't draw yet
  res@gsnFrame             =False    ; don't advance frame yet
  res@gsnAddCyclic         = False     ; plotted dataa are not cyclic

LeslieTse

dashuhecanghais 发表于 2023-7-3 10:10
看了一下，数据类型是type，不过数据不是全球的而是中国区的，这会不会有影响呢？

区域不存在影响吧，你这个是分布看起来是错乱的。type不是数据类型吧，用ncl_filedump看一下变量是不是short类型，如果是的话转成float类型。
通常变量都是float或者double类型，如下：
   variables:
      double time ( time )
         standard_name :        time
         long_name :        time
         bounds :        time_bnds
         units :        days since 1851-01-01 00:00:00
         calendar :        365_day
         axis :        T

      double time_bnds ( time, bnds )

      double lon ( lon )
         standard_name :        longitude
         long_name :        longitude
         units :        degrees_east
         axis :        X

      double lat ( lat )
         standard_name :        latitude
         long_name :        latitude
         units :        degrees_north
         axis :        Y

      float TREFHT ( time, lat, lon )
         long_name :        Reference height temperature
         units :        K
         cell_method :        time: mean