IDL提取多个.float高光谱数据平均反射率并存为CSV

1. 首先，需要将多个.float格式的高光谱数据读入IDL中。可以使用IDL的READ_FLOAT命令读取单个文件，或者使用FILES_GET命令读取多个文件（如一个文件夹下所有的.float文件），然后使用READ_FLOAT读取每个文件：

;假设需要读取的文件名存在一个字符串数组file_names中
for i=0, n_elements(file_names)-1 do begin
    file_id = openr(file_names[i], /get_lun)
    data = read_float(file_id, /endian)
    close, file_id
    ;在这里进行处理
endfor

2. 处理高光谱数据，提取每个波段的平均反射率。由于每个波段的波长可能有所差异，需要先读取波长信息并进行插值，以使所有波段的反射率计算使用相同的波长信息。假设读取到的波长数据存在于一个数组wavelength中，反射率数据存在于数组data中：

;插值波长信息
wave_ref = findgen(751) * 0.1 + 350.0    ;设定参考波长范围，最小波长350nm，间隔0.1nm
interpolate = interpolate(wave_ref, wavelength, data, /edge, /nan, /linear)
;计算平均反射率
mean_ref = total(interpolate, 1) / n_elements(data[*, 0])

3. 将计算得到的平均反射率写入CSV文件。可以使用IDL的WRITE_CSV命令将数组写入CSV文件。假设需要将平均反射率存在一个名为output.csv的文件中：

write_csv, mean_ref, output.csv, /double, header=['wavelength', 'mean_reflectance'], format='(%f, %.6f)'

以上代码片段可以整合到一个程序中，如下：

;假设需要读取的文件名存在一个字符串数组file_names中
;插值波长信息
wave_ref = findgen(751) * 0.1 + 350.0    ;设定参考波长范围，最小波长350nm，间隔0.1nm
;循环读入和处理每个文件
result = fltarr(751, n_elements(file_names))
for i=0, n_elements(file_names)-1 do begin
    file_id = openr(file_names[i], /get_lun)
    data = read_float(file_id, /endian)
    close, file_id
    ;插值到参考波长下
    interpolate = interpolate(wave_ref, wavelength, data, /edge, /nan, /linear)
    ;计算平均反射率
    mean_ref = total(interpolate, 1) / n_elements(data[*, 0])
    ;存储结果
    result[*,i] = mean_ref
endfor
;写入CSV文件
header_string = string(['wavelength', 'mean_reflectance'], /format)
format_string = '(%f, %.6f)'
write_csv, result, 'output.csv', /double, header=header_string, format=format_string

在程序中，result数组用于存储计算得到的平均反射率数据，每一列是一个文件对应的反射率数据。在循环中计算出每个文件的平均反射率后，将结果存储在result数组中。最后使用WRITE_CSV命令将result数组写入CSV文件。