What is in this document? Mathematical and algorithmic explanation of this process. I will put a naive Python implementation of this algorithm to make it more clear.

Abstract: A Toeplitz-block-Toeplitz (TBT) matrix is block Toeplitz with Toeplitz blocks. TBT systems of equations arise in 2D interpolation, 2D linear prediction and 2D least-squares deconvolution ...