Following a geometric paradigm, the estimation of a Toeplitz structured covariance matrix is considered. The estimator minimizes the distance from the Sample Covariance Matrix (SCM) while complying with some specific constraints modeling the covariance structure. The resulting constrained optimization problem is solved globally resorting to the Dykstra' projection framework. Each step of the procedure involves the solution of two convex sub-problems, whose minimizers are available in closed form. Simulation results related to typical radar environments highlight the effectiveness of the devised method.