This thesis focusses on the intersection 3D Boolean set operation on voxel data represented as surfaces, with other 3D surface based vector data, in GIS. The Application of this research makes the linkage with geohydrology possible: based on the vectorised GeoTOP (TNO) data the water infiltration and storage capacity of specific locations can be determined.