GIS数据处理-IDL区块读取写入TIFF文件

在GIS的数据预处理中经常用要读取大数据文件，通常情况下我会选择用IDL做结果试验，成功以后再用C++编码。

使用IDL读取TIFF文件，会遇到文件过大的情况，这时候如果全部读取进来很浪费存储空间和处理效率，通常情况可以选择区块读取TIFF文件。

下面是示例代码

;Name: IDL 读取 tiff 文件
;Time: 2014年3月1日 21:13:38
;Author: 尹全超
PRO readtif

tifarr = READ_TIFF('d:\aa.tif',SUB_RECT=[0,0,512,512])

img = image(tifarr,title="aa.tif")

WRITE_TIFF,'d:\bb.tif',tifarr

END

这里主要通过READ_TIFF的SUB_RECT设置读取区块大小，四个参数分别代表起始坐标和宽度、高度。

如果图像文件是多波段的，返回的就是一个3维的数组；如果是单波段的，返回就是一个二维的图像数组。

后面WRITE_TIFF是写入到一幅新的TIFF文件里面，但是通过本文的方法写入的TIFF文件无法用windows图片浏览器打开，采用photoshop可以打开。