Among the natural hazards, landslides are attracting more and more attention due to its increasing effect on economic and human losses. Landslide hazard zonation maps are considered as a prior and important base to identify the vulnerable zones for designing early warning systems and adequate mitigation measures in landslide-prone areas. 

Many techniques have been proposed in the literature for landslide hazard zonation (Hansen, 1984; Van Westen et. al.1997). These can be generally divided into two groups: direct or semi-direct hazard mapping in which the degree of hazard is determined by the mapping expert and, indirect hazard mapping in which either statistical or deterministic models are used to predict landslide-prone areas based on information obtained from the interrelation between landscape factors and the landslide distribution.

With the introduction of GIS, in particular indirect methods have gained enormous popularity due to the capacity of GIS to handle and analyze data with high spatial variability. In the statistical approach mathematical relationships between the observed landslide distribution and their controlling factors are made on the basis of polygons as mapping units. However, attitudes of geological structures are marked usually as linear or point measurements in factor maps. This holds e.g. for the influence of structural attitudes such as strike or dip directions and dip angles with respect to slope directions and slope angles, and influence of weak zones like fold axis, faults, joints. Even though these factors are considered very important for slope stability, they are often neglected or dealt frivolously in landslide hazard zonation analysis.

In this paper, a possible approach is discussed how GIS capabilities can be used efficiently to integrate the effect of structural attitudes for slope stability problems as an essential part of the final landslide hazard zonation mapping.