This paper develops and illustrates a new maximum-likelihood based method for the identification of Hammerstein-Wiener model structures. A central aspect is that a very general situation is considered wherein multivariable data, non-invertible Hammerstein and Wiener blocks, and coloured stochastic disturbances both before and after the Wiener non- linearity are all catered for.