文章详情

短信预约-IT技能 免费直播动态提醒

请输入下面的图形验证码

提交验证

短信预约提醒成功

数据预处理方法整理(数学建模)

2023-09-10 09:05

关注

这篇文章主要是整理了一些作者在各种建模比赛中遇到的数据预处理问题以及方法,主要针对excel或csv格式的数据,为后续进行机器学习或深度学习做前期准备

数据清洗

导入库和文件,这里使用的是绝对路径,可改为相对路径

import pandas as pdimport numpy as npimport matplotlib.pyplot as plt# 读取xls数据文件data = pd.read_csv(r'D:\1112222.csv',encoding='gbk')

传入的为csv格式的文件,如果是xlsx格式的文件,建议先使用excel另存为csv结尾格式的文件再进行操作,选择下图所示选项:

合并多个文件数据

可能同时存在多个csv文件,需要进行合并

import glob,osfilenames_in = r'D:\in'  # 输入文件的文件地址filenames_out = r'D:\inner'  # 新文件的地址path_in  = r'D:\in'file_names = os.listdir(path_in)file_paths = glob.glob(os.path.join(path_in,'*.csv'))print(file_paths) df1 = pd.DataFrame()for file in file_paths:    df2 = pd.read_csv(file,sep=',',header=None)    #df2=df2.iloc[:,2]  #只取第三列    df1 = pd.concat([df1, df2], axis=0)      #axis=0意思是纵向拼接,=1的时候是横向拼接print('dataframe的维度是:', df1.shape)#print(df1)# 输出数据到本地df1.to_csv(r'D:\inner\result.csv', index=False, sep=',')

填补缺失值

对于数字类型的数据,用均值填补空值;对于字符类型的数据,用众数填补空值

def fill_missing_values(df):    """用DataFrame中各列的均值或众数来填补空值"""    for column in df:        if df[column].dtype == np.number:  # 如果数据是数字类型            mean = df[column].mean()            df[column].fillna(mean, inplace=True)  # 用均值填补空值        else:  # 如果数据不是数字类型            mode = df[column].mode().iloc[0]  # 找到最频繁出现的项            df[column].fillna(mode, inplace=True)  # 用众数填补空值    return dfdata=fill_missing_values(data)

去除数据中的符号

对于有些数据,可能含有空格、中英文标点,需要进行去除

import stringfrom zhon.hanzi import punctuationpunctuation_string = string.punctuationfor i in punctuation_string:    data= data.replace(i, '')punctuation_str = punctuationfor i in punctuation_str:    data = data.replace(i, '')

当然,如果遇到这种数据,我会更建议使用excel自带的功能进行手动去除,比调代码更快(毕竟有时候可能调代码半天手动早就做好了)

如图,c列是一列数据前面有空格的数据,这时我们只需要在c列后插入新的一列,并在第一行中输入c列无空格的数据

输入好后,直接快捷键CTRL+E,即可自动将D列填充为去除空格后的数据(注意要确保此时其他数据没有别的空缺,否则会把别的空缺也自动填充上)

 四万个数据不到一秒就全部填充好了,速度还是比较快的。然后直接把c列删除即可。

如果是数据中包含特殊字符,想去除也可以用这个方法,只是对某几列做操作时,比调代码快。

去除冗余数据

有些数据,可能一整列都是同一个数据,没有变化,这些变量对机器学习没有帮助,但因为数据太多,不可能人工判断每列数据情况,需要代码调试进行去除

for col in data.columns:    # 如果这一列所有的值都相等    if data[col].nunique() == 1:        # 则删除这一列        data = data.drop(col, axis=1)

格式转换

将true和false类型数据转换为int型

data['11'] = data['11'].astype(int)

其他类型同理,将“int”改为想转换的类型即可

字母或字符串转换为对应数字

如果想将字母或字符串转换为数字在后续作为变量进行处理,建议直接在excel中快捷键CTRL+h进行转换,代码总是会出现各种错误。

合并某几列数据

如果给出的数据有年、月、日,需要将其合并成一列

data['timestamp'] = data['月'].astype(str) + '-' + data['日'].astype(str) + '-' + data['具体时间']data = data.drop(['月', '日', '具体时间'], axis=1)# 将时间戳列设置为索引#data.set_index('timestamp', inplace=True)data['time'] = pd.to_datetime('2023-' + data['timestamp'], format='%Y-%m-%d-%H:%M:%S')

对多列数据求平均值并合并为一列:

class_df = (data['ROLL_ATT1']+data['ROLL_ATT2'])/2data['ROLL_ATT1']=class_dfdata = data.drop(['ROLL_ATT2','MAGNETIC_HEADING'], axis=1)

数据可视化

可以将已处理好的数据进行可视化

以某个自变量为横坐标(如时间),其余为纵坐标画出折线图,画出所有变量随时间变化的折线图:

import matplotlib.pyplot as plt#画出所有变量随时间变化图像feature = data.columns[1:]for feas in feature:    plt.plot(data['time'], data[feas])    plt.xlabel('Time')    plt.ylabel(feas)    plt.show()

判断某个变量的分布情况可以绘制分布图

# 绘制分布图plt.hist(df2['train1'], bins=20)plt.xlabel('train1')plt.ylabel('Frequency')plt.title('Takeoff Weight Distribution')plt.show()

 散点图绘制

import seaborn as snssns.pairplot(data , hue ='label')plt.savefig(r"D:\pairplot001.png")

效果如下:

也可以画出相关性系数热力图

import seaborn as snssns.set(style="ticks")sns.heatmap(data.corr(), annot=True, cmap="YlGnBu");plt.savefig(r"D:\heatmap.png")

效果如下:

 

 针对机器学习及深度学习数据预处理

导入库(keras库的内容是深度学习才需要用到的,仅进行机器学习可以不导入)

import numpy as npimport pandas as pdimport kerasfrom keras.models import Sequentialfrom keras.wrappers.scikit_learn import KerasClassifierfrom keras.utils import np_utils, plot_modelimport matplotlib.pyplot as plfrom sklearn import metricsfrom sklearn.model_selection import cross_val_score, train_test_split, KFoldfrom sklearn.preprocessing import LabelEncoderfrom keras.layers import Dense, Dropout, Flatten, Conv1D, MaxPooling1Dfrom keras.models import model_from_jsonimport matplotlib.pyplot as pltfrom sklearn.metrics import confusion_matriximport seaborn as snsimport osfrom sklearn.preprocessing import StandardScalerimport itertools

首先是对数据集进行划分

data2=data.drop(['label'],axis=1)X = np.expand_dims(data2.astype(float), axis=2)Y = data['label']print(X.shape)print(Y.shape)x_train, x_test, y_train, y_test = train_test_split(X, Y, test_size=0.7, random_state=20)x_valid, x_test, y_test, y_valid=train_test_split(X, Y, test_size=0.5, random_state=20)

数据标准化

standard = StandardScaler()# 对训练集进行标准化,它会计算训练集的均值和标准差保存起来x_train = standard.fit_transform(x_train)# 使用标准化器在训练集上的均值和标准差,对测试集进行归一化x_test = standard.transform(x_test)

接下来都是深度学习需要进行的预处理操作,机器学习不用

深度学习需要对label进行onehot编码

from keras.utils import to_categoricaly_test = to_categorical(y_test)y_train = to_categorical(y_train)from keras import backend as KK.set_image_dim_ordering("tf")#one_hot编码转换def one_hot(Train_Y, Test_Y):    Train_Y = np.array(Train_Y).reshape([-1, 1])    Test_Y = np.array(Test_Y).reshape([-1, 1])    Encoder = preprocessing.OneHotEncoder()    Encoder.fit(Train_Y)    Train_Y = Encoder.transform(Train_Y).toarray()    Test_Y = Encoder.transform(Test_Y).toarray()    Train_Y = np.asarray(Train_Y, dtype=np.int32)    Test_Y = np.asarray(Test_Y, dtype=np.int32)    return Train_Y, Test_Yy_train, y_test = one_hot(y_train, y_test)

接下来需要根据训练参数进行调整,假设训练参数设置如下

# 训练参数batch_size = 128epochs = 40 #训练轮数num_classes = 6 #总共的训练类数length = 2048BatchNorm = False # 是否批量归一化number = 1000 # 每类样本的数量normal = False # 是否标准化

 重塑训练参数,否则传入模型时会出错

x_train=x_train.reshape((x_train.shape[0],x_train.shape[1],1))x_test = x_test.reshape((x_test.shape[0], x_test.shape[1],1))

需要根据betch_size的大小改变选择的训练数据数目,否则传入模型训练会出错

# # 改变dataset的大小,变成batch_size的倍数def change_dataset_size(x, y, batch_size):    length = len(x)    if (length % batch_size != 0):        remainder = length % batch_size        x = x[:(length - remainder)]        y = y[:(length - remainder)]    return x, yx_train,y_train=change_dataset_size(x_train,y_train,batch_size)x_valid, y_valid=change_dataset_size(x_valid, y_valid,batch_size)

这一步之后基本就可以直接传入模型进行训练了 

来源地址:https://blog.csdn.net/weixin_53488599/article/details/130553461

阅读原文内容投诉

免责声明:

① 本站未注明“稿件来源”的信息均来自网络整理。其文字、图片和音视频稿件的所属权归原作者所有。本站收集整理出于非商业性的教育和科研之目的,并不意味着本站赞同其观点或证实其内容的真实性。仅作为临时的测试数据,供内部测试之用。本站并未授权任何人以任何方式主动获取本站任何信息。

② 本站未注明“稿件来源”的临时测试数据将在测试完成后最终做删除处理。有问题或投稿请发送至: 邮箱/279061341@qq.com QQ/279061341

软考中级精品资料免费领

  • 历年真题答案解析
  • 备考技巧名师总结
  • 高频考点精准押题
  • 2024年上半年信息系统项目管理师第二批次真题及答案解析(完整版)

    难度     807人已做
    查看
  • 【考后总结】2024年5月26日信息系统项目管理师第2批次考情分析

    难度     351人已做
    查看
  • 【考后总结】2024年5月25日信息系统项目管理师第1批次考情分析

    难度     314人已做
    查看
  • 2024年上半年软考高项第一、二批次真题考点汇总(完整版)

    难度     433人已做
    查看
  • 2024年上半年系统架构设计师考试综合知识真题

    难度     221人已做
    查看

相关文章

发现更多好内容

猜你喜欢

AI推送时光机
位置:首页-资讯-后端开发
咦!没有更多了?去看看其它编程学习网 内容吧
首页课程
资料下载
问答资讯