Studies of Gamma-Ray Burst (GRB) properties, such as duration and spectral hardness, have found evidence for additional classes, beyond the short/hard and long/soft prototypes, using model-dependent methods. In this paper, a model-independent approach was used to analyse the gamma-ray light curves of large samples of GRBs detected by BATSE, Swift/BAT and Fermi/GBM. All the features were extracted from the GRB time profiles in four energy bands using the Stationary Wavelet Transform and Principal Component Analysis. t-distributed Stochastic Neighbourhood Embedding (t-SNE) visualisation of the features revealed two distinct groups of Swift/BAT bursts using the T100 interval with 64 ms resolution data. When the same analysis was applied to 4 ms resolution data, two groups were seen to emerge within the first second (T1) post-trigger. These two groups primarily consisted of short/hard (Group 1) and long/soft (Group 2) bursts, and were 95% consistent with the groups identified using the T100 64 ms resolution data. Kilonova candidates, arising from compact object mergers, were found to belong to Group 1, while those events with associated supernovae fell into Group 2. Differences in cumulative counts between the two groups in the first second, and in the minimum variability timescale, identifiable only with the 4 ms resolution data, may account for this result. Short GRBs have particular significance for multi-messenger science as a distinctive EM signature of a binary merger, which may be discovered by its gravitational wave emissions. Incorporating the T1 interval into classification algorithms may support the rapid classification of GRBs, allowing for an improved prioritisation of targets for follow-up observations.