IDL基础【影像读取-转存数组-输出】

IDL是一门非常适合高分辨率和大量图像数据处理的语言，在envi的安装过程中自带，通过学习IDL，可以实现十分简单的像元运算功能，方便快捷易上手。

学习所有计算语言，入门分为三个内容：数据输入 —— 数据处理 —— 数据输出

IDL也不例外，本文将采用IDL8.5版本，对采用IDL进行图像读取、转换为数组、图像输出三个部分的内容进行讲解

讲解图像读取之前，首先需要了解一下，任何一个IDL脚本都是从pro开始，end结束，所以这个开头和结束不能省略，并且，pro后的名称代表该脚本/函数名，所以文件也需要存储为对应的名字。

开始输入前，首先新建IDL脚本，构建程序框架并保存文件名为函数名：

pro my_program

end

一、IDL图像读取

compile_opt idl2 ;更改IDL的编译配置，推荐值为IDL2，参考：http://www.harrisgeospatial.com/docs/compile_opt.html

;如果IDL程序中调用了ENVI函数，那么在发布的时候需要加上以下语句：
ENVI,/restore_base_save_files
ENVI_BATCH_INIT

;选择图像文件
image_file=DIALOG_PICKFILE(PATH=’examples\data’, $
TITLE=’Select DICOM Patient File’, FILTER=’*.tif’) ;TITLE为数据选取窗口的标题，PATH为打开文件的默认目录，TILTER为过滤出的数据格式

;读取影像文件

e=envi()
raster=e.OpenRaster(image_file)

二、图像转换到数组

;获取文件ID

fid = ENVIRasterToFID(raster)
;获取对应ID的文件基本信息

MAP_INFO = ENVI_GET_MAP_INFO(fid=fid);获取投影信息
ENVI_FILE_QUERY, fid, ns=ns, nl=nl, nb=nb, dims=dim, data_type=data_type ;获取图像行列数、维度信息和数据类型信息

;获取数据

data = raster.getdata();这里data就是得到的图像数组

三、输出数据

ENVI_WRITE_ENVI_FILE, data, out_name=’c:\desktop’, /nocopy, $
ns=ns,nl=nl,nb=nb, offset=offset, bnames=bnames, map_info=map_info

以上三个部分内容，主要完成了遥感影像数据的读入并获取相关数据的信息、转换成可计算的数组以及进行带投影的输出，其中部分功能为ENVI5.0版本以后新更新的用法，而输出功能则是传统用法。

在学习过程中，个人觉得比较难把握的就是如何带上投影信息进行数据转存，目前这种输出方法只能输出为标准ENVI STANDARD格式，如何带投影输出为TIF等格式还需要进一步学习。