NCL 根据两期下垫面数据，画土地利用变化图

FrankieLJY

各位坛友，想用NCL做两期geogrid.nc文件中LU_INDEX（土地利用数据）的变化图，想通过两个nc文件来做差，如果值为0就是没有变化的地方，如果变化不为0，就设置为1，体现出哪些地方发生了变化。想请教大家如何实现？！
现在只能对一个nc文件的这个属性进行绘图：
begin
    Dir="./"
    filename="geo_em.d03"
    inputFile=addfile(Dir+filename+".nc", "r")
    type="pdf"
    wks=gsn_open_wks(type, filename)
    ;-------------------------------------------------------------
    res=True
    pltres=True
    mpres=True

    res@tiMainString="MyTest"
    res@cnFillOn=True
    res@cnFillMode="RasterFill"
    res@cnLinesOn=False
    res@cnFillPalette="NCV_roullet"
    res@cnMinLevelValF=1
    res@cnMaxLevelValF=24
    res@cnLevelSpacingF=1
    res@mpDataBaseVersion="HighRes"
    mpres@mpNationalLineColor= "Black"
    mpres@mpOutlineBoundarySets="AllBoundaries"
    ;-------------------------------------------------------------
    it=0
    land=wrf_user_getvar(inputFile, "LU_INDEX", it)
    land@lat2d=wrf_user_getvar(inputFile, "XLAT_M", it)
    land@lon2d=wrf_user_getvar(inputFile, "XLONG_M", it)
    printVarSummary(land)

    resParam=res
    resParam@cnFillOn=True
    counter=wrf_contour(inputFile, wks, land, resParam)
    plot=wrf_map_overlays(inputFile, wks, (/counter/), pltres, mpres)
end

FrankieLJY

因为对NCL不是很熟 以为可以直获取两个属性做差： land=wrf_user_getvar(inputFile, "LU_INDEX", it) -  land=wrf_user_getvar(inputFile2, "LU_INDEX", it),但是好像是不可以的

FrankieLJY

然后听说有个程序叫CDO，但是也是才上手还不太懂如何操作，只是把程序装好了，如果有需要CDO程序的，我也分享出来

jinljie

不知道楼主解决了没有呢？最近我也想到要做这些事情

放逐流年

用CDO好像确实可以，不过地理信息好像会丢失，楼主是否用NCL 解决了？

FrankieLJY

jinljie 发表于 2017-1-1 19:07
不知道楼主解决了没有呢？最近我也想到要做这些事情

现在还没有解决 还是用ArcGIS做的

FrankieLJY

放逐流年 发表于 2017-1-2 10:48
用CDO好像确实可以，不过地理信息好像会丢失，楼主是否用NCL 解决了？

还没有解决这个 CDO和NCL这两个都不太会用，目前是用ArcGIS做的，不知您会吗？

jinljie

我都不会呀

FrankieLJY

jinljie 发表于 2017-1-3 14:28
我都不会呀

ArcGIS是全图形化的 比较容易上手

jinljie

楼主，我又来找你了，能不能教一下arcgis做这个的方法