In a direction-finding process, high-resolution subspace-based algorithms are the most popular ones. It is well-known that their performance of direction of arrival estimation mainly depends on the accuracy of the signal subspace. However, the traditional methods of capturing the signal subspace do not mine the information hidden in the array output in depth, which may restrict their application to some extent. In this study, we elaborate on a novel scheme to extract the signal subspace through refinement of the correlation matrix of the array output. In the developed scheme, a collection of spatial temporal correlation matrices is firstly established. Then, we define a weighting vector for the correlation matrices, and take the weighted average of the correlation matrices as the covariance matrix of the array output. It is clear that this covariance matrix is more general than the traditional covariance matrix, and the signal subspace can be optimized through adjustment of the weighting vector. In this study, we present an optimal weighting vector by adopting the particle swarm optimization. Simulation results demonstrate that the proposed approach has better performance in root mean square error compared to the existing schemes.