getPointIndexInSlice.m
bf44e1088b42master

getPointIndexInSlice.m
View Options

			function output = getPointIndexInSlice(obj,point)
			% The current implementation actually returns the point in
			% the slice that is the closest to the given point, even if
			% the latter is clearly not in the slice (!).
			%
			% The former implementation has been kept but commented out.
			% It was close to the solution I (Matthieu Boubat) proposed
			% here https://ch.mathworks.com/matlabcentral/answers/524181-how-can-i-retrieve-x-y-coordinates-of-a-sliced-plane-through-a-3d-volume-using-the-obliqueslice-func.
			% This solution has been discarded because it failed sometime
			% (did not find point indices).
			%
			% This method might be refactored at some point.

			pointIndices = obj.getPointIndexInVolume(point);

			coordinatesError = (abs(obj.slicedX-pointIndices(1)) + abs(obj.slicedY-pointIndices(2)) + abs(obj.slicedZ-pointIndices(3)));
			[i,j] = ind2sub(size(obj.sliced),...
			find( coordinatesError == min(coordinatesError,[],'all') ));
			% [i,j] = ind2sub(size(obj.sliced),...
			% find( (fix(obj.slicedX) == pointIndices(1) \| ceil(obj.slicedX) == pointIndices(1)) &...
			% (fix(obj.slicedY) == pointIndices(2) \| ceil(obj.slicedY) == pointIndices(2)) &...
			% (fix(obj.slicedZ) == pointIndices(3) \| ceil(obj.slicedZ) == pointIndices(3)) ));
			assert(not(isempty([i,j])),'Point has not been found in slice');
			output = [i(1),j(1)];
			end