In Lijiang City, as a typical example, 93 soil samples were collected from the study area, and soil pH； organic matter； and heavy metals arsenic （As）, mercury （Hg）, copper （Cu）, zinc （Zn）, lead （Pb）, cadmium （Cd）, and chromium （Cr） were determined. We explored the sources of heavy metals in the study area by means of Positive Definite Matrix Factorization （PMF） modeling and analyzed the impact of influencing factors by combining seven heavy metals with 13 influencing factors in a GeoDetector. The results showed that the mean values of soil heavy metals ω（As）, ω（Hg）, ω（Cu）, ω（Zn）, ω（Pb）, ω（Cd）, and ω（Cr） in the study area were 17.55, 0.19, 86.75, 164.84, 28.95, 0.39, and 167.87 mg·kg-1, respectively, which were greater than the background values of soils in Yunnan Province （except for As and Pb）. Regarding spatial distribution, the high values of Cu and Cr content were mainly concentrated in Yulong Naxi Autonomous County； the high value areas of As, Hg, Pb, and Cd were mainly concentrated in Ninglang Yi Autonomous County； and the high value of Zn content was mainly concentrated in Huaping County. Correlation analysis and PMF modeling revealed that the main sources of heavy metals As and Hg in the study area were industrial sources, Zn was from transportation pollution sources, Cr and Cu were from natural sources, and Cd and Pb were from agricultural sources. Further, the factor detector of the GeoDetector found that soil pH and organic matter （OC） had strong explanatory power for the content of seven heavy metals, and the interaction detector found that the results following the interaction of different influencing factors were nonlinear enhancement or two-factor enhancement, in which the interaction of OC and pH was the dominant factor for the spatial differentiation of heavy metals. This provides an important scientific basis for the protection of the soil environmental health and sustainable development in Lijiang City.