ENVI批量输出ASCII文件

Question

ENVI将文件输出为其他格式的时候，不能批处理吗

一百多个遥感图像，想输出为ASCII格式，只能一个文件一个文件转吗，好麻烦

找了一下IDL调用ENVI的文章，没啥收获捏

求指导

Answer 1

是的，可以通过IDL脚本调用ENVI的批处理功能实现遥感图像的批量ASCII格式转换。以下是具体实现方法与步骤：

---

1. 核心实现原理
   ENVI的IDL二次开发接口支持通过无界面（Headless）模式启动批处理流程，结合ENVI_OPEN_RASTER读取栅格数据，再通过WRITE_ASCII或自定义数据导出逻辑将像素值写入文本文件。关键点包括：

• 无界面模式启动：使用ENVI(/HEADLESS)避免图形界面加载，提升处理效率。
• 数据分块读取：针对大文件，通过分块读取内存数据（如ENVI_GET_SLICE）避免内存溢出。
• 格式转换逻辑：将栅格数据矩阵转为ASCII文本格式，支持自定义分隔符（如逗号、空格）和精度控制。

---

1. 具体实现步骤

2. 1 初始化批处理环境
   idl
   ; 启动ENVI无界面模式
   e = ENVI(/HEADLESS)
   ; 设置输入输出路径
   inputDir = 'C:\遥感数据\输入'
   outputDir = 'C:\遥感数据\ASCII输出'
   files = FILE_SEARCH(inputDir + '/*.dat') ; 支持.dat、.img等格式

3. 2 遍历文件并处理
   idl
   FOREACH file, files DO BEGIN
   ; 打开栅格文件
   raster = e.OpenRaster(file)
   ; 获取数据信息
   ns = raster.NCOLUMNS
   nl = raster.NROWS
   nb = raster.NBANDS
   dataType = raster.DATA_TYPE
   ; 读取全幅数据（小文件适用）
   data = raster.GetData()
   ; 分块读取（大文件优化）
   ; slice = raster.GetSlice(BANDS=0, SUB_RECT=0,0,ns-1,nl-1)

   ; 构建输出文件名
   baseName = FILE_BASENAME(file, '.dat')
   outFile = outputDir + '' + baseName + '.asc'

   ; 写入ASCII文件头（可选）
   OPENW, lun, outFile, /GET
   PRINTF, lun, 'ncols ' + STRTRIM(ns,2)
   PRINTF, lun, 'nrows ' + STRTRIM(nl,2)
   PRINTF, lun, 'xllcorner 0.0' ; 根据实际坐标调整
   PRINTF, lun, 'yllcorner 0.0'
   PRINTF, lun, 'cellsize 30.0' ; 设置像元大小
   PRINTF, lun, 'NODATA_value -9999'

   ; 逐行写入数据
   FOR i=0, nl-1 DO BEGIN

    line = data*,i
    PRINTF, lun, STRJOIN(STRING(line, FORMAT='(F8.2)'), ' ')
   

   ENDFOR

   FREE_LUN, lun
   ; 释放内存
   raster.Close
   ENDFOREACH

---

1. 功能扩展与优化

2. 1 多波段处理
   若需将多波段数据合并为单个ASCII文件，可调整循环逻辑：
   idl
   FOR band=0, nb-1 DO BEGIN
   data = raster.GetData(BANDS=band)
   ; 按波段追加到文件或生成独立文件
   ENDFOR

3. 2 分块处理与内存优化
   对超大文件，使用分块读取（Tile-based Processing）：
   idl
   tileSize = 1000 ; 分块行数
   FOR j=0, nl-1, tileSize DO BEGIN
   subRect = 0, j, ns-1, MIN(j+tileSize-1, nl-1)
   slice = raster.GetSlice(BANDS=0, SUB_RECT=subRect)
   ; 逐块处理并写入
   ENDFOR

4. 3 格式精度控制
   通过FORMAT关键字指定数值精度（如保留两位小数）：
   idl
   PRINTF, lun, STRJOIN(STRING(line, FORMAT='(F8.2)'), ',')

---

1. 注意事项

• 数据类型兼容性：确保输入数据为浮点或整型，分类数据需特殊处理。
• 元数据保留：若需保留地理坐标信息，可通过ENVI_GET_MAP_INFO获取投影参数并写入头文件。
• 性能平衡：ASCII文件体积较大，建议压缩（如.zip）或优先考虑二进制格式（如CSV）。

---

1. 替代方案
   若需更高效的批量处理，可直接调用ENVI内置任务ENVI_SAVE_AS_ASCII_DOIT（需确认版本兼容性）：
   idl
   ENVI_DOIT, 'ENVI_SAVE_AS_ASCII_DOIT', FID=fid, DIMS=dims, POS=pos, OUT_NAME=outFile

---

通过上述方法，可灵活实现遥感数据的批量ASCII转换，满足自动化处理需求。实际应用中需根据数据规模和格式调整代码逻辑。