In the first step, all the geological data (observed, measured and interpreted) of the study area were organized in a database using a geographical information system. The standard tools of a computer-aided design (CAD) environment were used to build geological profiles using the observed data (outcropping contact lines) extracted from the database, as well as, to create the contact surfaces via interpolation using non-uniformal rational B-spline surfaces. The intersection lines between each of these surfaces and the digital elevation model yield the model contact lines (MCLs). Discrepancies between the MCLs and the corresponding outcropping contact lines were analysed to further constrain the surfaces of the model. Starting from these surfaces, a voxel-based solid model was built, for storing and displaying each rock formation.