请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4535|回复: 6

NCL画土地利用类别图

[复制链接]

新浪微博达人勋

发表于 2022-11-14 21:36:03 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 笑个YY 于 2022-11-15 18:18 编辑

最近NCL读取LU_INDEX变量画土地利用类别图
labelbar色标出来的是数字,想把它变成相应的土地类别,请问怎么做啊

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin
  a=addfile("geo_em.d03.nc","r")
  wks_type = "png"
  wks_type@wkWidth = 3000
  wks_type@wkHeight = 3000
  wks = gsn_open_wks(wks_type,"landuse_21")
  gsn_define_colormap(wks,"NCV_roullet")
   it=0
   land=wrf_user_getvar(a, "LU_INDEX", it)
   printVarSummary(land)
   land@lat2d=wrf_user_getvar(a, "XLAT_M", it)
   land@lon2d=wrf_user_getvar(a, "XLONG_M", it)
   printVarSummary(land)

   res=True
   pltres=True
   mpres=True
   res@cnFillOn=True
   res@cnFillMode="RasterFill"
   res@cnLinesOn=False
   res@mpDataBaseVersion="MediumRes"
   res@mpDataSetName ="Earth..4"
   res@cnMinLevelValF=1
   res@cnMaxLevelValF=21
   res@cnLevelSpacingF=1

    res@lbTitleOn =False
    res@lbOrientation="vertical"
    res@lbLabelOffsetF=0.007
    res@pmLabelBarOrthogonalPosF=-1
    res@pmLabelBarParallelPosF=1.2
    res@lbLabelAlignment="BoxCenters"                                                                                                                                   
    res@lbLabelStrings=(/"Evergreen Needleleaf Forest","Evergreen Broadleaf Forest","Deciduous Needleleaf Forest","Deciduous Broadleaf Forest","Mixed Forests","Closed Shrublands","Open Shrublands","Woody Savannas","Savannas","Grasslands","Permanent Wetlands","Croplands","Urban and Built-Up","Cropland/Natural Vegetation Mosic","Snow and Ice","Barren or Sparsely Vegetated","Water","Wooded Tundra","Mixed Tundra","Barren Tundra","Lake"/)
    mpres@mpNationalLineColor= "Black"
    mpres@mpOutlineBoundarySets="AllBoundaries"
   
    resParam=res
    resParam@cnFillOn=True
    counter=wrf_contour(a, wks, land, resParam)
    plot=wrf_map_overlays(a, wks, (/counter/), pltres, mpres)
end

landused_modis03.png
landuse_21.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2022-11-15 09:07:15 | 显示全部楼层
labelstring可以单独设置
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-11-15 11:06:50 来自手机 | 显示全部楼层
用lbLabelStrings这个属性设置每个对应位置的名称就行了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-11-15 18:19:14 | 显示全部楼层
Zhangyuanping 发表于 2022-11-15 11:06
用lbLabelStrings这个属性设置每个对应位置的名称就行了

您好,根据意见设置了但还是不太对,我把代码放出来了,能帮忙看看哪里有错吗
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-11-15 18:42:29 | 显示全部楼层
随缘 发表于 2022-11-15 09:07
labelstring可以单独设置

您好,根据意见修改了但还是不对,代码如上,请问哪里有问题呀
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-11-15 19:23:56 来自手机 | 显示全部楼层
笑个YY 发表于 2022-11-15 18:19
您好,根据意见设置了但还是不太对,我把代码放出来了,能帮忙看看哪里有错吗

具体提示错误是什么样呢,或者说出图什么样,要不你加上把lbLabelBarOn这个属性设置为True看看
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-11-15 19:34:01 | 显示全部楼层
Zhangyuanping 发表于 2022-11-15 19:23
具体提示错误是什么样呢,或者说出图什么样,要不你加上把lbLabelBarOn这个属性设置为True看看

不好意思,刚没看到你放的图,我刚试了一下那个属性是可以设置的,只是可能你要注意找到对应的位置,不然设置的不在显示的范围内可能就显示不了,我放个图片你参考,算了,刚来家园,不知道怎么回复图,res = True
res@gsnFrame=True
res@gsnDraw=True
res@lbLabelBarOn                     =True
res@cnFillPalette                    ="precip_diff_1lev"
res@cnLevelSelectionMode             ="ManualLevels"
res@cnMaxLevelValF                   =0.12
res@cnMinLevelValF                   =-0.12
res@cnLevelSpacingF                  =0.12
res@lbLabelStrings                   =(/"Negative","","Positive"/)
res@gsnLeftString                    ="500hPa omega"
res@gsnRightString                   ="Pa/s"
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表