Deep Learning algorithms, in particular neural networks have been steadily gaining popularity among the gravitational wave community for the last few years. The reliability and accuracy of Deep Learning approaches in gravitational wave detection, parameter estimation and glitch classification have already been proved and verified by several groups in recent years. In this paper, we report on the construction of a deep Artificial Neural Network (ANN) to localize simulated gravitational wave signals in the sky with high accuracy. We have modelled the sky as a sphere and have considered cases where the sphere is divided into 18, 50, 128, 1024, 2048 and 4096 sectors. The sky direction of the gravitational wave source is estimated by classifying the signal into one of these sectors based on it's right ascension and declination values for each of these cases. In order to do this, we have injected simulated binary black hole gravitational wave signals of component masses sampled uniformly between 30-80 M into Gaussian noise and used the whitened strain values to obtain the input features for training our ANN. We input features such as the delays in arrival times, phase differences and amplitude ratios at each of the three detectors Hanford, Livingston and Virgo, from the raw time-domain strain values as well as from analytical versions of these signals, obtained through Hilbert transformation. We show that our model is able to classify gravitational wave samples, not used in the training process, into their correct sectors with very high accuracy (> 90%) for coarse angular resolution using 18, 50 and 128 sectors. We also test our localization on test samples with injection parameters of the published LIGO binary black hole merger events GW150914, GW170818 and GW170823 for 1024, 2048 and 4096 sectors and compare the result with that from BAYESTAR and Parameter Estimation (PE). In addition, we report that the time taken by our model to localize one GW signal is around 0.018 secs on 14 Intel Xeon CPU cores. composition of these compact objects. [14][15][16][17][18][19].LIGO and Virgo began their third observation run (O3) in April, 2019 and have already detected several new GW signals. The technique used by search pipelines to identify and characterize GW signals from noise is matched filtering, which uses a bank of template waveforms of different component masses and spins [20][21][22][23][24][25][26][27]. These waveform models cover the inspiral, merger and ringdown phases of a compact binary coalescence and are generated by combining post-Newtonian theory [28][29][30][31], the effective-one-body formalism [32] and numerical relativity simulations [33]. While this approach has been very successful, the algorithms used are computationally expensive due to the fact that they probe a large parameter space with increasingly longerduration waveforms, as the low-frequency sensitivity of the interferometers improves. Also, future observation runs will involve more detectors running simultaneously over longer periods of time, th...