NCL绘制青藏高原的边界（已有边界经纬度坐标）

贝壳女孩

老师给了一个青藏高原边界的经纬度坐标的txt文件，读入文件之后不知道怎么画出边界轮廓，求各位大神解答，谢谢！下面是我的脚本和图：
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/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"
begin
minGlat=25.
maxGlat=40.
minGlon=74.
maxGlon=104.
;Three river source header"
slat=31.7
elat=36.2
slon=89.5
elon=102.2
slat2=31.7-1.
elat2=36.2-1.
slon2=89.5-1.
elon2=102.2-1.
f = addfile("$NCARG_ROOT/laiavgrcp85.nc","r")
lai1 = f->lai
wks = gsn_open_wks("png","laiavgrcp85")
res= True
res@gsnDraw            = False             ; don't draw yet
res@gsnFrame           = False             ; don't advance frame yet
res@cnLevelSpacingF = 0.5
res@gsnAddCyclic      = False
res@cnFillOn             = True           ; Turn off contour lines
res@cnLinesOn            = False           ; Turn off contour lines
res@cnFillPalette      =  "NOC_ndvi"
res@cnFillMode           = "RasterFill"
plot = gsn_csm_contour(wks,lai1(0,:,:),res)
cnres           = True
cnres@china     = False       ;draw china map or not
cnres@river     = True       ;draw changjiang&huanghe or not
cnres@province  = True       ;draw province boundary or not
cnres@nanhai    = False       ;draw nanhai or not
cnres@diqu      = False       ; draw diqujie or not

chinamap = add_china_map(wks,plot(0),cnres)
;TP boundary
coord = asciiread("$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/TP_boundary.txt",(/7300,2/),"float");读入文件下面该怎样绘图
;panel plot

draw(plot)
frame(wks)
end

一大碗年糕

试试看gsn_polygon_ndc绘制多边形
或者gsn_add_polygon绘制折线，这个的话应该要n+1个点，第一个点和第n+1个点一样形成闭合多边形
gsFillColor修改填充颜色

贝壳女孩

一大碗年糕 发表于 2020-4-10 19:50
试试看gsn_polygon_ndc绘制多边形
或者gsn_add_polygon绘制折线，这个的话应该要n+1个点，第一个点和第n+1 ...

好的，感谢！我去试试！