IDL程序设计笔记

1、目录Part one：文件读写操作 1.1、格式化输入与输出（read'print） 12、 ASCII文件读写2.3、二进制文件读写 5.4、图像格式文件读写 7.5、科学数据格式读写* 8.Part two： 图形绘制 1.1（一）plot过程绘制 11（二）plotg（）函数绘制12（三）散点图12（四）柱状图、条形图 15（五）绘制颜色条 colorbar18Part three: ENVI 二次开发201、自定义波段运算函数202、IDL中常见波段运算 213、ENVI classic创建自定义菜单 224、 ENVI classic+IDL 二次开发错误!未定义书签。Par

2、tone :文件读写操作1、格式化输入与输出（read'print）1）Read:该函数用于从键盘键入数据，默认数据类型为浮点型，若要输入其他类型 数据，需要先定义数据类型。Eg：IDL> read,口ecthp匸-* iWS® Ad.ar.alS : f请输入dateH : 123print / d旦匚乩123.000IDL> arr = 1 1IDL> read, arproir.pt = 】青f俞人占：rr值：* 请输入己工工值：the xainIDL> print, arrkisg the rain.2) Print:该函数用于将数据输出打印到

3、控制台，这里主要讲格式化控制输出format关键字。表01常用格式化控制符格式代码输出效果aN字符或者字符串按照 N个字符宽度输出fn.m浮点型数组按照N个字符宽度输出，小数点保留M位dn .mdouble型数组按照N个字符宽度输出，小数点保留M位en .m按照N个字符宽度的科学计数法输出，小数点后精确到M位Nx输出N个空格字符串/H直接引用字符串输出或者用H直接输出字符串c()用于输出日期数据Eg：IDL> print, ' LC7E CITY! I,fanrat=( (10)' LOVE CITY!IDL> print, 1 LC_?i 匚二】丫！ p 上口工盯

4、启匸二 $：吕)r LOVE CITY!IDL> print, 57B3, famiHt;=r i5 *67B9IDL> print, 6789, forrr.at-f (ilO) 1IDL7 print,520 > 520r fainaU = T (f" .3) F 520|lDL> vai = 7工print K vai , fannat=1 ( nalTae=*f iOS1 velue=00007IDL> vai - fLCVJ： CHENGDU!'IDt> print且：r” £Cir:niflt:=，(才仃比丄口已=*#

5、 了15) 1 value- LOVE 匚HENGDV!IDL> : Tf&58岀FtID1> print, format1rr2IT¥nF n匸HENGD口！)*LOVE CI T YCHEITSDU!IDL> print, forniac-t (1,LQVif 2Xf2x, rC-iENGDU! n) 1LOVE CITY CHENGDU!IDL> ;前系统吋间systiir.e 0IDL> pri ntr &ystime()Sun Jul 30 1:53:06 201?IDL> ;输出协调也齐时/UTC：IDL> pri

6、nt,systimeC/UTCSun Jul 30 11;5：15 2017IDL> ；输出秒制曲当前系统时间从、W0年1月1日算起IDLy pri ntf systitn« (11.50141Sae=+OCi92、ASCII文件读写常见ASCII文件：txt、matlab程序文件（*m）、c程序文件（*.c）、遥感影像头文件（如ENVI格式的头文件*.hdr）1）选择文件：fn = dialog_pick '选择 ASCII 文件：，get_path = work_dir）cd ,work_dir ;用于跳转到当前工作路径2）查询文件：（）函数Eg:查找当前工作路径下

7、所有的txt文件，并返回文件数目IDL> ad,'E:ENVIkIDL*IDL> file uujn=±i 丄已 search- txt1 f coEiT=f nums）IDL> print,fnams3IDL> print., file numzengsk*七x七 zengskOS rtxt IDL>*: （fname）函数用于查询文本文件的行数。3）打开文件：在IDL中读写ASCII码或者二进制文件，首先需要将文件与一个逻辑设备号 关联起来。语法：openr/openw/openu, lun, fname, /get_lunwidth=1O,

8、/append4）读写操作：IDL中利用readf函数读取文本文件、利用 printf函数写入ASCII文件 Readf, lu n, var01, var02 prin tf, lu n, var01, var025）文件关闭：在对文件操作完成之后，需要关闭文件的逻辑设备号，eg: freee_lun, lunEG:1 呷 J江）常规读取31-pro tes七 ASCIIa 唱fn - dialog pickfile1 SJfASCIIfr : '$5j gecath wark_dix)cd., work dir7varstrarr(5)openrrlunr fn,/get lun哲

9、readfr Lunr var10freelun11printsvar(012printrvar(113print,var(2)14print,var 15printrvar(4end'对于某种遥感数据，列数已知，行数若干，读取方法有两种 注循环读取数据pro test_A5CII±Ti = di al eg piekf lie (title= ': ' £tget path = woxk dixadtwozk dirmajcline = 1000trrp=iritarr (2)arr = intarr(2 rRAXlinelutl, / get l

10、tlfl，循环诱取每一行i = owhile eof (lun) rue 1 do begin readf lunr crnp air i = rrr.r丨=二丄endwhile data arr*,0;i-lprints r文件行数匸r r 1Xreelun,lun print,data erudM)利用"* line.国数获取文件行數pro testfn - dialog piakfile(title-'fifSASGIIft: r$ ,q e n_pd th = warkdii)cd,work dir|nl = file_lines (fn|data intarr(2r

11、nlopenr, lum f f n, j g已匸 lnreaafj Lur,freeun, Lunprintdftraend*如何获取某个文件你的列数：先读取一行数据，然后利用空格进行字符串拆分；最后利用n_eleme nts（计t mpArr-trgpli t(tmp.a(tmpArr)13算元素个数即为ns11I DL>IDL>丨DL>1DL>knL> :将字符串LOVE city! 1写入一个ASCH文件IDL> f ane=(ii a log pickfile (tile= f 文件保存为: ' r ge. path=work dirf f

12、 j.ltei= 1 * . txt' ) IDL> cd,vork_di rIDL> var =L0VE CITY!'IDL> openw,lunffn&me,/get_lunIDL> prirtf,lurifvarIDL> free_lun,lunIDL>对于某些遥感数据的头文件前面几行是说明文件，后一部分才是数据的，读取有两种方式：1）定义一个临时变量temp，将说明文字用该变量存储2）直接跳行：skip_lu n, lu n, 3, /lin es ;跳过文件的前三行3、二进制文件读写多数遥感数据的数据文件都是二进制文件。1)

13、 函数readu, lun, var01, var02 用于读取二进制文件* :对于某些数据文件开头是说明信息，读取时可以用point_lun, lun, position直接按字节跳过。Eg:读取IDL文件下的envi.img数据IDL> fnAmedialcgpiokrile(>IDL5> dit-i = bytarr ( G40r 400, & >IDL> openr, Lun, fnanier/get_lunIDA血匸nIDL> f ree_lcuif LunI DL> tielp, dataDATABYTE= Arxay40, 400

14、,IDL>2) 函数writeu, lun, var01, var02 用于写入二进制文件。Eg01:将上面读取的envi.img计算其NDVI并保存为二进制文件llDL> NIR=daza*F *F 3I七$1兴f *#2】：&L> NDVI-(float(1TIR)+IRIDL> help, NDVIKWIFLOAT= Array 640, 4D0IDL> o_fn = dialog_pickfile)openwr lunr o_fnr /GET_LU7II>L>- writau, lun, NDVItree l口口丄口11IDL>

15、|Eg02:将百个逋段數扳搜强段顺序保存対仔个浊亶的丈件PT O WTT 1 七0 RS;读取魏捉fn = dialoa_Dickf ile (title=ri?c-|xirfl- : 1 Sf get pa.th = work dixJg 严of cixdata = toytarr(640r4Q0f 6)npfinr, lunr fnf /ge*t_lunradu luiif d-a-afree lon lun，写文伴叵就一个艾件茗规fi3z = size(data)ns = sz 111 & nl sz (2) i nb = sz (3)prints ns f nl ni>匸

16、g = lndgen (IW) + JLfnarae - 'band1 4 string(rtr foxmat=,(i02)1)卡 1-datf;循秫保存文件for iO| nij-1 do bginopenwn Iunfxf /get_lunvri teu lun dataif roa luni inn ennrorJ"endEg03:;读取cha06下的风云02卫星数据pro read_AWXfile ='E:zengskIDLCha06dataFY2C TBB IR1 OTG 20061130 AOAD.AWX'ope nr , lu n, file,

17、/get_lu n;跳过前20个字节poin t_lu n,lu n, 20Headline =indgen ( 3);定义三个整型数据，分别读取记录长度、头文件记录数、数据记录数readu , lun, Headlineprint ,headli ne 0, headli ne1 ,headli ne2;定位到第58个字节，读取数据日期poin t_lu n, lu n, 58Begindate =indgen (5);记录开始日期readu , lun, BegindateEnddate =indgen (5);记录结束日期readu , lu n, En ddatelat_long =i

18、ndgen ( 4);记录网格左上和右上角的经纬度readu , lu n, lat_ longpri nt,Beg in datepr int ,En ddatepri nt,lat_ long;基于Headline信息数组，定义字节数组，其行列数分别为记录长度和数据记录数data = bytarr (Headl in e2,Headli ne0); 定位到数据部分，其头文件结束的字节位置为 “头文件记录数 * 记录长度 point_lun, lun, Headline0 * Headline1print ,Headline0 * Headline1; 读取数据部分readu , lun,

19、datahelp ,data; 显示数据window , 0, xsize =500 , ysize =500tv , congrid (data, 500 , 500 ) end4、图像格式文件读写IDL 自带了丰富的多种图像读写函数， 如 BMP，JPG，PNG, JPEG2000， DICOM ， TIFF 等相关函数： 文件查询 ： query_image（fname, dimensions=dimensions, chann els=nb）;cha nn els返回波段数、dime nsio ns返回图像的行列数读取图像：read_image （ fname ）;结果返回图像数组写入

20、图像：write_image,fname可带路径）, jpg/bmp/tif，data,/order ;order关键字用于设置图像的纵坐 标从上往下算起，默认为从下往上。*注：I、以上三个函数对于IDL支持的图像格式都适用。II、显示图像时：tv用法：tv, data, true = 1/2/3其中关键字true = 1表示 数据格式为（ 3, m,n）、 true = 2 表示数组格式为（ m, 3, n）、 true = 3 表示格式为 （m, n,3）III、tvrd:屏幕拷贝函数，返回当前直接图形窗口或设备的指定矩形部分 的内容。IV、对遥感图像而言，其数据存放方式有 BSQBIPB

21、IL，以一个 三维m列n行的遥感数据为例：BSQ的表示方式为m, n, 3、BIP表达方式为3, m ,n、BIL的表达方式为m, 3, n，常用的图像文件一般是BIP方式存储，遥感 文件一般是 BSQ。1)JPEG格式：读写 write_jpegEg01:IDL> jpgf lie *E'.IDL'11- Cha j6 - data' rose.'IDL> :查看图像倍亀IDL> result = query image (jpgfile, dimensiQus = dzLmensiQns, channel3 = rLb) 韦 Loaded

22、DLH： JPEG.idl> print,result1IDL> print, dimensions, nb22714? IDL 0- XIDL> ；用reaa_LiEage'LyIDL>= 丁已.-.1_:. -j-if： (jpgfile)IDL> hfllp,dataOlDA.TAO1BYTE- ArrayfS, 227r H9IDL> :利用r&ad_-peg函数i卖取IDL> read_jpeg, ipgfile, ggIDL> help, imgIMGBYTE= Array3, 227, 149IDL> ；显示图

23、借IDL> .indL-w 0r xsize 227f ysize - 149IDL> tv, uJf nr nt = 1IDL>Eg02:2) BMP格式：(不采用任何压缩)，数据存储量一般很大读：read_imge(fname) or read_bmp(fname, /rgb) 写: write_bmp, fname, data, /rgb3) TIFF 格式：TIFF可以存储多波段图像，还可以包含投影信息，如Iandsat-7就是直接用.tiff存储的。Eg01:IDL> riff = PE:VENVIklDLVeast.tifpI DL> res.11 =

24、 tjuery image (ti壬壬尸 dzLrDen5j.ons=diineTi3iQns )IDL> print j41354135IDL> data = read_tiff(tiff)IDL> h.elp, dataDATABYTE- Array3, 4133t 413S工DLA ;设魁羔键字匸侯耳2指点IK域囱像IDL> 3ui>_daca = read_tiff(匸丄f工尸uto_recr=(33,43,iodofiodoIDL> IielPtSITBDATA - BYTE= Array3, 10C0r 1000JTDt> ".!

25、 iido'Yf 0f xsLze = 1000r ysixe = 12 C 匚IDL>3Ub_d3LZaf true = 15、科学数据格式读写*IDL中支持的科学数据格式包括CDF、HDF、HDF5、HDF-EOS、NetCDF 等等。1)HDF4 文件：目前国内外的多种卫星传感器都是将 HDF 作为标准数据格式，包括 EOS/MODIS、EOS/OMI、HJ-01/HSI、FY-3/MERSI。EG01:；读取E:ENVIIDL 下的MODIS数据的经纬度以及其 1KM的反射率数据集 pro test_hdffname= dialog_pick = ' 选择数据文件

26、： ' ,get_path=cur_dir) cd , cur_dir hdf_id= hdf_sd_start (fname) ; 打开 hdf 文件，返回一个文件 idhdf_sd_ ,hdf_id,sd_nums,attribute;用于获取HDF文件的数据集数目和属性数目print ,sd_nums,attribute ; 读取经度lat_index= hdf_sd_nametoindex (hdf_id, 'Latitude' ) ; 根据数据集名称来获取数据 集的索引号lat_id= 的idhdf_sd_select(hdf_id,lat_index); 利

27、用索引号选择数据集，返回一个数据集hdf_sd_getinfo,lat_id,name=name,unit=unit; 查询已打开的数据集的基本信息print ,name, " " , unit hdf_sd_getdata ,lat_id,lat; 读取数据，参数lat 用于返回读取的数据hdf_sd_endaccess ,lat_id ; 关闭已经打开的数据集 help ,lat; 读取纬度lon_index= hdf_sd_nametoindex(hdf_id, 'Longitude' )lon_id= hdf_sd_select(hdf_id,lon

28、_index)hdf_sd_getinfo print ,name,I! I!,lon_id,name=name,unit=unit, unithdf_sd_getdata hdf_sd_endaccess,lon_id,lon,lon_idhelp ,lon;读取1KM反射率数据ref_index= hdf_sd_nametoindex(hdf_id, 'EV_1KM_RefSB' )ref_id= hdf_sd_select(hdf_id,ref_index)hdf_sd_getinfo,ref_id,name=name,unit=unitprint ,name, hdf_

29、sd_getdata, unit,ref_id,ref_valuehdf_sd_endaccess ,ref_id help ,ref_valuehdf_sd_end , hdf_id ;关闭打开的文件 end运行效果冬 CarrpiLed mod口1_已：TEST HDf.313日degreesLATFLOATLaajitude degreesLOWFLOATEV_lKM_R&fS2 neneREF_VJO.nEUIMTIDL>=Jlrray271, 40£=Array271, 4D|= Array13S, 2030, IS2) *h5、*hdf5、*he5 格式的读

30、取：J；下面介绍一个读取环境卫星*h5数据格式的例子pro read_hdf5fname= dialog_pick ='选择数据文件：',get_path=cur_dir)cd , cur_dir;打开一个图形用户界面查看h5文件result= h5_browser (fn ame);h5_message=h5_parse(fname);查询文件的基本信息，返回一个结构体变量;prin t,h5_messageh5_id= h5f_open (fname);打开一个h5文件，返回一个文件的idsdd= h5d_open (h5_id,'ImageData/BandDat

31、a');打开相应的数据集，数据集名称可带路径Bandata= h5d_read (sdd);读取数据;关闭文件h5d_close ,sd_idh5f_close ,h5_idhelp ,Bandata;将Ban data 转成BIP格式存储Bandata_BIP =transpose(Bandata, 2, 0, 1 )help ,Bandata_BIPwrite_tiff , 'E:ENVIIDLHJ-1.tif', Bandata_BIPend2）netCDF文件读取：n etCDF（ network Com mon Data Form,网络通用数据格式），常用于气

32、象科 学数据存储。;下面介绍一个读取ECMWF数据的例子，ECMWF为*.nc格式数据= ' 选择数据文件： ' ,get_path=cur_dir)pro test_nc fname= dialog_pick cd , cur_dir nc_id= ncdf_open (fname) ; 打开 netCDF 文件nc_message= ncdf_inquire (nc_id) ; 对打开的文件进行查询，返回文件基本信息 help ,nc_message; 变量查询; 后一个参数为变量索引号或者名称 ; 这里我读取的是它的第四个变量 'V10' var_mess

33、age= ncdf_varinq (nc_id, 4)help ,var_messagevar_id= ncdf_varid (nc_id, 'v10' ) ; 根据变量名称获取对应变量的索引号 print ,var_id;ncdf_varget 获取数据 ncdf_varget ,nc_id,var_id,value; 文件关闭ncdf_close , nc_idhelp ,value end4) Grid 格式读取：Part two ：图形绘制一) plot 过程绘制1、plot 过程： 1) plot, xdata, ydata, /nodata2)oplot, xdat

34、a, ydata ;用于在现有窗口中添加新曲线2、线形符号设置：1) 关键字psym修改数据点的符号，直方图线形psym=10, psym= -a表示除显示各数据点的符号外，还将各数据点连接起来。2) linestyle 关键字设置曲线线形， 0 为实线、 2为虚线、3 为点 划线。3、坐标轴设置：xytitle 设置标题； xycharsize 设置字体大小；xyrange 用于设置坐标轴范围， 注意加上 xystyle=1 来强行设定 坐标轴范围xyticks 和 xyminor 分别设置主刻度和最小刻度间隔数目 xytickname 设置刻度的名称4、颜色设置:color: 设置图像颜色

35、；一般用十六进制来表示颜色； 'FFFFFF'xl 表示白色；'000000'为黑色backgro und:设置背景颜色5、添加标注:xyouts过程用于在图像窗口中添加标注信息。示例: xyouts, 0.18, 0.80,RMSE'= , color= ' FFFFFF' xl, charsize=1.2,/normal6、图形保存为文件：一般方法为利用 TVRD 函数拷贝图形窗口的内容，在写入图像文件。img = tvrd(x, y, channel=value, true=123, /order)(二) plotg()函数绘制1

36、、 plot 函 数 ： curve=plot(x, y, /buffer, /current, dimensions=width, height, margin=num, title= ' ' , nam/oev=er'plot', w, indow_title=' ' , /nodata)*notice： /buffer 指将图形保存在缓存中/current 设置在当前窗口中绘制图形dime nsions设置窗口的大小title设置图形标题，name设置图像对象的名称*函数 plot 的返回结果为一个图像对象，具有多种对象和方法2、符号、线形

37、设置： symbol、linestyle、thick 用于改变曲线线宽3、坐标轴设置：见 plot 过程4、添加标注：函数text用于在图形窗口中添加标注信息，结果返回一个图形对象Eg: label=text(x, y, 标注内容 ' , target=curve, color='red', font_size=1.2)5、添加图例：函数lege nd用于在图形窗口中添加标注信息，该函数只适用于图形对象(三) 散点图Eg01 :读取文本文件中的数据并利用plot过程绘制相应的散点图 ; 利用过程 plot 绘制散点图pro plot_scatter.*读取数据*fnam

38、e= dialog_pick cd , work_dir选择数据:, get_path=work_dir)nl=(fname) print , nltemp_var=data= fltarr( 2, nl- 1)openr ,lun,fname,/get_lun result= fstat (lun) help ,result readf , lun, temp_var readf , lun, data free_lun, lunprint ,temp_varprint ,data.*绘制散点图*x = data0,* ; 观测数据y = data 1,* ; 估算数据plot , x, y

39、, psym=2 ,xrange=5, 15 ,yrange= 5, 15,$xtitle= 'Observed value', ytitle='Estimated value' , $background= 'FFFFFF'xl , color= 'oooooo'xl , $ xstyle= 1, ystyle=1, $charsize= 1.3 , /nodataoplot , x, y, psym=5, color= '000000'xl; 绘制基准线x2= 5, 15 & y2= 5, 15oplo

40、t , x2, y2, color='FF0000'xl.*计算MAE与 RMSE并在散点图中添加标注*MAE= mean ( abs (x-y) ;MAE 平均绝对误差RMSE= sqrt ( mean (x-yF2);均方根误差MAE_label= 'MAE=' +string (MAE, format= '(f5.2)')RMSE_label= 'RMSE=' +string (RMSE, format= '(f5.2)'); 添加标注xyouts , 6, 13.0 , MAE_label, color=

41、xyouts , 6, 12.3 , RMSE_label, color='000000'xl, charsize='000000'xl, charsize=1.2 , /data1.2 , /data; 保存散点图为图像文件，格式为 .pngimg= tvrd (true= 1);法o_fn= dialog_pick ' 文件保存为： ' ) ; 键入的文件名要包含后缀名 write_png , o_fn, img;法;write_png, 'E:ENVIscatter.png', img endEg02读取文本文件中的数据并利

42、用plot函数绘制相应的散点图;plot 函数绘制散点图pro plot_scatter02; 读取数据（title= ' 选择数据： ' ,get_path=work_dir） cd , work_dirvar=I!nl=(file)data= fltarr( 2, nl- 1)openr , lun, file, /get_lunreadf , lun, varreadf , lun, datafree_lun , lunprint , data; 绘图ob=data0, *Es=data1, *graphic01= plot (ob, Es, xrang= 7, 15 ,

43、yrang= 7, 15, xminor= 5, yminor= 5, $ ; 先数据title= 'scatter' , xtitle= 'Observed Value' , ytitle= 'Estimated Value' , $ ; 再坐标轴$ ; 然后symbol= 'X' , sym_size= 1.2 , linestyle= 'none' , color= 'black' 线形margin= 0.1 )x2= 7,15 & y2= 7, 15graphic02= plot (

44、x2, y2, linestyle= 1, /overplot); 计算平均绝对误差和均方根误差并添加标注 MAE= mean ( abs (Es-ob) ;RMSE= sqrt ( mean (Es-ob2)MAE_label= 'MAE=' +string(MAE, format='(f5.2)')RMSE_label= 'RMSE=' +string (RMSE, format= '(f5.2)')t01= text ( 0.15 , 0.8 , MAE_label, target=graphic01, font_size=

45、 12) ;fonts_size 用于设置标注的字体大小t01= text ( 0.15 , 0.75 , RMSE_label, target=graphic01, font_size=; 保存文件o_fn= dialog_pick ' 图形保存为： ' ) graphic01.save , o_fngraphic02.save , o_fnend(四) 柱状图、条形图EgO1: *利用bar_plot()函数绘制柱状图(不建议用直接图形法来绘制) *利用error_plot添加误差线;barplot 功能函数绘制柱状图pro test_histogramO2;# 读取文本文

46、件数据 #(title= ' 选择文本文件： ' , get_path=work_dir) cd , work_dirnl=(file)var= strarr ( 1,nl) openr , lun, file, /get_lun readf , lun, var free_lun , lun print , var ; 拆分area= strarr (nl) number= intarr (nl) for i= O, nl- 1dobeginareai= strmid (vari, O, 8)numberi= fix ( strmid (vari, 11, 3) endfor

47、 help , area print , area help , number print , number 3 ;# 画柱状图 #baselines= replicate ( 65 ,nl)His= barplot (number, $bottom_values=baselines, xticklen=O, xrange=- O.5 , 5.5 ,xtickname=area, $yrange= 65, 1OO , title= 'Histogram' , ytitle= 'Impact Facter' , $ width= O.5 , dimensions=

48、 9OO , 5OO , $bottom_color= 'EEEEEE'xl , fill_color= 'AOAOFF'xl , $ margin= O.2 );# 绘制误差线 #error= 1.7 ,2.7 ,2.5 ,2.4 ,1.2 ,1.0 error_lin es= errorplot(n umber, error, li nestyle=6,errorbar_capsize=0.2 , $errorbar_color='blue', /overplot);#文件保存为 jpg#o_fn= dialog_pick '图片保存

49、为：)+ '.jpg'注意后面加一个后缀名字符串更加方便His. save , o_fnerror_lines. save , o_fn end运行结果:chmodu kiihou站神時 yibmg ntnchonq 出冋営HliWgrfrnWB5M 門 益34 W8-E-wieaDle¥軒jiEg02: *在利用bar_plot()函数绘制图像时，添加/horizontal关键字用于将图像横 过来，即绘制条形图J;bar_plot函数绘制条形图 chartJpro test_chart;# 读取文本文件数据 #(title= '选择文本文件：',ge

50、t_path=work_dir)cd , work_dirnl=(file)var= strarr ( 1, nl)openr , lun, file, /get_lunreadf , lu n, varfree_lun , lu npri nt , var;拆分area= strarr (nl)n umber=in tarr(nl)for i= 0, nl- 1 dobeginareai= strmid (vari, 0,8)numberi= fix ( strmid (vari,11,3)en dforhelp , areapri nt , areahelp , numberpr int

51、, n umber3；#画柱状图 #baselines= replicate ( 65 ,nl)His= barplot (n umber, bottom_values=baseli nes, $yticklen=0, yrange=-0.5 , 5.5 , ytickname=area, $xrange= 65,100 , title= 'Histogram' , xtitle= 'Impact Facter', $width= 0.5 , dimensions=700 , 500 , $bottom_color= 'EEEEEE'xl , f

52、ill_color='AOAOFF'xl , $margin= 0.2 , /horizontal);#文件保存为 jpg#o_fn= dialog_pick '图片保存为：)+ '.jpg'注意后面加一个后缀名字符串更加方便 His. save , o_fnEnd运行结果：suinjngnanchongyibmgmianyatigluzhouchengdu6S 707580859095100Impact FacterRistoqrfim（五）绘制颜色条colorbarEg：1) 创建一个快速可视化窗口IDL> w = window (dimens

53、ions=600, 400)2) ;创建一个Colorbar，有自定义颜色表，不设置标注位置。IDL> cb = colorbar (POSITION = 0.05 ,0.85 , 0.7 , 0.9 , RGB_TABLE= 72)0501001502002503) ; 创建一个Colorbar具有自定义颜色表，设置标标注在顶部 texpos=1 。IDL> cb = colorbar (POSITION = 0.05 , 0.85 , 0.7 , 0.9 ,RGB_TABLE= 72,TEXTPOS= 1)4) ;创建Colorbar具有自定义颜色表，设置自定义范围自定义范围r

54、ange=*,*IDL> cb = colorbar (POSITION = 0.05 , 0.75 , 0.7 , 0.8 , RGB_TABLE =72 ,RANGE = 0,1)0002040608105) ; Colorbar具有自定义范围和自定义标签tick name= A array。IDL> cb = COLORBAR(POSITION= 0.05 , 0.55 , 0.7 , 0.6 , RGB_TABLE= 72 ,RANGE= 500 , 1000 ,TICKVALUES= 550 , 650 , 750 , 850 , 950 ,TICKNAME= '

55、A' , B , C , 'D' , 'E'A9C06) ;创建一个 colorbar,设置垂直色条orientation=1，范围 0-100range，锥形端tapar=1 (左图)。IDL> cb = COLORBAR(POSITION= 0.8 , 0.2 , 0.85 , 0.9 , TAPER= 1 , RGB_TABLE= 72 , RANGE= 0, 100 , ORIENTATION= 1)*;离散色彩条垂直方向；标注在右边一一textpos=1 ，添加轮廓线一一BORDER=1 (见右 图)。IDL> cb = COLORBAR(POSITION= 0.9 , 0.15 , 0.95 , 0.95 , RGB_TABLE=rgb,IDL> TICKNAME=tick name, BORDER=1 , TAPER= 1 , $> ORIENTATION= 1 , TEXTPOS= 1)100BO60