Check if there are repeated elements in a matrix

Question

Yesbol 2018-3-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/388589-check-if-there-are-repeated-elements-in-a-matrix

编辑： Guillaume 2018-3-16

coord.mat

Dear all,

I have a set of lines/polylines represented in a matrix coord 380x4 (attached with this post), first and second columns representing X and Y coordinates respectively, and 4th column representing ID of a line/polyline. I want to check the continuity of this set. To do that, I need to check any point within one line/polyline have another point with the same coordinates on another line/polyline. Now, let me explain what I mean. For example, if we refer to the matrix coord attached with this thread, first ten rows (1:10) have the ID of 42 meaning that this polyline consists of 10 points connected one after another. The next polyline is presented from rows 11 to 17 having ID of 43. Now we can see that

coord(10,1) == coord(11,1) & coord(10,2) == coord(11,2)

which means that these two polylines are connected. The program is expected to check every set of points within each line/polyline and identify where there is a gap if it exists. For example, if the polyline with ID of 56 did not have any points present in other polylines, it would mean that this polyline is disconnected. And, if the line/polyline is disconnected, I want to see which ones are disconnected by showing their indexes.

Thanks to the answer by Guillaume, I realised that the approach needs to be changed. So I am not happy if there are more than two connected poly(lines) are disconnected. What I want to check is the continuity between the first line and the rest of the lines.

Any ideas guys please? I am stuck now.

Thank you

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

请先登录，再进行评论。

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

Answer 1

Guillaume 2018-3-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/388589-check-if-there-are-repeated-elements-in-a-matrix#answer_310235

在 MATLAB Online 中打开

The following should work so long as for a single ID the same point is not repeated. Also it won't detect if you have two disconnected groups of connected lines. Finally, the points coordinates have to be exactly equal down to the last significant digits:

[~, ~, pointid] = unique(coords(:, [1 2]), 'rows');  %assign same id to identical points
idrepeat = accumarray(pointid, 1);   %compute histogram of each id
pointrepeat = idrepeat(pointid);     %how many times each point of coords is found in coords
[label, ~, labelid] = unique(coords(:, 4));
labelconnections = accumarray(labelid, pointrepeat, [], @max) %for each label get the max of pointrepeat. A label has at least that many connections

Any label for which labelconnections is 1 does not have any of its point found in another label.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Guillaume 2018-3-15

在 MATLAB Online 中打开

I assume when you mean continuity you don't just mean that the end points of one polyline connect to the end points of another as that's fairly easy to check. Instead polyline can connect through any point (i.e. they cross each other).

In that case, I would build a connectivity graph and use the graph connectivity tools:

[~, ~, pointid] = unique(coord(:, [1 2]), 'rows');  %assign same id to identical points
[label, ~, labelid] = unique(coord(:, 4));
adjmat = zeros(numel(label));
adjcols = accumarray(labelid, pointid, [], @(pts) {labelid(ismember(pointid, pts))});
adjrows = repelem(1:numel(label), cellfun(@numel, adjcols))';
adjmat(sub2ind(size(adjmat), adjrows, vertcat(adjcols{:}))) = 1;
adjmat = adjmat .* ~eye(size(adjmat));
g = graph(adjmat, compose('%d', label));

You can easily see the connectivity graph of your labels:

plot(g);

You can clearly see on that plot that with your example matrix lines 58 and 59 are connected to each and nothing else. To get the groups:

groups = accumarray(conncomp(g)', label, [], @(labels) {labels})

Yesbol 2018-3-16

Thanks, Guillaume. Yes, you're right, by continuity I mean that polyline can connect through any point. Is there any other way to check it? Because after all I want that the disconnected line or set of lines will connect to the closest point of a line(s) which is/are connected

Thank you