{"metadata":{"kernelspec":{"language":"python","display_name":"Python 3","name":"python3"},"language_info":{"pygments_lexer":"ipython3","nbconvert_exporter":"python","version":"3.6.4","file_extension":".py","codemirror_mode":{"name":"ipython","version":3},"name":"python","mimetype":"text/x-python"}},"nbformat_minor":4,"nbformat":4,"cells":[{"cell_type":"markdown","source":"## Import the libraries\n","metadata":{"execution":{"iopub.status.busy":"2022-02-11T15:15:42.09471Z","iopub.execute_input":"2022-02-11T15:15:42.095103Z","iopub.status.idle":"2022-02-11T15:15:42.099724Z","shell.execute_reply.started":"2022-02-11T15:15:42.095016Z","shell.execute_reply":"2022-02-11T15:15:42.098645Z"}}},{"cell_type":"code","source":"!pip install pyunpack\n!pip install patool\n!pip install py7zr\n!pip install sounddevice\n!pip install noisereduce\n!pip install librosa\n! pip install python_speech_features\n! pip install tensorflow==2.4\n! pip install malaya_speech\n! pip install webrtcvad","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:01:36.092384Z","iopub.execute_input":"2022-06-25T19:01:36.093181Z","iopub.status.idle":"2022-06-25T19:03:47.866198Z","shell.execute_reply.started":"2022-06-25T19:01:36.093074Z","shell.execute_reply":"2022-06-25T19:03:47.865231Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"! pip install --upgrade transformers","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:03:47.867722Z","iopub.execute_input":"2022-06-25T19:03:47.868043Z","iopub.status.idle":"2022-06-25T19:04:02.098887Z","shell.execute_reply.started":"2022-06-25T19:03:47.868001Z","shell.execute_reply":"2022-06-25T19:04:02.097975Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"pip install git+https://github.com/huggingface/transformers","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:04:02.101162Z","iopub.execute_input":"2022-06-25T19:04:02.101596Z","iopub.status.idle":"2022-06-25T19:04:43.863905Z","shell.execute_reply.started":"2022-06-25T19:04:02.101541Z","shell.execute_reply":"2022-06-25T19:04:43.862700Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"import numpy as np                        # linear algebra library\nimport pandas as pd                       # data frames processing\nimport matplotlib.pyplot as plt          # visualization library\nimport seaborn as sn                      # visualization library\n\n\n# audio processing library\nimport librosa                          \nimport IPython.display as ipd            \nfrom scipy.io import wavfile\nimport noisereduce as nr\nfrom malaya_speech import Pipeline\nimport malaya_speech\nfrom python_speech_features import mfcc\nfrom sklearn.preprocessing import LabelEncoder\n\n\n# unpacking the dataset\nfrom py7zr import unpack_7zarchive    \n\n#operating system libraries\nimport shutil\nimport warnings\nwarnings.filterwarnings(\"ignore\")\nimport tensorflow \nimport os","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:04:43.865833Z","iopub.execute_input":"2022-06-25T19:04:43.866174Z","iopub.status.idle":"2022-06-25T19:04:48.875889Z","shell.execute_reply.started":"2022-06-25T19:04:43.866124Z","shell.execute_reply":"2022-06-25T19:04:48.872057Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"# unpacking the dataset\nshutil.register_unpack_format('7zip', ['.7z'], unpack_7zarchive)\nshutil.unpack_archive('/kaggle/input/tensorflow-speech-recognition-challenge/train.7z', '/kaggle/working/tensorflow-speech-recognition-challenge/train/')","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:04:48.877386Z","iopub.execute_input":"2022-06-25T19:04:48.877785Z","iopub.status.idle":"2022-06-25T19:14:01.956790Z","shell.execute_reply.started":"2022-06-25T19:04:48.877749Z","shell.execute_reply":"2022-06-25T19:14:01.955909Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# <center> Implementing the Speech Recognition Model in Python\n**Dataset used for our Speech Recognition Project**\n    \nIt is a set of 10 numbers each is repeated 2000 times with different accents and different back ground conditions. TensorFlow recently released the Speech Commands Datasets. It includes 65,000 one-second long utterances of 30 short words, by thousands of different people. We’ll build a speech recognition system that understands simple spoken commands. <br>    \n    \n__You can download the dataset from__ [here](https://www.kaggle.com/c/tensorflow-speech-recognition-challenge).\n","metadata":{}},{"cell_type":"markdown","source":"**Data Exploration and Visualization**\n\nData Exploration and Visualization helps us to understand the data as well as pre-processing steps in a better way. \n\n","metadata":{}},{"cell_type":"code","source":"train_audio_path = '/kaggle/working/tensorflow-speech-recognition-challenge/train/train/audio/' #path of the training data","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:14:01.958163Z","iopub.execute_input":"2022-06-25T19:14:01.958495Z","iopub.status.idle":"2022-06-25T19:14:01.963095Z","shell.execute_reply.started":"2022-06-25T19:14:01.958461Z","shell.execute_reply":"2022-06-25T19:14:01.962008Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"**Accessing each file in data**","metadata":{}},{"cell_type":"code","source":"samples, sample_rate = librosa.load(train_audio_path+'on/5a3712c9_nohash_1.wav', sr = 16000)  # loading a sample to be explored carefully\nfig = plt.figure(figsize=(14, 8))\nax1 = fig.add_subplot(211)\nax1.set_title('Raw wave of ' + '../input/train/audio/on/0a7c2a8d_nohash_0.wav')\nax1.set_xlabel('time')\nax1.set_ylabel('Amplitude')\nax1.plot(np.linspace(0, sample_rate/len(samples), sample_rate), samples)","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:14:01.964444Z","iopub.execute_input":"2022-06-25T19:14:01.965095Z","iopub.status.idle":"2022-06-25T19:14:02.160723Z","shell.execute_reply.started":"2022-06-25T19:14:01.965057Z","shell.execute_reply":"2022-06-25T19:14:02.158788Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"fs=16000\nipd.Audio(samples, rate=fs)  #listen to audio file before noise reduction\nprint(\"Sample Rate:\",fs)\nsr=fs","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:14:02.163356Z","iopub.execute_input":"2022-06-25T19:14:02.163674Z","iopub.status.idle":"2022-06-25T19:14:02.169851Z","shell.execute_reply.started":"2022-06-25T19:14:02.163646Z","shell.execute_reply":"2022-06-25T19:14:02.168685Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"# Noise Reduction\ntime = np.linspace(0, len(samples - 1) / fs, len(samples - 1))\nreduced_noise1 = nr.reduce_noise(y=samples, sr=fs,stationary=True)\nplt.plot(time, reduced_noise1)  # plot in seconds\nplt.xlabel(\"Time [seconds]\")\nplt.ylabel(\"Voice amplitude\")\nplt.show()","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:14:02.172118Z","iopub.execute_input":"2022-06-25T19:14:02.172690Z","iopub.status.idle":"2022-06-25T19:14:02.890458Z","shell.execute_reply.started":"2022-06-25T19:14:02.172654Z","shell.execute_reply":"2022-06-25T19:14:02.889552Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"ipd.Audio(reduced_noise1, rate=sample_rate)  #listen to audio file after noise reduction","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:14:02.891812Z","iopub.execute_input":"2022-06-25T19:14:02.892289Z","iopub.status.idle":"2022-06-25T19:14:02.899582Z","shell.execute_reply.started":"2022-06-25T19:14:02.892248Z","shell.execute_reply":"2022-06-25T19:14:02.898811Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"#Silence Removal\nvad = malaya_speech.vad.webrtc()\ny=reduced_noise1\ny_= malaya_speech.resample(y, sr, 16000)\ny_ = malaya_speech.astype.float_to_int(y_)\nframes = malaya_speech.generator.frames(y, 30, sr)\nframes_ = list(malaya_speech.generator.frames(y_, 30, 16000, append_ending_trail = False))\nframes_webrtc = [(frames[no], vad(frame)) for no, frame in enumerate(frames_)]\ny_ = malaya_speech.combine.without_silent(frames_webrtc)\ny_","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:14:02.901089Z","iopub.execute_input":"2022-06-25T19:14:02.901626Z","iopub.status.idle":"2022-06-25T19:14:03.068450Z","shell.execute_reply.started":"2022-06-25T19:14:02.901581Z","shell.execute_reply":"2022-06-25T19:14:03.067616Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"ipd.Audio(y_, rate = sr )     #listen to audio file after noise reduction","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:14:03.069661Z","iopub.execute_input":"2022-06-25T19:14:03.070066Z","iopub.status.idle":"2022-06-25T19:14:03.079190Z","shell.execute_reply.started":"2022-06-25T19:14:03.070024Z","shell.execute_reply":"2022-06-25T19:14:03.078261Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"# padding signal with zeros\nzero = np.zeros((1*sr-y_.shape[0]))\nsignal = np.concatenate((y_,zero))\nsignal.shape\ntime = np.linspace(0, len(signal - 1) / fs, len(signal - 1))","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:14:03.080991Z","iopub.execute_input":"2022-06-25T19:14:03.081552Z","iopub.status.idle":"2022-06-25T19:14:03.086895Z","shell.execute_reply.started":"2022-06-25T19:14:03.081491Z","shell.execute_reply":"2022-06-25T19:14:03.086075Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"**Visualization of Audio signal in time series domain**\n\nNow, we’ll visualize the audio signal in the time series domain:","metadata":{}},{"cell_type":"code","source":"plt.plot(time,signal)","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:14:03.088192Z","iopub.execute_input":"2022-06-25T19:14:03.088744Z","iopub.status.idle":"2022-06-25T19:14:03.239153Z","shell.execute_reply.started":"2022-06-25T19:14:03.088709Z","shell.execute_reply":"2022-06-25T19:14:03.238302Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"labels=os.listdir(train_audio_path)     #Extracting labels to determine classes","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:14:03.240462Z","iopub.execute_input":"2022-06-25T19:14:03.241014Z","iopub.status.idle":"2022-06-25T19:14:03.246312Z","shell.execute_reply.started":"2022-06-25T19:14:03.240974Z","shell.execute_reply":"2022-06-25T19:14:03.245571Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"## **Distribution of the Data set**","metadata":{}},{"cell_type":"code","source":"#find count of each label and plot bar graph\nno_of_recordings=[]\nfor label in labels:\n    waves = [f for f in os.listdir(train_audio_path + '/'+ label) if f.endswith('.wav')]\n    no_of_recordings.append(len(waves))\n    \n#plot\nplt.figure(figsize=(30,5))\nindex = np.arange(len(labels))\nplt.bar(index, no_of_recordings)\nplt.xlabel('Commands', fontsize=12)\nplt.ylabel('No of recordings', fontsize=12)\nplt.xticks(index, labels, fontsize=15, rotation=60)\nplt.title('No. of recordings for each command')\nplt.show()","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:14:03.247758Z","iopub.execute_input":"2022-06-25T19:14:03.248185Z","iopub.status.idle":"2022-06-25T19:14:03.591038Z","shell.execute_reply.started":"2022-06-25T19:14:03.248149Z","shell.execute_reply":"2022-06-25T19:14:03.590096Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"## Words used","metadata":{}},{"cell_type":"code","source":"labels=[\"zero\",\"one\",\"two\",\"three\",\"four\",\"five\",\"six\",\"seven\",\"eight\",\"nine\"]","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:14:03.592426Z","iopub.execute_input":"2022-06-25T19:14:03.593050Z","iopub.status.idle":"2022-06-25T19:14:03.598211Z","shell.execute_reply.started":"2022-06-25T19:14:03.593005Z","shell.execute_reply":"2022-06-25T19:14:03.596882Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"## **Preprocessing the audio waves**\n\nlet us read the audio waves and use the below-preprocessing steps :\n\n* Noise Reduction\n* Silence Removal\n* Extracting MFCCs\n\nLet us define these preprocessing steps in the below code snippet:","metadata":{}},{"cell_type":"code","source":"sr=16000    # sample rate\nvad = malaya_speech.vad.webrtc()\nall_wave = []     #intitialize array to stack wave files of the whole data set in it \nall_label = []    #intitialize array to stack label of wave files of the whole data set in it \nfor label in labels:\n    print(label)\n    waves = [f for f in os.listdir(train_audio_path + '/'+ label) if f.endswith('.wav')] # access on each file\n    for wav in waves:\n        samples, sample_rate = librosa.load(train_audio_path + '/' + label + '/' + wav, sr = 16000)\n        samples = nr.reduce_noise(y=samples, sr=sr,stationary=True)  #noise reduction\n        y_= malaya_speech.resample(samples, sr, 16000)               # silence removal\n        y_ = malaya_speech.astype.float_to_int(y_)\n        frames = malaya_speech.generator.frames(samples, 30, sr)\n        frames_ = list(malaya_speech.generator.frames(y_, 30, 16000, append_ending_trail = False))\n        frames_webrtc = [(frames[no], vad(frame)) for no, frame in enumerate(frames_)]\n        y_ = malaya_speech.combine.without_silent(frames_webrtc)\n        zero = np.zeros(((1*sr+4000)-y_.shape[0]))                 \n        signal = np.concatenate((y_,zero))     # concatenation with zeros to adust length of the vector\n        all_wave.append(signal)     #append waves one by one \n        all_label.append(label)     #append corresponding label one by one","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:14:03.599733Z","iopub.execute_input":"2022-06-25T19:14:03.600219Z","iopub.status.idle":"2022-06-25T19:25:08.804837Z","shell.execute_reply.started":"2022-06-25T19:14:03.600177Z","shell.execute_reply":"2022-06-25T19:25:08.803993Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"print(\"shape of waves array\",np.array(all_wave).shape)\nprint(\"shape of labels array\",np.array(all_label).shape)\n\n#Inspecting random sample\ntime = np.linspace(0, len(signal - 1) / fs, len(signal - 1))\nplt.plot(time,np.array(all_wave)[2000,:])\nprint(np.array(all_label)[2000])\nipd.Audio(np.array(all_wave)[2000,:], rate = sr )","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:25:08.806081Z","iopub.execute_input":"2022-06-25T19:25:08.806406Z","iopub.status.idle":"2022-06-25T19:25:19.699167Z","shell.execute_reply.started":"2022-06-25T19:25:08.806372Z","shell.execute_reply":"2022-06-25T19:25:19.698387Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"all_mfcc=[]    #intitialize array to stack MFCCs of the whole data set in it\nfor wave in all_wave:\n    i=0\n    mfcc_feat = mfcc(wave , fs, winlen=256/fs, winstep=256/(2*fs), numcep=13, nfilt=26, nfft=256,\n                 lowfreq=0, highfreq=fs/2, preemph=0.97, ceplifter=22, appendEnergy=True, winfunc=np.hamming)\n    mfcc_feat= np.transpose(mfcc_feat)\n    all_mfcc.append(mfcc_feat)","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:25:19.700423Z","iopub.execute_input":"2022-06-25T19:25:19.700785Z","iopub.status.idle":"2022-06-25T19:26:31.892684Z","shell.execute_reply.started":"2022-06-25T19:25:19.700747Z","shell.execute_reply":"2022-06-25T19:26:31.891666Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"\nprint(\"shape of MFCCs\",np.array(all_mfcc).shape)\nprint(\"shape of Corresponding lables\",np.array(all_label).shape)\n\n# dimensions of the data\nd1=np.array(all_mfcc).shape[1]\nd2=np.array(all_mfcc).shape[2]\nd=d1*d2","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:26:31.897760Z","iopub.execute_input":"2022-06-25T19:26:31.900187Z","iopub.status.idle":"2022-06-25T19:26:32.417156Z","shell.execute_reply.started":"2022-06-25T19:26:31.900134Z","shell.execute_reply":"2022-06-25T19:26:32.416324Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"op_mfcc=np.array(all_mfcc)\n#op_mfcc=op_mfcc.reshape(np.array(all_mfcc).shape[0],-1)\nop_mfcc.shape","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:26:32.418382Z","iopub.execute_input":"2022-06-25T19:26:32.418729Z","iopub.status.idle":"2022-06-25T19:26:32.578362Z","shell.execute_reply.started":"2022-06-25T19:26:32.418695Z","shell.execute_reply":"2022-06-25T19:26:32.577492Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"# Adjusting shape of vector before feeding it to the model \nop_mfcc=op_mfcc.reshape(23666,13,156,-1)\nop_mfcc.shape","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:26:32.582644Z","iopub.execute_input":"2022-06-25T19:26:32.582897Z","iopub.status.idle":"2022-06-25T19:26:32.589180Z","shell.execute_reply.started":"2022-06-25T19:26:32.582871Z","shell.execute_reply":"2022-06-25T19:26:32.588380Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### One hot encoding","metadata":{}},{"cell_type":"code","source":"#all_label = all_label.tolist()\n\nle = LabelEncoder()\ny=le.fit_transform(all_label)\nclasses= list(le.classes_)","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:26:32.591286Z","iopub.execute_input":"2022-06-25T19:26:32.591852Z","iopub.status.idle":"2022-06-25T19:26:32.607398Z","shell.execute_reply.started":"2022-06-25T19:26:32.591816Z","shell.execute_reply":"2022-06-25T19:26:32.606633Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"## **Model based on CNN** ","metadata":{}},{"cell_type":"code","source":"! pip install --upgrade tensorflow\n! pip install --upgrade tensorflow-gpu\n! pip install keras==2.3.1","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:26:32.609328Z","iopub.execute_input":"2022-06-25T19:26:32.609989Z","iopub.status.idle":"2022-06-25T19:28:40.948852Z","shell.execute_reply.started":"2022-06-25T19:26:32.609947Z","shell.execute_reply":"2022-06-25T19:28:40.947905Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"# Importing model libraries from keras\nfrom keras.optimizers import SGD\nfrom keras.constraints import maxnorm\nfrom tensorflow.keras import Sequential\nfrom keras.callbacks import EarlyStopping, ModelCheckpoint\nfrom tensorflow.keras.layers import Dense, Dropout, Flatten, Conv2D, Input, MaxPooling2D\nfrom tensorflow.keras.layers import BatchNormalization","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:28:40.952983Z","iopub.execute_input":"2022-06-25T19:28:40.953514Z","iopub.status.idle":"2022-06-25T19:28:41.005826Z","shell.execute_reply.started":"2022-06-25T19:28:40.953470Z","shell.execute_reply":"2022-06-25T19:28:41.004754Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"y=tensorflow.keras.utils.to_categorical(y, num_classes=len(labels), dtype='float32')\ny.shape","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:28:41.007228Z","iopub.execute_input":"2022-06-25T19:28:41.007817Z","iopub.status.idle":"2022-06-25T19:28:41.017399Z","shell.execute_reply.started":"2022-06-25T19:28:41.007775Z","shell.execute_reply":"2022-06-25T19:28:41.016499Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"from sklearn.model_selection import train_test_split\nx_tr, x_val, y_tr, y_val= train_test_split(op_mfcc,np.array(y),stratify=y,test_size = 0.2,random_state=777,shuffle=True)","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:28:41.018780Z","iopub.execute_input":"2022-06-25T19:28:41.019419Z","iopub.status.idle":"2022-06-25T19:28:41.791852Z","shell.execute_reply.started":"2022-06-25T19:28:41.019379Z","shell.execute_reply":"2022-06-25T19:28:41.790963Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"# Inspect shapes of each set\nprint(\"shape of training samples\",x_tr.shape)\nprint(\"shape of training labels\",y_tr.shape)\nprint(\"shape of test samples\",x_val.shape)\nprint(\"shape of test labels\",y_val.shape)","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:28:41.793223Z","iopub.execute_input":"2022-06-25T19:28:41.793575Z","iopub.status.idle":"2022-06-25T19:28:41.800858Z","shell.execute_reply.started":"2022-06-25T19:28:41.793522Z","shell.execute_reply":"2022-06-25T19:28:41.800036Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### **Model Architecture**","metadata":{}},{"cell_type":"code","source":"\n#Model Architecture\nmodel = Sequential()\nmodel.add(Conv2D(32, kernel_size=(4, 4), activation='relu', input_shape=(d1,d2,1)))\nmodel.add(BatchNormalization())\n\nmodel.add(Conv2D(48, kernel_size=(3, 3), activation='relu'))\nmodel.add(BatchNormalization())\n\nmodel.add(Conv2D(120, kernel_size=(3, 3), activation='relu'))\nmodel.add(BatchNormalization())\n\nmodel.add(MaxPooling2D(pool_size=(2, 2)))\nmodel.add(Dropout(0.5))\n\nmodel.add(Flatten())\n\nmodel.add(Dense(128, activation='relu'))\nmodel.add(BatchNormalization())\nmodel.add(Dropout(0.5))\nmodel.add(Dense(64, activation='relu'))\nmodel.add(BatchNormalization())\nmodel.add(Dropout(0.5))\nmodel.add(Dense(len(classes), activation='softmax'))","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:28:41.802089Z","iopub.execute_input":"2022-06-25T19:28:41.802649Z","iopub.status.idle":"2022-06-25T19:28:45.119178Z","shell.execute_reply.started":"2022-06-25T19:28:41.802608Z","shell.execute_reply":"2022-06-25T19:28:45.118343Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"tensorflow.keras.utils.plot_model(model, 'model.png',show_shapes=True)","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:28:45.120464Z","iopub.execute_input":"2022-06-25T19:28:45.120813Z","iopub.status.idle":"2022-06-25T19:28:46.363837Z","shell.execute_reply.started":"2022-06-25T19:28:45.120779Z","shell.execute_reply":"2022-06-25T19:28:46.362755Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"model.compile(loss='categorical_crossentropy',optimizer='adamax',metrics=['accuracy'])\n","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:28:46.365188Z","iopub.execute_input":"2022-06-25T19:28:46.365500Z","iopub.status.idle":"2022-06-25T19:28:46.411579Z","shell.execute_reply.started":"2022-06-25T19:28:46.365468Z","shell.execute_reply":"2022-06-25T19:28:46.410656Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"es = EarlyStopping(monitor='val_loss', mode='min', verbose=1, patience=10, min_delta=0.0001) \nmc = ModelCheckpoint('best_model.hdf5', monitor='val_acc', verbose=1, save_best_only=True, mode='max')","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:28:46.415480Z","iopub.execute_input":"2022-06-25T19:28:46.417449Z","iopub.status.idle":"2022-06-25T19:28:46.430639Z","shell.execute_reply.started":"2022-06-25T19:28:46.417412Z","shell.execute_reply":"2022-06-25T19:28:46.423140Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"history=model.fit(x_tr, y_tr,validation_data=(x_val,y_val), epochs=300, batch_size=64)","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:28:46.434717Z","iopub.execute_input":"2022-06-25T19:28:46.437684Z","iopub.status.idle":"2022-06-25T19:46:41.171671Z","shell.execute_reply.started":"2022-06-25T19:28:46.437647Z","shell.execute_reply":"2022-06-25T19:46:41.170772Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"## Model Evaluation","metadata":{}},{"cell_type":"code","source":"train_score = model.evaluate(x_tr, y_tr, batch_size=12)\nprint(train_score)\n\nprint('----------------Training Complete-----------------')\n\ntest_score = model.evaluate(x_val, y_val, batch_size = 12)\nprint(test_score)","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:46:41.173385Z","iopub.execute_input":"2022-06-25T19:46:41.173824Z","iopub.status.idle":"2022-06-25T19:46:47.029116Z","shell.execute_reply.started":"2022-06-25T19:46:41.173777Z","shell.execute_reply":"2022-06-25T19:46:47.028230Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"history.history.keys()","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:46:47.030456Z","iopub.execute_input":"2022-06-25T19:46:47.030822Z","iopub.status.idle":"2022-06-25T19:46:47.038402Z","shell.execute_reply.started":"2022-06-25T19:46:47.030785Z","shell.execute_reply":"2022-06-25T19:46:47.037488Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"plt.plot(history.history['loss'], label='train')  # losses learning curve of training set.\nplt.plot(history.history['val_loss'], label='test') # losses learning curve of validation set.\nplt.legend()\nplt.title(\"losses learning curves\")\nplt.show()","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:46:47.039876Z","iopub.execute_input":"2022-06-25T19:46:47.040453Z","iopub.status.idle":"2022-06-25T19:46:47.213172Z","shell.execute_reply.started":"2022-06-25T19:46:47.040417Z","shell.execute_reply":"2022-06-25T19:46:47.212354Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"plt.plot(history.history['accuracy'])      # Accuracy learning curve of training set.\nplt.plot(history.history['val_accuracy'])  # Accuracy learning curve of validation set.\nplt.title('model accuracy')\nplt.ylabel('accuracy')\nplt.xlabel('epoch')\nplt.legend(['train', 'val'], loc='upper left')\nplt.title(\"Accuracy learning curves\")\nplt.show()","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:46:47.214501Z","iopub.execute_input":"2022-06-25T19:46:47.215032Z","iopub.status.idle":"2022-06-25T19:46:47.362129Z","shell.execute_reply.started":"2022-06-25T19:46:47.214995Z","shell.execute_reply":"2022-06-25T19:46:47.361383Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"y_predict=model.predict(x_val)\nconf_mat=tensorflow.math.confusion_matrix(np.argmax(y_val,axis=1) , np.argmax(y_predict,axis=1))","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:46:47.363187Z","iopub.execute_input":"2022-06-25T19:46:47.363492Z","iopub.status.idle":"2022-06-25T19:46:48.550604Z","shell.execute_reply.started":"2022-06-25T19:46:47.363467Z","shell.execute_reply":"2022-06-25T19:46:48.549603Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"df_cm = pd.DataFrame(np.array(conf_mat), index = [i for i in classes],\n                  columns = [i for i in classes])\nplt.figure(figsize = (13,7))\nax = sn.heatmap(df_cm, annot=True)\nplt.title(\"Confusion Matrix\", fontsize=20)\nplt.ylabel(\"True Class\"     , fontsize=20)\nplt.xlabel(\"Predicted Class\", fontsize=20)\nplt.show()","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:46:48.552397Z","iopub.execute_input":"2022-06-25T19:46:48.552967Z","iopub.status.idle":"2022-06-25T19:46:49.405459Z","shell.execute_reply.started":"2022-06-25T19:46:48.552928Z","shell.execute_reply":"2022-06-25T19:46:49.404532Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"## loading model to be tested locally.","metadata":{}},{"cell_type":"code","source":"from keras.models import load_model\nmodel.save(\"numbers_cnn.h5\")    ","metadata":{"execution":{"iopub.status.busy":"2022-06-25T19:46:49.510606Z","iopub.status.idle":"2022-06-25T19:46:49.511316Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"","metadata":{"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"","metadata":{},"execution_count":null,"outputs":[]}]}