2013-12-11 78 views
0

我使用了三角测量的曲面近似,其在这里复制的码MATLAB ERROR IN PINV COMMAND

[totalTris,three] = size(tri); 
[totalPoints,two] = size(registeredPts); 

% 1. Find the 3 equations for each vertex, and 
%  place them in c_equations matrix; 
% c_equations = [A for vertex 1; 
%    A for vertex 2; ... 
%    A for vertex totalPoints] 
% c(point,row,:) gives one row from an A matrix 
Btotal = zeros(3,totalPoints); 
c_equations = zeros(3*totalPoints,3,9); 
for pointNum = 1:totalPoints 
    % B = [pixVal; x gradient; y gradient] at this vertex 
    z = pixelVals(pointNum); 
    B = [z; vGradientVecs(pointNum,1); vGradientVecs(pointNum,2)]; 

    % Compile all B matrices into a vector 
    Btotal(:,pointNum) = B; 

    % B = Ac 
    x = registeredPts(pointNum,1); 
    y = registeredPts(pointNum,2); 
    A = [1 x y x^2 y^2 x^3  (x^2)*y x*(y^2) y^3; ... 
     0 1 0 2*x 0 3*(x^2) 2*x*y y^2  0; ... 
     0 0 1 0 2*y 0  x^2  2*x*y 3*(y^2)]; 

    % Compile all A matrices into a vector 
    c_equations(pointNum,1,:) = A(1,:); 
    c_equations(pointNum,2,:) = A(2,:); 
    c_equations(pointNum,3,:) = A(3,:); 
end 

% 2. Find the c values for each triangle patch 
c = zeros(totalTris,9); 
c9 = zeros(9,9); 
for triNum = 1:totalTris 
    p1 = tri(triNum,1); 
    p2 = tri(triNum,2); 
    p3 = tri(triNum,3); 

    B9 = [Btotal(:,p1); Btotal(:,p2); Btotal(:,p3)]; 
    c9 = [c_equations(p1,1,:); c_equations(p1,2,:); c_equations(p1,3,:); ... 
      c_equations(p2,1,:); c_equations(p2,2,:); c_equations(p2,3,:); ... 
      c_equations(p3,1,:); c_equations(p3,2,:); c_equations(p3,3,:)]; 
    c(triNum,:) = pinv(c9)*B9; %linsolve(c9,B9); 
end 

end 

这是一个更大的码的一部分,当我运行的代码它给出了下面的错误MATLAB

Error using svd 
Input must be 2-D. 

Error in pinv (line 29) 
    [U,S,V] = svd(A,0); 
Error in findBPolyCoefficients (line 50) 
    c(triNum,:) = pinv(c9)*B9; 

根据我的搜索我有PINV的正确的语法,“输入必须是2D”的错误超出了我的理解,如果任何一个可以电话我的方式来解决这个问题

+0

'c9'的大小是多少?在MATLAB中执行'size(c9)'并将其粘贴到此处。 –

+0

我做了大小(c9)并得到了9 1 9。我必须通过这个c9,我从中了解到的是这9个矩阵的9行和1列不知道我是否正确 – Jav

回答

2

阿科根据你给出的信息,看起来c9是一个平凡的3-D矩阵,即它的第二维是1。这样的尺寸称为单体尺寸。 Squeeze是用于折叠此类尺寸的命令。你应该简单地做,

c9=squeeze(c9); %before doing pinv 

这将使c9 2 d矩阵,然后可以计算出它的SVD。

+0

感谢你让我的一天: )在这样做之后得到代码运行 – Jav

+1

@Jav只需确认'squeezing'得到的'c9'就能给你想要的最终结果。它可能已经解决了错误,但你仍然需要检查这是否是正确的做法。 –

+0

我无法得到正确的结果为ma代码,最后的形象是给问题,我认为这是因为它,但是当我会找到实际的原因,我会更新 – Jav