How to remove outliers in a 3D surface

Question

HAMID 2021-11-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1580479-how-to-remove-outliers-in-a-3d-surface

评论： Chris 2021-11-7

I have a 3D surface (see below) construced using surf command based on excel data:

T = readtable('https://in.mathworks.com/matlabcentral/answers/uploaded_files/791474/excel-surface.xlsx')

x = T.(1) ;

y = T.(2) ;

z = T.(3) ;

xi = linspace(min(x),max(x));

yi = linspace(min(y),max(y)) ;

[X,Y] = meshgrid(xi,yi) ;

Z = griddata(x,y,z,X,Y) ;

surf(X,Y,Z)

shading interp

colorbar

The resulting surface is as follows:

However, I know that the surface should be smooth without these "peaks" that appear randomly from above or below.

Is there a way to remove these outliers and interpolate to construct a smooth surface? (NOTE, I AM NEW TO MATLAB).

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Matt J 2021-11-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1580479-how-to-remove-outliers-in-a-3d-surface#answer_825309

You could try rmoutliers().

5 个评论
显示 3更早的评论隐藏 3更早的评论

Matt J 2021-11-6

I would have to see what you tried in order to be able to guess what's not working.

HAMID 2021-11-7

I tried fill outlier instead of rmoutliers(), but I am obtaining this:

https://www.mathworks.com/matlabcentral/answers/1580479-how-to-remove-outliers-in-a-3d-surface#comment_1819454

请先登录，再进行评论。

Answer 2

Chris 2021-11-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1580479-how-to-remove-outliers-in-a-3d-surface#answer_825329

编辑：Chris 2021-11-6

在 MATLAB Online 中打开

I would recommend trying the "clean outlier data" task in the live editor, after sorting the data:

T = sortrows(T,1); % Sort by the first column ('x')

Open a New Live Script, and select the task.

Select the table and input variable z.

I can get decent results with spline interpolation, a centered moving median, and a low threshold. But I'm not sure how, computationally, you would create a perfectly smooth surface that retains the contours you would want to see, as there are quite a few clumped outliers.

Perhaps repeated sorting by x and y, with outlier cleaning in between, could work.

2 个评论
显示无隐藏无

HAMID 2021-11-7

clear; clc;

% Read table:

T = readtable('https://in.mathworks.com/matlabcentral/answers/uploaded_files/791474/excel-surface.xlsx')

% Fill outliers:

T = sortrows(T,1);

T.Z = filloutliers(T.Z, "spline", "movmedian", 4);

T = sortrows(T,2);

T.Z = filloutliers(T.Z, "spline", "movmedian", 4);

% Plot surface:

x = T.(1) ;

y = T.(2) ;

z = T.(3) ;

xi = linspace(min(x),max(x));

yi = linspace(min(y),max(y)) ;

[X,Y] = meshgrid(xi,yi) ;

Z = griddata(x,y,z,X,Y) ;

surf(X,Y,Z)

shading interp

colorbar

I tried sorting multiple times but it is not working. I now have a new weird outlier.

Before:

After:

Chris 2021-11-7

在 MATLAB Online 中打开

clear; clc;
% Read table:
T = readtable('https://in.mathworks.com/matlabcentral/answers/uploaded_files/791474/excel-surface.xlsx');
% Fill outliers:
for idx = 1:2
    for jdx = 8:-2:2
        T = sortrows(T,1);
        T.Z = filloutliers(T.Z, "spline", "movmedian", jdx);
        T = sortrows(T,2);
        T.Z = filloutliers(T.Z, "spline", "movmedian", jdx);
    end
end
% Plot surface:
x = T.(1) ;
y = T.(2) ;
z = T.(3) ;
xi = linspace(min(x),max(x));
yi = linspace(min(y),max(y)) ;
[X,Y] = meshgrid(xi,yi) ;
Z = griddata(x,y,z,X,Y) ;
surf(X,Y,Z)
% shading interp
colorbar

Getting pretty close.

With default shading, you'll see there are still a few bumps. But you can also see interesting features. Repeating the process, and changing the window size, continues to smooth the surface.

请先登录，再进行评论。