FFilter.m
33bba04159d6master

FFilter.m
View Options

			% collect domains to filter, then apply
			classdef FFilter
			properties
			domains {} % cell of domains to verify
			% each element of the cell is composed by
			% an array [dom, flag] where:
			% dom = domain of filtering
			% flag = invert if false (apply domain filtering for all points outside the domain instead of inside)
			end
			methods
			function obj=FFilter(domains)
			K = length(domains);
			obj.domains = cell(K);
			% verify domain integrity
			for k=1:K
			if length(domains{k}) == 2
			obj.domains{k} = {domains{k}{1}, domains{k}{2}};
			else
			warning('Invalid input, skipping it')
			obj.domains = obj.domains{1:end-1};
			end
			end
			end
			% add a specific domain with the corresponding flag
			function add(obj,dom,flag)
			obj.domains = {obj.domains, [dom, flag]};
			end

			% img = a discmap object
			function out=apply(obj, img)
			[Ni, Nj] = size(img); % get image boundaries
			map = FDiscMap(img); % create a discrete map of the image
			out = zeros(Ni,Nj);
			minval = min(min(img));
			for i=1:Ni
			for j=1:Nj
			[x,y] = map.point(i,j); % get corresponding point
			if obj.mustfilter(x,y)
			out(i,j) = minval; % apply filter, minval supports log scaling
			else
			out(i,j) = img(i,j); % don't apply
			end
			end
			end
			end
			% find whether the point is inside at least one of the specific domains
			function match=mustfilter(obj,x,y)
			match = false;
			K = length(obj.domains);
			k = 1;
			while k <= K && ~match
			dom = obj.domains{k}{1};
			flag = obj.domains{k}{2};
			% oscure if
			match = match \| ~xor(dom.inside(x,y), flag);
			k = k + 1;
			end
			end
			end
			end