In this paper, a curve-fitting channel estimation method is proposed for orthogonal frequency division multiplexing (OFDM) system in a time-varying frequency-selective fading channel. The method can greatly improve channel state information (CSI) estimation accuracy by performing smoothing and interpolation through consecutive curve-fitting processes in both time domain and frequency domain. It first evaluates least-squares (LS) estimates using pilot symbols and then the estimates are approximated to a polynomial with proper degree in the LS error sense, starting from one preferred domain in which pilots we densely distributed. Smoothing, interpolation, and prediction are performed subsequently to obtain CSI estimates for data transmission. The channel estimation processes are completed by smoothing and interpolating CSI estimates in the other domain once again using the channel estimates obtained in one domain. The performance of proposed method is influenced heavily on the time variation and frequency selectivity of channel and pilot arrangement. Hence, a proper degree of polynomial and an optimum approximation interval according to various system and channel conditions are required for curve-fitting. From extensive simulation results in various channel environments, we see that the proposed method performs better than the conventional methods including the optimal Wiener filtering method, in terms of the mean square error (MSE) and bit error rate (BER).