Soil property maps are essential resources for agricultural land use. However, soil properties mapping is costly and time-consuming, especially in the regions with complicated topographic conditions. This study was conducted in a hilly region of Central Vietnam with the following objectives: (i) to evaluate the best environmental variables to estimate soil organic carbon (SOC), total nitrogen (TN), and soil reaction (pH) with a regression kriging (RK) model, and (ii) to compare the accuracy of the ordinary kriging (OK) and RK methods. SOC, TN, and soil pH data were measured at 155 locations within the research area with a sampling grid of 2 km × 2 km for a soil layer from 0 to 30 cm depth. From these samples, 117 were used for interpolation, and the 38 randomly remaining samples were used for evaluating accuracy. The chosen environmental variables are land use type (LUT), topographic wetness index (TWI), and transformed soil adjusted vegetation index (TSAVI). The results indicate that the LUT variable is more effective than TWI and TSAVI for determining TN and pH when using the RK method, with a variance of 7.00% and 18.40%, respectively. In contrast, a combination of the LUT and TWI variables is the best for SOC mapping with the RK method, with a variance of 14.98%. The OK method seemed more accurate than the RK method for SOC mapping by 3.33% and for TN mapping by 10% but the RK method was found more precise than the OK method for soil pH mapping by 1.81%. Further selection of auxiliary variables and higher sampling density should be considered to improve the accuracy of the RK method.