3D Intersection operations for voxel data represented as surfaces in GIS