Earthquakes have caused considerable economic damage and loss of lives. The ability to determine earthquake-prone areas has utmost importance to crisis management. Given the gravity of the issue, this research proposes to incorporate geographic information system-based fuzzy-analytic network process (FAHP) and fuzzy-analytic hierarchy process (FANP) methods to determine earthquake-prone areas in the southwest of Iran. In addition, this paper is to utilize the feature selection algorithm (FSA) to select the most relevant parameters influencing earthquakes and prepare seismic activity maps. The final maps obtained using the relevant data were compared with maps obtained using the entire data. For this aim, a fuzzy membership function was initially employed to generate a fuzzy map for different layers, after which the analytic network process and analytic hierarchy process methods were utilized for assigning proper weights to each layer and eventually procuring earthquake maps. The input information for the mentioned methods consisted of landscape, slope, elevation, lithology, land use, soil, road, and fault maps. According to the FSA approach, the most relevant parameters include distance to fault (DTF), lithology, slope, and landscape. The results of FANP and FAHP methods utilizing the receiver-operating characteristic curve showed that the FANP method has a higher accuracy for determining earthquake-prone areas. As the areas located in the northern parts of the region are more at risk. In addition, Finally, the results of earthquake maps obtained using solely the relevant data from FSA were compared with the maps obtained using the entire data. The results of FSA showed that the Best first and Greedy-Stepwise models have high accuracy (R = 0.99) and the most important data for determining earthquake-prone areas using these models are DTF, lithology, slope, aspect. Therefore, by saving time and money, it is possible to prepare maps of earthquake-prone areas with high accuracy.