2

我想将 3D 表面（CylCoors 300000x3）的纹理投影到 2D 平面（Image380x360）中。为此，我采用 Z ( UniqueZ=unique(CylCoors(:,3))) 和 Theta ( UniqueTheta=unique(CylCoors(:,1))) 中的每个唯一值，并投影所有纹理值 ( PointValues 300000x1 )，它们都像这样相遇：

Image=zeros(max(UniqueH),max(UniqueTheta)); %380x360

tic

HMat=bsxfun(@eq,CylCoors(:,3),UniqueH');  % 300000x380
ThetaMat=bsxfun(@eq,CylCoors(:,1),UniqueTheta'); %300000x360


for ii=1:length(UniqueH)          % Sloooow and not nice :(
    for jj=1:length(UniqueTheta)

        Image(ii,jj)=sum(PointValues.*...
            HMat(:,ii).*ThetaMat(:,jj))/...
            sum(HMat(:,ii).*ThetaMat(:,jj));

    end
end

toc

这是一个带有修剪变量的示例：

CylCoorsSample = [263.0000  184.2586   10.0000    
                  264.0000  183.0417   10.0000    
                  264.0000  182.1572   10.0000    
                  82.0000   157.4746   11.0000    
                  80.0000   158.2348   11.0000    
                  86.0000   157.3507   11.0000    
                  84.0000   157.7633   11.0000]   

PointValuesSample = [0.4745
                     0.5098
                     0.5020
                     0.4784
                     0.4510
                     0.4431
                     0.5804]


UniqueTheta = [80              
               82              
               84
               86
               263
               264]

UniqueH =[10
          11]

ThetaMat =                             HMat =   

 0     0     0     0     1     0       1   0
 0     0     0     0     0     1       1   0
 0     0     0     0     0     1       1   0
 0     1     0     0     0     0       0   1
 1     0     0     0     0     0       0   1
 0     0     0     1     0     0       0   1
 0     0     1     0     0     0       0   1 

Image =        % size: length(UniqueH)xlength(UniqueTheta)

   NaN       NaN       NaN       NaN    0.4745    0.5059
0.4510    0.4784    0.5804    0.4431       NaN       NaN

有没有办法使用矢量化for循环accumarray？有些东西告诉我情况就是这样，但我发现这个函数的多维子索引相当混乱。

也欢迎任何其他矢量化解决方案（使用 、 和我假设的智能组合reshape）bsxfun，sum但我真的很想accumarray更好地理解这种用法​​。

4

1 回答 1

2

鉴于：

CylCoorsSample = [263.0000  184.2586   10.0000    
                  264.0000  183.0417   10.0000    
                  264.0000  182.1572   10.0000    
                  82.0000   157.4746   11.0000    
                  80.0000   158.2348   11.0000    
                  86.0000   157.3507   11.0000    
                  84.0000   157.7633   11.0000]   

PointValuesSample = [0.4745
                     0.5098
                     0.5020
                     0.4784
                     0.4510
                     0.4431
                     0.5804]

您可以使用accumarray如下（使用第三个输出unique来生成所需的subs输入）：

[UniqueTheta, ~, subsTheta]=unique(CylCoorsSample(:,1))
[UniqueH,~,subsH]=unique(CylCoorsSample(:,3))
sz = [numel(UniqueH), numel(UniqueTheta)]

Image = accumarray([subsH, subsTheta], PointValuesSample, sz, @mean, NaN)
于 2015-06-10T13:17:33.350 回答