Log In Sign Up
Models:
Tobias Harvey
tobiasharvey/
Lund-and-Heart-Segmentati
0
Downloads: 1
[1fedde]: / Thoracic Organs Segmentation code / MAP / fillingProc.m

Download this file

104 lines (81 with data), 3.3 kB 

 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
%function [finalImg]=fillAorta(i,j,aortaValue,filteredImg)

%function finalImg=fillAorta(i,j,filteredImg,inputImg)

function finalImg = fillingProc(centre_i,centre_j,filteredImg);





%maxPixelValue=250;

setValue=1;%max(max(filteredImg));

finalImg=zeros(size(filteredImg));

 

%threshold choice

threshold=1;%max(max(filteredImg));%mean(mean(filteredImg));



centre_i=centre_i;

centre_j=centre_j;



stack=-ones(255^2,2);

stack_index=1;

stack(stack_index,1)=centre_i;

stack(stack_index,2)=centre_j;

i=centre_i;j=centre_j;

while(stack_index>0 )%&& j+1<size(filteredImg,1) && j-1>0 && i+1<size(filteredImg,2) && i-1>0)

    %oldValue=finalImg(pixel_i,pixel_j);

    i=stack(stack_index,1);

    j=stack(stack_index,2);

    stack_index=stack_index-1;

    finalImg(i,j)=setValue;



    %if( checkNeighbour(filteredImg(i,j),filteredImg(i,j+1))==1 && finalImg(i,j+1)~=setValue )

    if(j+1<=size(filteredImg,2))

        if( finalImg(i,j+1)~=setValue)

        

                   if( checkNeighbour(filteredImg(i,j),filteredImg(i,j+1),threshold)==1 )

                         stack_index=stack_index+1;

                         stack(stack_index,1)=i;

                         stack(stack_index,2)=j+1;

                   else

                         finalImg(i,j+1)=setValue;

                   end



        end

    end

    

    %if( checkNeighbour(filteredImg(i,j),filteredImg(i,j-1))==1 && finalImg(i,j-1)~=setValue )

    if (j-1>0)

         if( finalImg(i,j-1)~=setValue)

        

                    if( checkNeighbour(filteredImg(i,j),filteredImg(i,j-1),threshold)==1 )

                         stack_index=stack_index+1;

                         stack(stack_index,1)=i;

                         stack(stack_index,2)=j-1;

                         %fillAorta(i,j-1);

                    else

                        finalImg(i,j-1)=setValue;

                    end

        end

    end



    %if( checkNeighbour(filteredImg(i,j),filteredImg(i+1,j))==1 && finalImg(i+1,j)~=setValue )

    if (i+1<=size(filteredImg,1))

        if( finalImg(i+1,j)~=setValue )

        

                    if( checkNeighbour(filteredImg(i,j),filteredImg(i+1,j),threshold)==1 )

                        stack_index=stack_index+1;

                        stack(stack_index,1)=i+1;

                        stack(stack_index,2)=j;

                        %fillAorta(i+1,j)

                    else

                        finalImg(i+1,j)=setValue;

                    end

      

        end

    end



    %if( checkNeighbour(filteredImg(i,j),filteredImg(i-1,j))==1 && finalImg(i-1,j)~=setValue )

    if (i-1>0)

    if( finalImg(i-1,j)~=setValue )

        

                if( checkNeighbour(filteredImg(i,j),filteredImg(i-1,j),threshold)==1 )

                       stack_index=stack_index+1;

                       stack(stack_index,1)=i-1;

                       stack(stack_index,2)=j;

                       %fillAorta(i-1,j)

                else

                       finalImg(i-1,j)=setValue;

                end

    end

      

    end

  

end



function isValid=checkNeighbour(curPixel,nextPixel,threshold)

isValid=1;

if (curPixel==threshold) 

    isValid=0;

end

if(abs(nextPixel-curPixel)>0) 

    isValid=0;

end