In this paper, we propose the improved initial image estimation method for a fast fractal image decoding. When the correlation between a domain and a range is given as the linear equation, the value of initial image estimation using the conventional method is the intersection between its linear equation and y=x. If the gradient of linear equation is large, that the difference of the value between each adjacent pixels is large, the conventional method has disadvantage which has the impossibility of exact estimation. The method of the proposed initial image estimation performs well by two steps. he first step can improve the disadvantage of the conventional method. The second step upgrades the range value which was found previous step by referring information of its domain. Though the computational complexity for the initial image estimation increses slightly, the total computational complexity decreases by 30% than that of the conventional method because of diminishing in the number of iterations.