填色图 马赛克图 值不相同

Gabriella

大家好，我用ncl画了一个马赛克图，如图一，值都是对的；
然后把  res@cnFillMode ="RasterFill" 注释掉，其他不变，出来的填色图图2与图1相差很多，值都不同，实在不知道问题在哪里，求助大家，谢谢。
脚本直接贴了画图部分，前面数据处理无误，算出来score（60,60），画score（60,60）的分布；

脚本如下：

图一 rasterfill

图2

;--------画肉质期X5 雨量分布----------;
    wks=gsn_open_wks("png", "F:/data/product_CQ/HUAJIAO_2019/caishou/X6_prep_2019_larger")
     ;gsn_define_colormap(wks, "grads_rainbow")
    res                         = True
    res@gsnMaximize             =True              
    res@gsnDraw                 =False               
    res@gsnFrame                =False
    res@gsnAddCyclic            =False   ; regional data, don't add pt

    res@cnFillOn              = True
    res@cnLinesOn             = True
    res@cnLineLabelsOn        = True
    ;res@cnFillMode            ="RasterFill"     图2 只注释掉了这一行

    res@cnLevelSelectionMode = "ExplicitLevels" ; set manual contour levels
    res@cnLevels             =  (/0.6,0.8,1.0/)              ; set min contour level
    res@cnFillColors         =  (/"orange","yellow","blue","green"/)              ; set max contour level
    res@lbOrientation        = "vertical"              

    res@mpMaxLatF  =30
    res@mpMinLatF  =28.0333
    res@mpMaxLonF  =106.9667
    res@mpMinLonF  =105
     res@tmXBMode             ="Explicit"
     res@tmXBValues           =(/105,105.5,106.0,106.5,106.9/)
     res@tmXBLabels           =(/"105","105.5","106.0","106.5","106.9~S~o~N~E"/)
     res@tmYLMode             ="Explicit"
     res@tmYLValues           =(/28.1,28.5,29,29.5,30/)
     res@tmYLLabels           =(/"28.1","28.5","29","29.5","30~S~o~N~N"/)

    res@mpDataBaseVersion           = "MediumRes"      ; use finer database


     ;---------setting---------------;
    res@tmXBLabelDeltaF = -0.5                       ;使其更靠近X轴
    res@tmYLLabelDeltaF = -0.5
    res@tmXBMajorLengthF = 0.01                      ;刻度线长短
    res@tmYLMajorLengthF = 0.01            
    res@tmXBMajorThicknessF = 5.0             ;刻度线粗细
    res@tmYLMajorThicknessF = 5.0
    res@tmXBMinorOn =False
    res@tmYLMinorOn =False
    res@tmXBLabelFontHeightF = 0.02                 ;字体大小
    res@tmYLLabelFontHeightF = 0.02
    res@tmXBLabelFontThicknessF = 3
    res@tmYLLabelFontThicknessF = 3
    res@tmBorderThicknessF     = 5.0   
     plot=gsn_csm_contour_map(wks, score, res)

      lon_add=(/105.8167, 106.6334 ,106.6334, 105.8167, 105.8167/)
    lat_add=(/28.4667, 28.4667, 29.4667, 29.4667, 28.4667 /)
    ;add 区域
    lnres  =True
    lnres@gsLineColor ="red"
    lnres@gsLineThicknessF=10.0
    lnres@gsLineDashPattern =1
    Polygon=gsn_add_polyline(wks,plot,lon_add,lat_add,lnres)


    draw(plot)
    frame(wks)

千里草

你这原始数据是格点数据吧？配完经纬度后没有插值是吗？所以导致以栅格出图，它是每个格点上的真实值，如果去掉那句话，ncl会自动用周边格点数据插值，我认为是这样

Gabriella

千里草 发表于 2019-8-27 11:10
你这原始数据是格点数据吧？配完经纬度后没有插值是吗？所以导致以栅格出图，它是每个格点上的真实值，如果 ...

谢谢您 我刚试了一下改cnlevel
把原来的
res@cnLevelSelectionMode = "ExplicitLevels" ; set manual contour levels
res@cnLevels             =  (/0.6,0.8,1.0/)              ; set min contour level
res@cnFillColors         =  (/"orange","yellow","blue","green"/)              ; set max contour level
res@lbOrientation        = "vertical"
改成
colors=(/"orange","yellow","blue","green"/)
    res@cnFillPalette        = colors
    res@cnLevelSelectionMode ="ManualLevels"
     res@cnMinLevelValF  =0.6
     res@cnMaxLevelCount =1.0
     res@cnLevelSpacingF =0.2
    res@lbOrientation        = "vertical"  
是能出图的，但是labelbar的最大值不是1 是0.975 不知道怎么改；

Gabriella

通过改原来的画图等值线设置属性    ; res@cnLevelSelectionMode = "ExplicitLevels" ; set manual contour levels    ; res@cnLevels             =  (/0.6,0.8,1.0/)              ; set min contour level
    ; res@cnFillColors         =  (/"orange","yellow","blue","green"/)              ; set max contour level
    ; res@lbOrientation        = "vertical"              
改为：
    colors=(/"orange","yellow","blue","green"/)
    res@cnFillPalette        = colors
    res@cnLevelSelectionMode ="ManualLevels"
     res@cnMinLevelValF  =0.6
     res@cnMaxLevelCount =1.0
     res@cnLevelSpacingF =0.2

     ;laberbar设置
    res@lbLabelsOn =True
    res@lbOrientation        = "vertical"              
    res@lbLabelStrings       = (/"0.6","0.8","1.0"/)

还有一个问题就是黄色绿色的边界总是有蓝色的描边，按理说不应该有，这又是为啥呢？

kaiofchina

{:eb302:}

王小婷enjoy

会不会是填色的时候，插值间隔设置太大了呢？

你是最美七月天

学习一下！学一下！！

云云云云云

学习学习，谢谢分享