NCL绘制单省图

紫菊花

begin
fname ="/mnt/e/huitu/data/td.txt"
data=asciiread(fname,(/119,4/),"float")   ;119*4的矩阵
lat=data(:,0)  ;0表示第一列
lon=data(:,1)
td=data(:,2)
time=data(:,3)
nrow =119              ;共有119个时次的数据
wks=gsn_open_wks ("x11","409_508td_station")
;--------------------------提取内蒙省---------------------------------
res=True
res@gsnMaximize=True
res@gsnFrame = False
res@gsnDraw =False
res@mpOutlineOn=True
res@mpMinLatF  =35
res@mpMaxLatF  =55
res@mpMinLonF  =95
res@mpMaxLonF  =130
res@pmTickMarkDisplayMode ="Always"   ;坐标轴标签上需要加°
res@mpDataBaseVersion="Ncarg4_1"
res@mpDataSetName  ="Earth..4"
res@mpAreaMaskingOn =True
res@mpMaskAreaSpecifiers =(/"China: Nei Mongol"/)
res@mpOceanFillColor =0       ;让海洋、陆地、湖泊等不填色
res@mpInlandWaterFillColor =0
res@mpLandFillColor =0
res@mpOutlineSpecifiers ="China: Nei Mongol"
;res@mpGeophysicalLineThicknessF = 2.0    ;加粗边界线
res@mpOutlineBoundarySets ="NoBoundaries"   ;不绘制设定范围的其它地图
;res@mpGridAndLimbOn           = False   ;不绘制地图投影边缘线
map=gsn_csm_map(wks,res)
;----------------添加散点--------------------------------
mkres=True
mkres@gsMarkerIndex =16
mkres@gsMarkerSizeF=0.005
mkres@gsMarkerColor="black"
do i=0, nrow-1  ;nrow =119
   str=unique_string("string")
   map@$str$ = gsn_add_polymarker(wks,map,lon(i),lat(i),mkres)
end do
  draw(map)
  frame(wks)     ; advance the frame
end

ccznizuibang

感谢楼主      

whhh昊

感谢楼主 学习到了