Restoring high value medicinal and industrial plant species is quite important in terms of economical and social point of view. Nowadays, progresses in statistical methods and geographic information techniques facilitate habitat potential mapping using more advanced modeling approaches. Species distribution model can be defined as a statistical - analytical algorithm which predicts species distribution from field observations and available environmental maps. This study aimed to generate potential distribution map of Astragalus verus Olivier, a valuable medicinal plants in Fereydunshahr rangelands in Isfahan province using Logistic regression (LR) and Ecological-niche factor analysis (ENFA) methods. Soil maps were obtained from 70 soil profiles and climatic maps were produced from 10 climatic stations data. Geo-statistical methods were used to convert the obtained data to digital maps with 90 by 90 meter resolutions. Potential distribution map of the species in the study area was produced using presence-absence data of 100 sites using a stratified- random sampling approach. For each sampling site, presence and absence of the species as well as the environmental variables were recorded. The relationships of the plant species distribution and environmental variables were studied using LR and ENFA methods. In the next step, habitat potential maps of the species were produced using these two methods. ENFA as a profile model is one of the most frequently used techniques to generate species distribution map which uses species “presence-only” data. This method integrates geographic information system and factor analysis to produce habitat distribution map from the species occurrence and environmental data. According to the results, the environmental factors influencing determination of the ecological niche of the species in order of importance were soil electrical conductivity, Caco 3 content, pH, slope degree and monthly mean temperature. The produced species distribution model had high accuracy indicated by calculated Kappa coefficient and ROC area under curve plots being 0.72 and 0.85 respectively. According to ENFA results, 10000 hectare or 10% of the study area was recognized as suitable area for the species establishment. LR from discrimination group models is another modeling approach that enables to produce distribution maps using presence and absence data. According to LR results, annual mean precipitation, soil clay content, mean temperature of warmest season and slope degrees were the most important environmental factors influencing the distribution of the species. According to this model analysis, 21000 hectare or 21 percent of the whole study area had excellent habitat condition for the species. The produced species distribution model had acceptable accuracy as calculated Kappa coefficient and ROC area under curve plots were 0.75 and 0.5 respectively. Analysis of sampling sizes on the accuracy of the studied modeling approaches revealed that ROC area under curve plots had less variability to the Kappa coefficient and is more appropriate method to analyze the models outputs accuracy. According to produced potential maps of the species and response curve results, appropriate Astragalus verus habitats characterized by moderate to relatively heavy soil texture (44 to 48 percent clay content), moderate slope (10 to 20 percent), low electrical conductivity ( less than 0.3 dS/m), saturated water content of 45 to 47 % and pH equivalent to 7.5. The studied modeling approaches can be used broadly to allocate appropriate areas for rehabilitation and protection of valuable native plant species such as Astragalus verus . Key words: Species distribution model, Astragalus verus Oliver, Ecological-niche factor analysis, Logistic regression, Geo-statistic, Geographic information system, Fereydounshahr, Isfahan Province.