--- a
+++ b/preprocessData.m
@@ -0,0 +1,69 @@
+function preprocessData(filename_id, rhythm, filename)
+% add preprocessing steps
+    % Baseline wanderer removal
+    % Noise removal
+    % Make every 12 lead data as clean as possible
+    % create a datastore based on the rhythm label on data
+    fs = 500;
+    table = readtable(filename_id);
+    [b,a] = butter(5, 45/fs); %for filtering of the power line noises frequncy at 60 Hz
+    [d,c] = butter(5, 10/fs, 'high');  %for filtering of baseline wander 
+    processed_data = zeros(size(table));
+    % doing of butterworth filtering from both side to remove the phase
+    % shift problem 
+    for i = 1:12
+        data = table2array(table(:,i));
+        data_1 = filter(b, a, data);
+        data_2 = filter(b, a, flip(data_1,1));
+        data_2 = flip(data_2,1);
+        data_3 = filter(d, c, data_2);
+        data_4 = filter(d, c, flip(data_3,1));
+        data_4 = flip(data_4,1);
+        mask = smooth(smooth(data_4,300),10);
+        data_4 = data_4 - mask;
+        processed_data(:,i) = data_4;
+    end
+    
+    rhythm = cell2mat(rhythm);
+    if rhythm == "SB"
+        filepath = "C:\Users\tripats\Documents\Biomedical Signal Analysis\Grad Single Project\ECG\SinusBradycardia";
+    elseif rhythm == "SR"
+        filepath = "C:\Users\tripats\Documents\Biomedical Signal Analysis\Grad Single Project\ECG\SinusRhythm";
+    elseif rhythm == "AFIB"
+        filepath = "C:\Users\tripats\Documents\Biomedical Signal Analysis\Grad Single Project\ECG\AtrialFibrilation";
+    elseif rhythm == "ST"
+        filepath = "C:\Users\tripats\Documents\Biomedical Signal Analysis\Grad Single Project\ECG\SinusTachycardia";
+    elseif rhythm == "AF"
+        filepath = "C:\Users\tripats\Documents\Biomedical Signal Analysis\Grad Single Project\ECG\AtrialFlutter";
+    elseif rhythm == "SA"
+        filepath = "C:\Users\tripats\Documents\Biomedical Signal Analysis\Grad Single Project\ECG\SinusIrregularity";
+    elseif rhythm == "SVT"
+        filepath = "C:\Users\tripats\Documents\Biomedical Signal Analysis\Grad Single Project\ECG\SupraventricularTachycardia";
+    elseif rhythm == "AT"
+        filepath = "C:\Users\tripats\Documents\Biomedical Signal Analysis\Grad Single Project\ECG\AtrialTachycardia";   
+    elseif rhythm == "AVNRT"
+        filepath = "C:\Users\tripats\Documents\Biomedical Signal Analysis\Grad Single Project\ECG\AtrioventricularNodeReentrantTachycardia";
+    elseif rhythm == "AVRT"
+        filepath = "C:\Users\tripats\Documents\Biomedical Signal Analysis\Grad Single Project\ECG\AtrioventricularReentrantTachycardia";
+    elseif rhythm == "SAAWR"
+        filepath = "C:\Users\tripats\Documents\Biomedical Signal Analysis\Grad Single Project\ECG\SinusAtriumtoAtrialWanderingRhythm";
+    end
+
+    EKG = processed_data;
+    Leadc = {'I','II','III','aV_R','aV_L','aV_F','V_1','V_2','V_3','V_4','V_5','V_6'};
+    figu = figure('Visible','off');
+    e = 1;
+    for f = 1:4
+        for g = 1:4:9
+            sbpt = (f-1)+g;
+            subplot(3,4,sbpt)
+            plot(EKG(1:1500,e))
+            title(Leadc{e})
+            grid
+            e = e+1;
+        end
+    end
+    
+    figu.Position = [0 0 512 512];
+    saveas(figu, fullfile(strcat(filepath,'\',filename,'.png')))
+end