01 实验数据

诸多.float后缀的影像文件(但以ENVI默认格式存储)

02 实验思路

迭代循环所有影像文件所在的文件夹， 获取每一个float后缀的影像文件，并对每一个影像文件进行循环，获取循环文件的每一个波段影像的中值，最后将其输出为csv文件（每一列表示一个影像文件各个波段的中值数）

03 部分函数说明

Result = SIZE( Expression [, / L64 ] [, / DIMENSIONS | , / FILE_LUN | , / FILE_OFFSET | , / N_DIMENSIONS | , / N_ELEMENTS | , / SNAME , | , / STRUCTURE | , / TNAME | , / TYPE ] )

解释：expresion表示传入的数组，size函数依次返回该数组的维度、列数、行数、数组元素的数据类型(返回类型代码)、数组元素的总个数；（以一维数组形式返回）

（除了L64之外，上述所有关键字参数用于确定SIZE函数的返回值，并且它们是互斥的——最多只能指定上述关键字之一）

部分参数说明：

DIMENSIONS

将此关键字设置为返回Expression的维度。如果Expression是标量，则结果是包含0的标量。对于数组、列表和哈希表，结果是包含数组或列表维度的数组。如果Expression中的元素数量需要，结果为32位整数，如果需要，为64位整数。设置L64以强制在所有情况下返回64位整数。如果Expression未定义，则IDL报告八个维度。（例如是3列4行的数组，那么返回3 4）

L64

默认情况下，SIZE的结果在可能的情况下为32位整数，如果Expression中的元素数量需要，则为64位整数。设置L64以强制在所有情况下返回64位整数。除了影响默认结果外，L64还影响DIMENSIONS、FILE_OFFSET、N_ELEMENTS和STRUCTURE关键字的输出。

N_DIMENSIONS

将此关键字设置为返回Expression中的维数（如果它是数组或列表）。如果Expression是标量或未定义的，则返回0。

N_ELEMENTS

将此关键字设置为返回Expression中的数据元素数量。设置此关键字相当于使用N_ELEMENTS函数。如果Expression中的元素数量需要，则结果将为32位整数，如果需要，则为64位整数。设置L64以强制在所有情况下返回64位整数。如果Expression未定义，则返回0。

其中，数据类型代码的指代如下：

Result = MEDIAN(Array [, Width] [, /DOUBLE] [, DIMENSION=value] [, /EVEN])

其中，参数的含义为：

- Array：要计算中位数的数组。

- Width：指定中位数计算的窗口大小。默认情况下，中位数计算的窗口大小为1，即只考虑每个元素本身的值。如果指定了Width参数，则中位数计算的窗口大小为2*Width+1。例如，如果Width=2，则计算的窗口大小为5，包括当前元素及其相邻的4个元素。如果Width=0，则计算的窗口大小为1，即只考虑当前元素本身的值。

- /DOUBLE：指定计算中间结果时使用双精度浮点数。默认情况下，MEDIAN函数使用单精度浮点数计算中间结果。

- DIMENSION=value：指定在哪个维度上计算中位数。默认情况下，MEDIAN函数在整个数组上计算中位数。如果指定了DIMENSION参数，则在指定的维度上计算中位数。例如，如果DIMENSION=1，则在数组的第一维(行)上计算中位数。

- /EVEN：指定在计算中位数时，如果窗口大小为偶数，则取中间两个数的平均值作为中位数。默认情况下，如果窗口大小为偶数，则取第二个中间数作为中位数。

04 代码部分

pro batch_extract_median
  ; 二次开发接口
  compile_opt IDL2
  e = envi(/headless) ; Launch the application
  ; 所有需要提取的文件的所在目录
  in_path = "C:\Users\HeZehuang\Desktop\66"
  ; 获取in_path目录下所有满足后缀为.float的文件
  files_path = file_search(in_path, "*.float", count=files_amount)
  ; 创建存储所有影像中值的csv文件
  csvfile_path = "C:\Users\HeZehuang\Desktop\median.csv"
  ; 存储所有影像中值的矩阵
  box_median = fltarr(files_amount, 128)
  ; 进入循环对每一个文件做批量处理
  for file_i = 0, files_amount - 1 do begin
    ; 获取当前循环下的文件路径
    file_path = files_path[file_i];
    file_name = FILE_BASENAME(file_path, ".float")
    ; 获取该文件数据
    dat = e.OpenRaster(file_path)
    dat_data = dat.GetData()
    ; 获取该文件的行列数以及波段数，这里需要注意，这里size函数的输出结果与前文的讲解不一致，可能
    ; 与该文件是ENVI默认格式有关，一般读取tiff等格式时就是正常的，所以大家最好打开ENVI看一下行列
    ; 以及波段数在数组的哪个索引上
    dat_size = size(dat_data)
    dat_col = dat_size[1]
    dat_row = dat_size[3]
    dat_channel = dat_size[2]
    ; 对当前循环的文件再次循环，找到每一个波段的中值
    for channel_i = 0, dat_channel - 1 do begin
      ; 获取当前波段影像的中值
      data = dat_data[channel_i, *, *]
      channel_median = median(data[where(data ne 0.0)], /EVEN)
      ; 存储当前循环文件下的第channel_i波段的中值
      box_median[file_i, channel_i] = channel_median
    endfor
    print, file_name, "提取中值完毕"
  endfor
  ; 将保存在矩阵中的中值全部存储为csv格式
  write_csv, csvfile_path, box_median
end