Iterative Tensor Factorization Based On Krylov Subspace-Type Methods With Applications To Image Processing