A standard use of triangulation in GIS is to model terrain surface using TIN. In many simulation models of physical phenomena, triangulation is often used to depict the entire spatial domain, which may include buildings, landmarks and other surface objects in addition to the terrain surface. Creating a seamless surface of complex building structures together with the terrain is challenging and existing approaches are laborious, time-consuming and error-prone. We propose an efficient and robust procedure using computational geometry techniques to derive triangulated building surfaces from 2D polygon data with a height attribute. We also propose a new method to merge the resultant building surfaces with the triangulated terrain surface to produce a seamless surface for the entire study area. Using Oklahoma City data, we demonstrate the proposed method. The resultant surface is used as the input data for a simulated transport and dispersion event in Oklahoma City. The proposed method can produce the seamless surface data to be used for various types of physical models in a fraction of the time required by previous methods
