We report new photometric observations of the transiting exoplanetary system HAT-P-20, obtained using CCD cameras at Yunnan Observatories and Ho Koon Nature Education cum Astronomical Centre, China, from 2010 to 2013, and Observatori Ca l'Ou, Sant Marti Sesgueioles, Spain, from 2013 to 2015. The observed data are corrected for systematic errors according to the coarse de-correlation and SYSREM algorithms, so as to enhance the signal of the transit events. In order to consistently model the star spots and transits of this exoplanetary system, we develop a highly efficient tool STMT based on the analytic models of Mandel & Agol and Montalto et al. The physical parameters of HAT-P-20 are refined by homogeneously analyzing our new data, the radial velocity data, and the earlier photometric data in the literature with the Markov chain Monte Carlo technique. New radii and masses of both host star and planet are larger than those in the discovery paper due to the discrepancy of the radius among K-dwarfs between predicted values by standard stellar models and empirical calibration from observations. Through the analysis of all available mid-transit times calculated with the normal model and spotted model, we conclude that the periodic transit timing variations in these transit events revealed by employing the normal model are probably induced by spot crossing events. From the analysis of the distribution of occulted spots by HAT-P-20b, we constrain the misaligned architecture between the planetary orbit and the spin of the host star.