|
Angelika |
Forum-Newbie
|
|
Beiträge: 1
|
|
|
|
Anmeldedatum: 15.06.18
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 15.06.2018, 09:39
Titel: Region Growing
|
|
|
|
|
Hallo meine Lieben,
ich muss ein Region Growing programmieren, was mir im Idealfall die Leber als Region erkennt. Ich hab einen aktuellen Algorithmus programmiert, jedoch muss das Kriterium (jetzt Q) optimiert werden, denn er nimmt immer zu viel mit. Hat einer/eine von euch eine Idee oder Tipps, was ich anders machen und verbessern kann?
Wäre wirklich super nett !
DANKE!
I=imread('17.bmp');
imshow(I);
%Mausklick
[y,x] =getpts;
x=uint16(x);
y=uint16(y);
%definiere Saatpunkt
S=[x y];
%ermittle Bildgröße
[height ,width,farbtiefe]=size(I);
% erstelle ein leeres Bild, indem später alle Saatpunkte abgespeichert
% werden
Iopenseeds=zeros(height,width);
%übertrage den ersten Saatpunkt in dieses Bild
Iopenseeds(S(1,1), S(1,2))=1;
imtool(Iopenseeds);
Idone=Iopenseeds;
%erstelle eine Variable, die immer dann 1 ist, wenn die Region weiter
%wächst, bzw 0 ist, wenn sie nicht wächst
regiongrows=1;
%Grauwertschwelle
T=50;
while regiongrows==1
regiongrows=0;
for ix=2:width-1
for iy=2:height-1
if Iopenseeds (iy,ix)>0
%Untersuche Pixel in 8er-Nachbarschaft
Q=abs (int8(I(iy-1,ix))-int8(I(iy,ix)));
if (Q<T) && Idone(iy-1,ix) ==0
Iopenseeds(iy-1,ix)=1;
Idone(iy-1,ix)=255;
regiongrows=1;
end
Q=abs (int8(I(iy+1,ix))-int8(I(iy,ix)));
if (Q<T) && Idone(iy+1,ix) ==0
Iopenseeds(iy+1,ix)=1;
Idone(iy+1,ix)=255;
regiongrows=1;
end
Q=abs (int8(I(iy,ix-1))-int8(I(iy,ix)));
if (Q<T) && Idone(iy,ix-1) ==0
Iopenseeds(iy,ix-1)=1;
Idone(iy,ix-1)=255;
regiongrows=1;
end
Q=abs (int8(I(iy,ix+1))-int8(I(iy,ix)));
if (Q<T) && Idone(iy,ix+1) ==0
Iopenseeds(iy,ix+1)=1;
Idone(iy,ix+1)=255;
regiongrows=1;
end
Q=abs (int8(I(iy-1,ix+1))-int8(I(iy,ix)));
if (Q<T) && Idone(iy-1,ix+1) ==0
Iopenseeds(iy-1,ix+1)=1;
Idone(iy-1,ix+1)=255;
regiongrows=1;
end
Q=abs (int8(I(iy-1,ix-1))-int8(I(iy,ix)));
if (Q<T) && Idone(iy-1,ix-1) ==0
Iopenseeds(iy-1,ix-1)=1;
Idone(iy-1,ix-1)=255;
regiongrows=1;
end
Q=abs (int8(I(iy+1,ix+1))-int8(I(iy,ix)));
if (Q<T) && Idone(iy+1,ix+1) ==0
Iopenseeds(iy+1,ix+1)=1;
Idone(iy+1,ix+1)=255;
regiongrows=1;
end
Q=abs (int8(I(iy+1,ix-1))-int8(I(iy,ix)));
if (Q<T) && Idone(iy+1,ix-1) ==0
Iopenseeds(iy+1,ix-1)=1;
Idone(iy+1,ix-1)=255;
regiongrows=1;
end
Iopenseeds(iy,ix)=0;
end
end
end
end
imtool(Idone);
%markierten Bereich rot anzeigen
Irot=double(zeros(height,width,farbtiefe));
for iy=1:height
for ix=1:width
gw1=Idone(iy,ix);
if gw1>50
Irot(iy,ix,1)=255;
end
end
end
imtool(Irot);
|
|
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
|
|
Impressum
| Nutzungsbedingungen
| Datenschutz
| FAQ
| RSS
Hosted by:
Copyright © 2007 - 2024
goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
|
|