-
Notifications
You must be signed in to change notification settings - Fork 1
/
RegionGrowing.m
126 lines (104 loc) · 5.43 KB
/
RegionGrowing.m
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
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
function J = RegionGrowing(I, x, y, Threshold)
% Copyright (C) 2013 Heidelberg University
%
% Developed at CKM (Computerunterstützte Klinische Medizin),
% Medical Faculty Mannheim, Heidelberg University, Mannheim,
% Germany
%
%
% LICENCE
%
% CKM PhM Software Library, Release 1.0 (c) 2013, Heidelberg
% University (the "Software")
%
% The Software remains the property of Heidelberg University ("the
% University").
%
% The Software is distributed "AS IS" under this Licence solely for
% non-commercial use in the hope that it will be useful, but in order
% that the University as a charitable foundation protects its assets for
% the benefit of its educational and research purposes, the University
% makes clear that no condition is made or to be implied, nor is any
% warranty given or to be implied, as to the accuracy of the Software,
% or that it will be suitable for any particular purpose or for use
% under any specific conditions. Furthermore, the University disclaims
% all responsibility for the use which is made of the Software. It
% further disclaims any liability for the outcomes arising from using
% the Software.
%
% The Licensee agrees to indemnify the University and hold the
% University harmless from and against any and all claims, damages and
% liabilities asserted by third parties (including claims for
% negligence) which arise directly or indirectly from the use of the
% Software or the sale of any products based on the Software.
%
% No part of the Software may be reproduced, modified, transmitted or
% transferred in any form or by any means, electronic or mechanical,
% without the express permission of the University. The permission of
% the University is not required if the said reproduction, modification,
% transmission or transference is done without financial return, the
% conditions of this Licence are imposed upon the receiver of the
% product, and all original and amended source code is included in any
% transmitted product. You may be held legally responsible for any
% copyright infringement that is caused or encouraged by your failure to
% abide by these terms and conditions.
%
% You are not permitted under this Licence to use this Software
% commercially. Use for which any financial return is received shall be
% defined as commercial use, and includes (1) integration of all or part
% of the source code or the Software into a product for sale or license
% by or on behalf of Licensee to third parties or (2) use of the
% Software or any derivative of it for research with the final aim of
% developing software products for sale or license to a third party or
% (3) use of the Software or any derivative of it for research with the
% final aim of developing non-software products for sale or license to a
% third party, or (4) use of the Software to provide any service to an
% external organisation for which payment is received. If you are
% interested in using the Software commercially, please contact
% Prof. Dr. Lothar Schad (lothar.schad@medma.uni-heidelberg.de).
%
%Method Description:
%-------------------
%perform region growing based on seed points x and y and a certain
%threshold
if numel(x) ~= numel(y)
Error('WARNING: numel(x) ~= numel(y)')
else
[SizeY, SizeX] = size(I);
ListX = nan(1, numel(I));
ListY = nan(1, numel(I));
ListX(1, 1 : numel(x)) = x;
ListY(1, 1 : numel(x)) = y;
Counter = 1;
NumOfPixels = numel(x);
J = zeros(SizeY, SizeX);
NeighborsX = [0, 0, 1, -1];
NeighborsY = [1, -1,0 0];
while (~isnan(ListX(Counter)))
CurrentX = ListX(Counter);
CurrentY = ListY(Counter);
Counter = Counter + 1;
for NeighborIndex = 1 : 4
%Define Neightbor
NeighborX = CurrentX + NeighborsX(NeighborIndex);
NeighborY = CurrentY + NeighborsY(NeighborIndex);
%Check if Neighbor is in Image and not handled yet
if NeighborY > 1 && NeighborY <= SizeY && ...
NeighborX > 1 && NeighborX <= SizeX
if I(NeighborY, NeighborX) >= Threshold && J(NeighborY, NeighborX) == 0
J(NeighborY, NeighborX) = 1;
NumOfPixels = NumOfPixels + 1;
ListX(NumOfPixels) = NeighborX;
ListY(NumOfPixels) = NeighborY;
elseif J(NeighborY, NeighborX) == 0
J(NeighborY, NeighborX) = 2;
end
%End Valid Neighbor
end
%End Neighbor Index
end
%End While
end
%End Error Handling
end
end