Optical flow computation methods are traditionally classified to two categories: local and global. Several previous works have investigated the combination of them by exploiting their complementary effects. Unlike previous works, this paper is motivated by the common weakness of the two schemes and proposes linking local and global computation by subspace regularization. We will show that the subspace regularization can effectively reduce the error caused by ill-conditioned local computation systems, while inhibiting error diffusion suffered by the global regularization. Experimental results on benchmark sequences validate the enhanced performance of the proposed combination scheme.