Multivariable canonical forms for model reduction of 2-D discrete time systems