using aerial photos and field study, then landslide distribution map was produced. Effective factors on landslide occurrence such as: slope, aspect, elevation, lithology, land use, distance to fault, distance to road and distance to drainage network were digitized in GIS environment as shape digital layers and then used for analysis of fuzzy sets theory. Fuzzy analysis was accomplished after determining membership function for each effective parameter in landslide occurred using Matlab 7.1 software. Then the results representing hazard degree for each pixel in case study map was transported in ILWIS software and landslide hazard zonation was produced. The results showed that identified landslides were located in the class very high (29.42%), high (32.19%), moderate (30.36%) and low (8.04%) in susceptibility zones.