ncl 中遇到的等值线平滑问题

Nalanalanaa

1. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
2. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"  
   
3. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
4. load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"
   
5. begin
   
6. fili0 = "D:/data/topograph.bin"
   
7. setfileoption("bin","ReadByteOrder","BigEndian")
   
8. topo =fbindirread(fili0,0,(/1080,2160/),"float")
   
9. 
   
10. lat = fspan(-90,90,1080)
    
11. lon = fspan(0,360,2160)
    
12. lat@units = "degrees_north"
    
13. lon@units = "degrees_east"
    
14. topo@units = "m"
    
15. topo@long_name = "topograph"
    
16. topo!0 = "latitude"
    
17. topo!1 = "longitude"
    
18. topo&latitude = lat
    
19. topo&longitude = lon
    
20. 
    
21. ;打开作图空间
    
22. wks  = gsn_open_wks("png","topo")
    
23. 
    
24. dres = True
    
25. ;dres@mpMinLatF = 34
    
26. ;dres@mpMaxLatF = 39
    
27. ;dres@mpMinLonF = 106
    
28. ;dres@mpMaxLonF = 112
    
29. dres@tmXTOn = False
    
30. dres@tmYROn = False
    
31. ;dres@tmXBLabelsOn = False
    
32. ;dres@tmYLLabelsOn = False
    
33. ;dres@gsnMajorLatSpacing = 2
    
34. ;dres@gsnMajorLonSpacing = 2
    
35. dres@gsnAddCyclic = False
    
36. dres@gsnMaximize = True
    
37. dres@cnFillOn = False
    
38. ;dres@cnFillPalette = "prcp_1"
    
39. dres@cnLinesOn = True
    
40. dres@cnLineLabelsOn = True
    
41. ;dres@cnLevelSelectionMode = "ExplicitLevels"
    
42. ;dres@cnLevels = (/100,500,1000,1500,2000,2500,3000,3500/)
    
43. dres@cnLevelSelectionMode = "ManualLevels"
    
44. dres@cnMinLevelValF = 0
    
45. dres@cnMaxLevelValF = 3500
    
46. dres@cnLevelSpacingF = 200
    
47. ;dres@cnFillMode = "RasterFill"
    
48. ;dres@cnRasterSmoothingOn = True
    
49. dres@lbLabelAutoStride = True
    
50. ;dres@lbLabelBarOn = False
    
51. dres@gsnLeftString  = " "
    
52. ;dres@gsnRightString  = " "
    
53. 
    
54. plot = gsn_csm_contour(wks,topo({35:39},{106:112}),dres)
    
55. 
    
56. end

复制代码

应该是平滑出问题了，试过几种方法都没做出来。有没有大神指点一下

某某蒙

我也遇到了这个问题，请问楼主解决了吗？

Nalanalanaa

某某蒙 发表于 2020-2-12 23:07
我也遇到了这个问题，请问楼主解决了吗？

解决了，

1. tres = True
   
2. tres@tmXTOn = False
   
3. tres@tmYROn = False
   
4. ;tres@tmXBLabelsOn = False
   
5. ;tres@tmYLLabelsOn = False
   
6. ;tres@gsnMajorLatSpacing = 2
   
7. ;tres@gsnMajorLonSpacing = 2
   
8. tres@gsnAddCyclic = False
   
9. ;tres@gsnMaximize = True
   
10. tres@cnFillOn = False
    
11. ;tres@cnFillPalette = "prcp_1"
    
12. tres@cnLinesOn = True
    
13. tres@cnLineLabelsOn = True
    
14. tres@cnInfoLabelOn = False
    
15. tres@cnLineLabelPlacementMode = "computed"
    
16. tres@cnLineThicknessF = 1
    
17. tres@cnLineLabelFontHeightF = 0.005
    
18. tres@cnLineLabelInterval = 1
    
19. tres@cnLineDashSegLenF = 0.3  
    
20. tres@cnLineColor = "SlateGray"
    
21. tres@cnLineLabelFontColor = "black"
    
22. tres@cnLineDashPattern = 7
    
23. ;tres@cnLevelSelectionMode = "ExplicitLevels"
    
24. ;tres@cnLevels = (/100,500,1000,1500,2000,2500,3000,3500/)
    
25. tres@cnLevelSelectionMode = "ManualLevels"
    
26. tres@cnMinLevelValF = 500
    
27. tres@cnMaxLevelValF = 3500
    
28. tres@cnLevelSpacingF = 300
    
29. ;tres@cnFillMode = "RasterFill"
    
30. ;tres@cnRasterSmoothingOn = True
    
31. tres@cnSmoothingDistanceF=0.01
    
32. tres@cnSmoothingOn = True
    
33. tres@cnSmoothingTensionF = 1
    
34. tres@lbLabelAutoStride = True
    
35. ;tres@lbLabelBarOn = False
    
36. tres@gsnLeftString  = " "
    
37. ;tres@gsnRightString  = " "

复制代码

调整了平滑设置

某某蒙

Nalanalanaa 发表于 2020-2-14 16:37
解决了，
调整了平滑设置

谢谢！！！

李头大

jiangbing

一维数据可以平滑吗?