the clean seismic data which is 1*2500 were taken as the training labels of the resnet18 in matlab

main.m
% Load data
data = load('seismic_data.mat');
seismic_data = data.seismic_data;

% Preprocessing
seismic_data = detrend(seismic_data);
seismic_data = seismic_data - mean(seismic_data);

% Bandpass filter
fs = 2500; % Sampling frequency
f1 = 1;    % Lower cutoff frequency
f2 = 50;   % Upper cutoff frequency
[b, a] = butter(4, [f1, f2]/(fs/2), 'bandpass');
seismic_data_filt = filtfilt(b, a, seismic_data);

% Noise reduction
seismic_data_denoised = wdenoise(seismic_data_filt);

% Normalization
seismic_data_clean = normalize(seismic_data_denoised);
550 chars
21 lines