To exploit the advantageous properties of isogeometric analysis (IGA) in a scan-based setting, it is important to extract smooth geometric domain from scan data (e.g., voxel data). IGA-suitable domains can be constructed by convoluting grayscale using B-splines. A negative side-effect this convolution technique is, however, that induce topological changes process smoothing when features with si...