The creation of earthquake hazard maps requires various datasets with selected attenuation relations. Based on the selected attenuation relation, the calculation time varies from half an hour to a couple of days. The length of time needed to create an earthquake hazard map also depends on the resolution of the resulting map. The time gets longer as the resolution of the resulting earthquake hazard map gets higher. The basic form of an attenuation relation requires complex calculation algorithms including geospatial information related to the region of interest. Nowadays, next-generation attenuation (NGA) models are introduced to generate more realistic earthquake hazard maps. However, the more complex the attenuation relation is, the longer time will be required to create a hazard map. This paper offers a new method to create high-resolution earthquake hazard maps, faster than using traditional attenuation relation methods, by using an analytic hierarchy process of spatial multi-criteria decision analysis and geographic information systems. This method has been generated and tested for the city of Istanbul. The resulting maps are compared with the earthquake hazard maps created for the city of Istanbul by using the NGA model of Boore and Atkinson (in Boore-Atkinson NGA ground motion relations for the geometric mean horizontal component of peak and spectral ground motion parameters (trans: Engineering Co, University of California B). Pacific Earthquake Engineering Research Center 2007). A second output of this paper is a map of the elements at risk (EaR) for the population and buildings of Istanbul, and the introduction of a new approach of net elements at risk (NEaR).