HOUGH : Punkte aus dyn. Array löschen, Lokale Maxima finden

[bernhard_LA]

unter http://sourceforge.net/projects/houghtransforma/ gibt es ein Demoprojekt um
Kreise und Linien in einem Bild zu erkennen.

Jeder Kreismittelpunkt wird als lokales Maximum einer Funktion Z=f(x,y) dargestellt. ( aHoughResult : THoughResult)
Wenn ich an einer Position ein Maximum gefunden habe kann im Abstand R um diesen Punkt kein weiterer Kreis liegen. Das kleinere Maximum innerhalb meines Radius R um einen Maxima- Punkt entsteht zB. durch Rauschen im Bild oder anderes Effekte.

Ich muss also in einem dyn. array alle Punkte überprüfen ob sie näher als ein gegebener Abstands-Wert MinPeakDistance beisammenliegen , den Punkt mit dem größten Z Wert beihalten und alle Punkte mit kleineren Z Werten dann löschen.

Sieht jemand warum der Code unten nicht dieses Problem löst ?

Delphi-Quellcode:

///
///   exchange two array elements
///
procedure ExchangeArrayElements (var AHoughResultsArray: THoughFinal;
  const AIndex, BIndex: integer);
var   aPoint3D : THPoint3D;
begin
      aPoint3D := AHoughResultsArray[AIndex];

      AHoughResultsArray[AIndex] := AHoughResultsArray[BIndex];

      AHoughResultsArray[BIndex] := aPoint3D;
end;

///
/// aHoughFinal -> T3DPoint (x,y,z)
///      algo: Für alle Array Elemente ....
///              solange Entfernung zwischen (P1-P2) < Min. Entferung
///                     lösche Punkt mit kleinerem z Wert
///
procedure HoughResultToParameterDynamic(aHoughResult: THoughResult; Range: Real;
  var aHoughFinal: THoughFinal; MinPeakDistance: integer; var deleteElementCount : Integer);  overload;

var
  len_Hough: integer;
  i, j, iter: integer;
  Distance: integer;
  rescan_from_i_position : Boolean;
begin

  deleteElementCount := 0;

  /// extract from Hough Accu results only the peaks
  HoughResultToParameterDynamic(aHoughResult, Range, aHoughFinal);

  iter := 0;

  ///  lösche nahe beieinander liegende Maxima Peaks
  repeat

  i := Iter;

    repeat

    inc(i);

    rescan_from_i_position := false;

      /// wie weit sind 2 Punkte entfernt ?
      Distance := abs(aHoughFinal[iter].X - aHoughFinal[i].X) +
        abs(aHoughFinal[iter].Y - aHoughFinal[i].Y);

      ///  zunahe ??  -> dann lösche Punkt mit kleinerem Z Wert
      if (Distance <  MinPeakDistance) then
        if (aHoughFinal[iter].z > aHoughFinal[i].z) then
        begin
          DeleteArrayElement(aHoughFinal, i);

          inc(deleteElementCount);

        end
        else
        begin
          rescan_from_i_position:=true;

          ExchangeArrayElements(aHoughFinal, i, iter);

          DeleteArrayElement(aHoughFinal, i);

          inc(deleteElementCount);
        end;

      len_Hough := length(aHoughFinal);

    until ((i > len_Hough - 1) or (rescan_from_i_position=true));

    if not rescan_from_i_position then Inc(iter);

  until (iter > len_Hough - 2);

end;

[gammatester]

Nur eine Vermutung: Du redest zwar von Kreisen und Radien, benutzt aber die L1-Norm statt der Euklidischen. Wenn Du keine Wurzeln ziehen willst, kannst Du mit den Euklidischen Abstandsquadraten arbeiten. Ob das allerdings Dein Problem löst, weiß ich nicht.

Gruß Gammatester

[bernhard_LA]

wir haben ein paar hundert tausend Kreise - Euklid spart hat die Wurzelberechnung geht schnell und ist hinreichend genau für uns

[gammatester]

Häh? Du benutzt doch wie schon bemerkt überhaupt keine Euklidische Norm. Beispiel:

Code:

1. Punkt (4,4): Deine L1-Norm = 8, L2-Norm (Euklid) = sqrt(32) = 5.66
2. Punkt (1,6): Deine L1-Norm = 7, L2-Norm (Euklid) = sqrt(37) = 6.08

Das heißt, in Deiner Metrik ist der erste Punkt weiter vom Ursprung entfernt als der zweite, euklidisch ist es umgekehrt.

[bernhard_LA]

ja, was wir machen nennt sich doch http://en.wiktionary.org/wiki/Manhattan_distance , aber das Thema wie wir den Abstand berechne ist nicht das Problem ,
ich habe Probeme struktueriert alle Punkte mit Nebenmaxima aus dem Array zu löschen ....

[gammatester]

Vielleicht solltest Du dann mal verraten, was nicht funktioniert. Meine Bemerkung bezog sich darauf, das eventuell falsche Punkte gelöscht werden, weil die Abstandsberechnung nicht angemessen ist.

[Amateurprofi]

@bernhard_LA:
Du rufst in der Prozedur "HoughResultToParameterDynamic" eine überladene Prozedur von "HoughResultToParameterDynamic" und eine Prozedur "DeleteArrayElement" auf.

Ich würde ja gern versuchen dir zu helfen, jedoch ohne die Source Codes dieser Prozeduren ist das nicht möglich, denn das Problem könnte in einer dieser Prozeduren liegen.