IDL中contour和colorbar的使用技巧

/使用IDL繪制等值線和填充對應(yīng)的顏色【學(xué)習(xí)內(nèi)容】：等值線和畫法和顏色的填充主要程序：ＩDL8.０使用IＤL程序的圖形系統(tǒng)對數(shù)據(jù)極性成像化分析,研究ＩDL的畫圖方法和比較畫圖效果，得出比較滿意的圖畫.程序源代碼如下:proｔest_4；filｅ='M’+[’Ｘ'，'Y’，'Z＇，＇F＇,＇Ｉ’，'D'，’H＇］+’＿＇+＇2０1０＿0.0'＋'.dat’ｆile＝＇ＭX＿20１0＿0．0.ｄaｔ'fｏri=0，6dｏｂegin（fiｌe［i]）name=strmｉd（fiｌe[i］,0，7）ｄ1=ｄｅviｃe，ｄecomｐose=１xx=-89+findgen（１7９)*（１８0./１７9)ｙｙ＝—18０+findgen(３6１）ｎｌeｖels＝2１sｔep=(max（d1）－ｍin(d1)）/ｎleveｌｓvａls=findgｅn(nlｅvels)*ｓteｐ+mｉｎ（d1)ｖalｓ=[miｎ(ｄ1）—１,ｖals,max(d１)+１］xnaｍe＝[’－１８0’,＇—13５'，'—90’，＇—4５’，'0'，'４5＇,’９０’,’１35’,'1８０'］yｎamｅ=[’-90’，’—60＇,'-30'，＇０',’30','6０’，＇9０＇]p1=ｃoｎtoｕr(d1,yｙ，xｘ，ＲGB_tａble＝25，＄c_valuｅ=vals,＄xstyｌｅ=1，ysｔyle=1，$;N_leｖels＝N_eｌemeｎｔs(vａls)，$；設(shè)置填充層數(shù)pｏsitｉon=［0。10，0。20，0.8５，０。75］，$；設(shè)置圖像顯示位置,正交坐標(biāo)系【0,０，１,１】xrange=［-180,１80］,yrange=[—９0，9０］，＄;用來局部放大數(shù)據(jù)空間，即只顯示給定范圍的圖像planar=１，＄；設(shè)置顯示模式,0為三維，1為平面view_tiｔle=naｍe，＄xtｉtｌe＝’Ｌogitude’，$xminor=０，xticｋnamｅ＝xｎame,ｘtickleｎ=０。０5，＄ｙmｉnｏr=0,ytｉｃknａme=yname，ｙｔicklen＝０。05，＄ytitｌe＝'Lａt(yī)ｉtuｄe’，＄；設(shè)置標(biāo)題／fｉll）;填充模式ｃbaｒ=COＬORBAR(ｔｒaget＝ｐ1,ＯRIENTATIOＮ=1,＄fｏｎｔ_sｉze=7.5，POSＩTION=［0。86,0.20，0.9１，0。75］，tｅxtpos＝１）p2＝contouｒ（d1，yｙ,xx，c_value=vaｌs,/fｏｌlｏw,＄ｃ＿label_ｓhow=1,ｆont_sｉze=7．5,pｌanar=1，/ｏｖerplot）ｅndforend使用該程序可以對地球主磁場的7要素進(jìn)行繪圖處理,圖像如下：下面對其中的幾條主要的指令進(jìn)行說明：step=（max(d１)—min（d１））/nlｅveｌsvａlｓ=finｄgen(nｌｅvels）＊steｐ+min（ｄ1)vaｌs=［mｉn(ｄ１）—1，ｖａｌｓ,max（ｄ1）＋1]以上三條命令是設(shè)置等值線的線條數(shù)目和對應(yīng)的顏色數(shù)目,為了得到合理的等值線和條數(shù)和使用顏色填充是，對應(yīng)的顏色和等值線想吻合，在每兩條等值線之間都有對應(yīng)的顏色填充。自動(dòng)設(shè)置兩條等值線的間隔為ｓtep＝（maｘ（d１)－mｉn(ｄ1）)/nleｖels,這樣的話可以確保等值線的合理性和下一步設(shè)置等值線的數(shù)值是的對應(yīng)關(guān)系,vals=finｄgｅn(nｌｅvｅls）＊sｔep＋min(d1)，其中vals的數(shù)值和原始數(shù)據(jù)相對應(yīng)。由于IＤL填充顏色是從最小值開始填充,那樣的話在最大值和背景顏色之間就會(huì)出現(xiàn)一個(gè)顏色“空洞"，如:顯然這個(gè)很慢滿足要就要求,為了將所有等值線之間填充上顏色（最大值和背景色之間）,需要改變填充的最值,即vａlｓ=fｉndgen(nlｅveｌs)*steｐ+ｍin（ｄ１）的需要做出調(diào)整，經(jīng)過不斷的嘗試，選取vａls=［min（d1）-1，ｖａlｓ,ｍaｘ（d1）+1］,在原始的基礎(chǔ)上多加兩層等值線，確保數(shù)據(jù)范圍都在填充顏色范圍內(nèi)。結(jié)果如上圖（右),顯然重新設(shè)置后的填充效果要好很多。在contouｒ中使用RGＢ_TABLE＝2５，導(dǎo)入所需的顏色列表，c_valｕｅ=vaｌｓ用來設(shè)置繪制的等值線，ｘminor=0,xtickｎamｅ=ｘnamｅ,ｘtｉcklen=0．05，$ymiｎｏr=０，yticｋnａme＝y(tǒng)ｎame,yｔｉcklen=０．0５,＄用來設(shè)置坐標(biāo)軸的顯示信息，/fiｌl表明為顏色填充，若取消該關(guān)鍵字則不對圖像進(jìn)行顏色點(diǎn)沖，效果如下:ｃbａr=ＣOLORBAＲ(traｇet=p1，ORＩENTATIＯN=1,＄fｏnt_sizｅ=7。5,POSＩTION=[0。８6，0．2０,0。91，0．75］,ｔextpos＝1)用來設(shè)置色標(biāo)，ｏrｉentaｔion＝1為垂直方向,默認(rèn)值為0（水平方向）ｐ2=conｔｏuｒ（ｄ１,yy，xx，ｃ_ｖaｌue=vals,/fｏlｌoｗ,$c_label_shｏw=１，c_ｌｉnｅｓｔyｌe＝2,font＿size=7.5，planar＝１,/ｏverplot)使用／oｖerpｌoｔ的contｏur命令在顏色填充圖上繪制等值線，最終效果件上圖（右).當(dāng)然也可以只顯示其中的幾條等值線,用一下代碼取代相應(yīng)位置，p２=cｏntour(d1，yｙ，xx,c＿vaｌue=vals＊2,/foｌlow，＄c_lａbｅl_sｈｏw＝１,c＿ｌｉnestylｅ=２,fｏnｔ＿s