padarray矩阵的高效快速方式

时间:2017-05-24 06:04:17

标签: matlab performance zero-pad

我尝试了padarray超过1000张图片。但是,当我计算代码时,此特定行需要最长的时间才能完成

I=abs(padarray(I, [2, 2], 'replicate', 'both'));

主要是因为padarray算法的第35行(内部探查器): images \ private \ padarray_algo

b = a(aIdx{:});

有什么方法可以提高效率?也许用另一种方法?谢谢!

1 个答案:

答案 0 :(得分:2)

您可以使用repmatmatrix concatenation获得相同的结果:

r=@repmat;
pad=@(I,d)[r(I(1),d)    r(I(1,:),d(1),1)   r(I(1,end),d)
       r(I(:,1),1,d(2)) I                  r(I(:,end),1,d(2))
       r(I(end,1),d)    r(I(end,:),d(1),1) r(I(end),d)];

用法:

pad(I,[2 2])

如果所有图像的大小相同,您可以创建图像的线性索引矩阵并对其应用padarray,然后使用填充索引数组填充图像:

%create matrix of indices
Idx = reshape(1:numel(I),size(I));
%pad the index
Idx_padded = padarray(Idx, [2, 2], 'replicate', 'both'); 
%use the padded index to pad images
result = I(Idx_padded);
result2 = I2(Idx_padded);