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 BiFeO 3 films with 71°and 109°domain wall, and, tetragonal BiFeO 3 , Pb(Zr 0.2 Ti 0.8 )O 3 , and BaTiO 3 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.