Add files via upload

Author: dineshramsai

DSP-Project-Final/18BEC0042-Subsystem-1/Subsystem - 1/Confusion Chart.PNG

@@ -0,0 +1,64 @@
+
+%18BEC0042 - T.DINESH RAM SAI
+
+function ecg2cwtscg(ecgdata,cwtfb,ecgtype)
+nos=10; %number of signals pieces. I.e, there are 65536 samples. We cut them into 500 sample blocks, i.e each block contains 500 samples each. We take first 10 blocks
+nol=500; %signal length
+colormap=jet(128); %Defining the colormap for scalogram images
+if strcmp(ecgtype, 'ARR') %equivalent to ecgtype=='ARR'
+ folder=strcat('C:\Users\DINESH RAM SAI\Desktop\Subsystem-1\ecgdataset\arr\arr'); %Destination Folder for saving Images
+ findx=0; %Required Variable Declaration for using in loops
+ for i=1:30 
+ indx=0; %Required Variable Declaration for using in loops
+ for k=1:nos
+ ecgsignal=ecgdata(i,indx+1:indx+nol); %Extracting the ECG Signal from the ecgdata matrix. It equivaletly extracts 500 sample blocks
+ cfs = abs(cwtfb.wt(ecgsignal)); %Finding the wavelet coefficients using filter 'fb'. abs() is used because coefficients are complex.
+ im = ind2rgb(im2uint8(rescale(cfs)),colormap); %by this command, we are converting the images into scalogram images. 
+ %First rescaling the coefficients from 0 to 1.Then converting them into unsigned interger 8 so that it becomes image
+ %and converting this index mode into rgb color mode 
+ filenameindex=findx+k; %Required Variable Declaration for using in naming the images
+ filename=strcat(folder,sprintf('%d.jpg',filenameindex)); %Naming the images. D is the loop execution number. Hence for first time execution, it produces 1.jpg
+ imwrite(imresize(im,[227 227]),filename); %Stored in the current directory using imwrite. Before storing, it changes image dimensions to 227x227 px which is applicable to AlexNet
+ indx=indx+nol; %Incrementing the value of appropriate variable for running the next iteration
+ end
+ findx=findx+nos; %Incrementing the value of appropriate variable for running the next iteration
+ end
+elseif strcmp(ecgtype, 'CHF') %Procedure is same, except the ECG Datatype is CHF
+ folder=strcat('C:\Users\DINESH RAM SAI\Desktop\Subsystem-1\ecgdataset\chf\chf');
+ findx=0;
+ for i=1:30
+ indx=0;
+ for k=1:nos
+ ecgsignal=ecgdata(i,indx+1:indx+nol);
+ cfs = abs(cwtfb.wt(ecgsignal));
+ im = ind2rgb(im2uint8(rescale(cfs)),colormap);
+ filenameindex=findx+k;
+ filename=strcat(folder,sprintf('%d.jpg',filenameindex));
+ imwrite(imresize(im,[227 227]),filename);
+ indx=indx+nol;
+ end
+ findx=findx+nos;
+ end
+elseif strcmp(ecgtype, 'NSR')%Procedure is same, except the ECG Datatype is NSR
+ folder=strcat('C:\Users\DINESH RAM SAI\Desktop\Subsystem-1\ecgdataset\nsr\nsr');
+ findx=0;
+ for i=1:30
+ indx=0;
+ for k=1:nos
+ ecgsignal=ecgdata(i,indx+1:indx+nol);
+ cfs = abs(cwtfb.wt(ecgsignal));
+ im = ind2rgb(im2uint8(rescale(cfs)),colormap);
+ filenameindex=findx+k;
+ filename=strcat(folder,sprintf('%d.jpg',filenameindex));
+ imwrite(imresize(im,[227 227]),filename);
+ indx=indx+nol;
+ end
+ findx=findx+nos;
+ end
+end
+
+
+
+
+
+

DSP-Project-Final/18BEC0042-Subsystem-1/Subsystem - 1/ecg2cwtscg.m

@@ -0,0 +1,36 @@
+
+%18BEC0042 - T.DINESH RAM SAI
+
+clc; clear all; close all;
+% create CWT Image Database from ECG Signals
+% Loading ECG Data
+%load(fullfile('C:','Users','DINESH RAM SAI','Documents','VIT', '5th Semester - Fall Sem 2020 - 2021','Digital Signal Processing (ECE2006)','Project','Project Files','1 - ECG Signal Classification','ECGData.mat'));
+load('ECGData.mat'); % Loading the ECG Data
+data=ECGData.Data; %Getting Database from the loaded Data. It is a 162x65536 matrix which means it consists of 162 ECG Signals each with sample size of 65536
+labels=ECGData.Labels; %Getting Labels from the loaded Data. Name of the Signals Can be get from Labels
+
+ARR=data(1:30,:); %Taken first 30 Recordings of each ECG type
+CHF=data(97:126,:);
+NSR=data(127:156,:);
+signallength=500; %Limiting the Sample Size of Signal to 500 instead of 65536
+
+%Defining filterbanks for CWT with amor wavelet and 12 filters per octave
+fb=cwtfilterbank('SignalLength',signallength,'Wavelet','amor','VoicesPerOctave',12); %used to define cwt coefficients.
+
+%Making Folders
+
+mkdir('ecgdataset'); %MainFolder
+mkdir('ecgdataset\arr'); %SubFolder
+mkdir('ecgdataset\chf');
+mkdir('ecgdataset\nsr');
+
+
+
+ecgtype={'ARR','CHF','NSR'};
+
+%Function to convert ECG to Image
+ecg2cwtscg(ARR,fb,ecgtype(1)); %User Defined function to convert ecg signal into scalogram image using cwt
+ecg2cwtscg(CHF,fb,ecgtype(2));
+ecg2cwtscg(NSR,fb,ecgtype(3));
+
+