Abstarct: At the preliminary prospecting and exploration stages, the areas that must be paid attention to in more detailed studies are identified by mineral potential mapping (MPM). The final goal of this study is prepare the MPM of iron ore mineralization in Esfordi 1:100000 sheet. The Esfordi 1:100,000 Sheet is located in theBafgh-Posht_e_Badam district of the central Iran structural zone. Iron oxide-apatite, lead and zinc mineralizations have been observed in different areas of this region. The simultaneous utilization of the Prediction-Area plot and DEMATEL, ANP, VIKOR, and TOPSIS methods in fuzzy space has resulted in reliable results. The P-A plot and fuzzy ANP methods were used to weighting and integrated the exploratory laxyers. In the P-A plot method, the weight of each laxyer was calculated using simultaneous analyses of the value-area fractal curve and the use of 22 Fe-bearing occurrences in the study region. The host rock laxyer has the most weight with the intersection point of 78%. Causal and effect relationships of the criteria were assessed by fuzzy DEMATEL method. Eight exploratory laxyers were extracted from four various types of data and utilized; thus, anomaly laxyers of Fe2O3 and TiO2 were prepared from the stream sediments data, magnetic body laxyer from airborne geophysics, host rock laxyer from geology map, the gossan zone, iron oxide, phyllic alteration and lineament laxyers were prepared from satellite images. The use of three types of satellite imagery, including Sentinel-2, Landsat-8 and ASTER, caused that the results are reliable and powerful. Meanwhile, host rock laxyer was identified as the most effective criteria. The results from DEMATEL method were integrated with fuzzy ANP method and mineral potential model was prepared. Several areas in the resulting map were introduced as zones with high potential of mineralization that five areas were selected then the fuzzy VIKOR and the fuzzy TOPSIS methods were employed to prioritize them. In field surveys, regions with high potential of iron ore mineralization were examined. The sample analysis results of each region largely confirm the findings of the used methods. In the meantime, some parts of the southeast of Esfordi area have been introduced as new zones with high potential of iron ore mineralization. The results indicate that rhyolite, rhyodacite and calcareous dolomites are mostly units associated with mineralization.