Parameter estimation of multi-dimensional hidden Markov models - a scalable approach