The identification of non-Gaussian signatures in cosmic microwave background (CMB) temperature maps is one of the main cosmological challenges today. We propose and investigate altenative methods to analyse CMB maps. Using the technique of constrained randomisation we construct surrogate maps which mimic both the power spectrum and the amplitude distribution of simulated CMB maps containing non-Gaussian signals. Analysing the maps with weighted scaling indices and Minkowski functionals yield in both cases statistically significant identification of the primordial non-Gaussianities. We demonstrate that the method is very robust with respect to noise. We also show that Minkowski functionals are able to account for non-linearities at higher noise level when applied in combination with surrogates than when only applied to noise added CMB maps and phase randomis