Gamma-Ray Bursts (GRBs) have been conventionally bifurcated into two distinct categories: ``short'' and ``long'' with durations less than and greater than two seconds respectively. However, there is a lot of literature (although with conflicting results) regarding the existence of a third intermediate class. To investigate this issue, we extend a recent study (arXiv:1612.08235) on classification of GRBs to two dimensions by incorporating the GRB hardness in addition to the observed durations. We carry out this unified analysis on GRB datasets from four detectors, viz. BATSE, RHESSI, Swift (observed and intrinsic frame), and Fermi-GBM. We consider the duration and hardness features in log-scale for each of these datasets and determine the best-fit parameters using Gaussian Mixture Model. This is followed by information theoretic criterion (AIC and BIC) to determine if a three-component fit is favored compared to a two-component one or vice-versa. For BATSE, we find that both AIC and BIC show preference for three components with decisive significance. For Fermi and RHESSI, both AIC and BIC show preference for two components, although the significance is marginal from AIC, but decisive using BIC. For Swift dataset in both the observed and rest frame, we find that three components are favored according to AIC with decisive significance, and two are preferred with BIC with marginal to strong significance.