Real-Time Interference Artifacts Suppression in Array of ToF Sensors †

Time of Flight (ToF) sensors are the source of various errors, including the multi-camera interference artifact caused by the parallel scanning mode of the sensors. This paper presents the novel Importance Map Based Median filtration algorithm for interference artifacts suppression, as the potential 3D filtration method. The approach is based on the processing of multiple depth frames, using the extraction of the interference region and application of the interpolation. Considering the limitations and good functionalities of proposed algorithm, the combination with some standard methods was suggested. Performance of the algorithm was evaluated on the dataset consisting of the real-world objects with different texture and morphology against popular filtering methods based on neural networks and statistics.


Introduction
Due to the demand in the last decades, there are several types of depth cameras available on the market. They use the different operation principles of depth sensors, that allow reconstructing the geometric properties of the environment. Often in addition to the depth of the scene, also the radiometric information about the scene is obtained by a coupled RGB sensor [1]. The popular principle for determining the depths of scene-points is called Time of Flight (ToF). With the advantage of non-invasive nature and eye-safety of ToF cameras, medical research and imaging is the potential area of use.
However, besides all advantages, ToF cameras are subject to a large variety of measurement error sources. The number of investigations of these errors were reported and have shown that they are caused by sensor parameters and properties or environment configuration [2]. Essential ToF camera errors include multi-path (MPI) and multi-camera interference (MCI).
In our work, we deal with the analysis of MCI occurs due to the reception of a reflected modulated signal generated from another ToF camera. We also propose the filtering method of extreme values generated by the interference of multiple ToF cameras. This method is designed for object scanning applications where the cameras are directed to one common point. In comparison with other hardware-based interference suppression methods, the proposed algorithm is universal and does not require physical interventions, such as the frequency modulation of the IR light source. The disadvantage of a software-based solution is that there is a possible loss of non-corrupted data. Therefore, we aim to maximize the amount of true interference data removed while retaining as much non-corrupted data as possible. The filtering algorithm is based on the combination of temporal filtering and interpolation of areas damaged by interference. It is possible to recover the 3D model of an object from its multi-view images [3]. Our algorithm is able to reconstruct the damaged data by interpolation from one view. The algorithm should filter the data precisely, which means that the useful points should be preserved. The performance of the algorithm is evaluated on the interference-corrupted point cloud model via comparison with other popular algorithms (Statistical Outlier Removal, Radius Outlier Removal, and PointCleanNet).

Related Work
Multi-path interference is one of the basic ToF camera errors. The assumption that the light returned to each pixel of the sensor comes only from a single position in the scene is not practically possible and correct. The reason is that the light can travel multiple paths. The well known is the scene based MPI, where the light scattering of multiple points in the scene arrives in the same pixel of the sensor is caused by inhomogeneities of the scene; or intra-camera multi-path interference caused by the light refraction and reflection of an imaging lens and aperture. This interference is described in detail in References [4][5][6][7].
MCI occurs in multi-camera systems where the cameras operate with the light signal of the same modulation frequency and are directed to a common point. This interference is due to the reception of a reflected modulated signal generated from another ToF camera. The depth map error generated by MCI is much more significant than MPI [4].
One of the still remaining problems of ToF cameras are motion artifacts in dynamic scenes, that may lead to false values in the measured object. Creating the most accurate models of dynamic objects such as the pediatric patient is dependent on the scanning time. With a longer time needed for scanning, the motion artifacts and motion errors are introduced into the resulting model [8]. In general, motion artifacts can be minimized with a shorter image acquisition time. Placing the demands on the fast image acquisition excludes the usage of single-camera scanning for such an application with children. A possible solution is to use a multi-camera system that works in parallel mode and collects the data from multiple viewing angles. Although, if we want to use the multi-camera system, it is necessary to reduce the impact of MCI on the generated 3D models.
One concept of the multi-camera system uses a sequential mode where the cameras work at certain synchronized intervals (in one time only one camera is active). It is the opposite method of parallel cooperation of cameras. However, this method is not suitable for the scanning of dynamic objects. Another option is to use cameras that offer multiple modulation frequencies of emitted light. Many recent works dealing with the problem of multi-camera interference state that ToF-cameras must not be operated using the same modulation frequency as a solution [9]. They rely on orthogonal functions, such as sinusoids of different frequency modulation. Although it does not solve the problem if we are limited by the camera's availability or price. There is also software complication with the connection of different types of cameras in one application [10][11][12].
Using the least-squares estimation of non-interference signal (LSENS), it is possible to restore depth and amplitude information from serious interference. This method is based on the analysis of the statistical properties of signal interference [13]. However, the interfering signal must have a positive and negative fluctuation of the sinusoidal function. In situations where the cameras are directed against each other, the fluctuation of the signal may be extremely high. In this case, the LSENS method will not be effective. The method for MCI error reduction that was proposed by Lee [14] is a Stochastic Exposure Coding (SEC), which includes dividing a camera's integration time into multiple slots and switching the camera off and on stochastically during each slot. However, the experiment was not carried out with dynamic objects.

Time of Flight Camera Measurement Model
ToF cameras are optical sensors that provide information about the depth of the scene. They contain an active light source that generates an amplitude modulated signal. The signal may have a continuous or impulse character. Most of the ToF cameras emit amplitude modulated continuous wave (AMCW) with a frequency near IR for illumination of the scene [15]. Depth measurement is based on the amplitude measurement of the phase shift of the transmitted and received modulated signal ( Figure 1). The depth information for each pixel can be calculated by the synchronous demodulation of the received modulated light in the detector. The demodulation can be performed by interleaving with the original modulated signal. This process is called cross-correlation. In general, the correlation function c(τ) is defined as follows: where s(t) is the received optical signal and the g(t) is transmitted (original) signal. Using specific functions for ToF cameras the equations are: where a is modulation amplitude and ϕ is a phase shift. This function is calculated for four different ωt arguments that are shifted from 0 by 90 • . The received signal is mostly superimposed on a background image, which requires adding the offset b to the correlation function: With these four selected points it is possible to calculate the correlation function and determine the phase ϕ and amplitude a of s(t): The depth d is calculated by the following equation: where c is the speed of light and f is the IR modulation frequency [16].

Multi-Camera Interference Model
In a multi-view ToF camera system each camera uses the same modulation frequency and IR wavelength, so the received signals interact with one another. In the experimental topology of sensors, we use three ToF sensors S1-S3 that are static and placed in the scanning cabin. This system is shown in Figure 2, which describes the interrelationships between the individual cameras and the object. Each camera generates a signal g(t) and receives the reflected signal s(t) which represents the connection of the generated sub-signals from all cameras S. These signals can be described as follows: Each sub-signal can be written as: s xy (t) = b xy + a xy · cos(ωt − ϕ xy ). (11) The x index represents the target camera and index y source camera. Next equations serve as an example of the computational model for one signal s 1 . For rest of cameras the mathematical model is derived the same way: After 4-phase correlation of interfering signals, Equations (7) and (8) transform to the following form:φ 1 = atan a 11 sin(τ 11 ) + a 12 sin(τ 12 ) + a 13 sin(τ 13 ) a 11 cos(τ 11 ) + a 12 cos(τ 12 ) + a 13 sin(τ 13 ) , where:ā = a 11 a 12 [sin(ϕ 11 + ϕ 12 ) + cos(ϕ 11 + ϕ 12 )]+ a 11 a 13 [sin(ϕ 11 + ϕ 13 ) + cos(ϕ 11 + ϕ 13 )]+ a 12 a 13 [sin(ϕ 12 + ϕ 13 ) + cos(ϕ 12 + ϕ 13 )]. (15) This mixing of signals from different cameras causes significant measurement errors and therefore the output depth map contains artifacts [13].

Multi-Camera System Interference Analysis
Equations (13) and (14) are valid when each received signal s x (t) is the combination of the sub-signals of all the cameras S1-S3, where x is going from 1 to 3. The interactions of the cameras that are located according to Figure 2 depends on the location of the object in the scene being scanned. If the object is large enough, it can prevent interference with the s 32 and s 23 interference sub-signals. Also, if the object is too distant, the S1 interference between the other cameras disappears. The shape of the object also affects interference. When scanning dynamic objects, the effect of interference on output images may change, so image filtering must be responsive to changes (adaptive).
ToF cameras provide RGB, IR and depth images. The depth image encodes the distance of the scanned surface from the sensor to the grayscale image. Thus, the depth image contains 3D information in the 2D plane of the image. The value of individual pixels (32-bit float) represents absolute distance. Having this information, it is possible to reproject the scanned scene, where internal and external parameters of RGB and IR cameras must be known.
Interference affects only the IR image (not RGB) and subsequently, interference artifact occurs in the depth image. In Figure 3 it is possible to see IR and depth image with and without interference. In the depth image, interference is manifested as extreme values considering values from the nearest surrounding. This fact causes some deformities in the 3D reconstructed model. Figure 4 shows the difference between 3D models and points clouds reconstructed from depth images with and without interference. Figure 3. Images from Microsoft Kinect v2 sensor: (a) Infrared (IR) frame without interference.

Methodology
There are several overview articles dealing with filtration and restoration techniques used in 3D imaging and scanning. From all of them, Reference [17] seems to be very useful, because it contains a review of approximately 40 methods and categorizes them into several classes. Many of these methods are applicable directly to 3D point clouds.
In our case, we considered the interference as a specific place in the signal, where the points are extremely distant from the real surface. Since outliers appear as isolated points with extreme values, interference areas look like specific continuous clusters, as we can see later in Figure 15. The interference suppression means removing these interference points while retaining the useful points. In this paper, we propose several approaches described in the following section.

Experimental Setup
To verify the functionality of the selected filtration methods, we performed an experiment under real conditions on a static object. Three ToF cameras of type Microsoft Kinect v2 were used in the experiment. These cameras were placed in the scanning cabin, the design of which is shown in Figure 5d. In the center of the cabin, the model object is located in the range from 60 cm to 80 cm from the sensors. The object was represented by a plastic head with a relatively complex shape. The size of the object corresponds to normal dimensions of a child aged 5-10 years. The RGB images received from each camera are shown in Figure 5a, the depth frames of the scanned scene without interference (sequentially working cameras) for the individual cameras are shown in Figure 5b. These depth maps were used as reference data with no interference artifacts. Figure 5c is a reconstructed 3D model of a static object. To objectively evaluate the performance of interference suppression we need to know about every point in the dataset, whether the point is true interference or not. The data distorted by interference were obtained in the parallel scanning mode of cameras. For obtaining ground truth information we used the spatial filtering of interference point clouds against the reference. The reference mesh, shown in Figure 5c, was obtained from the reference point cloud. The overall process of method evaluation is shown in Figure 6. The steps in the process are described as follows. • Buffering-capturing the several point clouds from a particular view and merging to one cloud. • Filtering-applying one of the filtering methods described in Section 5.2.

•
Comparing with reference -computation of Hausdorff distances against reference mesh for a particular view. • Spatial filtering-statistical decision, whether the point is interference or not based on its Hausdorff distance from reference. As interference, we consider the statistical outliers from the normal distribution.

•
Analysis-the evaluation of results using the known metrics from ROC analysis.
In this experimental part, we use the following metrics to evaluate the performance of particular methods.

•
Accuracy reflects the overall algorithm performance, so it also takes into account the true negative predictions (the point is not interference). The accuracy is high when the number of true positive and true negative predictions is large.

•
Recall reflects relatively, how many true interfering points we have removed, without the knowledge about false positives. The recall is high, when the number of true predictions (true interference points filtered) is large, no matter how many useful points we lose.

•
Precision reflects relatively the number of true predictions with respect to the number of false-positive predictions. • F1 score is a balance value between recall and precision. This means that the F1 score is not distorted by a large number of true negative predictions and may be considered as decisive.

•
Hausdorff distance measures how far two point sets are from each other. We use the mean value of distances from a point in the cloud to the closest surface in the reference mesh.

Filtering Methods
Actually there are several approaches in the filtering of 3D data. Many of them are based on statistical properties between the given point and its surroundings and there are some based on artificial intelligence. In the following part, we describe our method called Importance Map Based Median (IMBM) filter and also several other standard methods for further result comparison (Statistical Outlier Removal (SOR), Radius Outlier Removal (ROR), PointCleanNet).

SOR Filter
Statistical Outlier Removal is the outlier filtering method based on the statistical analysis on each point's neighborhood. An outlier can be classified as a misplaced or isolated point or set of points in the point cloud. In the SOR method, we establish the number of points, that will be considered as neighbors (user-defined number of k nearest neighbors [18]) and then the mean distance of each point to its neighbors is calculated. The points whose mean distances do not meet the criteria and are outside an interval which is defined by the global distances mean and standard deviation are considered as outliers and trimmed from the dataset [19]. The SOR is a widely used and efficient method implemented as a part of the PCL library, although there is a time processing limitation in large 3D datasets [20].

ROR Filter
Radius Outlier Removal is a quite simple statistical filter. If values of a spherical radius rad and a minimum number of neighbors k in the point cloud are given, the filter removes the points which have less than k points in a sphere of radius rad centered in the point [21].

POINTCLEANNET
PointCleanNet is a two-stage cleaning deep-learning algorithm that removes outliers and reduces noise in unordered point clouds. The architecture is based on PCPNet (a deep-learning-based approach for estimating local 3D properties in point clouds [22]) with few modifications. In the first step, the outliers are removed. The network uses a novel loss function that removes the outliers without knowledge of the noise characteristics or the information about the surface. In the second step, the correction vectors for the remaining points are estimated. As a deep learning-based method, it takes advantage of high accuracy but the training can be time-consuming and resources-intensive [23]. In the experimental part, we use only the pre-trained model for the removal of the outliers.

Importance Map Based Median Filter
For the suppression and filtering of interference artifacts in the multi-camera system, we proposed an algorithm (see Figure 7). The algorithm is based on interference region extraction and interpolation in a series of images obtained from the camera. The algorithm was proposed mainly for the application of scanning parts of the human body but is not strictly intended only for this application field. The important prerequisite for filtering and interpolation is to have a sufficient number of subsequent depth frames in the buffer. The basic condition is that objects must have the same or very close position in each frame in buffer. This condition is met by fact, that object is static. In the case of the object which can make a move, we have to minimize the buffer size to avoid possible motion artifacts in the resulting depth map. Firstly, we must set the adequate number of frames for buffer. This number is discussed in experimental results. Each frame going into the buffer is thresholded using threshold values determined by the topology and properties of our system. If d(x;y) is a single pixel of the depth map, the threshold is performed using the following equation: where d(x;y) is 32-bit float number and T LOW , T H IGH are lower and upper depth threshold values. In our case they are set to 80 and 1300 (these numbers represent real distance in mm). The scanning cabin is shown in Figure 5d and the dimension of the bottom is 1520 × 1520 mm. After depth thresholding, the background is removed and images include only the object of interest (Figure 8). This simple step is important and effective because it often removes areas with interference in which pixels tend to take extreme values (overflowing 32-bit stack variable). In the next step, the thresholded image with a minimum number of non-zero values is selected as a reference image in the buffer. It is because this image with the highest probability contains a minimal number of interfering pixels. The binary map is created from this reference image and each depth image in the buffer is masked with it. In this step, input images are divided into three-pixel regions-background (also including extreme interference regions), interference region pixels (also including some object pixels) and object pixels (also including interference pixels). Interference regions are dilated by structuring element 1 × 1 square because extreme interference regions are mostly surrounded with lower intensity of interference regions. In this way, the sum of interference pixels is decreased but the number of object pixels in the interference map encreases.
These morphologically processed frames in the buffer are used to compute the median frame and importance map. The median frame pixels are computed as medians of corresponding pixels through all buffer frames, where all zero values are excluded from the median computation. If all the values in number series are zeros, the resulting value is zero. This step is analogous to temporal averaging filtering of image data and also reduces noise processes [24].
Each median value (pixel in median frame) is annotated with information about the number of non-zero values used for median computation. This additive information (2D array) is called importance map (Figure 9a).
Using importance map we have knowledge about the significance of each pixel in the median frame. The small importance value of a given pixel means a high probability of interference occurring at a particular point. For this reason, it seems to be useful to remove such areas by threshold settings. Pixels in median frame (Figure 9b) with selected importance higher than k · N are preserved, other pixels are set to zero. The k represents the importance map threshold and N is the size of the buffer. The coefficient k can take a range of values from 0 to 1.
The next phase in IMBM filter is an interpolation. Data removed from the median frame could be reconstructed by using bilinear interpolation as an example. The effect and influence of different interpolation are discussed in Reference [25].
Bilinear interpolation can be replaced by another interpolation method (e.g., bicubic, biquadratic...). In the aspect of computation, bilinear transformation is a good compromise between the nearest neighbor method and other more complex interpolation methods. The limitation of interpolation occurs if there is a depth map with a relatively big area of holes. The limitation of interpolation is also the holes situated in the place of a mouth, eye corners, and so forth.
The used interpolation method is adapted for scanning the application of the human head and face. We must avoid and restrict the false reconstruction in the mentioned areas (mouth, eyes...). Therefore we do not perform the interpolation over the data where the depth significantly varies, over the areas between small isolated particles or large areas affected by interference. Interpolation can be disabled in IMBM algorithm or can be replaced by other interpolation methods.
Using the median frame we can subsequently repair input depth maps with interference. In the repairing process, the pixels of depth maps are compared with pixels of the median frame. If the absolute difference value is greater than the selected difference threshold (this value is discussed in the experimental part in the Section 6.4), the pixel in original depth map is replaced by median value (or zero value, depending on the mode of IMBM method). The value of the difference threshold depends on camera noise intensity for the defined distance. If this value is too low, the big number of non-interfering pixels will be replaced by the median value. On the other side, the filter could preserve a large amount of interference in the original image. The filter output can be in the form of the median frame or in the form of non-interfering input images.
Generally, we can say, that applying the proposed algorithm improves depth image (Figure 9c) by reducing flying pixels, noise processes and filling the holes in the resulting 3D model.

Experimental Results
For a suitable comparison, we need to set up every filtering method for the best performance. Therefore, in the first step, we tested multiple combinations of filter parameters and get the best options for particular filters. As a standard and decisive metric, we consider the F1 score. In this case, the high accuracy may be confusing because of a large number of true negative predictions (number of not filtered points), that distort the results. Also, the low Hausdorff distance in many cases does not mean high filter performance. This means, that using Hausdorff distance, we are not able to see the point, where the filter removes points that do not suffer from interference. The example of the difference in evaluation using mean Hausdorff distance and F1 Score is shown in Figure 10. Measurements of processing-time consumption were performed on Ubuntu OS with processor Intel(R) Core(TM) i5-6440HQ CPU @ 2.60GHz and 16 GB (Zilina, Slovakia) operating memory.
(a) (b) (c) Figure 10. The difference in evaluation using Hausdorff distance and F1 Score. (a) Point cloud with interference points. The mean Hausdorf distance is large. (b) Precisely filtered point cloud. The mean Hausdorff distance is low and F1 score is large. (c) In this case, large number of useful points were lost. The mean Hausdorff distance is still low but also F1 score is low.

SOR Filtration
The SOR filter allows us to adjust the sdm-parameter (the standard deviation multiplier threshold) and k-parameter (the number of points to use for mean distance estimation). Following color maps ( Figure 11) shows, how the performance of filtration evaluated as the F1 score is changed depending on parameter selection and buffer size.
Based on the analysis in Figure 11 we obtained the best setup for the SOR filter. The performance seems to have an increasing trend up to the buffer size 5. The ideal combination of parameters for buffer size 5 is sdm 1.5 and k 100. At this setup, the algorithm achieves the best F1 score of 0.68, as seen in Table 1.

ROR Filtration
The ROR filter allows us to adjust the radius (the sphere radius used for determining the k-nearest neighbors) and minimal pts-parameter (the minimum number of neighbors that a point needs to have in the given search radius in order to be considered an inlier). Following color maps shows, how the performance of filtration evaluated as the F1 score is changed depending on parameter selection and buffer size.
Similar to the previous case, the analysis in Figure 12 shows the best setup for the ROR filter. The performance of filtration is highest for buffer size 5. The ideal combination of parameters is rad 0.004 and pts 10. At this setup, the filter achieves the best F1 score of 0.63, as seen in Table 2.

PointCleanNet Filtration
By using the pre-trained PointCleanNet model for outliers removal, there are no parameters to set up. Therefore, the filtration performance, in this case, is evaluated as depends only on buffer size. The computation time of the filtration is in the order of minutes, so it is irrelevant for the real-time application.
The performance has an increasing trend with buffer size. As seen in Table 3, the best result obtained based on the F1 score is only 0.16 for buffer size 6.

IMBM Filtration
The IMBM filter allows us to adjust the depth threshold (removing the background and extreme interference pixels), importance map threshold (removing regions with low significance) and difference threshold (setting the depth camera noise threshold). Optimal setting of difference threshold for buffer size 2 is shown in Figure 13a and for buffer size 6 in Figure 13b. Based on the analysis of these figures, the best difference threshold for buffer size 2 is value 13 and for buffer size 6 is value 8. This parameter was statically set to 10 in the next section of the experiment. The proposed filtering method achieves the highest F1 score for buffer size 5. In this point, the improvement in mean Hausdorff distance is 12% (mean distance for the input noised point cloud is 0.001255 m).

Final Comparison
After getting the best setups for particular filtering methods we show the comprehensive comparison of their performance. Figure 14 shows the visualization of Hausdorff distances of the outputs from particular filters for buffer size 5. The buffer size 5 for this comparison was selected based on the best results from Table 4 (except PointCleanNet, which achieves the best results using buffer size 6).  The interfered areas contained most complex details (eyes, chin, nose,...) and in output from SOR and ROR filters, they are lost (holes in models). The output from our IMBM filter recovers corrupted points in the interference area, as seen in Figure 14.
The dataset of several real-world objects was used for the final evaluation of filtering methods. The dataset was collected by authors and contains 23 scans of objects (the ball, the plastic bottle...) captured from various views. The following Table 5 shows the medians of each evaluated metric in the whole dataset. All the results are for buffer size 5. In many cases, none of the filtering methods is capable of sufficient removal of the outliers as good as interference artifacts. As seen in Figure 15, there remains an amount of interfering points after SOR filtering. On the other hand, the IMBM filtered cloud contains outliers, which reduce the overall F1 score. Therefore, the best results according to the F1 score were achieved by using both filtering methods applied subsequently, as seen in Table 5.

Discussion
The performance of the proposed algorithm was evaluated in comparison to known methods for point cloud filtering based on how precise they remove the interference artifacts. However, the SOR, ROR and the PointCleanNet filters are developed mainly for denoising and outliers removal [20,21,23,26]. In many cases, the interference artifact has specific statistical attributes and we cannot consider them as outliers. In this paper, we evaluated their performance on the task, for which they are not primarily designed. Nevertheless, based on Hausdorff distance and F1 score, the SOR and ROR filters achieve comparable results with our proposed method for the filtration of interference artifacts. As figured above, our filtering method does not remove the outlying points sufficiently. Therefore, we propose the usage of our algorithm together with some of the outliers filtering methods. Based on results in this research we recommend using filtering methods for specific applications by following their functionalities described in Table 6. The problematic parts for SOR and ROR filters are corners and edges, that are rounded on the output. As described in Reference [18], the SOR filter perceives points at borders as outlying points. These filters also remove the small areas (clusters of points that are not interference points) [27,28]. As described in Reference [20], SOR filtering is not entirely appropriate for real-time performance.
The IMBM filter can be divided into the filtration part and the interpolation part (which repairs corrupted data). When focusing only on the filtration phase (common for all discussed methods) our method takes the shortest time. With approx. 200 ms for 6-frame buffer computation, we are theoretically able to process 5 full buffers while framing ratio of the sensor can be set to 20-30 fps. In the case of static objects or quasi-static objects we assume this framing frequency as sufficient.
Considering the limitations and also good functionalities of individual methods, the suitable solution is the combination of the IMBM method and SOR or ROR. When known filtration methods (SOR, ROR, PCN) removes outliers, it can partially also suppress the interference. On the other hand, our IMBM method removes interference but also it can partially remove outliers. These combinations caused increase in the time consumption and rounding of the corners. However, if the removal of interference and also denoising is requested we accept the cost of limitations mentioned above. The IMBM algorithm we propose does not only aim to remove the interference but also provides the interpolation of missing points. According to this fact, the processing time is increased. The probability that the filter will detect and remove the interference rises with the higher variance in the positions of the interference in between the buffer frames.

Conclusions
As mentioned in the introduction, the main objective of our research is the analysis and the suppression of multi-camera interference artifacts. The scanning system demands several ToF cameras in its spatial topology, which results in the multi-camera interference error.
In this paper, we describe the theoretical analysis of interference artifacts and propose the novel algorithm for multi-camera interference suppression, as the innovation and extending to standard methods. There are several known outlier removal and denoising methods which are not primarily designed for interference suppression. We conclude that interference and outliers differ in character and we cannot consider them as the same type of data corruption. Therefore, we proposed IMBM algorithm, which achieves the best performance in combination with some of the outliers filtering methods. The proposed combination leads to suppression of noisy artifacts as well as interference, what extends the functionality of known methods. Unlike the other filtering methods, IMBM also provides the interpolation of removed points. Nevertheless, it is still suitable for real-time applications. The specific medical applications focusing on the 3D scanning of head and face are potential areas of usage for the proposed algorithm. In the future work, the research may potentially focus on optimization of the script for the multi-thread processing.