分辨率为12.5km的nc数据应如何绘图

爱睡觉的青团

从NSIDC下载的冰龄数据，是12.5km的，且经纬度是二维的，若插值到1°×1°的数据和原数据相比会有些许差距，尤其是靠近海岸线边缘没有值，想请教一下大家ncl怎么绘制分辨率为12.5km或者25km的数据，尤其是以极地为中心的，照着官网的例子画了但是底图和填色不对应。希望万能的朋友们能给出一个具体参考。下面第一张是插值绘制的图，第二张是原数据给的例图，同一天的。

爱睡觉的青团

害 突然就搞懂了
begin
    filename = "iceage_nh_12.5km_19840101_19841231_v4.1.nc"
    f = addfile(filename, "r")
    ice_age0 = f->age_of_sea_ice(0,:,:) ;1984年第一周的数据
    ice_age = byte2flt(ice_age0)
    ice_age@lat2d = f->latitude ;加上这两句就可以了
    ice_age@lon2d = f->longitude

    wks = gsn_open_wks("png","nic_ims")
    res = True
    res@cnFillOn = True
    res@cnLinesOn = False
    res@cnLineLabelsOn = False
    res@gsnPolar = "NH"                       
    res@mpMinLatF  = 50

    plot = gsn_csm_contour_map_polar(wks,ice_age,res)

end

爱睡觉的青团

同样的方法画海冰运动的数据files = "icemotion_daily_nh_25km_19781101_19781231_v4.1.nc"
就出不来图，只有地图
warning:ContourPlotInitialize: no valid values in scalar field; ContourPlot not possible:[errno=1101]

爱睡觉的青团

爱睡觉的青团 发表于 2024-11-2 13:37
同样的方法画海冰运动的数据files = "icemotion_daily_nh_25km_19781101_19781231_v4.1.nc"
就出不来图， ...

搞清楚了 是我读取的这个数据刚好全部都是缺测值
建议先加上下面这句话检查一下数据是否有值
print("Is viab all missing ? " + all(ismissing(viab)))