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如果有问题，欢迎一起探讨.
<p>炒茄子</p>
