about summary refs log tree commit diff
path: root/sourcecodes/bnt-master/KPMtools/rectintSparseC.m
blob: 401d7dc3fc34d884b4451dd683b38f7e67594494 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
function [overlap, normoverlap] = rectintSparseC(A,B)
%
% function [area, normarea] = rectintSparseC(A,B)
% A(i,:) = [x y w h]
% B(j,:) = [x y w h]
% out(i,j) = area of intersection
%
% Same as built-in rectint, but uses less memory.
% Also, returns area of overlap normalized by area of patch.
% See rectintSparse

if isempty(A) | isempty(B)
  overlap = [];
  normoverlap = [];
  return;
end

leftA = A(:,1);
bottomA = A(:,2);
rightA = leftA + A(:,3);
topA = bottomA + A(:,4);

leftB = B(:,1)';
bottomB = B(:,2)';
rightB = leftB + B(:,3)';
topB = bottomB + B(:,4)';

numRectA = size(A,1);
numRectB = size(B,1);

verbose = 0;
[overlap, normoverlap] = rectintSparseLoopC(leftA, rightA, topA, bottomA, leftB, rightB, topB, bottomB, verbose);