IDL程序设计技术方案笔记.docxVIP

目录Part one：文件读写操作11、格式化输入与输出（read\print）12、ASCII文件读写23、二进制文件读写54、图像格式文件读写75、科学数据格式读写***8Part two：图形绘制11（一）plot过程绘制11（二）plotg()函数绘制12（三）散点图12（四）柱状图、条形图15（五）绘制颜色条colorbar18Part three: ENVI二次开发201、自定义波段运算函数202、IDL中常见波段运算213、ENVI classic 创建自定义菜单224、ENVI classic+IDL 二次开发23Partone：文件读写操作1、格式化输入与输出（read\print）Read:该函数用于从键盘键入数据，默认数据类型为浮点型，若要输入其他类型数据，需要先定义数据类型。Eg:Print：该函数用于将数据输出打印到控制台，这里主要讲格式化控制输出format关键字。表01 常用格式化控制符格式代码输出效果aN字符或者字符串按照N个字符宽度输出fn.m浮点型数组按照N个字符宽度输出，小数点保留M位dn.mdouble型数组按照N个字符宽度输出，小数点保留M位en.m按照N个字符宽度的科学计数法输出，小数点后精确到M位Nx输出N个空格字符串/H直接引用字符串输出或者用H直接输出字符串c()用于输出日期数据Eg：2、ASCII文件读写常见ASCII文件：txt、matlab程序文件(*.m)、c程序文件（*.c）、遥感影像头文件（如ENVI格式的头文件*.hdr）选择文件：fn = dialog_pickfile(title=选择ASCII文件：,get_path = work_dir)cd,work_dir 。用于跳转到当前工作路径查询文件：file_search()函数Eg:查找当前工作路径下所有的txt文件，并返回文件数目*: File_lines(fname)函数用于查询文本文件的行数。打开文件：在IDL中读写ASCII码或者二进制文件，首先需要将文件与一个逻辑设备号关联起来。语法：openr/openw/openu, lun, fname, /get_lun, width=10,/append读写操作：IDL中利用readf函数读取文本文件、利用printf函数写入ASCII文件①Readf, lun, var01, var02….②printf, lun, var01, var02….文件关闭：在对文件操作完成之后，需要关闭文件的逻辑设备号，eg: freee_lun，lunEG:*如何获取某个文件你的列数:先读取一行数据，然后利用空格进行字符串拆分；最后利用n_elements()计算元素个数即为ns对于某些遥感数据的头文件前面几行是说明文件，后一部分才是数据的，读取有两种方式：定义一个临时变量temp，将说明文字用该变量存储直接跳行：skip_lun, lun, 3, /lines 。跳过文件的前三行3、二进制文件读写多数遥感数据的数据文件都是二进制文件。函数readu, lun, var01, var02….用于读取二进制文件*：对于某些数据文件开头是说明信息，读取时可以用point_lun, lun, position直接按字节跳过。Eg:读取IDL文件下的envi.img数据函数writeu, lun, var01, var02…用于写入二进制文件。Eg01:将上面读取的envi.img计算其NDVI并保存为二进制文件Eg02:Eg03:。读取cha06下的风云02卫星数据proread_AWX file = E:\zengsk\IDL\Cha06\data\FY2C_TBB_IR1_OTGAOAD.AWXopenr, lun, file, /get_lun。跳过前20个字节point_lun,lun,20 Headline = indgen(3)。定义三个整型数据，分别读取记录长度、头文件记录数、数据记录数readu, lun, Headlineprint,headline[0], headline[1],headline[2]。定位到第58个字节，读取数据日期point_lun, lun, 58Begindate = indgen(5)。记录开始日期readu, lun, BegindateEnddate = indgen(5)。记录结束日期readu, lun, Enddatelat_long = indgen(4)。记录网格左上和右上角的经纬度readu, lun, lat_longprint,Begindateprint,Enddateprint,lat_long。基于Headline信息数组，定义字节数组，其行列数分别为记录长度和数据记录数 data