ncl画散点图

无敌大壮壮

1. f=readAsciiTable("20210720120000.txt",5,"float",1);读取文本文件，这种读取方式可以跳过表头，后面那个的意思是表头占的行数
   
2. data=f(:,:);把数据取出来
   
3. lat=data(:,1)
   
4. lon=data(:,2)
   
5. pre1=data(:,4)
   
6. 
   
7. lat@_FillValue=999990
   
8. lon@_FillValue=999990
   
9. pre1@_FillValue=999990
   
10. pre1(ind(ismissing(pre1))) = 0 ; ind将a中为缺省值的数组索引去除，并将其中对应缺省值的元素赋值为0
    
11. 
    
12. ; lat1 = where(lat.lt. 31.and.lat.gt. 37,lat@_FillValue,lat)
    
13. ; lon1 = where(lon.lt.110.and.lon.gt.117,lon@_FillValue,lon)
    
14. ;;;d = ind(pre1.le.50.and.pre1.gt.0)   ;;;;返回满足这个条件中pre1的索引
    
15. lat1 = lat(ind(lat.lt. 36.and.lat.gt. 31));;;;这个范围的经纬
    
16. lon1 = lat(ind(lon.lt.117.and.lon.gt.110))
    
17. 
    
18. A=new(6,float)
    
19. A=(/0,50,100,250,500,800/)
    
20. 
    
21. sizes  = ispan(30,55,5)/5000.    ; 0.00025 to 0.00055
    
22. colors = (/"orange","green","red","yellow","purple"/);,"blue","red","brown","red2","skyblue"/)
    
23. labels  = (/"0~F25~~~50","50~F25~~~100","100~F25~~~250","250~F25~~~500","500~F25~~~800"/)
    
24. xpos3    = (/  0.18,   0.33,  0.48,    0.63,  0.78/)
    
25. xpos2    = (/  0.15,   0.30,  0.45,    0.60,  0.75/)
    
26. markers1 = (/     16,      16,     16,       16,     16/)
    
27. 
    
28. 
    
29.    wks  = gsn_open_wks("pdf","Scatter")
    
30.     gsn_define_colormap(wks, "WhViBlGrYeOrRe")
    
31. 
    
32. ; 填色
    
33. ;************************************************
    
34.    res              = True
    
35.    res@gsnDraw      = False
    
36.    res@gsnFrame     = False
    
37.    res@tiMainString = "~F25~Obeserve Precipitation"
    
38.    res@tiMainFontHeightF    = 0.022   ;set main title fontsize
    
39.    
    
40.    res@mpMinLatF = 31
    
41.    res@mpMaxLatF = 37
    
42.    res@mpMinLonF = 110
    
43.    res@mpMaxLonF = 117
    
44.    res@mpFillOn    = False
    
45.    res@mpFillColor = False
    
46.    res@tmXBTickSpacingF = 1   ;坐标轴标注间距
    
47.    res@tmYLTickSpacingF = 1
    
48.    res@tmXBLabelFontHeightF = 0.016
    
49.    res@tmYLLabelFontHeightF = 0.016
    
50. 
    
51. 
    
52.     res@mpOutlineOn                  = True  ; Use outlines from shapefile
    
53.     res@mpOutlineBoundarySets        = "Geophysical"
    
54.     res@mpNationalLineColor          = "black"
    
55.     res@mpProvincialLineColor        = "black"
    
56.     res@mpGeophysicalLineColor       = "black"
    
57.     res@mpNationalLineThicknessF     = 2
    
58.     res@mpProvincialLineThicknessF   = 2.5
    
59.     res@mpGeophysicalLineThicknessF  = 2
    
60.     res@mpOutlineSpecifiers          = (/"China","China:States"/)
    
61.     res@mpDataBaseVersion            = "MediumRes"
    
62.     res@mpDataSetName                = "Earth..4"
    
63.         
    
64.         ;;;;色标设置
    
65.    ; res@lbLabelBarOn = True    ;显隐开关
    
66.    ; res@lbLabelFontHeightF = 0.25
    
67.    ; res@cnFillPalette="wh-bl-gr-ye-re"
    
68.    ; res@cnLabelBarEndStyle   = "ExcludeOuterBoxes" ; 色标两端标注风格，ExcludeOuterBoxes是去掉色标两端色块，IcludeMinMaxLabels是色标两端显示最大最小值，默认是色标两端无标注
    
69.    ; res@pmLabelBarWidthF     = 0.8      ;色标宽度
    
70.    ; res@pmLabelBarHeightF        = 0.08   ; 色标高度
    
71.    ; res@pmLabelBarOrthogonalPosF = 0.12  ;色标垂直偏移坐标，默认为0.02，正向上，负向下
    
72.    ; res@cnFillOpacityF          = 0.5  ;不透明度，默认为1
    
73.    res@pmTickMarkDisplayMode= "Always"
    
74.     res@lgOrientation          = "horizontal"     ; Default is vertical
    
75.   res@pmLegendWidthF         = 0.75             ; Make it wider
    
76.   res@pmLegendOrthogonalPosF = -0.1             ; Move it up slightly
    
77.   res@lgPerimOn              = False            ; Turn off the perimeter box
    
78. 
    
79. 
    
80.    ;;色标名字设定
    
81.    ; res@lbTitleOn            =  True
    
82.    ; res@lbTitleString        =  "mm"
    
83.    ; res@lbLabelFontHeightF   = 0.02;;;字体大小
    
84.    ; res@lbTitlePosition      = "Bottom"
    
85.            ; res@cnLevelSelectionMode =  "ExplicitLevels";等值线手动挡开始
    
86.         ; res@cnLevels             =  (/0,50,100,250,500,800/)        
    
87. 
    
88.    map = gsn_csm_map(wks,res)
    
89.    
    
90.    
    
91. 
    
92.    dum_fill   = new(5,graphic)
    
93.    mkres               = True
    
94.    mkres@gsMarkerThicknessF = 2.0      ; Twice as thick
    
95.    do i=0,4  ;;;6个区间
    
96.        ii = ind(pre1.le.A(i+1).and.pre1.gt.A(i))
    
97.        mkres@gsMarkerIndex = markers1(i)          ; Filled dots
    
98.        mkres@gsMarkerSizeF = sizes(i)
    
99.        mkres@gsMarkerColor = colors(i)
    
100. 
     
101.            txres               = True
     
102.            txres@txFontHeightF = 0.016             ;图例上文本的大小
     
103.        txres@txJust        = "CenterLeft"      ;文本放在图例的那个位置
     
104. 
     
105.            
     
106.            ;;;;;其实没什么用，操作两小时，回到中华人民共和国成立前，只需要对地图哪儿设置就行了。。。。但是我还是要保留他
     
107.            x=new(dimsizes(lon(ii)),float)
     
108.            x=lon(ii)
     
109.            y=new(dimsizes(lat(ii)),float)
     
110.            y=lat(ii)
     
111.            
     
112.            dd=ind(lon(ii).gt.110.0 .and. lon(ii).lt.116.5)
     
113.            bb=ind(lat(ii).lt.36.5 .and. lat(ii).gt.31.5)
     
114.                         ;print(x(dd))
     
115.                         ;print(in(lon(ii)))
     
116.                         ; if(all(lat(ii).lt.37.and.lat(ii).gt.31))
     
117.                 dum_fill(i)=gsn_add_polymarker(wks,map,x(dd),y(bb),mkres)
     
118.                 gsn_polymarker_ndc(wks,           xpos2(i),0.15,mkres)
     
119.                 gsn_text_ndc      (wks,labels(i),xpos3(i),0.15,txres)
     
120.                
     
121.                 delete(ii)
     
122.                 delete(dd)
     
123.                 delete(x)
     
124.                 delete(bb)
     
125.                 delete(y)
     
126.    end do
     
127.    
     
128. 
     
129.    
     
130.     draw(map)
     
131.     frame(wks)   
     
132. 
     

复制代码

ahtoffff

楼主  请问，当文件中缺测值不止一个的时候该怎么设置呀

森下下士

ahtoffff 发表于 2025-1-6 09:19
楼主  请问，当文件中缺测值不止一个的时候该怎么设置呀

把其他缺测值都换成一个缺测值即可

ahtoffff

森下下士 发表于 2025-1-6 09:56
把其他缺测值都换成一个缺测值即可

太聪明了 点醒我了 我怎么没想到从原数据下手。。。。 感谢！！