Prev: want to some code
Next: Legendre coefficients
From: Manal on 18 Jun 2010 10:09
"someone" <someone(a)somewhere.net> wrote in message <hvftm5$f59$1(a)fred.mathworks.com>...
>
> % To find the index (idx) of values that have the minimum distance:
> mind = min(dt(:,1));
> idx = find(dt(:,1) - mind)

idx = find(dt(:,1) - mind)
this returns a vector of values
while it should be two values ( x values and y values ) !!!!!!
From: Manal on 18 Jun 2010 10:50
% Reorder the Matrix s
n=x-1;
u=1;
for k=1:n
A(u,1)=s(k,1);
A(u,2)=s(k,2);
s(k,:)=0;
i=1;
dt=0;
for o=1:n
xt=s(o,2);
yt=s(o,1);
dt(i,1) = sqrt((xt-A(u,2))^2 + (yt-A(u,1))^2);
dt(i,2) = A(u,1);
dt(i,3) = A(u,2);
i=i+1;
end
u=u+1;
mind = min(dt(:,1));

idx = find(dt(:,1)== mind);
A(u,1)=dt(idx,2);
A(u,2)=dt(idx,3);
u=u+1;

idx2 = find(s(:,1)== dt(idx,2));
s(idx2,:)=0;
end
disp(A);
disp('------------------------------------------');
disp(s);


this is the code that I have implemented so far
but I've got an errors in these lines

A(u,1)=dt(idx,2);
A(u,2)=dt(idx,3);

idx2 = find(s(:,1)== dt(idx,2));

any help ?? :(
From: someone on 18 Jun 2010 15:12
"Manal " <be.san(a)live.com> wrote in message <hvg12v$o3o$1(a)fred.mathworks.com>...
> % Reorder the Matrix s
> n=x-1;
> u=1;
> for k=1:n
> A(u,1)=s(k,1);
> A(u,2)=s(k,2);
> s(k,:)=0;
> i=1;
> dt=0;
> for o=1:n
> xt=s(o,2);
> yt=s(o,1);
> dt(i,1) = sqrt((xt-A(u,2))^2 + (yt-A(u,1))^2);
> dt(i,2) = A(u,1);
> dt(i,3) = A(u,2);
> i=i+1;
> end
> u=u+1;
> mind = min(dt(:,1));
>
> idx = find(dt(:,1)== mind);
> A(u,1)=dt(idx,2);
> A(u,2)=dt(idx,3);
> u=u+1;
>
> idx2 = find(s(:,1)== dt(idx,2));
> s(idx2,:)=0;
> end
> disp(A);
> disp('------------------------------------------');
> disp(s);
>
>
> this is the code that I have implemented so far
> but I've got an errors in these lines
>
> A(u,1)=dt(idx,2);
> A(u,2)=dt(idx,3);
>
> idx2 = find(s(:,1)== dt(idx,2));
>
> any help ?? :(

doc find

% Basically, (if i remember right),
% find returns the index of a vector whose value is zero.
% so, try replacing
idx2 = find(s(:,1)== dt(idx,2));
% with something like
idx2 = find(s(:,1) - dt(idx,2));
From: us on 18 Jun 2010 15:56
"Manal " <be.san(a)live.com> wrote in message <hve70o$2dg$1(a)fred.mathworks.com>...
> Hi all
> I have complicated case and I can't deal with
> hope some could give suggestions
>
> http://img687.imageshack.us/img687/7439/screenshot20100617at230.png
>
>
> in the picture above the border's values are stored in Matrix B ( the index x and y for each point)
>
>
> it is stored ordered by the value over the axis y
>
>
> I want to change the matrix order to be start with one point draw all the next points and end with the same point ( like drawing a circle )
>
> its complicated :'(
>
> any ideas ???
>
>
> thanks in advance

a hint:
- this FEX submission will help you to achieve what you want...
- look at the accompanying published M-file (towards the end)...

http://www.mathworks.com/matlabcentral/fileexchange/6760

us
From: ImageAnalyst on 18 Jun 2010 18:14
On Jun 18, 6:53 am, "Manal " <be....(a)live.com> wrote:

>
> thanks for your replay
>
> the centroid is not a problem since its drawn in another stage
> the problem is only in the boundary
>
> I was thinking to solve it from earlier stage as you have suggested but I did not found a solution
> this boundary has been extracted from the image using Sobel edge detection
> B = edge(im1,'sobel'); %finding edges
>
> after that I extracted the none zero values from the B to be the boundary points
> if B(i,j)>0
>     s(x,1)=i;
>     s(x,2)=j;
>     x=x+1;
> end
>
> I know that the problem start from here
> but I don't have any idea how can I solve it
---------------------------------------------------------------------------------------

Aha! Now the problem has suddenly become much clearer. You started
out with a bad image to begin with - the one that edge gave you. What
you should have done was to use fspecial() and imfilter() instead of
edge() and then pick a more appropriate threshold than edge() does.
Then use imclose() to do a morphological closing on your edge blobs to
close up the gaps. Then use bwboundaries. Like I do in this demo to
find boundaries around that difficult cell in the standard cell.tif
image:

clc; % Clear the command window.
close all; % Close all figures (except those of imtool.)
imtool close all; % Close all imtool figures.
clear; % Erase all existing variables.
workspace; % Make sure the workspace panel is showing.
fontSize = 20;

% Change the current folder to the folder of this m-file.
if(~isdeployed)
cd(fileparts(which(mfilename)));
end

% Read in a standard MATLAB gray scale demo image.
folder = 'C:\Program Files\MATLAB\R2010a\toolbox\images\imdemos';
baseFileName = 'cell.tif';
fullFileName = fullfile(folder, baseFileName);
grayImage = imread(fullFileName);
% Get the dimensions of the image. numberOfColorBands should be = 1.

% Display the original gray scale image.
subplot(2, 2, 1);
imshow(grayImage, []);
title('Original Grayscale Image', 'FontSize', fontSize);
set(gcf, 'Position', get(0,'Screensize')); % Enlarge figure to full
screen.

% Do a Sobel filter - a real one, not the thresholded one edge() does.
H = fspecial('sobel');
verticalEdges = imfilter(grayImage, H, 'replicate');
horizontalEdges = imfilter(grayImage, H', 'replicate');
edgeImage = max(verticalEdges, horizontalEdges);
subplot(2,2,2);
imshow(edgeImage, []);
title('Edge Image', 'FontSize', fontSize);

% Make the edges more "solid" by closing them
se = strel('disk',10);
closedImage = imclose(edgeImage, se);
subplot(2,2,3);
imshow(closedImage, []);
title('Closed Image', 'FontSize', fontSize);

binaryImage = closedImage > 50;
subplot(2,2,4 );
imshow(binaryImage, []);
title('Thresholded Image', 'FontSize', fontSize);

% Find boundaries
boundaries = bwboundaries(binaryImage);

% Done! Now, plot the boundaries over the original image
numberOfBoundaries = size(boundaries);
subplot(2,2,1);
hold on;
for k = 1 : numberOfBoundaries
thisBoundary = boundaries{k};
plot(thisBoundary(:,2), thisBoundary(:,1), 'g', 'LineWidth', 2);
end
hold off;

Don't be put off by the above - it's mostly just code to initialize
and display things - the key code (calls to fspecial, imfilter,
imclose, thresholding, and bwboundaries) is really only 5 lines. Try
it with your image. If it doesn't work, then post your original image
somewhere.

It's much better to start off doing things in the right way than to
try some programmatical gynmnastics to fix up a bad result.


First  |  Prev  |  Next  |  Last
Pages: 1 2 3 4 5 6
Prev: want to some code
Next: Legendre coefficients