An image degradation process is considered as equivalent to a linear transformation of original image matrix processed by transfer function and noise, while the image restoration process is equivalent to trying to get the original image using the least squares method. When the transfer function is separable, the problem is transformed into finding the inverse matrix of a Toeplitz matrix. The simulation results verify the validity of the method and analyze its numerical stability.