鱼C论坛

 找回密码
 立即注册
查看: 720|回复: 3

InceptionTime

[复制链接]
发表于 2023-10-12 16:18:11 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能^_^

您需要 登录 才可以下载或查看,没有账号?立即注册

x
现在我有一个InceptionTime程序,如何使用我的数据运行这个程序.train数据集作为训练集训练模型,数据的第一列是类别,第二列到最后一列是光谱数据,一共有8个类别
# resnet model
import keras
import numpy as np
import time

from utils.utils import save_logs
from utils.utils import calculate_metrics
from utils.utils import save_test_duration


class Classifier_INCEPTION:

    def __init__(self, output_directory, input_shape, nb_classes, verbose=False, build=True, batch_size=64,
                 nb_filters=32, use_residual=True, use_bottleneck=True, depth=6, kernel_size=41, nb_epochs=1500):

        self.output_directory = output_directory

        self.nb_filters = nb_filters
        self.use_residual = use_residual
        self.use_bottleneck = use_bottleneck
        self.depth = depth
        self.kernel_size = kernel_size - 1
        self.callbacks = None
        self.batch_size = batch_size
        self.bottleneck_size = 32
        self.nb_epochs = nb_epochs

        if build == True:
            self.model = self.build_model(input_shape, nb_classes)
            if (verbose == True):
                self.model.summary()
            self.verbose = verbose
            self.model.save_weights(self.output_directory + 'model_init.hdf5')

    def _inception_module(self, input_tensor, stride=1, activation='linear'):

        if self.use_bottleneck and int(input_tensor.shape[-1]) > 1:
            input_inception = keras.layers.Conv1D(filters=self.bottleneck_size, kernel_size=1,
                                                  padding='same', activation=activation, use_bias=False)(input_tensor)
        else:
            input_inception = input_tensor

        # kernel_size_s = [3, 5, 8, 11, 17]
        kernel_size_s = [self.kernel_size // (2 ** i) for i in range(3)]

        conv_list = []

        for i in range(len(kernel_size_s)):
            conv_list.append(keras.layers.Conv1D(filters=self.nb_filters, kernel_size=kernel_size_s[i],
                                                 strides=stride, padding='same', activation=activation, use_bias=False)(
                input_inception))

        max_pool_1 = keras.layers.MaxPool1D(pool_size=3, strides=stride, padding='same')(input_tensor)

        conv_6 = keras.layers.Conv1D(filters=self.nb_filters, kernel_size=1,
                                     padding='same', activation=activation, use_bias=False)(max_pool_1)

        conv_list.append(conv_6)

        x = keras.layers.Concatenate(axis=2)(conv_list)
        x = keras.layers.BatchNormalization()(x)
        x = keras.layers.Activation(activation='relu')(x)
        return x

    def _shortcut_layer(self, input_tensor, out_tensor):
        shortcut_y = keras.layers.Conv1D(filters=int(out_tensor.shape[-1]), kernel_size=1,
                                         padding='same', use_bias=False)(input_tensor)
        shortcut_y = keras.layers.normalization.BatchNormalization()(shortcut_y)

        x = keras.layers.Add()([shortcut_y, out_tensor])
        x = keras.layers.Activation('relu')(x)
        return x

    def build_model(self, input_shape, nb_classes):
        input_layer = keras.layers.Input(input_shape)

        x = input_layer
        input_res = input_layer

        for d in range(self.depth):

            x = self._inception_module(x)

            if self.use_residual and d % 3 == 2:
                x = self._shortcut_layer(input_res, x)
                input_res = x

        gap_layer = keras.layers.GlobalAveragePooling1D()(x)

        output_layer = keras.layers.Dense(nb_classes, activation='softmax')(gap_layer)

        model = keras.models.Model(inputs=input_layer, outputs=output_layer)

        model.compile(loss='categorical_crossentropy', optimizer=keras.optimizers.Adam(),
                      metrics=['accuracy'])

        reduce_lr = keras.callbacks.ReduceLROnPlateau(monitor='loss', factor=0.5, patience=50,
                                                      min_lr=0.0001)

        file_path = self.output_directory + 'best_model.hdf5'

        model_checkpoint = keras.callbacks.ModelCheckpoint(filepath=file_path, monitor='loss',
                                                           save_best_only=True)

        self.callbacks = [reduce_lr, model_checkpoint]

        return model

    def fit(self, x_train, y_train, x_val, y_val, y_true, plot_test_acc=False):
        if len(keras.backend.tensorflow_backend._get_available_gpus()) == 0:
            print('error no gpu')
            exit()
        # x_val and y_val are only used to monitor the test loss and NOT for training

        if self.batch_size is None:
            mini_batch_size = int(min(x_train.shape[0] / 10, 16))
        else:
            mini_batch_size = self.batch_size

        start_time = time.time()

        if plot_test_acc:

            hist = self.model.fit(x_train, y_train, batch_size=mini_batch_size, epochs=self.nb_epochs,
                                  verbose=self.verbose, validation_data=(x_val, y_val), callbacks=self.callbacks)
        else:

            hist = self.model.fit(x_train, y_train, batch_size=mini_batch_size, epochs=self.nb_epochs,
                                  verbose=self.verbose, callbacks=self.callbacks)

        duration = time.time() - start_time

        self.model.save(self.output_directory + 'last_model.hdf5')

        y_pred = self.predict(x_val, y_true, x_train, y_train, y_val,
                              return_df_metrics=False)

        # save predictions
        np.save(self.output_directory + 'y_pred.npy', y_pred)

        # convert the predicted from binary to integer
        y_pred = np.argmax(y_pred, axis=1)

        df_metrics = save_logs(self.output_directory, hist, y_pred, y_true, duration,
                               plot_test_acc=plot_test_acc)

        keras.backend.clear_session()

        return df_metrics

    def predict(self, x_test, y_true, x_train, y_train, y_test, return_df_metrics=True):
        start_time = time.time()
        model_path = self.output_directory + 'best_model.hdf5'
        model = keras.models.load_model(model_path)
        y_pred = model.predict(x_test, batch_size=self.batch_size)
        if return_df_metrics:
            y_pred = np.argmax(y_pred, axis=1)
            df_metrics = calculate_metrics(y_true, y_pred, 0.0)
            return df_metrics
        else:
            test_duration = time.time() - start_time
            save_test_duration(self.output_directory + 'test_duration.csv', test_duration)
            return y_pred
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复

使用道具 举报

发表于 2023-10-12 16:18:58 | 显示全部楼层
要使用InceptionTime程序运行数据,需要将数据集准备成合适的格式并调用fit方法。首先需要将数据集拆分成训练集和验证集,然后将它们都转换成numpy数组的形式。假设训练集和验证集分别为x_train和x_val,对应的标签为y_train和y_val,代码如下:
import numpy as np

# 加载数据集,假设数据集已经处理好并保存在numpy数组中
data = np.load('data.npy')
x = data[:, 1:]  # 特征数据
y = data[:, 0]   # 标签数据

# 将数据集拆分成训练集和验证集
val_split = 0.1  # 验证集比例
val_size = int(len(x) * val_split)
x_val = x[:val_size]
y_val = y[:val_size]
x_train = x[val_size:]
y_train = y[val_size:]

# 调用fit方法进行训练
from Classifier_INCEPTION import Classifier_INCEPTION

output_directory = 'output/'  # 输出目录
input_shape = (x_train.shape[1], 1)  # 输入形状
nb_classes = 8  # 类别数

classifier = Classifier_INCEPTION(output_directory, input_shape, nb_classes)
classifier.fit(x_train, y_train, x_val, y_val, y_true=y_val)

在调用fit方法之前,需要先实例化一个Classifier_INCEPTION对象,并传入输出目录、输入形状和类别数等参数。然后就可以调用fit方法进行训练了。在训练过程中,程序会自动保存最佳模型,并将训练日志保存到输出目录中。训练完成后,可以调用predict方法对测试集进行预测,代码如下:
# 加载测试集数据
x_test = np.load('x_test.npy')
y_test = np.load('y_test.npy')

# 调用predict方法进行预测
df_metrics = classifier.predict(x_test, y_true=y_test)
print(df_metrics)

在调用predict方法之前,需要先加载测试集数据,并将其转换成numpy数组的形式。然后就可以调用predict方法进行预测了。predict方法会返回一个DataFrame对象,包含了预测结果的各种指标。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2023-10-12 16:19:39 | 显示全部楼层
你的InceptionTime模型代码已经准备好了,现在你需要将你的数据集加载到模型中进行训练。以下是一种可能的方法:

首先,你需要将你的数据集读入Python。假设你的数据集是CSV格式,你可以使用pandas库来读取它。然后,你需要将数据集分割为特征和标签,然后将标签进行one-hot编码,以便于模型的训练。
import pandas as pd
from keras.utils import to_categorical

# 读取数据集
data = pd.read_csv('your_dataset.csv')

# 分割特征和标签
X = data.iloc[:, 1:].values
y = data.iloc[:, 0].values

# 将标签进行one-hot编码
y = to_categorical(y, num_classes=8)

接下来,你需要将数据集分割为训练集和验证集。你可以使用sklearn库的train_test_split函数来实现这一点。
from sklearn.model_selection import train_test_split

# 分割数据集
X_train, X_val, y_train, y_val = train_test_split(X, y, test_size=0.2, random_state=42)

然后,你可以创建你的InceptionTime模型,并使用你的数据来训练它。
# 创建模型
model = Classifier_INCEPTION(output_directory='output_directory', input_shape=X_train.shape[1:], nb_classes=8)

# 训练模型
model.fit(X_train, y_train, X_val, y_val)

注意,你需要将上述代码中的'your_dataset.csv'和'output_directory'替换为你的数据集的实际路径和你希望保存模型的实际路径。

此外,你可能需要根据你的数据集和你的需求调整数据预处理和模型训练的代码。例如,你可能需要调整train_test_split函数的test_size参数来改变训练集和验证集的大小,或者调整Classifier_INCEPTION类的参数来改变模型的配置。
请注意,本回答来自GPT-4
如果问题已经解决,请设置最佳答案
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2023-10-12 16:20:00 | 显示全部楼层
首先,你需要准备好你的训练数据集。数据集应该是一个二维数组,其中第一列是类别,第二列到最后一列是光谱数据。你提到共有8个类别。

接下来,在你的代码中找到以下几行:
model = Classifier_INCEPTION(output_directory, input_shape, nb_classes)model.fit(x_train, y_train, x_val, y_val, y_true)
将上面的代码替换为以下内容:
from sklearn.model_selection import train_test_splitfrom sklearn.preprocessing import LabelBinarizer# 准备数据X = your_data[:, 1:]  # 提取光谱数据y = your_data[:, 0]   # 提取类别# 使用LabelBinarizer对类别进行one-hot编码lb = LabelBinarizer()y = lb.fit_transform(y)# 划分训练集和验证集x_train, x_val, y_train, y_val = train_test_split(X, y, test_size=0.2, random_state=42)# 创建并训练模型model = Classifier_INCEPTION(output_directory, input_shape, nb_classes)model.fit(x_train, y_train, x_val, y_val, y_true)
以上代码中,我们使用sklearn库中的train_test_split函数将数据划分为训练集和验证集,并使用LabelBinarizer对类别进行one-hot编码。然后,我们使用划分好的数据集来训练模型。

记得根据你的实际数据调整相关参数,例如test_size可以调整验证集大小。

希望对你有所帮助!如有任何进一步的问题,请随时提问。

球一个最佳答案谢谢啦!这对我非常重要!
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2024-9-21 11:12

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表