请问文本格式的反射率怎么用NCL画出来————问题已解决

积极歪歪

txt格式： LNG(°)      LAT(°)      V(dBZ)
122.366973     37.146756         N/A
122.356447     37.143512       -28.5
122.345920     37.140268         4.5
122.335393     37.137024         8.5
122.324867     37.133780         4.0
122.366912     37.146885         N/A
122.356324     37.143770       -28.0
122.345736     37.140656       -13.0
122.335149     37.137541         0.0
122.324561     37.134426         9.5
122.366851     37.147021         N/A
122.356201     37.144043       -25.5
122.345552     37.141064        -3.0
122.334902     37.138086         7.5
122.324253     37.135107         9.0
122.366792     37.147159         N/A
122.356084     37.144317       -24.5
122.345375     37.141476        -1.0
122.334667     37.138634         8.0
122.323959     37.135793         7.5
122.366738     37.147290         N/A
122.355977     37.144580       -22.0
122.345215     37.141870         0.0
122.334453     37.139159         9.0
122.323692     37.136449        11.0
122.366685     37.147428         N/A


我编的NCL:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"

begin

wks = gsn_open_wks("png","201801100003_R_CAPPI_1.5km")     ; Open a workstation
;-----------------------------设置色标-------------------------

gsn_define_colormap(wks,"default") ;


;-----------------读取数据-----------------------------
diri = "./"
  fili = "201801100003_R_CAPPI_1.5km.txt"
  nlvl = 58288
  ncol = 3
  TestData = asciiread (diri+fili , (/nlvl,ncol/), "float")

lon = TestData (:,0)  
lat  = TestData (:,1)   
dbz = TestData (:,2)  

;---Create lat and lon coordinate variables
lon!0          = "lon"
lon@long_name  = "lon"
lon@units      = "degrees-east"
lon&lon        = lon
lat!0          = "lat"
lat@long_name  = "lat"
lat@units      = "degrees_north"
lat&lat        = lat

;-------------------------------------参数设置-----------------------------------------




mpres              = True
mpres@gsnMaximize   = True                                       ;图片最大化设置？
mpres@gsnDraw       = False                                         ;暂不绘制
mpres@gsnFrame      = False                                         ;暂不翻页
mpres@mpMinLonF    = 117.8
mpres@mpMaxLonF    =  123.0
mpres@mpMinLatF    =  35.0
mpres@mpMaxLatF    =  39.0
mpres@mpOutlineOn                 = True                          ;绘制边界               
mpres@mpOutlineSpecifiers         = (/"China:States"/)    ;绘制的边界名称
mpres@mpOutlineBoundarySets       = "NoBoundaries"  ;绘制国界线，但中国国界线（"National"）有误
mpres@mpDataBaseVersion           = "MediumRes"       ;中等分辨率地图数据库
mpres@mpDataSetName               = "Earth..4"              ;该数据集中包含了中国省界
;mpres@mpAreaMaskingOn             = True                    ;使能填充覆盖
;mpres@mpLandFillColor             = "white"                    ;填色陆地的颜色
;mpres@mpInlandWaterFillColor      = "white"                ;填色内陆水体的颜色
;mpres@mpOceanFillColor            = "white"                   ;填色海洋的颜色
mpres@mpFillDrawOrder             = "PostDraw"              ;
mpres@mpGeophysicalLineThicknessF = 2.                           ;地图线条粗细
map = gsn_csm_map(wks,mpres)


;-----------------------主要部分------------------------------------------------

draw(map)
res              = True
res@cnFillOn       =True
plot =gsn_csm_contour_map(wks,dbz,res)

end


初学者，对错误实在无从下手。我查了很多资料，应该读出数据后，定义经纬度，然后插值。定义经纬度老是出错。。。。










问题已解决，这个雷达数据本身就很密集了，所以没有插值，直接画出来，效果与原图相似。

星海

用asciiread命令读入数据，然后正常绘图就行

星海

读取一维的经纬度和数据后，创建一个新的二维经纬度网格，用obj_anal_ic_deprecated命令把一维数据插值到二维网格中绘图，具体用法见 http://www.ncl.ucar.edu/Document ... ic_deprecated.shtml

积极歪歪

问题已经解决，脚本在附件中。可以正确出图。