We present a new method of interpolation for the pixel brightness estimation in astronomical images. Our new method is simple and easily implementable. We show the comparison of this method with the widely used linear interpolation and other interpolation algorithms using 1000 astronomical images obtained from the Sloan Digital Sky Survey. The comparison shows that our method improves bad pixel brightness estimation with four times lower mean error than the presently most popular linear interpolation, and performs better than any other examined method. The presented idea is flexible and can be also applied to present and future interpolation methods. The proposed method is especially useful for large sky surveys image reduction but can be also applied to single image correction.