We extend thresholding methods for numerical realization of mean curvature flow on obstacles to the anisotropic setting where interfacial energy depends on the orientation of the interface. This type of schemes treats the interface implicitly, which supports natural implementation of topology changes, such as merging and splitting, and makes the approach attractive for applications in material science. The main tool in the new scheme are convolution kernels developed in previous studies that approximate...