使用MATLAB查找坐标列表中每个坐标的第一个和最后一个索引

作者: 雨儿
发布时间: 2024-12-15 11:21:13 (2月前)
转自：

4 条回复

0#
回复此人
甲基蓝 | 2019-08-31 10-32

<div class =“post-text”itemprop =“text”> <P> 假设坐标是正整数，特别是当坐标范围很小时，你可以使用 <a href="https://www.mathworks.com/help/matlab/ref/accumarray.html" rel="nofollow noreferrer"> <code> accumarray </code> </A> ： </p> <pre> <code> x1 = min(x); y1 = min(y); x2 = max(x); y2 = max(y); sz = [y2-y1+1, x2-x1+1]; ind = sub2ind(sz,y-y1+1, x-x1+1); ifirst = accumarray(ind(:), 1:numel(ind), [], @min); ilast = accumarray(ind(:), 1:numel(ind), [], @max); Mfirst = ifirst(ind); Mlast = ilast(ind); </code> </pre> <P> 对于更高的范围，您可以使用稀疏选项： </p> <pre> <code> ifirst = accumarray(ind(:), 1:numel(ind), [], @min,[],1); ilast = accumarray(ind(:), 1:numel(ind), [], @max,[],1); </code> </pre> </DIV>

编辑
1#
回复此人
早岁那知世事艰 | 2019-08-31 10-32

<div class =“post-text”itemprop =“text”> <P> 该 <a href="https://www.mathworks.com/help/matlab/ref/unique.html" rel="noreferrer"> <code> unique </code> </A> 功能可以做到这一点： </p> <pre> <code> [~, b, c] = unique([x(:) y(:)], 'rows', 'first'); first = b(c).'; [~, b, c] = unique([x(:) y(:)], 'rows', 'last'); last = b(c).'; </code> </pre> </DIV>

编辑
2#
回复此人
薄情 | 2019-08-31 10-32

<div class =“post-text”itemprop =“text”> <P> 如果你有30M点，那么你可能没有足够的内存用于这种方法......但对于较小的阵列来说它很快 </p> <pre> <code> x = [1 3 7 1 3]; y = [5 1 6 5 1]; xy = cat( 3, x, y ); chk = all( xy == permute( xy, [2 1 3] ), 3 ); [~,first] = max( chk ); [~,last] = max( flipud( chk ) ); last = size(chk,1) - last + 1; </code> </pre> <HR /> <P> 的<strong> 编辑 </强> 你也可以这样做 <code> findgroups </code> ，循环遍历唯一的坐标而不是每个坐标，以便有一个可能更短的循环...... </p> <pre> <code> x = [1 3 7 1 3]; y = [5 1 6 5 1]; g = findgroups( x, y ); first = zeros( size( x ) ); last = first; for ii = 1:max(g) idx = (ii==g); first( idx ) = find( idx, 1, 'first' ); last( idx ) = find( idx, 1, 'last' ); end </code> </pre> <HR /> <P> 的<strong> EDIT2 </强> 我认为这些相对于其他答案都相当慢...... </p> </DIV>

编辑

登录后才能参与评论