NCL利用站点值画图

叁玖壹

虽然有大佬有发教程但对于新手来说还是很难，自己用的脚本贴出来供大家参考

HDdd

你最好复制脚本 别人下载txtx还要扣除积分

叁玖壹

begin
data = asciiread("f:/htsj/ps.txt",(/687,4/),"float")
lon=data(:,1)
lat=data(:,2)
ps=data(:,3)
olon=new(127,"float")
olat=new(81,"float")
data_p=new((/81,127/),"float")
do i=0,63*2
olon(i)=72+i*0.5
end do
do l=0,40*2
olat(l)=17+l*0.5
end do

lon!0="lon"
lon@long_name="lonitude"
lon@units="degrees_east"
lon&lon=lon
lat!0="lat"
lat@long_name="latitude"
lat@units="degrees_north"
lat&lat=lat

olon!0="lon"
olon@long_name="lonitude"
olon@units="degrees_east"
olon&lon=olon
olat!0="lat"
olat@long_name="latitude"
olat@units="degrees_north"
olat&lat=olat
rscan=(/10,5,3/)
data_p=obj_anal_ic_deprecated(lon,lat,ps,olon,olat,rscan,False)

wks=gsn_open_wks("png","pinshu")
res=True
res@gsnAddCyclic=False
res@gsnMaximize=True
res@mpDataBaseVersion="MediumRes"
res@mpDataSetName="Earth..4"
res@mpOutlineOn=True
res@mpFillOn = True
res@mpOutlineSpecifiers=(/"China:states","Taiwan"/)
res@mpGeophysicalLineThicknessF=2.0
res@mpNationalLineThicknessF=2.0
res@mpMinLatF=15.0
res@mpMaxLatF=55.0
res@mpMinLonF=105
res@mpMaxLonF=135.0
res@cnLinesOn=True
res@cnFillOn                        = True
res@cnLineDrawOrder       = "PreDraw"
res@cnLineLabelsOn=True
;res@cnLevelSelectionMode="AutomaticLevels"
plot=gsn_csm_contour_map(wks,data_p,res)
end

sxc123

1111111111111

sxc123

1111111111111111

猪八戒背西瓜

谢谢谢谢谢谢感谢分享

xenos_an

能不能附张预览图

叁玖壹

xenos_an 发表于 2020-5-14 19:48
能不能附张预览图

这个是新改的
begin
data = asciiread("f:/jg/winter/R95.txt",(/663,4/),"float")

lon=data(:,1)
lat=data(:,2)
ps=data(:,3)/10
olon=new(127,"float")
olat=new(81,"float")
data_p=new((/81,127/),"float")
do i=0,63*2
olon(i)=72+i*0.5
end do
do l=0,40*2
olat(l)=17+l*0.5
end do

lon!0="lon"
lon@long_name="lonitude"
lon@units="degrees_east"
lon&lon=lon
lat!0="lat"
lat@long_name="latitude"
lat@units="degrees_north"
lat&lat=lat

olon!0="lon"
olon@long_name="lonitude"
olon@units="degrees_east"
olon&lon=olon
olat!0="lat"
olat@long_name="latitude"
olat@units="degrees_north"
olat&lat=olat
rscan=(/10,5,3/)
data_p=obj_anal_ic_deprecated(lon,lat,ps,olon,olat,rscan,False)

wks=gsn_open_wks("png","R95")
res=True
res@gsnAddCyclic=False
res@gsnMaximize=True
res@mpDataBaseVersion="MediumRes"
res@mpDataSetName="Earth..4"
res@mpOutlineOn=True
res@mpFillOn = True
res@mpOutlineSpecifiers=(/"China:states","Taiwan"/)
res@mpGeophysicalLineThicknessF=2.0
res@mpNationalLineThicknessF=2.0
res@mpMinLatF=15.0
res@mpMaxLatF=55.0
res@mpMinLonF=105
res@mpMaxLonF=135.0
  res@cnFillDrawOrder  = "PreDraw"   
   res@mpAreaMaskingOn = True                        
         res@mpMaskAreaSpecifiers = (/"China:states","Taiwan"/)
res@mpOceanFillColor = 0
             res@mpLandFillColor = 0         
          res@mpInlandWaterFillColor = 0  
res@cnLinesOn=False
res@cnFillOn                        = True
res@cnLineLabelsOn=False
res@cnLevelSelectionMode =  "ManualLevels"   
  res@cnMinLevelValF       = 0
  res@cnMaxLevelValF       =  65
  res@cnLevelSpacingF      =  5
res@gsnLeftString = "winter"

;res@cnLevelSelectionMode="AutomaticLevels"
plot=gsn_csm_contour_map(wks,data_p,res)
end

xenos_an

此帖仅作者可见

叁玖壹

xenos_an 发表于 2020-5-15 14:59
此帖仅作者可见

你再试试