Indexing for partial matrix

4 次查看（过去 30 天）
Stephen Thompson
Stephen Thompson 2020-7-17
评论： Stephen Thompson 2020-7-21
I have a 140 x 2000 matrix. It contains a blob of area (I know the linear indices) that has known peaks. I want to find the maximum peaks in part of the blob. Specifically I want to find the maximum peak in rows 1-77 and the maximum in rows 78-140 of the blob. There may be other blobs with other peaks, so it has to be by blob area.
Should I generate the linear indices for each section (e.g. all columns, rows 1:77) and then find which are common with the blob indices? How would I do that?

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

回答（1 个）

Matt J
Matt J 2020-7-17
编辑：Matt J 2020-7-17
[upperBlob,lowerBlob]=deal(yourMatrix);
upperBlob(78:end,:)=-inf;
lowerBlob(1:77,:)=-inf;
[imax1,jmax1]=find(upperBlob==max(upperBlob(:))); %upper blob max
[imax2,jmax2]=find(lowerBlob==max(lowerBlob(:))); %lower blob max
  3 个评论
显示 1更早的评论
Matt J
Matt J 2020-7-21
编辑：Matt J 2020-7-21
I don't see why it wouldn't have worked. The area outside your blobs appears to be zero, while the area inside the blobs appears to be greater than zero, so the maxima I've computed above would inherently have to fall within the blobs. Regardless, you could exclude the background with a few additional lines:
A=-inf(size(yourImage));
A(blobIndices)=yourImage(blobIndices):
[upperBlob,lowerBlob]=deal(A);
upperBlob(78:end,:)=-inf;
lowerBlob(1:77,:)=-inf;
[imax1,jmax1]=find(upperBlob==max(upperBlob(:))); %upper blob max
[imax2,jmax2]=find(lowerBlob==max(lowerBlob(:))); %lower blob max
Stephen Thompson
Stephen Thompson 2020-7-21
It is my fault for not explaining the problem well enough. Another example with 2 blobs - but I still would want them classed by row (corresponds for frequency but axes left generic) and by blob.
There could be multiple blobs with peaks above and below the row cut line, but I want to separate them into maxima by blob and by dichotomous row grouping, So it does have to be separated on a per-blob basis. Now that I can find the subscripts (ind2sub) for the blobs I can solve the problem. Thank you though for your solution.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Matrix Indexing 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by