本文介绍基于Python语言,对大量不同的Excel文件加以跨文件、逐单元格平均值计算的方法。
首先,我们来明确一下本文的具体需求。现有一个文件夹,其中有如下所示的大量Excel文件,我们这里就以.csv
文件为例来介绍。其中,每一个.csv
文件的名称都是如下图所示的Ref_XXX_Y.csv
格式的,其中XXX
表示三个字母,后面的Y
则表示若干位数字。
对于其中的每一个.csv
文件,都有着如下图所示的数据格式。
我们现在的需求是,希望对于每一个名称为Ref_GRA_Y.csv
格式的.csv
文件,求取其中每一个单元格在所有文件中数据的平均值。例如,对于上图中DOY
为1
的blue
这个单元格,那么求出来的平均值就是在全部名称为Ref_GRA_Y.csv
格式的.csv
文件之中,DOY
为1
且列名为blue
的单元格的平均值。此外,如果像上图一样,出现了部分单元格数值为0
的情况,表明在当前文件夹下,这个单元格是没有数据的,因此需要在计算的时候舍去(并且取平均值时候的分母也要减小1
)。
知道了需求,我们就可以开始代码的书写。其中,本文用到的具体代码如下所示。此外,本文实现的需求也和我们之前的文章基于Python读取多个Excel文件数据并跨越不同xlsx表格文件计算平均值(https://blog.csdn.net/zhebushibiaoshifu/article/details/115533619)有些类似,大家如果有需要,也可以参考之前的这一篇文章。
# -*- coding: utf-8 -*-
"""
Created on Fri Oct 6 13:07:48 2023@author: fkxxgis
"""import os
import glob
import pandas as pdfolder_path = "E:/04_Reconstruction/02_Data/01_RGBNINDVI_History"
output_path = "E:/04_Reconstruction/02_Data"
file_pattern = "Ref_GRA_*.csv"file_paths = glob.glob(os.path.join(folder_path, file_pattern))combined_data = pd.DataFrame()for file_path in file_paths:df = pd.read_csv(file_path)df_filtered = df[df != 0]combined_data = pd.concat([combined_data, df_filtered])average_values = combined_data.groupby('DOY').mean()output_file = "04_Data_YearAverage.csv"
average_values.to_csv(os.path.join(output_path,output_file), index=True)
其中,上述代码的具体介绍如下。
首先,我们导入必要的库——os
库用于文件路径操作,glob
库用于文件匹配,pandas
库用于数据处理和分析。同时,我们定义文件夹路径folder_path
,代表存储.csv
文件的文件夹路径;定义输出路径output_path
,代表保存结果文件的路径;定义文件匹配模式file_pattern
,用于匹配需要处理的.csv
文件的文件名模式。
随后,我们使用glob.glob()
函数结合文件夹路径和文件匹配模式,获取满足条件的.csv
文件的路径列表,存储在file_paths
变量中。创建一个空的数据框combined_data
,用于存储所有文件的数据。
接下来,我们使用一个循环,遍历file_paths
列表中的每个文件路径。对于每个文件路径,使用pd.read_csv()
函数加载.csv
文件,并将其存储在名为df
的数据框中。其次,使用条件筛选语句df[df != 0]
排除值为0
的数据,并将结果存储在名为df_filtered
的数据框中。紧接着,将当前文件的数据框df_filtered
合并到总数据框combined_data
中,这一步骤使用pd.concat()
函数实现。
完成所有文件的处理后,使用combined_data.groupby('DOY').mean()
计算所有文件的平均值,按照DOY
列进行分组并求平均值。随后,定义输出文件名output_file
,代表保存平均值结果的文件名。
最后,使用os.path.join()
函数结合输出路径和输出文件名,生成保存路径,并使用average_values.to_csv()
函数将平均值数据框average_values
保存为一个新的.csv
文件,指定index=True
以包含索引列。
运行上述代码,我们即可得到结果文件。如下图所示,可以看到结果文件中,已经是计算之后的平均值结果了。
至此,大功告成。
欢迎关注:疯狂学习GIS