The spatial geostatistical characteristics of the 5 heavy metals (Cr, Cu, Ni, Zn and  Cd) of the top 106 soil samples on Yancheng coast were studied with the aid of GS + software. With ARCGIS 9.2, the Hakanson Risks Index (HRI). A computational approach, including NDVI, determining Euclidean Distance to the main anthropogenic impact sources, and land use reclassification, was proposed to detect the ecological hotspots. The conclusions can be drawn. Metals (Cr, Cu and Cd) had strong spatial autocorrelations, while Ni had modest spatial autocorrelations, and the spatial autocorrelation of Zn could not be modeled.