An estimation method for correction filters for wavefront aberration correction in ultrasound imaging is given. The aberration is introduced due to spatial variations of the ultrasound propagation velocity in heterogeneous tissue. The correction filters are found as eigenvectors of the cross-correlation matrix of the element signals. In particular, the eigenvectors with a small spatial linear/plane phase component is sought. Detailed methods for calculating the eigenvectors are given.