We present the combined NuSTATR and XMM-Newton 0.6-79 keV spectral analysis of a Seyfert 2 galaxy, NGC 1358, which we selected as a candidate Compton thick (CT-) active galactic nucleus (AGN) on the basis of previous Swift/BAT and Chandra studies. According to our analysis, NGC 1358 is confirmed to be a CT-AGN using physical motivated models, at >3 $\sigma$ confidence level. Our best-fit shows that the column density along the 'line-of-sight' of the obscuring material surrounding the accreting super-massive black hole is N$\rm _H$ = [1.96--2.80] $\times$ 10$^{24}$ cm$^{-2}$. The high-quality data from NuSTAR gives the best constraints on the spectral shape above $\sim$10 keV to date on NGC 1358. Moreover, by combining NuSTAR and XMM-Newton data, we find that the obscuring torus has a low covering factor ($f_c$ <0.17), and the obscuring material is distributed in clumps, rather than uniformly. We also derive an estimate of NGC 1358's Eddington ratio, finding it to be $\lambda_{\rm Edd}$ $\sim$$4.7_{-0.3}^{+0.3}$ $\times$ 10$^{-2}$, which is in acceptable agreement with previous measurements. Finally, we find no evidence of short-term variability, over a $\sim$100 ks time-span, in terms of both 'line-of-sight' column density and flux.