Domain orientations and their volume ratios in ferroelectrics are recognized as a compelling topic recently for domain switching dynamics and domain stability in devices application. Here, an optimized second harmonic generation method has been explored for ferroelectric domain characterization. Combing a unique theoretical model with azimuth-polarization-dependent second harmonic generation response, the complex domain components and their distributions can be rigidly determined in ferroelectric thin films. Using the proposed model, the domain structures of rhombohedral BiFeO3 films with 71° and 109° domain wall, and, tetragonal BiFeO3, Pb(Zr0.2Ti0.8)O3, and BaTiO3 ferroelectric thin films are analyzed and the corresponding polarization variants are determined. This work could provide a powerful and all-optical method to track and evaluate the evolution of ferroelectric domains in the ferroelectric-based devices.