Abstract: In this paper, we propose a novel channel estimation (CE) algorithm for orthogonal frequency division multiplexing (OFDM) systems that exploits the unique characteristics of the delay-Doppler (DD) domain channel. Specifically, the time-frequency (TF) domain input-output relationship (IOR) is derived in a compact form by focusing solely on the non-zero elements of the TF domain channel matrix. Based on this compact IOR, a coarse TF domain CE is first performed using a linear minimum mean square error estimator. Then, the resultant TF domain estimates are transformed to the DD domain through a unitary transformation for further refinement. We reveal that the effective DD domain channel matrix can be viewed as an aggregation of multiple DD domain channel responses with different phase shifts. This allows us to devise a threshold-based estimation for DD domain channel parameters with high accuracy. The estimated DD domain channel parameters are then applied to form a refined estimate of TF domain channel. Our numerical results demonstrate that the proposed method can achieve substantial performance gains over conventional OFDM channel estimation techniques under the same pilot deployment.