finding submatrix and index

4 次查看（过去 30 天）
Mohammad Golam Kibria
Hi I have the following matrix as two small sample matrix. I want to fill the matrix value zero with 1 between two diagonal 1.
I1=
0 0 0
0 0 1
0 0 0
1 0 0
I2=
1 0 0
0 0 0
0 0 1
0 0 0
I want to have my output as follows:
I1=
0 0 0
0 0 1
0 1 0
1 0 0
I2=
1 0 0
0 1 0
0 0 1
0 0 0
Thanks in advanced.

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

回答（2 个）

Andrei Bobrov
Andrei Bobrov 2013-5-20
编辑：Andrei Bobrov 2013-5-20
I1out = imdilate(I1,flipud(eye(min(size(I1)))));
I2out = imdilate(I2,eye(min(size(I1))));
ADD after Mohammad Golam Kibria's comment
s = size(I1);
d = 1-size(I1):size(I1,2)-1;
x = spdiags(I1(end:-1:1,:),d);
x1 = cumsum(x)&cumsum(x(end:-1:1,:));
I1out = full(spdiags(flipud(x1),d,s(1),s(2)));
  2 个评论
Mohammad Golam Kibria
Mohammad Golam Kibria 2013-5-20
编辑：Mohammad Golam Kibria 2013-5-20
This does not work well. for example, if I1 =
0 0 1 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
1 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 1 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
1 0 0 0 0 0 0 0 0 0
then I1out = imdilate(I1,flipud(eye(min(size(I1))))); creates the following output:
-Inf 0 0 1 0 0 0 0 0 0
0 0 1 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0 0 0
1 0 0 0 0 0 0 1 0 0
0 0 0 0 0 0 1 0 0 0
0 0 0 0 0 1 0 0 0 0
0 0 0 0 1 0 0 0 0 0
0 0 0 1 0 0 0 0 0 0
0 0 1 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0 0 0
but I need the output as follows:
0 0 1 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0 0 0
1 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 1 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0 0 0
1 0 0 0 0 0 0 0 0 0
could u please help me any more.
Image Analyst
Image Analyst 2013-5-20
You would have to tell it which pairs to connect. Once you specify that, you can use my answer.

请先登录，再进行评论。

Image Analyst
Image Analyst 2013-5-20
You need to use Bresenham's Line Algorithm: http://en.wikipedia.org/wiki/Bresenham%27s_line_algorithm. If you have the Image Processing Toolbox, you can do that with imline(). Let me know if you want a demo of imline.
  2 个评论
Mohammad Golam Kibria
Mohammad Golam Kibria 2013-5-21
yes, I need a demo. then it will be clear to me.
Image Analyst
Image Analyst 2013-5-21
Here's the demo:
% Demo to write an ellipse and a line into the overlay of an image,
% and then to burn those overlays into the image.
%----- Initializing steps -----
% Clean up
clc;
clear all;
close all;
workspace; % Display the workspace panel.
hasIPT = license('test', 'image_toolbox');
if ~hasIPT
% User does not have the toolbox installed.
message = sprintf('Sorry, but you do not seem to have the Image Processing Toolbox.\nDo you want to try to continue anyway?');
reply = questdlg(message, 'Toolbox missing', 'Yes', 'No', 'Yes');
if strcmpi(reply, 'No')
% User said No, so exit.
return;
end
end
% Display images to prepare for the demo.
monochromeImage = imread('pout.tif');
subplot(2, 4, 1);
imshow(monochromeImage);
title('Original Image');
subplot(2, 4, 2);
imshow(monochromeImage);
title('Original Image with ellipse in overlay');
subplot(2, 4, 5);
imshow(monochromeImage);
title('Original Image');
subplot(2, 4, 6);
imshow(monochromeImage);
title('Original Image with line in overlay');
set(gcf, 'units','normalized','outerposition',[0 0 1 1]); % Maximize figure.
set(gcf,'name','Image Analysis Demo','numbertitle','off')
%----- Burn ellipse into image -----
% Create elliptical mask, h, as an ROI object over the second image.
subplot(2, 4, 2);
hEllipse = imellipse(gca,[10 10 50 150]); % Second argument defines ellipse shape and position.
% Create a binary image ("mask") from the ROI object.
binaryImage = hEllipse.createMask();
% Display the ellipse mask.
subplot(2, 4, 3);
imshow(binaryImage);
title('Binary mask of the ellipse');
% Let's try to add some text. (Doesn't work)
% hText = text(50, 100, 'Line of Text');
% textMask = hText.createMask();
% binaryImage = binaryImage & textMask;
% imshow(binaryImage);
% Burn ellipse into image by setting it to 255 wherever the mask is true.
monochromeImage(binaryImage) = 255;
% Display the image with the "burned in" ellipse.
subplot(2, 4, 4);
imshow(monochromeImage);
title('New image with ellipse burned into image');
%----- Burn line into image -----
burnedImage = imread('pout.tif');
% Create line mask, h, as an ROI object over the second image in the bottom row.
subplot(2, 4, 6);
hLine = imline(gca,[10 100],[10 100]); % Second argument defines line endpoints.
% Create a binary image ("mask") from the ROI object.
binaryImage2 = hLine.createMask();
% Display the line mask.
subplot(2, 4, 7);
imshow(binaryImage2);
title('Binary mask of the line');
% Burn line into image by setting it to 255 wherever the mask is true.
burnedImage(binaryImage2) = 255;
% Display the image with the "burned in" line.
subplot(2, 4, 8);
imshow(burnedImage);
title('New image with line burned into image');

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Read, Write, and Modify Image 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by