Summary: Because of the long propagation distances, Computational Aeroacoustics schemes must propagate the waves at the correct wave speeds and lower the isotropy error as much as possible. The spatial differencing schemes are most frequently analyzed and optimized for one-dimensional test cases. Therefore, in multidimensional problems such optimized schemes may not have isotropic behavior. In this work, optimized finite difference schemes for multidimensional Computational Aeroacoustics are derived which are designed to have improved isotropy compared to existing schemes. The derivation is performed based on both Taylor series expansion and Fourier analysis. Various explicit centered finite difference schemes and the associated boundary stencils have been derived and analyzed. The isotropy corrector factor, a parameter of the schemes, can be determined by minimizing the integrated error between the phase or group velocities on different spatial directions. The order of accuracy of the optimized schemes is the same as that of the classical schemes, the advantage being in reducing the isotropy error. The present schemes are restricted to equally-spaced Cartesian grids, so the generalized curvilinear transformation method and Cartesian grid methods are good candidates. The optimized schemes are tested by solving various multidimensional problems of Aeroacoustics.