这里将对ENVI-IDL二次开发程序的一个通用流程做一个总结。
1.首先是文件打开和数据读取:
文件打开work_dir = dialog_pickfile(title='选择路径',/directory) ;手动选择文件路径 fns = file_search('*.dat', count=fnums) ;搜寻路径中后缀为...的文件集合,接着可以用fn[i]进行遍历操作 result = envi_pickfile(/multiple_file,title='select images') ;打开图形对话框选择文件名 filename = '......' ;直接读取路径操作 envi_open_data_file,fname,r_fid=fid,/png ;在ENVI中打开ENVI支持的外部格式文件,并返回fid号filename = filepath('a.xml', root_dir='D:\', subdirectory=['temp']) ;root_dir文件根路径,subdirectory子路径数据读取result = envi_get_data(fid=file id, dims=array,pos=long integer) ;读取一个波段数据result = envi_get_slice(fid=file id, line=integer, pos=array,xs=value,xe=value,[/bil][,bip]) ;从打开的ENVI文件中读取一行的数img = read_image(fn) ;读取图像,结果为图像数组data = e.OpenRaster(fn) ;要打开envi使用,来读取数据
2.其次是数据查询:
envi_file_query, fid, ns=ns, nl=nl, nb=nb, dims=dims ;对于打开的ENVI文件进行查询,获取行列数、波段数、头文件偏移、数据类型、数据存放顺序、定标系数等信息 envi_get_map_info(fid=file id) ;返回投影坐标信息
3.接着是数据读取写入(二进制文件/ASCII码...):
openr, lun, fname, /get_lun ; 打开文件读操作 openw, lun, fname, /get_lun ; 写操作 openu, lun ,fname, /get_lun ; 读写操作 readf, lun, var ;读取ASCII码 readu, lun, var ;读取二进制文件 printf, lun, var ;写入ASCII码 writeu, lun, var ;写入二进制文件 result = read_csv(fname,count=n1,header=header) ;读取csv文件
4.文件保存:
envi_write_envi_file, data, out_name=variable, ns=ns, nl=nl, nb=nb ;文件保存 data.export, 'c:/1.tif' ;文件保存
5.异常处理:
Catch, errorStatus IF (errorStatus NE 0) THEN BEGIN Catch, /CANCEL err_msgs = [err_msgs, file+'--'+!ERROR_STATE.MSG] MESSAGE, /reset CONTINUE ENDIF
6.文件路径、名称更改:
FILE_BASENAME(file) ;获取文件名 FILE_DIRNAME(file) ;获取文件夹路径 file.replace('.tif','_ndvi.tif') ;更换文件名后缀,进行更改
7. 字符串操作方法:
除了上述基本操作以外,还有投影信息的变换、矢量数据读取等方法,下次进行介绍。
来源:https://www.cnblogs.com/ljwgis/p/12250039.html