IDL中contour和colorbar的使用技巧

/使用IDL绘制等值线和填充对应的颜色【学习内容】：等值线和画法和颜色的填充主要程序：ＩDL8.０使用IＤL程序的图形系统对数据极性成像化分析,研究Ｉ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)，$；设置填充层数pｏsitｉon=［0。10，0。20，0.8５，０。75］，$；设置图像显示位置,正交坐标系【0,０，１,１】xrange=［-180,１80］,yrange=[—９0，9０］，＄;用来局部放大数据空间，即只显示给定范围的图像planar=１，＄；设置显示模式,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ｉtuｄe’，＄；设置标题／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要素进行绘图处理,图像如下：下面对其中的几条主要的指令进行说明：step=（max(d１)—min（d１））/nlｅveｌsvａlｓ=finｄgen(nｌｅvels）＊steｐ+min（ｄ1)vaｌs=［mｉn(ｄ１）—1，ｖａｌｓ,max（ｄ1）＋1]以上三条命令是设置等值线的线条数目和对应的颜色数目,为了得到合理的等值线和条数和使用颜色填充是，对应的颜色和等值线想吻合，在每两条等值线之间都有对应的颜色填充。自动设置两条等值线的间隔为ｓtep＝（maｘ（d１)－mｉn(ｄ1）)/nleｖels,这样的话可以确保等值线的合理性和下一步设置等值线的数值是的对应关系,vals=finｄgｅn(nｌｅvｅls）＊sｔep＋min(d1)，其中vals的数值和原始数据相对应。由于IＤL填充颜色是从最小值开始填充,那样的话在最大值和背景颜色之间就会出现一个颜色“空洞"，如:显然这个很慢满足要就要求,为了将所有等值线之间填充上颜色（最大值和背景色之间）,需要改变填充的最值,即vａlｓ=fｉndgen(nlｅveｌs)*steｐ+ｍin（ｄ１）的需要做出调整，经过不断的尝试，选取vａls=［min（d1）-1，ｖａlｓ,ｍaｘ（d1）+1］,在原始的基础上多加两层等值线，确保数据范围都在填充颜色范围内。结果如上图（右),显然重新设置后的填充效果要好很多。在contouｒ中使用RGＢ_TABLE＝2５，导入所需的颜色列表，c_valｕｅ=vaｌｓ用来设置绘制的等值线，ｘminor=0,xtickｎamｅ=ｘnamｅ,ｘtｉcklen=0．05，$ymiｎｏr=０，yticｋnａme＝yｎame,yｔｉcklen=０．0５,＄用来设置坐标轴的显示信息，/fiｌl表明为颜色填充，若取消该关键字则不对图像进行颜色点冲，效果如下:ｃ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)用来设置色标，ｏrｉentaｔion＝1为垂直方向,默认值为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命令在颜色填充图上绘制等值线，最终效果件上图（右).当然也可以只显示其中的几条等值线,用一下代码取代相应位置，p２=cｏntour(d1，yｙ，xx,c＿vaｌue=vals＊2,/foｌlow，＄c_lａbｅl_sｈｏw＝１,c＿ｌｉnestylｅ=２,fｏnｔ＿s