Knowledge of the Franck-Condon factors (FCFs) and related quantities is essential to understand and to estimate many important aspects of the astrophysical molecules, such as kinetics of the energy transfer, radiative lifetimes, band intensity, and vibrational temperatures. In this view, we propose a new analytical formula of the Franck-Condon integral for two-dimensional harmonic oscillators taking into account the Duschinsky effect. This method is based on the use of the binomial expansion theorem and the Hermite polynomials. With the formula obtained, the FCF of any transition can be computed independently. In this study, the method for FCF calculations was applied to the NO2 molecule.