Hough transformation for Kreise

[bernhard_LA]

sieht jemand einen Bug in dieser Version der Hough Transformation für Kreise ?
Das thema Hough und Lines gabe es bereits hier : http://www.delphipraxis.net/156893-h...formation.html

Delphi-Quellcode:

///
///   Hough transformation for line detection
///
///
///   AnalysisBitmap : TBitMap;      ->  the image  for hough tranmsformation
///   aHoughResult :  THoughResult   ->  the result of the Hough transformation
///   r  :  Integer;                 ->  the search radius
///
///


procedure Hough_CircleDetection ( AnalysisBitmap : TBitMap; var aHoughResult :  THoughResult ;  r  :  Integer  );
var x,y, a,b    : integer;
    ImageWidth  : integer;
    ImageHeight : Integer;
    max_d       : Integer;
    help        : Real;

begin


   ///  size of hough array
   ImageWidth := AnalysisBitmap.Width;
   ImageHeight:= AnalysisBitmap.Height;



                            //  a        //  b
   SetLength(aHoughResult,ImageWidth, ImageHeight );

   // For all rows in image:
   for y:=0 to AnalysisBitmap.Height-1 do
   begin

   // For all pixel in one row :
   for x:=0 to AnalysisBitmap.Width-1 do
   begin

      // Is there a point  ?
      if IsPixel(x,y, AnalysisBitmap, 128 ) then
      begin
           //   circle
           //   r *r =  ( x-a) *(x-a)  + (y-b) *( y-b)
           for a:=0 to ImageWidth do
           begin
                // r * r   := (x-a) * (x-a)  + (y-b) * (y-b);
                // solve equation and get y ...
                help :=   r * r - (x-a) * (x-a);

                if ( help >= 0 ) then
                    begin

                    b :=  Round(  sqrt (help) + y ) ;

                    Inc(aHoughResult [a, b]);

                    end;



           end;
      end;
   end;
  end;


end;

[Corpsman]

Aus welcher Quelle hast du das denn umgesetzt ?

Die Houghtransformation so aus dem Kopf wird wohl keiner können ..

[bernhard_LA]

Kreise analog zu Linien behandelt, habe keinen Source Code incl. Akkumulator gefunden - vielleicht hat jemand schon mal das ganze in Delphi implementiert

die Kreisgleichung : r *r = (x-a) *(x-a) + (y-b) *(y-b)

über alles (x,y) falls ein Pixel gefunden ->
gibt es einen Kreis mit Mittelpunlkt (a,b) und Radius r, fix der die Gleichung löst ....


PS: den Ganzen Code habe ich unter Sourge forge eingestellt

[Corpsman]

Ich meinte nicht irgend einen Source,

Sondern eher die Mathematische Abhandlung darüber, dass es bei den Kreisen Prinzipiell das selbe ist wie bei den Geraden war mir bewusst.

Nun hast du mich auf jeden Fall neugierig genug gemacht, dass ich mich auch mal daran versuche, ich werde dann berichten..

[bernhard_LA]

unter http://users.ecs.soton.ac.uk/msn/boo.../houghCircles/

gibt eine Lösung in Java mit einem 1D acc -Array

Math. wurde hier x = a + r * sin q und y=b + y * sin q gelöst,
eigentlich sollte mein Code das gleiche erledigen nur habe ich mich von sin q und cos q befreit....

[Corpsman]

So, ich habe mal das Java Ding nach Lazarus "portiert", da ich ja kein Delphi mehr habe.

Ich habe aber extra keine LCL spezifischen Sachen benutzt.

Du müsstest relativ einfach mein "gewerksel" nach Delphi portieren können.

Im Zip ebenfalls die Ausgabe des Akkumulators für "basic.bmp"

Die ganzen 1D Sachen sind nur weil Java ( wie alle C ähnlichen Dinge ) so bescheiden mit 2D-Arrays umgeht. Die Portierung in 2D-Arrays sollte aber kein Problem für dich darstellen.