2023年IDL程序设计笔记.pdf

《2023年IDL程序设计笔记.pdf》由会员分享，可在线阅读，更多相关《2023年IDL程序设计笔记.pdf（32页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、目录Part one：文件读写操作.11、格式化输入与输出(readprint).12、ASCH文件读写.33、二进制文件读写.64、图像格式文件读写.95、科学数据格式读写*.11Part two：图形绘制.15()plot过程绘制.15(二)plotg()函数绘制.16(三)散点图.17(四)柱状图、条形图.21(五)绘制颜色条 colorbar.26Part three:ENVI 二次开发.291、自定义波段运算函数.292、IDL中常见波段运算.303、ENVI classic创建自定义菜单.314、ENVI classic+IDL二次开发.错误!未定义书签。Par t on e:文献

2、读写操作1、格式化输入与输出(reach p r int)1)R e a d:该函数用于从键盘键入数据，默认数据类型为浮点型，若要输入其他类型数据，需要先定义数据类型。Eg:IDL re a d,d a ta,prom pt=请生俞入d a ta值：1请愉入d a ta值：123IDL p r i n t,d a ta123.000IDL a r r =,IDL re a d,a r r,p ro叩 七=请谕入a r工值：请输入a r r值：k is s th e ra inIDL p rin ts a r rk is s th e ra in2)Print：该函数用于将数据输出打印到控制台,

3、这里重要讲格式化控制输出fo rm at关键字。表0 1 常用格式化控制符格式代码输出效果aN字符或者字符串按照N 个字符宽度输出fn.m浮点型数组按照N 个字符宽度输出，小数点保存M 位dn.mdou b Ie型数组按照N 个字符宽度输出，小数点保存M 位e n.m按照N 个字符宽度的科学计数法输出，小数点后精确到M位Nx输出N 个空格字符串/H直接引用字符串输出或者用H 直接输出字符串c()用于输出日期数据gE g:IDL p r in七,LOVE CITY!.fo rm a t*(a lO)(LOVE CITY!IDL p r i n t,*LOVE CITY!,form at=(a 1

4、 5)1LOVE CITY!IDL p r i n t,6789rform at=*(i5)6789IDL p r in t,6789,f o r m a l (ilO),6789IDL p r i n t,5 2 0.520z form at=(7.3)520.520IDL v ar=7IDL p r in ts v a rrform at=(nv alu e=n,i0 5),value=00007IDL v ar=LOVE CHENGDU!,IDL p r in ts v a rrform at=(nv alu e=n,a!5),value=LOVE CHENGDU!IDL；空格输出NxI

5、DL p r in t,fo rm at=,(LOVE1*,nCITYnr CHENGDU!w)LOVECITYCHENGDU!IDL p r in t,f o r m a t1(nLOVEnF 2xr CITY,2xr CHENGDU!0)LOVE CITY CHENGDU!IDL;输出当前系统时间sy stim e()IDL p rin ts s y s tim e()Sun J u l 30 19:53:06 2017IDL；输出协调世界时/UTCIDL p rin ts s y s tim e(/UTC)Sun J u l 30 11:54:15 2017IDL；输出秒制由当前系统时间

6、，从1970年1月1日算起IDL p r in t,sy stim e(1)1.5014158e+0092、ASCII文献读写常 见A S C II文献:tx t、mat 1 ab程序文献(*.m)、c程 序 文 献(*、遥感影像头文献(如ENVI格式的头文献*.hd r)1)选择文献：f n=d i a l o g _ p i c k f i l e (t i t le=选择AS C I 工文献：，g e t _ p at h =wo r k _ d i r)cd,w o r k _ d i r ;用于跳转到当前工作途径2)查询文献:file_se a rch()函数Eg:查找当前工作途径下

7、所有的txt文献，并返回文献数目IDL c d,E:E N V IID L IDL file_ n u m=file _ s e a rc h(1*.t x t ,count=fnum s)IDL p rin ts fnums3IDL p rin ts file_numze n g sk.tx r zengskO l,tx t zengsk02,tx tIDL*：Fi 1 e_ lines(f n a m e)函数用于查询文本文献的行数。3)打开文献：在I D L中读写ASC I I码或者二进制文献，一方面需要将文献与一个逻辑设备号关联起来。语法:openr/op e n w/o penu,1

8、 u n,fname,/getu n,w i d t h=1 0,/a p p e nd4)读写操作：I D L中运用read f函数读取文本文献、运用p r intf函数写入A SC II文献 Readf,1 u n,var 0 1,v a r02.(2)p r i n t f,1 un,v a r 01,v a r 02.5)文献关闭：在对文献操作完毕之后，需要关闭文献的逻辑设备号,e g:fre e eun,l uEG:2;1)常规读取3-pro test_ASCIIfn=dialog_pickfile 选择 ASCII 文件：$,get_path=work_dir)6 cdzwork_

9、dirvar=strarr(5)openrzlun,tn,/get_lun9 readfzlunf var10 free_lunr lun11 printsvar(0)12 print,var(1)13 prints var(2)14 print,var(3)15 print,var(4)16 end;对于某种遥感数据，列数已知，行数若干，读取方法有两种;1)循环读取数据p r o t e s t _ AS C I If n=d i al o g_ p i c k f i l e(t i t l e=，选择 AS C I I 文件：$,ge t _ p at h =w o r k _ d i

10、r)c drw o r k _ d i rm ax l i ne =1 0 0 0r m p=i nt ar r(2)ar r =i nt ar r(2 fm ax l i ne)o p e nrr l u n,f n,/ge t _ _ l x i n;循环读取每一行i =0w h i l e e o f(l u n)ne 1 d o b e gi nr e ad f r lun,t m par r *r i=七 即i =i +1e nd w h i l ed at a=ar r*r 0:i-1 p r i nt,，文件行数：f r e e _ l u nr l u np r i nt,d

11、at ae nd；#*#*#*#*#*#;2)利用f i l e _l i n e s函数获取文件行数p r o t e s t _A SC I If n =d i al o g _p i c k i 1e ，选择A SC工工文件：$,g e t _p at h =w o r k _d i r)c dr w o r k _d i r|n l =f i l e _l i n e s(f n)d at a=i n t ar r(2r n l)o p e n r,l u n,f n,/g e t _l i mr e ad frl u n,d at af r e e _l u nzl u np r i

12、 n t,d at ae n d*如何获取某个文献你的列数:先读取一行数据，然后运用空格进行字符串拆分;最后运用n_elem e nt s()计算元素个数即为n sI D L 511I D L I D L I D L-X C y -6-7-8-9 1012t m p A r r=s i t r s p l i t (t m pz*)0 3=3 2 5 5 0 (t m p A r r)n s13IDL；将字符串OVE CITY!,写入一个ASCI工文件IDL fnam 8=dialog_jpickile 文件保存为：1,g et_ p ath=w o rk _ d irz filte r=*.

13、zx z )IDL c d,w ork_dirIDL v a r=LOVE CITY!,IDL openw,lu n,fn am e,/g e j l u nIDL p r i n t frlu nr v arIDL fr e e _ lu n,lu nIDL对于某些遥感数据的头文献前面几行是说明文献，后一部分才是数据的，读取有两种方式：1)定义一个临时变量t e m p,将说明文字用该变量存储2)直接跳行:sk i pun,1 un,3,/l i n e s；跳过文献的前三行3、二进制文献读写多数遥感数据的数据文献都是二进制文献。1)函数readu,lun,varOl,v a r 0 2.用

14、于读取二进制文献*：对于某些数据文献开头是说明信息，读取时可以用P o in t J u n,lun,pos i t ion直接按字节跳过。Eg:读取IDL文献下的envi.img数据IDL fnaine=dialog_pick.f i l e ()IDL d a ta =b y t a r r(640r 400r 6)IDL o p en rzlu nz fnam er/g e t_ ltmIDL re a d urlu n,d a taIDL fre e _ lu nzlu nIDL h e lpz d a taDATA BYTE=A rray 6 4 0r 400,6IDL2)函数wri

15、t e u,lun,varOl,v ar02用于写入二进制文献。Eg 0 1 ：将上面读取的envi.img计算其ND V I并保存为二进制文献IDL NIR=data*r*r3IDL IR=data*r*r 2IDL NDVI=(float(NIR)-IR)/(float(NIR)+IR)IDL help,NDVINDVI FLOAT=Array640,400IDL o_fn=dialog_pickfile()IDL openw,lun,o_fnr/GET_LUNIDL writeur lun z NDVIIDL free_lunrlunIDL|Eg02：；将6个波段数据按波段顺序保存为6个

16、独立的文件pro w rite RS;读取数据fn=d ia lo g _ p ic k f i l e (七，选择文件：1$tg et_ p ath =w ork_dir)c d,w ork_dird ata=b y ta r r(6 4 0,400r 6)openr,lu n,fn,/g e t_ lu nreadu,lu n,d atafre e lu n,lun；写文件;生成一个文件名数组sz-s i z e(data)ns=s z(1)&n l=sz(2)&nb=s z(3)p r i n t,n s,n l,n bf t =in d g en(nb)+1fname-band+s t

17、 r in g(f tr form at=*(i0 2),)+.d a t1:循环保存文件fo r i=0,nb-1 do b eginopenw,lu n,fnante i ,/g e t_ lu nw rite uz lu n,d a ta*,*rifre e lu n,lunendforj-endEgO3:;读取cha06下的风云02卫星数据p ro read_AWxf i l e =E:zen g sk ID L C h a0 6 d ataF Y 2 C _ T B B _ lR l_ O T G _ 2 0 2 3 1130_AOAD.AWXopenr,1 un,f i l e,/

18、get_ 1 un;跳过前20个字节poin t _1 un,1 u n,20H ead lin e=indgen(3);定义三个整型数据，分别读取记录长度、头文献记录数、数据记录数re a du,1 un,H eadlin ep rin t,he a d li n e 0 z head 1 i ne 1,headl i ne 2;定位到第58个字节，读取数据日期p o i n t _1 u n,1 un,5 8Be g inda t e=indgen(5)；记录开始 口 期re a d u,lu n,B egindateE nddate=indgen(5);记录结束日期readu,1 u n

19、,Enddate1 a t _ Ion g=indgen(4)：记录网格左上和右上角的经纬度re a d u,1 u n z 1 a t _1 o n gp r in t,B eg i ndatep r in t,Enddatep rin t,1 at_ 1 ong;基于Headl in e信息数组，定义字节数组，其行列数分别为记录长度和数据记录数d a ta =b y t a r r (Headline 2 A He a d lin e 0 )；定位到数据部分，其头文献结束的字节位置为、头文献记录数,记录长度p o in t_ l u n,lu n,H eadl i ne 0*H e a d

20、 li n e 1 p r i n t,H eadlin e o*Hea d 1 i ne 1;读取数据部分r e a du,1 unz datah e lp,data;显示数据w indow,0,x si z e=50 0,y s ize=50 0tvz c on g r i d(data,500,50 0)end4、图像格式文献读写M DL自带了丰富的多种图像读写函数，如BMP,JPG,PNG,JPEG2023,DICOM,TIF F 等相关函数:文献查询：query_ ima g e(f n a me,dime n sio n s=d i me ns i o ns,chann e Is=

21、n b);ch a n nels 返回波段数、d im e n sion s返回图像的行列数。读取图像：read_ i m a ge(fname);结 果 返 回 图 像 数 组 写入图像：write_ i ma g e,fname(可带途径)，jpg/bm p/tiff,d ata,/o rd e r;ord e r关键字用于设立图像的纵坐标从上往下算起，默认为从下往上。*注：I、以上三个函数对于IDL支持的图像格式都合用。I I、显示图像时:t v用法:tv,data,tr u e=1/2/3;其中关键字true=1表达数据格式为(3,m,n)、true=2表达数组格式为(m,3,n)、t

22、 r ue=3表达格式为(m,n,3)III、tvrd:屏幕拷贝函数，返回当前直接图形窗口或设备的指定矩形部分的内容。IV、对遥感图像而言，其数据存放方式有BSQBIPBIL,以一个三维m列n行的遥感数据为例：B SQ的表达方式为 m,n,3、B IP表达方式为 3,m,n、B IL的表达方式为 m,3,n,常用的图像文献一般是B IP方式存储，遥感文献一般是BSQ。1)JPEG 格式：读 readjp e g 写 w r itej p egEgO 1:IDL jp g f ile =*E:z en g slcID L C h a0 6 d ataro se.jp g*IDL;查看图像信息ID

23、L r e s u lt=query_im age(jp g file,dim ensions=dim ensions,channels=nb)%Loaded DLM:JPEG.IDL p r in t,r e s u lt1IDL p r in t,dim ensions,nb227 1493IDL;用读取IDL dataO l=read_iniage(jp g file)IDL h elp,d a taO lDAT AOI BYTE=A rray 3rIDL;利用read_jpeg函数读取IDL read _ jp eg,jp g f ile,xmgIDL h elp,imgIMG _ B

24、YTE=A rray 3,IDL;显示图像IDL window,0,x siz e =227,y siz e=IDL tv,img,tru e =1IDLEg0 2：2)BMP格式：(不采用任何压缩)，数据存储量一般很大读：re a d_im g e(f n ame)or re a d _bmp(fn a m e,/r g b)写:write_ bm p,fname,d a t a,/r gb3)TIFF 格式：T IF F 可以存储多波段图像，还可以包含投影信息，如 lands a t-7 就是直接用.t iff存储的。E gOl:IDL t i f f =E:E N V IID L e a

25、 3 t.tifIDL r e s u lt=query image(t i f f z dim ensions=diinensions)IDL p rin ts dim ensions4135 4135IDL d a ta =r e a d _ t i f f(t i f f)IDL h elp r d a taDATA BYTE=A rray(3f 4135,4135IDL;设置关键字su b _ recc读取指点区域自像IDL sub_data=r e a d _ t i f f(t i f fr su b _ re c t=99r 49,1000r1000)IDL h e lp,sub

26、_dataSUB_DATA BYTE=A rra y 3,1000,1000IDL w indow,0rx siz e =1000,y siz e =1000IDL t vr sub d a taz tru e =15、科学数据格式读写*I D L 中支持的科学数据格式涉及CDF、HDF、H D F5、HDF-EOS、NetCDF 等等。1)HDF4 文献：目前国内外的多种卫星传感器都是将H D F 作为标准数据格式，涉 及 EOS/M ODIS、EOS/OML HJ-01/H SL F Y-3/MERS I EG01:f;读取E：EN V I I D L下的M O DIS数据的经纬度以及其1

27、K M的反射率数据集p r o t e s t _hdff n am e=d i a 1 o g _ p i c k f ile (t i t le=选择数据文献：,ge t _p at h=c ur _dir)c d,c u r _d i rh d f_i d=h d f_ sd _ s t a r t (fn a m e)；打开hd f 文献，返回一个文献 i dhd f _ sd _ f i 1 e in f o,hd f _id,s d _ n ums,at t r ib ut e;用于获取 H D F 文献的数据集数目和属性数目p r in t,s d_n um s,at t r i

28、 b u t e;读取经度1 at _ i n de x =h d f_sd_n a me to in d ex(h dfi d,L at it ud e 1)；根据数据集名称来获取数据集的索引号la t_ id=h d f_ s d _ s e l e c t (h d f_ i d,1 ajin dex);运用索弓 I 号选择数据集，返回一个数据集的idhdf_ s d_get i n f o ,1 a t _id,n ame=n a m ez un i t =u n it：杳询已打开的数据集的基本信息p r in t,n ame,*”,u n i th d f _ s d _ g et

29、dat a,lat _id,lat;读取数据,参数1 at用于返回读取的数据h df_s d_en dac c es s ,1 a t _ i d;关闭已经打开的数据集he 1 p,lat；读取纬度lo n _in dex=hdf_s d_n am e t o in dex (hdf_id,L o n git u de*)1 o n _i d=h d f_ s d_ s e 1 ec t (hdf _id,1 o n _ in de x )h d f_s d _ g e t i n f o z Io n _ id,n ame=n a me,u n it=un i tp r in t ,n a

30、m e,，u n ith df_ s d_ get d at a,1 o n i d,Io nhdf_ s d_en d ac c es s ,1 o n _ idhel p ,lo n;读取IK M反射率数据r ef_ in dex=hdf_ s d_n amet o in d e x (hdf_idz EV _1 K M _Ref S B)r ef_id=h df_ s d_ s e lec t (h df_id,r ef_in de x )hdf_ s d_g e t i n f o,r e f _ i d,n am e=n am e,u n it =u n itp r in t,n a

31、 me,u n i thd f_s d_g e t da 七a,r e f id,r ef_valu eh d f_s d_e n da c c e s s,r ef_ i dh elp,r e f _ v al u eh d f _ s d_ en d,hdf_ i d;关闭打开的文献en d运营效果%Coniled module:TEST_HDF.31 58Latitude degreesLAT FLOAT=Array271r 406Longitude degreesLON FLOAT=Array271F 406EV_lKM_RefSB noneREF_VALUE UINT=Array13

32、54r 2030,15IDL2)*.h5、*.hdf5、*.he5 格式的读取:f;下面介绍一个读取环境卫星火.h 5数据格式的例子p ro r ead_h d f5f n a m e=d i al o g_p ic k file(t it le=，选择数据文献：，g e t _p a t h=c ur _dir)cd,c ur _d i r；打开一个图形用户界面查看h5文献r es ult=h5 _ b r ows e r(f n ame)；h 5 _ m e s s age=h5 _p a r s e(f n a m e);查询文献的基本信息，返回一个结构体变量；p r i n t z h

33、 5 _ m e s s a g eh 5 _ i d=h5f_open(fn a m e)；打开一个h5文献，返回一个文献的idsd_id=h5d_op e n(h 5 _id,1 ImageD a t a/B a n d D a t a，)；打开相应的数据集，数据集名称可带途径B a n d a ta=h5d_r e ad(s d_id);读取数据;关闭文献h5d_c 1 os e z s d _idh5f_clos e,h5 _ i dh e lp,Ban da t at a转成B工P格式存储B a n d a t a_BI P =t r an s p o s e(Ban d a t a

34、,2,0,1 )h e l p.Ban d a t a_BIPw rite _ tiff,E：EN V I lDL H j-l.t i f ,Ban d a t a_B I Pen d2)netC D F文献读取：net C DF(network C o mmo n Data F o rm,网络通用数据格式)，常用于气象科学数据存储。;下面介绍一个读取E CMW F数据的例子，E C M W F为*.n c格式数据pro t e st_n cfnam e=dial o g _ p ic k fil e(t i t 1 e=选择数据文献：，get _ p at h=cur _di r )cd,c

35、 u r _ dirn c _id=n c d f_ope n(fn ame);打开 n et CDF文献n c _mes s age=n c d f _ in q u ir e (n c _ i d)；对打开的文献进行查询，返回I文献基本信息h el p,n c _ mes s age;变量查询;后一个参数为变量索引号或者名称；这里我读取的是它的第四个变量1 0 var _mes s a g e=n cdf_ v a rin q(n c _idz 4)h el p z v a r _me s s agev ar _ i d=ncdf_ v a r id (n c _ id,vlO 1;根据变

36、量名称获取相应变量的索5号p r in t z v a r id;ncdf_varge t获取数据n cdf_ v arget,nc_id,v a r _ i dz v a lu e;文献关闭ncdf_close,n c _ idh elp,valueend4)Grid格式读取：Part two：图形绘制(一)p lo t过程绘制1、p 1 ot 过程：1 )p lot,x d at a,y data,/nodata2)oplot,x d a ta,ydata;用于在现有窗口中添加新曲线2、线形符号设立：1 )关键字p s y m修改数据点的符号，直方图线形p s ym=10,psym=-a

37、表达除显示各数据点的符号外，还将各数据点连接起来。2)line style关键字设立曲线线形，0 为实线、2 为虚线、3 为点划线。3 坐标轴设立：x yt i t le设立标题;x y c h ars i ze设立字体大小；xyr a nge用于设立坐标轴范围，注意加上x y s tyle=l来强行设定坐标轴范围xytick s 和 xymino r 分别设立主刻度和最小刻度间隔数目xyt i cknam e 设立刻度的名称4、颜色设立：color：设 立 图 像 颜 色；一 般 用 十 六 进 制 来 表 达 颜色；TFFFFFxliififerOOOOO 0 1b a c kground

38、：设立背景颜色5、添加标注：x you t s 过程用于在图像窗口中添加标注信息。示例：xy o u t s,0.1 8,0.80,RMSE=L c o lor=TF FFFFx 1 ,charsi z e=1.2,/normal6、图形保存为文献：一般方法为运用TVRD函数拷贝图形窗口的内容，在写入图像文献。img=t v r d(x,y,c hannel=value,tr u e=1 2 3,/orde r)（二）plotg（通数绘制1 pl o t 函数:curv e=p 1 o t（x,y,/buf f er,/curren t,dimensions=width,height,ma r

39、 gin=num,t i t 1 e=*n a m e=/o verp 1 ot,wi n dow_ titl e=，,/n o data）o *not i ce:/b u ffer指将图形保存在缓存中g o /c urre n t 设立在当前窗口中绘制图形gg d i me n sio n s 设立窗口的大小。t itle设立图形标题,n a me设立图像对象的名称*函 数 p l o t 的返回结果为一个图像对象，具有多种对象和方法2、符号、线形设立：s ymb o 15 lines t yle t h i ck用于改变曲线线宽3、坐标轴设立:见plo t 过程4、添加标注：。函 数 te

40、x t 用于在图形窗口中添加标注信息,结果返回一个图形对象 Eg：lab e l=t e xt(x,y,标注内容，t a r g et=c u r v e,color=,r edf ont_ size=l.2)5、添加图例：“函 数 lege n d 用于在图形窗口中添加标注信息，该函数只合用于图形对象(三)散点图Eg 0 1:读取文本文献中的数据并运用p 1 o t 过程绘制相应的散点图;运用过程Plot绘制散点图pro plot_s c atter *读取数据*fn a m e=di a lo g _p i c kf i le(ti t le=选择数据：get_path=work_dir)

41、cd,w ork_di rnl=f il e _ lines(f n a m e)pr i nt,nltem p _var=data=f 1 tarr(2,nl-1)op e nr,1 u n,f n a m e /g et _ l u nre s u l t =f st a t (1 u n)he 1 p,res u 1 treadf,1 u n z t em p _ va rr e adf,l u nz d a t afree_lu n,1 u npr i n tz t ern p _va rp ri n tz d a t a；*绘制散点图*x=da t a 0z*；观测数据y=d a t

42、 a 1,*;估算数据p 1 o tz x,y r p sym=2,x r a n g e=5,1 5 ,yra n g e=5,1 5 ,$x t i t l e=1 O b s e rv e d va 1 u e y t i t l e=1 E st i m a t e d v a l u e ,$b a c k g ro u n d=*F F F F F F *x lz c o l o r=o o o o o o xl,$x s t yl e=1 ,y s t yl e=l,$c h a rsi ze=l.3,/n o da t aop 1 ot,X,V，p sym=5,c ol o r

43、=,0 0 0 0 0 0 1 xl;绘制基准线x2=5,1 5&y 2=5,1 5 oplot,x 2fy 2 ,c o l o r=*F F 0 0 0 0 x 1；*计算M A E与R M S E并在散点图中添加标注*M A E=m e a n (a b s(x-y)；M AE平均绝对没菜R MS E =s q r t (me a n(x y)A2);均方根误差MA E _ 1 a b el=,MA E=,4-s t rin g (MA E,f o r m a t=f(f 5 .2 )R MS E _ l a b el=*R M S E=+s t r ing(RMSE,f o rm a

44、t=,(f 5.2)1);添加标注x y o u ts,6,13.0 z MA E _ l a b el,c o 1 o r=O O O O O O l,c h a rsi z e=1 .2,/da t ax yo u ts,6,1 2.3,R MS E _ l a b e 1 ,c o 1 o r=1 0 0 0000 1 x lf c h a rsi ze=l.2,/da t a;保存散点图为图像文献，格式为.p n gim g=tvrd(t ru e=l);法一：o _ f n=d i a l o g _ p ic k fi 1 e(t i t l e=，文献保存为：1);键入的文献名要

45、包含后缀名write_pngz o _ f nz i m g;法二;wr i t e_ p n g,1 E ：E N V I sc a t t er.p n g 1 r i m ge ndEg02:读取文本文献中的数据并运用P lot函数绘制相应的散点图；p l o t函数绘制散点图p r o p lo t_ sc a tte r0 2;读取数据f i 1 e=d ia l o g_p i c k f ile (t i t l e=选择数据：，g e t _p a t h=wo r k _ di r)c d,wo rk di rv a r=nl=fi 1 e_ 1 in e s(f i l e

46、)da t a=f l t arr(2,n l-1)ope n rz l u n,f i l e,/g e t _ l u nreadf,l u n ,v a rre a d f z l u n,da t afree_lun,l u np r int,d a t a;绘图o b=d a t a 0,*E s=da t a 1 ,*g r a p h i c 0 l=p l o t (o b,E s r x r a n g=7,1 5 ,yra n g=7,1 5 ,xm i n o r=5z ym i n o r=5,$；先数据t i t l e=sc a t t e r 1,xt i t l

47、 e=1 O b s erved V a 1 u e*,y t i t l e=E st i m a t e d V a i u e*,$；再坐标轴sym b o l=,X *,s y m _ s i ze=l .2,l i n es t yl e=1 n o n e 1 r c o 1 o r=*b l a ck ,$;然后线形m a rg i n=0.1)X2 =7,1 5&y2=7,1 5 g ra p h i c 0 2=p l o t(x2,y 2 r 1 i n e s t y l e=l,/o v e rp 1 o t);计算平均绝对误差和均方根误差并添加标注M A E=mea

48、n(abs(E s-o b )；RMSE=sq r t (me a n(Es-o b)人2)MA E _ 1 a b el=*MA E=*+st r i n g(MA Ez f o rm a t =*(f 5 .2)RM SE_lab e 1 =,R MS E =,+s tring(R M SE,f o rm a t=,(f 5.2)t0 1=te x t(0.15z 0.8,M A E _l a b el,t a rg et=g ra p h i c 0 1,f o n t s iz e=1 2)；f o n t s_ s i ze用于设立标注的字体大小t 0 1 =t ex t(0.1 5

49、,0.7 5,R MS E _ 1 a b el,t a rg et=g ra p h i c 0 1,f o n t _ s ize=12);保存文献o_fn=dialog_pic k f i l e (title i图形保存为:*)g r a p h i c 0 1.s a v ez o _f ng r a p h i c 0 2.s a v e,o _fne n d(四)柱状图、条形图EgOl:*运用bar_plot()函数绘制柱状图(不建议用直接图形法来绘制)*运用erro r _ p lo t添加误差线;b a rp l o t功能函数绘制柱状图p ro t es t _ h i s

50、 t o g ram02;#读取文本文献数据#f i l e=d i a l o g _ p ic k ile (t i t l e=，选择文本文献：,g e t _p a t h=wo rk _ di r)cd,wo r k _ di rn l =f i l e 1 i n e s(f i 1 e)v ar=s t r a r r (1,n 1 )op e n r,lu nz f ile,/g e t_ lu nread f，lun,varfre e _ lu nz 1 u np r i n t,var;拆分ar e a=stra r r(n l)n um b e r=i n ta rr(n