Many popular direction-of-arrival (DOA) estimators rely on the fact that the array response vector of the array is Vandermonde, for example, that of a uniform linear array (ULA). Array interpolation is a preprocessing technique to transform the array response vector of a planar array of arbitrary geometry to that of a ULA over an angular sector. While good approximation within the target sector is attained in the existing array interpolation approaches, the response of the interpolated array in the out-of-sector region is at best partially controlled. Accordingly, out-of-sector signals, especially those highly correlated with the in-sector signals, can degrade significantly the performance of DOA estimators (e.g., MUSIC with spatial smoothing) that rely on the Vandermonde form to work correctly. In this paper, we propose an improved array interpolation approach that takes into account the array response over the full azimuth. We present also numerical examples to demonstrate the shortcomings of the existing approaches and the effectiveness of our proposal.