Kreis im Viereck bestimmen < Matlab < Mathe-Software < Mathe < Vorhilfe
|
| Status: |
(Frage) beantwortet  | | Datum: | 17:44 Mo 03.12.2007 | | Autor: | chebsam |
| Aufgabe | | [Dateianhang nicht öffentlich] |
Hallo,
mit einem Tool kann ich mir 3d-Volumen in matlab angucken. Habe mit imread('a.jpg') die Datei im Anhang geladen und eine 172x210x3 Matrix erhalten. Mit dem Tool gucke ich mir aber einzelne Schichten an.
Dann kann ich mit der Maus die linke obere und rechte untere Ecke des Kreises (eigentlich Kugel aber wegen der schichtweisen Betrachtung eben Kreis) anklicken und erhalte die Koordinaten der Punkte. Jetzt will ich die Punkte des Kreises speichern und am Ende soll der Mittelpunkt bestimmt werden. Dazu muss ich aber irgendwie die Informationen der Eckpunkte benutzen...wie weiß ich aber nicht.
Pseudocode:
for x_oben bis x_unten
for y_oben bis y_unten
vektor=find(imread('a.jpg')>wert);
Mittelpunkt=[mean(x), mean(y)];
end;
end;
Bin für andere Ideen dankbar.
LG
chebsam
Dateianhänge:
Anhang Nr. 1 (Typ: jpg) [nicht öffentlich]
|
|
| |
|
Hallo,
zuerst muss ich einen Kommentar loswerden:
Ich habe noch nie einen Kreis mit Ecken gesehen.
So, jetzt zum Thema:
Geht es darum, einen Kreis mit der Maus anzuklicken und dann die weißen Pixel darin zu zählen?
Oder sollen die Koordinaten der Mausklicks gemittelt werden?
Soll also nur zweimal geklickt werden oder zweimal am Anfang und dann beliebig oft?
Wozu sind die beiden Klicks? Zum Markieren des Kreises? Reicht da nicht ein Klick?
Deine Beschreibung und dein Pseudocode klingen irgendwie unterschiedlich. Vielleicht kannst du das ja nochmal genauer erklären.
Gruß
Martin
|
|
|
| |
|
| Status: |
(Frage) beantwortet | | Datum: | 19:57 Mo 03.12.2007 | | Autor: | chebsam |
Hallo,
also zunächst muss ich zum Thema sagen, dass es nicht exakt ein Kreis ist, sondern es tritt auch Streuung auf!
Mittlerweile hab ich die z-Koordinate (z=3), bei dem ich die Punktwokle bzw. "Kreis" mithilfe der beiden Klicks bestimmt.
Diese habe ich mir mit
a=imread('a.jpg');
schicht=squeeze(a(:,:,3));
imagesc(schicht);
angeguckt.
Beim ersten Klick lege ich die obere linke Ecke fest und beim zweiten die rehcte untere, d.h. ich muss irgendwie die x und y Koordinaten als Anfangs- bzw. Endpunkte nehmen für meine ROI(region of interest). Danach könnten wir schwellwertsegmentieren(find oder so) und per Mean die Koordinaten des Miitelpunkts ausgeben!?
Ich muss zweimal klicken, weil es mehrere "Kreise" pro Schicht gibt und ich die Koordinaten der Mittelpunkte nach und nach bestimmen möchte.
Ja, richtig...mit den beiden Klicks markiere ich den Kreis.
LG
chebsam
|
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 22:52 Mo 03.12.2007 | | Autor: | chebsam |
Hallo,
habe mir gerade eine Notlösung gebastelt. Problem ist aber, dass der Mittelpunkt nicht exakt ausgegeben wird. Der ist leicht nach links verschoben...hab da keine Erklärung für!?
LG
chebsam
Dateianhänge:
Anhang Nr. 1 (Typ: m) [nicht öffentlich]
|
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig | | Datum: | 08:30 Di 04.12.2007 | | Autor: | Martin243 |
Hallo,
das Problem ist, dass eine quadratische Gleichung in den meisten Fällen zwei Lösungen hat, wenn denn überhaupt welche. Da du aber immer nur die linke auswertest, bekommt dein Schwerpunkt Schlagseite und deine Maxima bilden einen nach links zeigenden Tropfen.
Ich habe das mal behoben (und leicht beschleunigt). Hier die korrigierte Fassung des betreffenden Teils:
%Transformation in den Parameterraum
for i=1:Anzahl_Einsen
a=round([ x(i)-sqrt(radius_quad-(y(i)-(1:c_y)).^2);
x(i)+sqrt(radius_quad-(y(i)-(1:c_y)).^2)]);
a_real = radius >= abs(y(i)-(1:c_y));
j = find(all([a_real; a(1,:)>0]));
houghmatrix(j, a(1,j))=houghmatrix(j,a(1,j))+1;
a_real = radius > abs(y(i)-(1:c_y)); %andere Bedingung wegen doppelter Lösungen
j = find(all([a_real; a(2,:)<=c_x]));
houghmatrix(j, a(2,j))=houghmatrix(j,a(2,j))+1;
end;
|
|
|
| |
|
Hallo,
schau dir doch mal meine ![[a]](/images/paperclip.gif "Dateianhänge") Bastelei an. Hier warte ich auf zwei Mausklicks, in deren Mittelpunkt iregendein Punkt innerhalb aus dem Kreis liegen muss. Danach breitet sich mein Template aus, bis der mit coverage angegebene Anteil eines Kreisinhalts überdeckt ist. Damit legst du in etwa fest, wie stark Ausreißer mit einbezogen werden.
Sobald der Kreis komplett überdeckt wurde (bis dahin einfach immer wieder RETURN in der Konsole drücken), kann man sich das aktuelle Template nehmen und dem entsprechenden Bildausschnitt elementweise multiplizieren. Dann kann man den Mittelwert der Koordinaten der Punkte mit 1 berechnen.
Hier ist das etwas modifizierte Bild:
[Dateianhang nicht öffentlich]
Hier die Überdeckung mit coverage=1:
[Dateianhang nicht öffentlich]
Hier die Überdeckung mit coverage=0.8;
[Dateianhang nicht öffentlich]
Vielleicht kannst du damit etwas anfangen...
Gruß
Martin
Dateianhänge:
Anhang Nr. 1 (Typ: m) [nicht öffentlich]
Anhang Nr. 2 (Typ: jpg) [nicht öffentlich]
Anhang Nr. 3 (Typ: jpg) [nicht öffentlich]
Anhang Nr. 4 (Typ: jpg) [nicht öffentlich]
|
|
|
|
