关于NCL绘制AWX格式数据的问题

zly4814624

是这样的，本人现需要用风云2E的AWX格式数据进行NCL绘图，因为NCL无法直接读取AWX格式的数据，所以先用fortran转换成grd格式的文件在进行读取，在进行NCL编译的过程也用到了国家卫星网的对照经纬度，但是。画出来的图非常奇怪，我查了半天也不知道是哪里出了问题，烦请各位大神帮帮看一下是什么问题，以下是我的图和脚本：

;**********************************************************
load "$WORKDIR/include/library.ncl"
;***********************************************************
begin

        f1    = "/home/suzhou/AWX/Z_SATE_C_BAWX_20160830045501_P_FY2E_SEC_IR1_LCN_20160830_0430.grd"
        f=fbindirread(f1,0,(/512,512/),"integer")
        tbb=int2flt(f(:,:))      ; [512] x [512]

        f2    = "/home/suzhou/AWX/land.dat"
        nlon     = 512
        nlat     = 512
        latlon2d = fbindirread(f2,-1,(/512,512,2/),"float")

        lon2d  = latlon2d(:,:,0)
        lat2d  = latlon2d(:,:,1)
        ; printVarSummary(lon2d)     ; [512] x [512]
        ; printVarSummary(lat2d)     ; [512] x [512]

        wks = gsn_open_wks("png","/home/suzhou/AWX")
        lat2d@units               = "degrees_north"
        lon2d@units               = "degrees_east"

;---- Add meta data for TBB.
        tbb@lat2d                 = lat2d
        tbb@lon2d                 = lon2d
        tbb@units                 = "degree Kelvin"
        tbb@long_name             = "Temperature of Bright Blackbody"
        tbb@coordinates           = "lat2d lon2d"
        tbb@_FillValue            = -999.

        tbb = tbb-273.15        ;(convert degree Kelvin to degree Celsius)
        tbb@units = "degree Celsius"

;---- Begin to plot.
        wks_type = "png"
        wks_type@wkWidth  = 1000          ;工作台宽度
        wks_type@wkHeight = 1000          ;工作台高度

        res = True

        res@gsnDraw                                 = False
        res@gsnFrame                                = False
        res@gsnMaximize                             = True
        res@gsnPaperOrientation                     = "portrait"
        res@gsnAddCyclic                            = False         ; regional data
        ;tfDoNDCOverlay should be set False (default) when lat2d and lon2d are used.
        ;res@tfDoNDCOverlay                         = True                        
        res@pmTickMarkDisplayMode                   = "Always"

        res@mpMinLonF                               = 70
        res@mpMaxLonF                               = 140
        res@mpMinLatF                               = 0
        res@mpMaxLatF                               = 55

        ;---- Map Set
        res@mpDataSetName                           = "Earth..4"
        res@mpDataBaseVersion                       = "MediumRes"
        res@mpOutlineSpecifiers                     = (/"China","China:Provinces"/)
        res@mpOutlineBoundarySets                   = "NoBoundaries"
        res@mpNationalLineColor                     = "black"
        res@mpProvincialLineColor                   = "black"
        res@mpGeophysicalLineColor                  = "black"
        res@mpNationalLineThicknessF                = 3
        res@mpProvincialLineThicknessF              = 3
        res@mpGeophysicalLineThicknessF             = 3

        res@trGridType                              = "TriangularMesh"
        res@cnFillOn                                = True
        res@cnFillMode                              = "RasterFill"
        res@cnLinesOn                               = False

        res@gsnLeftString                           = "TBB"
        res@gsnRightString                          = "~S~o~N~C"

        ;---- Color set
        ; cmap  = read_colormap_file("precip3_16lev")
        ; cmap1 = cmap(4::,:)
        ; cmap1(0,:) = cmap(0,:)
        ; cmap1 = cmap1(::-1,:)
        ; res@cnFillPalette = cmap1

        ; res@cnLevelSelectionMode = "ExplicitLevels"
        ; res@cnLevels = (/-30,-40,-50,-60,-70,-80/)
        ; res@lbOrientation = "vertical"

        plot = gsn_csm_contour_map(wks,tbb,res)

        draw(plot)
        frame(wks)

end

雪颜

很好很有艺术感

zly4814624

雪颜 发表于 2016-9-19 14:58
很好很有艺术感

我应该转专业的。

盈星

经纬度与像素点不匹配，少了投影经纬度配置。

weidinghe

太需要了，真是好人