Hello All,
I want to filter some seismic data I have like an array (27*4352) elements and I want to filter this data with Pass-Band filter so only Frequency between 100-2000 Hz will access I was reading a lot of codes and functions like Butterworth or filter(filter) that could help but I'm not able to perform the filter untill now? any suggestions to do that!
I use trying something but not sure!
%addpath('~/segymat1.6')
% read the data set in mtx files
filename = ('shot_1.mtx')
seismogram = readSeismogram(filename);
DT= 0.00003125;
T=1*DT:DT:size(seismogram,2)*DT;
Fs = 1/DT;
T = 0.135;
ns=round(T/DT);
% design filter
[b_pass, a_pass] = butter(2,[100,2000]/(Fs/2), 'bandpass');
%freqz(b_pass, a_pass)
plot(filtfilt(b_pass, a_pass, seismogram));