GtPhase: added method to get neighbors of a grain.

Signed-off-by: Nicola Vigano <nicola.vigano@esrf.fr>

git-svn-id: https://svn.code.sf.net/p/dct/code/trunk@569 4c865b51-4357-4376-afb4-474e03ccb993

zUtil_DataStructures/GtPhase.m

@@ -436,6 +436,22 @@ classdef GtPhase < handle
                              & obj.bboxExtremes(:,6) >= point(3) ...
                             );
         end
+
+        function IDs = getNeighbours(obj, grainID)
+            numGrains = length(obj.grains);
+            replOnes = ones(numGrains, 1);
+
+            halfLengths = obj.boundingBox(:, 4:6) / 2;
+            grainHalfLengths = halfLengths(grainID, :);
+
+            maxDistances = halfLengths + grainHalfLengths(replOnes, :);
+
+            centersInPixels = obj.boundingBox(:, 1:3) + halfLengths;
+            grainCenter = centersInPixels(grainID, :);
+            relativeCenters = abs(centersInPixels - grainCenter(replOnes, :));
+
+            IDs = setdiff(find(all(relativeCenters < maxDistances, 2)), grainID);
+        end
         
         function idspot = buildIdSpot(obj, grainid)
             difspotID = obj.summary{grainid}.difspotID'; % row vector