How to add values from vertcat to a script

Question

Tipu Sultan 2019-4-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/458689-how-to-add-values-from-vertcat-to-a-script

编辑： Stephen23 2019-4-29

I have script which calls a function and the return values of the function are stored in this way:

M = vertcat(C{:});

C is 1*99 cell and each cell has 4 values.Now I want to store individual value of 1st column to an array 2nd column to another array and so on.So is it possible to do so?

11 个评论
显示 9更早的评论隐藏 9更早的评论

Tipu Sultan 2019-4-29

在 MATLAB Online 中打开

I am calling this function as below.

D = 'C:\Users\Tipu_Standard\Desktop\My_project_cropped_images1';
S = dir(fullfile(D,'Original*.jpg')); % pattern to match filenames.
N = numel(S);
C = cell(1,N);
for k = 1:N
    F = fullfile(D,S(k).name);
    I = imread(F);
    C{k} = fun1(I);
end
M = vertcat(C{:});

My function is below.

function [z,distance,angle] = fun1(I)
    B=rgb2gray(I);
    A=adapthisteq(B);
    % Display the original image
    figure,imshow(A);
    distance =0;
    angle=0;
    rmin=7;
    rmax=22;
    % Find all the circles with radius >= 8 and radius <= 30
    [centersdark,radiidark] = imfindcircles(A,[rmin rmax],'Objectpolarity','dark','Sensitivity',0.913);
    k=numel(radiidark);
    viscircles(centersdark,radiidark ,'EdgeColor','b');          
    if(k==0)
        z=0;
    else
        z=1;
        p=2*radiidark;
        f=2.1;
        w=62;
        h=144;
        s=2.88;
        distance=(f*w*h)/(p*s);
        a=(192/2)-centersdark(1); 
        b=144-centersdark(2);
        angle=atan2(a,b)*(180/pi);
    end
    %tri=[z,distance,angle];
end

Tipu Sultan 2019-4-29

在 MATLAB Online 中打开

Thanks @stephen.

I changed the code as follows:

D = 'C:\Users\Tipu_Standard\Desktop\My_project_cropped_images1';
S = dir(fullfile(D,'Original*.jpg')); % pattern to match filenames.
N = numel(S);
C = cell(1,N);
for k = 1:N
    F = fullfile(D,S(k).name);
    I = imread(F);
    %C{k} = fun1(I);
    Zvec = zeros(1,20);
    Dvec = zeros(1,20);
    Avec = zeros(1,20);
    [Zvec(k),Dvec(k),Avec(k)] = fun1(I);
end
M = vertcat(C{:});

Command window gives me following error:

In an assignment  A(I) = B, the number of elements in B and I must be the same.
Error in funCall (line 12)
    [Zvec(k),Dvec(k),Avec(k)] = fun1(I);

Is there any problem that I am not able to locate!

Stephen23 2019-4-29

编辑：Stephen23 2019-4-29

在 MATLAB Online 中打开

Move the preallocation before the loop (just like the preallocation documentation shows):

D = 'C:\Users\Tipu_Standard\Desktop\My_project_cropped_images1';
S = dir(fullfile(D,'Original*.jpg')); % pattern to match filenames.
N = numel(S);
Zvec = zeros(1,N);
Dvec = zeros(1,N);
Avec = zeros(1,N);
for k = 1:N
    F = fullfile(D,S(k).name);
    I = imread(F);
    [Zvec(k),Dvec(k),Avec(k)] = fun1(I);
end

You will also need to adjust the sizes of the output arrays AND the indexing to suit the sizes of the output arguments: if they are non-scalar then you will continue to get errors.

I do not have your data, you will have to check this yourself. After that, adjust the sizes of the output arrays. If the sizes change on between iterations, use cell arrays for the appropriate array/s (in fact if you just want to get your code working, use cell arrays for all three outputs).

请先登录，再进行评论。

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

Answer 1

madhan ravi 2019-4-26

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/458689-how-to-add-values-from-vertcat-to-a-script#answer_372406

在 MATLAB Online 中打开

Wanted = num2cell(M,1)
[a,b,c,d]=deal(Wanted{:})

3 个评论
显示 1更早的评论隐藏 1更早的评论

Tipu Sultan 2019-4-26

编辑：Tipu Sultan 2019-4-28

在 MATLAB Online 中打开

This is my first script:

D = 'C:\Users\Tipu_Standard\Desktop\My_project_cropped_images';
S = dir(fullfile(D,'Original*.jpg')); % pattern to match filenames.
N = numel(S);
C = cell(1,N);
for k = 1:N
    F = fullfile(D,S(k).name);
    I = imread(F);
    C{k} = fun1(I);
end
N = vertcat(C{:});

now the function 'fun1(I)' is like this:

function tri=fun1(I)    
    B=rgb2gray(I);
    A=adapthisteq(B);
    % Display the original image
    figure,imshow(A);
    distance =0;
    angle=0;
    rmin=8;
    rmax=30;
    % Find all the circles with radius >= 8 and radius <= 30
    [centersdark,radiidark] = imfindcircles(A,[rmin rmax],'Objectpolarity','dark','Sensitivity',0.913);
    k=numel(radiidark);
    viscircles(centersdark,radiidark ,'EdgeColor','b');          
    if(k==0)
        z=0;
    else
        z=1;
        p=2*radiidark;
        f=2.1;
        w=62;
        h=144;
        s=2.88;
        distance=(f*w*h)/(p*s);
        a=(192/2)-centersdark(1); 
        b=144-centersdark(2);
        angle=atan2(a,b)*(180/pi);
    end
    tri=[z,distance,angle];
end

And now I want to send the value of 'N' of the first script above to a new script which is as follows:

S = [100 0 0 0 0 0;
    0 100 0 0 0 0;
    0 0 100 0 0 0;
    0 0 0 100 0 0;
    0 0 0 0 100 0;
    0 0 0 0 0 100];
prev_S = S;
Big_lambda = eye(2);
a=0;
b=0;
c=0;
p=0;
q=0;
s=0;
est_vec=[a ; b; c; p; q; s];
theta = [45 60 90 ] ;
t= [ 10 20 30 ];
     
r= [ 10 20 30];
Wanted = num2cell(M,1)
[r,b,distance,angle]=deal(Wanted{:})
x = [ r; theta ; t];
 time = 40;
for i=1:3
     dif_x = [(-cos(theta(i))) (r(i).*sin(theta(i))) (2*a.*t(i)+b);(-sin(theta(i))) (-r(i).*cos(theta(i))) (2*p.*t(i)+q)];    
     
     W = dif_x(i)*Big_lambda*dif_x(i)';
     
     M = [t(i)^2 t(i) 1 0 0 0;
         0 0 0 t(i)^2 t(i) 1];
     
     y = M * est_vec + dif_x * x(:,i);
     
     K = prev_S*M'/((W + M*prev_S*M'));
     
     S = prev_S;
          
     est_vec_new = est_vec +  K*(y-M*est_vec);
     
     cond = abs(est_vec_new - est_vec);
     
     if  cond < 0.003
         break 
     end
     
     est_vec = est_vec_new;
     
     S_new = (eye(6) - K*M)*S;
     
     S = S_new; 
     meas_equa1 = a*time.^2+b*time+c; % measurement euation
     meas_equa2 = p*time.^2+q*time+r; % measurement equation
 end

Now I want to store the values of 'z,distance,angle' which are derived from the function fun1(I) to the last script array t,r,theta respectively i.e t=z,r=distance,theta=angle. So is it possible to do that?

Any help will be greatly appreciated.

madhan ravi 2019-4-26

That seems like a different question , I answered to the question you originally asked. It’s better you post it as a separate question.

Tipu Sultan 2019-4-28

@madhan I posted the question but it was closed due to duplicate!

So please help me with the problem.

请先登录，再进行评论。

How to add values from vertcat to a script

11 个评论
显示 9更早的评论隐藏 9更早的评论

回答（1 个）

3 个评论
显示 1更早的评论隐藏 1更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

How to add values from vertcat to a script

11 个评论 显示 9更早的评论隐藏 9更早的评论

回答（1 个）

3 个评论 显示 1更早的评论隐藏 1更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

11 个评论
显示 9更早的评论隐藏 9更早的评论

3 个评论
显示 1更早的评论隐藏 1更早的评论