How to remove multiple outlier data in a rectangular maze?

Question

Atanu 2022-6-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1738060-how-to-remove-multiple-outlier-data-in-a-rectangular-maze

评论： Atanu 2022-6-16

采纳的回答： Image Analyst

Maze4.mat

在 MATLAB Online 中打开

I have this rodent x,y coordinates as described in the attached picture. There are some outliers in the data, I need to remove them. I have to dynamically define a zone based on the extreme coordinates in the quardrants.

For example, let's say I need to remove the outlier data circled in red. The datapoint is in Maze4. I have attached the data for Maze4. I want to remove the bins where histcounts2 is < 2. I also need the 'xcoordinates2' and 'ycoordinates2' array after cleaning the outliers. I tried this so far.

h4 = histogram2(Maze4.xcoordinates2, Maze4.ycoordinates2, ...
            nbins,'DisplayStyle','tile','ShowEmptyBins','on');
counts4 = histcounts2(Maze4.xcoordinates2, Maze4.ycoordinates2, 25);
index4 = h4.Values(counts4<2);

But the index4 gives me 1-D array. How do I solve this?

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

请先登录，再进行评论。

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

Answer 1

Image Analyst 2022-6-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1738060-how-to-remove-multiple-outlier-data-in-a-rectangular-maze#answer_983200

在 MATLAB Online 中打开

First try to avoid the problem by not letting your rodents escape from your maze! 🐭🐹🐁🐀🤣

Then if you know your x coordinates must be between -80 and 80, use masking to extract only those good indexes:

% Find out which indexes are outside the -80 to 80 range.
mask = abs(x <= 80);
% Extract only those good indexes.
x = x(mask);
y = y(mask);

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

Image Analyst 2022-6-14

在 MATLAB Online 中打开

How about this, where you take all x data between the 5% and 95% points?

s = load('maze4.mat')
% Extract the table
t = s.Maze4
xOriginal = t.xcoordinates2;
yOriginal = t.ycoordinates2;
subplot(2, 2, 1)
plot(xOriginal, yOriginal, 'r.', 'MarkerSize', 8)
grid on
title('Original Data')
subplot(2, 2, 2)
[counts, edges] = histcounts(xOriginal, 100)
bar(edges(1:end-1), counts, 1)
grid on;
title('Histogram of X values')
theCDF = rescale(cumsum(counts) / sum(counts), 0, 100);
subplot(2, 2, 4)
plot(edges(1:end-1), theCDF, 'b-', 'LineWidth', 2)
title('CDF of Histogram')
grid on;
% Find the index of the 5% and 95% points
index1 = find(theCDF > 5, 1, 'first')
x1 = edges(index1)
xline(edges(index1), 'LineWidth', 2, 'Color', 'r')
% Find the index of the 5% and 95% points
index2 = find(theCDF > 95, 1, 'first')
x2 = edges(index2)
xline(edges(index2), 'LineWidth', 2, 'Color', 'r')
% Get a mask for the x values we want to exclude
indexesToKeep = t.xcoordinates2 >= x1 & t.xcoordinates2 <= x2;
xCleaned = xOriginal(indexesToKeep);
yCleaned = yOriginal(indexesToKeep);
subplot(2, 2, 3)
plot(xCleaned, yCleaned, 'r.', 'MarkerSize', 8)
grid on
title('Cleaned Data')