In the Hamersley Ranges of Western Australia, highly consistent marker shales are used to identify the stratigraphic location of the banded iron formation (BIF) or BIF-hosted ore. These marker shales produce distinctive signatures in the natural gamma downhole logs, which are currently processed by slow manual interpretation. The location of the marker shales is important when distinguishing lithology for geological and ore body mapping. Previously, two methods of automatically identifying these shales in the natural gamma logs have been proposed, using either Gaussian processes (GP) or dynamic time warping (DTW). There are advantages and disadvantages with each method and a combined result is desirable. However, the outputs from each technique are very different and cannot be directly combined. This paper proposes a method of standardizing the results to the same scale and combining them using a weighted average. The method is tested using two natural gamma signatures from a deposit that contains the Marra Mamba sequence of the Hamersley province. The results showed that the combined method outperforms both individual methods. The GP has accuracies of 87.0 to 92.0 %, increasing to 92.9 to 95.7 % when the method is certain. The DTW has accuracies of 90.0 to 90.6 %, increasing to 94.0 to 95.4 % when the method is certain. In comparison, the combined method has accuracies of 90.7 to 93.6 %, increasing to 94.4 to 95.9 % when this method is certain. The combined method balances the strengths of the individual methods, identifying a larger number of positive signatures than the GP, while retaining a higher classification accuracy for these signatures than the DTW.