Restoration of blurred images is a topic of continuous research today. Satellite images are affected by atmospheric turbulence and noise. Blurred images are represented by the matrix equation Ax=b, where A is the blur matrix, x is the restored image and b is the distorted blurred image. Our research objective is to find blur operators that can be represented by hierarchical structured matrices or HSM. Once blur operators have been found, they will be rewritten in hierarchical structured matrices. Iterative techniques are sometimes better to use, due to the high cost with direct inverse methods. These structured matrices are used as an approximation and basis for preconditioned iterative methods.