Soil salinity is one of the most severe environmental problems worldwide. It is necessary to develop a soil-salinity-estimation model to project the spatial distribution of soil salinity. The aims of this study were to use remote sensed images and digital elevation model(DEM) to develop quantitative models for estimating soil salinity and to investigate the influence of vegetation on soil salinity estimation. Digital bands of Landsat Thematic Mapper(TM) images, vegetation indices, and terrain indices were selected as predictive variables for the estimation. The generalized additive model(GAM) was used to analyze the quantitative relationship between soil salt content, spectral properties, and terrain indices. Akaike’s information criterion(AIC) was used to select relevant predictive variables for fitted GAMs. A correlation analysis and root mean square error between predicted and observed soil salt contents were used to validate the fitted GAMs. A high ratio of explained deviance suggests that an integrated approach using spectral and terrain indices with GAM was practical and efficient for estimating soil salinity. The performance of the fitted GAMs varied with changes in vegetation cover.Salinity in sparsely vegetated areas was estimated better than in densely vegetated areas. Visible red and near-infrared bands, and the second and third components of the tasseled cap transformation were the most important spectral variables for the estimation. Variable combinations in the fitted GAMs and their contribution varied with changes in vegetation cover. The contribution of terrain indices was smaller than that of spectral indices, possibly due to the low spatial resolution of DEM. This research may provide some beneficial references for regional soil salinity estimation. Soil is one of the most severe environmental problems worldwide. It is necessary to develop a soil-salinity-estimation model to project the spatial distribution of soil salinity. The aims of this study were to use remote sensed images and digital elevation model (DEM ) to develop quantitative models for estimating soil salinity and to investigate the influence of vegetation on soil salinity estimation. Digital bands of Landsat Thematic Mapper (TM) images, vegetation indices, and terrain indices were selected as predictive variables for the estimation. The generalized additive model (GAM) was used to analyze the quantitative relationship between soil salt content, spectral properties, and terrain indices. Akaike’s information criterion (AIC) was used to select relevant predictive variables for fitted GAMs. A correlation analysis and root mean square error between predicted and observed soil salt contents were used to validate the fitted GAMs. A high ratio of explained deviance sugge sts that an integrated approach using spectral and terrain indices with GAM was practical and efficient for estimating soil salinity. The performance of the fitted GAMs varied with changes in vegetation cover. Salinity in sparsely vegetated areas was estimated better than in densely vegetated areas. Visible red and near-infrared bands, and the second and third components of the tasseled cap transformation were the most important spectral variables for the estimation. Variable combinations in the fitted GAMs and their contribution varied with changes in vegetation cover. The contribution of terrain indices was smaller than that of spectral indices, possibly due to the low spatial resolution of DEM. This research may provide some positive references for regional soil salinity estimation.