Abstract: Accurate monitoring of farmland environment and soil pollution is of great significance for food security and ecologically sustainable development in recent years. Taking Xiong'an New Area as the research area, this study aims to inverse the heavy metal content in farmland. Spectral response data of heavy metal content was extracted from multi-source remote sensing data Sentinel-2 (L2A), Zhuhai-1 (OHS), and field measurement. Pollution levels of heavy metal were evaluated in farmland soil using the one-factor and Nemero indexes. The heavy metal elements were screened to cause soil pollution in farmland, such as copper (Cu), zinc (Zn) and plumbago (Pb). Partial least squares regression (PLSR) was integrated with principal component and multiple linear regression. The quadratic inversion model was then established for the contents of the three heavy metal elements. The quantitative inversion of the heavy metal contents was realized in a wide range of soil. Sentinel-2 was used to extract seven vegetation indices, while Zhuhai-1 was to extract the original spectral reflectance of the sample points, as well as four transformed spectral reflectance of sample points. Pearson's correlation coefficient was used to analyze the correlation between spectral reflectance and vegetation index with the content of three heavy metal elements. A comparison was also made on the correlation coefficients, sensitive bands and vegetation indices. The spectral response indicators were modeled as the independent variables, while the measured parameters of the soil heavy metal contents were the dependent variables. The models were finally evaluated using the coefficient of determination (R²), mean absolute error (MAE), and root mean squared error (RMSE). The results showed that excellent performance was achieved in the overall inversion accuracy of the three models. The R², RMSE, and MAE values of the Pb content inversion model were 0.490, 4.66 mg/kg, and 1.92 mg/kg, respectively. The inversion model was obtained for Cu with the R² of 0.491, the RMSE of 16.85 mg/kg, and the MAE of 3.69 mg/kg. While, the inversion model of Zn was achieved with the R² of 0.664, the RMSE of 20.63 mg/kg, and the MAE of 9.36 mg/kg. The spatial distribution of soil heavy metal content in farmland was mapped, according to the constructed inversion model of soil heavy metal content. The visualization of the inversion was conducive to the distribution of regional soil heavy metals from a more intuitive perspective. Among them, the Pb element content in soil in most areas was within the risk screening standard; Cu element exceeded the risk screening value of soil pollution in the southwestern and western parts of the test area; The western and southwestern parts of the test area were more seriously contaminated by Zn element, indicating the sporadic distributions in some farmlands in the southeastern part. The Cu and Zn contents in other areas were within the national risk control value of soil pollution. Therefore, the multi-source remote sensing data and spectral response can be expected to synergistically invert the Cu, Zn and Pb content of soil heavy metals. High feasibility and accuracy can also be obtained to investigate soil heavy metal pollution. At the same time, the findings can provide a theoretical basis for monitoring the heavy metal content of soil in large areas. A new idea can also be offered to construct the inverse model of heavy metal content in farmland soil.