How can I keep overlapping regions that I see in Imfuse and remove everything else in the image?

3 次查看（过去 30 天）

AAS 2020-7-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/560432-how-can-i-keep-overlapping-regions-that-i-see-in-imfuse-and-remove-everything-else-in-the-image

编辑： AAS 2020-7-13

I used region props and have converted frames in a video binary images. How do I keep only the overlapping regions and remove everythign else? I would like to store the residual as frames to a video

2 个评论
显示无隐藏无

jonas 2020-7-6

Upload example images with your code?

AAS 2020-7-6

编辑：Image Analyst 2020-7-6

The right most panel is the imfused image of the first and second panel. I would like to keep only the regions that have an overlap and discard everything else as shown in the next image i.e just keep the blue circled regions .

请先登录，再进行评论。

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

采纳的回答

Image Analyst 2020-7-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/560432-how-can-i-keep-overlapping-regions-that-i-see-in-imfuse-and-remove-everything-else-in-the-image#answer_461945

This can be done in one line of code with imreconstruct().

J = imreconstruct(marker,mask) performs morphological reconstruction of the image marker under the image mask, and returns the reconstruction in J. The elements of marker must be less than or equal to the corresponding elements of mask. If the values in marker are greater than corresponding elements inmask, then imreconstruct clips the values to the mask level before starting the procedure.

14 个评论
显示 12更早的评论隐藏 12更早的评论

Image Analyst 2020-7-7

在 MATLAB Online 中打开

Try this:

% Initialization steps.  Brute force cleanup of everything currently existing to start with a clean slate.
clc;    % Clear the command window.
fprintf('Beginning to run %s.m ...\n', mfilename);
close all;  % Close all figures (except those of imtool.)
clear;  % Erase all existing variables. Or clearvars if you want.
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 20;
pre = imread('pre.png');
post = imread('post.png');
% Convert to logical
pre = logical(pre(:,:,1));
post = logical(post(:,:,1));
% Show image.
subplot(3, 2, 1);
imshow(pre);
axis('on', 'image');
impixelinfo;
title('Pre', 'FontSize', fontSize);
% Show image.
subplot(3, 2, 2);
imshow(post);
axis('on', 'image');
impixelinfo;
title('post', 'FontSize', fontSize);
% Find common/overlap areas
overlapped = pre & post;
% Show image.
subplot(3, 2, 3);
imshow(overlapped);
axis('on', 'image');
impixelinfo;
title('overlapped', 'FontSize', fontSize);
fused = imfuse(pre,post,'falsecolor','Scaling','joint','ColorChannels',[1 2 0]);
subplot(3, 2, 4);
imshow(fused)
axis('on', 'image');
impixelinfo;
title('fused', 'FontSize', fontSize);
% Make the union of images.
both = pre | post;
% Show image.
subplot(3, 2, 5);
imshow(both);
axis('on', 'image');
impixelinfo;
title('Both -- ORed, Union of images', 'FontSize', fontSize);
% Make output image.
output = imreconstruct(overlapped, both);
% Show image.
subplot(3, 2, 6);
imshow(output);
axis('on', 'image');
impixelinfo;
title('Output', 'FontSize', fontSize);

In the lower images, I zoomed in on one particular blob so you can see that it's getting all of the pre and post together in one blob. In the output there are only the 4 small blobs (and of course the huge white frame that you included) so it looks pretty similar to the "overlapped" image, just with slightly bigger blobs.

Image Analyst 2020-7-13

在 MATLAB Online 中打开

Sure. Just simply use pre in imreconstruct() instead of both:

% Initialization steps.  Brute force cleanup of everything currently existing to start with a clean slate.
clc;    % Clear the command window.
fprintf('Beginning to run %s.m ...\n', mfilename);
close all;  % Close all figures (except those of imtool.)
clear;  % Erase all existing variables. Or clearvars if you want.
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 20;
pre = imread('pre.png');
post = imread('post.png');
% Convert to logical
pre = logical(pre(:,:,1));
post = logical(post(:,:,1));
% Show image.
subplot(3, 2, 1);
imshow(pre);
axis('on', 'image');
impixelinfo;
title('Pre', 'FontSize', fontSize);
% Show image.
subplot(3, 2, 2);
imshow(post);
axis('on', 'image');
impixelinfo;
title('post', 'FontSize', fontSize);
% Find common/overlap areas
overlapped = pre & post;
% Show image.
subplot(3, 2, 3);
imshow(overlapped);
axis('on', 'image');
impixelinfo;
title('overlapped', 'FontSize', fontSize);
fused = imfuse(pre,post,'falsecolor','Scaling','joint','ColorChannels',[1 2 0]);
subplot(3, 2, 4);
imshow(fused)
axis('on', 'image');
impixelinfo;
title('fused', 'FontSize', fontSize);
% Make output image.
output = imreconstruct(overlapped, pre);
% Show image.
subplot(3, 2, 6);
imshow(output);
axis('on', 'image');
impixelinfo;
title('Output', 'FontSize', fontSize);