Overlay operation is a critical and powerful tool in GIS that superimposes spatial and attribute information from various thematic map layers to produce new information. Overlay operations facilitate spatial analysis and modeling processes when being used with other spatial operations (e.g. buffer, dissolve, merge) to solve real-world problems. For both vector and raster data models, the input layers need to be spatially aligned precisely with each other to ensure a correct overlay operation. In general, vector overlay is geometrically and computationally complex. Some most used vector overlay operations include intersection, union, erase, and clip. Raster overlay combines multiple raster layers cell by cell through Boolean, arithmetic, or comparison operators. This article provides an overview of the fundamentals of overlay operations, how they are implemented in vector and raster data, and how suitability analysis is conducted.