1
如何使由gsn_csm_contour_map創建的繪圖的邊緣平滑而不是使用那些藍色方塊?
我已經嘗試了不同的填充模式,沒有一個幫助。 在ncl中是否有一個活動輪廓函數?是可以創建一個空值和非空值之間的輪廓?
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"
;======================================================================
; The main code
;======================================================================
begin
;---Read desired data
sfile = addfile(f,"r")
var = sfile->var1
;---For zooming in on map
minlat = foo
maxlat = bar
minlon = foo1
maxlon = bar1
;---Get dimentions
dims = dimsizes(var)
nlev = dims(0)
time = var&time
date_str_i = getDate(time)
;---Set some resources
res = True
[email protected] = True
;[email protected] = "RasterFill"
;[email protected] =True
[email protected] = False
[email protected] = True
[email protected] = "ManualLevels"
[email protected] = -100
[email protected] = 3000
[email protected] = 200 ; 300 ; 50 ; 150
[email protected] = minlat
[email protected] = maxlat
[email protected] = minlon
[email protected] = maxlon
[email protected] = "HighRes"
[email protected] = True
[email protected] = 0.005
[email protected] = 0.001
[email protected] = (minlon+maxlon)*0.5
[email protected] = (minlat+maxlat)*0.5
[email protected] = "Always"
[email protected] = 0.01
[email protected] = False ; this is regional data
;---Loop across each level and plot to a different PNG file every time
do n=4,nlev-1
wks_type = "png"
[email protected] = 2000
[email protected] = 2000
wks = gsn_open_wks(wks_type,fname(0)+sprinti("%03i",n)) ;
[email protected] = "Time:" + date_str_i(n)
[email protected] = 0.010
plot = gsn_csm_contour_map_ce(wks,var(n,:,:),res)
delete(wks)
end do
end
但我有點懷疑它是代碼相關。
你應該補充一點,你已經習慣了得到這樣的結果你的問題的代碼。 –
謝謝,我會添加它。 –