ellipt_ifilter.m
33bba04159d6master

ellipt_ifilter.m
View Options

			function out=ellipt_ifilter(L, a, b, C)
			[Ni, Nj] = size(L);
			out = zeros(Ni, Nj);
			Na = int32((Nj-1) * a); % x
			Nb = int32((Ni-1) * b); % y
			C = [int32((Nj-1) * C(1)), int32((Ni-1) * C(2))]; % center
			% ellipse parametrisation
			ellipt = @(x) int32(Nb * sqrt(1.0 - (double(x - C(2))/double(Na))^2));
			xmin = C(2) - Na + 1;
			xmax = C(2) + Na + 1;
			% loop through x
			for i=xmin:xmax
			range = ellipt(double(i-1));
			ymin = C(1) - range + 1;
			ymax = C(1) + range + 1;
			for j=ymin:ymax
			out(j,i) = L(j,i);
			end
			end
			end