Discerning multiple order parameters with interpretable machines: towards the automation of phase classification