: Based on regional geochemical prospecting using heavy mineral stream sampling in Hamadan geological rectangle area (scale1:100000), a number of promising prospects for both metallic and non metallic ores were introduced. Following this stage and through applying more detailed litho-geochemical and geophysical studies on the proposed area, four anomalous zones were suggested for further exploration in more details. The current study were defined in the same line of previous research to delineate more favorable areas in two better suited anomalies in a GIS environment. After applying geometric, radiometric corrections the data were smoothed and heavily vegetated area were removed for further preprocessing and enhancements routines. . In order to map geochmically anomalous zones a number of preprocessing methods such as sorting, validating, estimating censored data and statistical analysis and visualizations through histogram charting were applied on rock and soil geochemical sample values. Combined multi-element maps were produced using multivariate partial correlation, clustering and factor analysis carried on elemental enrichment indexes. Three different type of mineralization were determined based on type and intensity of multi elemental distributions called gold, Cu-Ni- and Sb-Pb. The geophysical data including DC resistivity and magnetometric were modeled using 3D algorithms and the surficial physical properties were extracted from them through removing topographic data and imported in GIS environment as evidential layers. Using mineral potential mapping methods such as index overlaying, weight of evidence and logistic regressio techniques three different favorability maps for Au, Cu-Ni- and Sb-Pb bearing zones were produced. The inputted data layers for GIS integration were composed of lithologig maps (Diroite, Granite, Gabbro and Aplite), alteration maps (OH , Fe and Alunite bearing formations), multiplicative geochemical elemts maps from rock samples (Au,S,Ag,Sb,Cd,Zn,Pb, Cu.Ni.Co) and soil samples (Sb.Pb.Cd.Zn,,Cu.Ni.Cr.Co), geophysical evidential layers (true surface susceptibility and resitivity distributions) plus faulting related lineaments. In index overlaying method the weights were determined using our knowledge on the type of different mineralization. In assessing the final weight of evidence results the independency of different layers were evaluated via conditional independency tests. For logistic regression method the final coefficients for each evidential layer were calculated and used in preparing the final favorability map. Overall evaluation of the maps produced from each method were through cross checking maps showed acceptable agreement and consistency.