Made PixelReplicate able to deal with 3D images

SRC/MATLAB/IndToCoord.m

+
+% coords = IndToCoord(arraySize, arrayIdx)
+%   This function will take a list of indices and convert them into a
+%   coordinate with the same dimentionality of the array size in.
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+% Copyright (c) 2016, Drexel University
+% All rights reserved.
+% 
+% Redistribution and use in source and binary forms, with or without
+% modification, are permitted provided that the following conditions are met:
+% 
+% * Redistributions of source code must retain the above copyright notice, this
+%   list of conditions and the following disclaimer.
+% 
+% * Redistributions in binary form must reproduce the above copyright notice,
+%   this list of conditions and the following disclaimer in the documentation
+%   and/or other materials provided with the distribution.
+% 
+% * Neither the name of PixelRep nor the names of its
+%   contributors may be used to endorse or promote products derived from
+%   this software without specific prior written permission.
+% 
+% THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+% AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+% IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
+% DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE
+% FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+% DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
+% SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
+% CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
+% OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+% OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+function coords = IndToCoord(arraySize, arrayIdx)
+    coords = zeros(length(arrayIdx),length(arraySize));
+
+    linSize = [1 cumprod(arraySize)];
+    partialIdx = arrayIdx;
+    for i = length(arraySize):-1:1
+        r = rem(partialIdx-1, linSize(i)) + 1;
+        q = floor((partialIdx-r) / linSize(i)) + 1;
+
+        coords(:,i) = q;
+        partialIdx = r;
+    end
+end
+

SRC/MATLAB/PixelReplicate.m

@@ -42,11 +42,9 @@
 function ptsReplicate = PixelReplicate(im, label)
 
 if nargin>1
-    [r c] = find(im==label);
-    pts = [c,r];
+    ind = find(im==label);
 else
-    [r c] = find(im);
-    pts = [c,r];
+    ind = find(im);
     % NOTE - this is turned off for speed -- but be sure to 
 %     [L num]=bwlabel(im);
 %     if num>1
@@ -54,26 +52,27 @@ else
 %     end
 end    
     
-bwd=DepthImageFromPoints(pts,im);
+bwd = DepthImageFromPoints(ind,im);
 
-ptsReplicate=zeros(sum(round(bwd(:))),2);
+indsReplicate = []; 
 idx = find(bwd);
-nInsert=1;
 
 % note - you could use a parfor here for speed with large components
 for i=1:length(idx)
-    [r c]=ind2sub(size(bwd),idx(i));
     nrep = round(bwd(idx(i)));
     % speed up by using pre-allocated ptsReplicate
-    ptsReplicate(nInsert:nInsert+nrep-1,:) = repmat([r,c],nrep,1);
-    nInsert=nInsert+nrep;
+    indsReplicate = [indsReplicate;repmat(idx(i),nrep,1)];
 end
 
-function [bwd bw] = DepthImageFromPoints(pts,im)
+ptsReplicate = SwapXY_RC(IndToCoord(size(im),indsReplicate));
+end
+
+function [bwd, bw] = DepthImageFromPoints(ind,im)
 
 bw = false(size(im));
-pts=round(pts);
-pts=unique(pts,'rows');
-pts = sub2ind(size(bw),pts(:,1),pts(:,2));
-bw(pts)=1;
+ind = unique(ind);
+
+bw(ind) = true;
 bwd = bwdist(~bw);
+end
+

SRC/MATLAB/SwapXY_RC.m

+% swapCoords = SwapXY_RC(coords)
+%   This will swap the first two colums and is useful when changing between
+%   row and column centric coordinate systems (e.g. (row,col) -> (x,y)).
+
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+% Copyright (c) 2016, Drexel University
+% All rights reserved.
+% 
+% Redistribution and use in source and binary forms, with or without
+% modification, are permitted provided that the following conditions are met:
+% 
+% * Redistributions of source code must retain the above copyright notice, this
+%   list of conditions and the following disclaimer.
+% 
+% * Redistributions in binary form must reproduce the above copyright notice,
+%   this list of conditions and the following disclaimer in the documentation
+%   and/or other materials provided with the distribution.
+% 
+% * Neither the name of PixelRep nor the names of its
+%   contributors may be used to endorse or promote products derived from
+%   this software without specific prior written permission.
+% 
+% THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+% AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+% IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
+% DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE
+% FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+% DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
+% SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
+% CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
+% OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+% OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+function swapCoords = SwapXY_RC(coords)
+    swapCoords = coords;
+    
+    if ( isvector(coords) )
+        swapCoords([1 2]) = swapCoords([2 1]);
+    elseif ( ismatrix(coords) )
+        swapCoords(:,[1 2]) = swapCoords(:,[2 1]);
+    end
+end