Cross in 2D

For reference see incontrastsRef.html.

function q = cross2D(X1,X2,varargin)

seti = varargin{end};
R = seti.rCD/2;

w = 1/8*R; % width
l = 7/8*R; % length
w = w/2;
l = l/2;

q1 = (-w < X1) & (X1 < w) & (-l < X2) & (X2 < l); % vertical line
q2 = (-w < X2) & (X2 < w) & (-l < X1) & (X1 < l); % horizontal line
q = q1 | q2; % to avoid q1+q2 (because some elements would have contrast 2...)

q = double(q);

end