]> Creatis software - CreaPhase.git/blob - octave_packages/geometry-1.5.0/geom2d/intersectBoxes.m
Add a useful package (from Source forge) for octave
[CreaPhase.git] / octave_packages / geometry-1.5.0 / geom2d / intersectBoxes.m
1 %% Copyright (c) 2011, INRA
2 %% 2010-2011, David Legland <david.legland@grignon.inra.fr>
3 %% 2011 Adapted to Octave by Juan Pablo Carbajal <carbajal@ifi.uzh.ch>
4 %%
5 %% All rights reserved.
6 %% (simplified BSD License)
7 %%
8 %% Redistribution and use in source and binary forms, with or without
9 %% modification, are permitted provided that the following conditions are met:
10 %%
11 %% 1. Redistributions of source code must retain the above copyright notice, this
12 %%    list of conditions and the following disclaimer.
13 %%     
14 %% 2. Redistributions in binary form must reproduce the above copyright notice, 
15 %%    this list of conditions and the following disclaimer in the documentation
16 %%    and/or other materials provided with the distribution.
17 %%
18 %% THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
19 %% AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
20 %% IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
21 %% ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
22 %% LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR 
23 %% CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
24 %% SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS 
25 %% INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
26 %% CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
27 %% ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
28 %% POSSIBILITY OF SUCH DAMAGE.
29 %%
30 %% The views and conclusions contained in the software and documentation are
31 %% those of the authors and should not be interpreted as representing official
32 %% policies, either expressed or implied, of copyright holder.
33
34 %% -*- texinfo -*-
35 %% @deftypefn {Function File} {@var{box} =} intersectBoxes (@var{box1}, @var{box2})
36 %% Intersection of two bounding boxes.
37 %% 
38 %%   Example
39 %%
40 %% @example
41 %%   box1 = [5 20 5 30];
42 %%   box2 = [0 15 0 15];
43 %%   intersectBoxes(box1, box2)
44 %%   ans = 
45 %%       5 15 5 15
46 %% @end example
47 %%
48 %% @seealso{boxes2d, drawBox, mergeBoxes}
49 %% @end deftypefn
50
51 function bb = intersectBoxes(box1, box2)
52
53   % unify sizes of data
54   if size(box1,1) == 1
55       box1 = repmat(box1, size(box2,1), 1);
56   elseif size(box2, 1) == 1
57       box2 = repmat(box2, size(box1,1), 1);
58   elseif size(box1,1) != size(box2,1)
59       error('geom2d:Error',"Bad size for inputs.\n");
60   end
61
62   % compute extreme coords
63   mini = min(box1(:,[2 4]), box2(:,[2 4]));
64   maxi = max(box1(:,[1 3]), box2(:,[1 3]));
65
66   % concatenate result into a new box structure
67   bb = [maxi(:,1) mini(:,1) maxi(:,2) mini(:,2)];
68   
69 endfunction
70
71 %!test
72 %! box1 = [5 20 10 25];
73 %! box2 = [0 15 15 20];
74 %! res  = [5 15 15 20];
75 %! bb = intersectBoxes(box1, box2);
76 %! assert (res, bb, 1e-6);