Due to multipath delay spread and relatively high sampling rate in OFDM systems, the channel estimation is formulated as a sparse recovery problem, where a hybrid compressed sensing algorithm as subspace orthogonal matching pursuit (SOMP) is proposed. SOMP first identifies the channel sparsity and then iteratively refines the sparse recovery result, which essentially combines the advantages of orthogonal matching pursuit (OMP) and subspace pursuit (SP). Since SOMP still belongs to greedy algorithms, its computational complexity is in the same order as OMP. With frequency orthogonal random pilot placement, the technique is also extend to MIMO OFDM systems. Simulation results based on 3GPP spatial channel model (SCM) demonstrate that SOMP performs better than OMP, SP and interpolated least square (LS) in terms of normalized mean square error (NMSE).