This paper proposes a realization method of multidimensional digital filters that has high modularity, regularity and parallelism enjoying the attributes for efficient VLSI implementation. The method shows that multidimensional transfer functions can be treated as two-dimensional transfer functions modifying the decomposition method of multidimensional transfer functions proposed by Venetsanopoulos etal, and then be displayed by multiplications and additions of one-dimensional transfer functions by applying the griangular decomposition theorem to the coefficient matrices of the two-dimensional transfer functions.