matlab/tools/ROIselectfull.m


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18

function V_out = ROIselectfull(input, ROI)

	s1 = size(input,1);
	s2 = size(input,2);
	[x y] = meshgrid(-(s2-1)/2:(s2-1)/2,(s1-1)/2:-1:-(s1-1)/2);
	A = Afstand(x,y,0,0);

	V_out = zeros(size(input));
	for slice = 1:size(input,3);
		V = input(:,:,slice);
		V(A > ROI/2) = 0;
		V_out(:,:,slice) = V; 
	end
end

function A = Afstand(x1,y1,x2,y2)
	A = sqrt((x1-x2).^2+(y1-y2).^2);
end