Saint-Sauvier !

unit Unit1;

interface

uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, ExtCtrls, StdCtrls, Grids;

type
  TForm1 = class(TForm)
    Image1: TImage;
    Button1: TButton;
    Grille: TStringGrid;
    Memo1: TMemo;
    procedure Button1Click(Sender: TObject);
  private

    { Déclarations privées }
    Function getpixel(ax,ay : integer): TColor;
  public
    { Déclarations publiques }
  end;

var

  Form1: TForm1;

implementation

{$R *.dfm}

type
  TRGBArray = ARRAY[0..0] OF TRGBTriple;   // élément de bitmap (API windows)

  pRGBArray = ^TRGBArray;     // type pointeur vers tableau 3 octets 24 bits

var
            // attention, pas de bitmap plus haut que 2048 lignes
              Lignes : array[0..2047] of pRGBArray;


function distance(a,b:TPoint):real; {carré de la distance AB}
begin
     result:=sqr(a.x-b.x)+sqr(a.y-b.y);
end;

//  Accès rapide à un pixel
Function Tform1.getpixel(ax,ay : integer): TColor;

const
   Filler : array[0..2] of byte = (0,0,0);   // noir par défaut

var interm:TRGBTRiple;  R,G,B : integer;
begin
  if (ax >= 0) AND (ay >= 0) AND             // limites

    (ax < Image1.Picture.Bitmap.width-1)  AND (ay < Image1.Picture.Bitmap.height-1) then
  interm := Lignes[ay,ax]
  else

  interm := Trgbtriple(Filler);  // transtypage des 3 bytes

  R := interm.RGBtRed;  
  G := interm.RGBtGreen;
  B := interm.RGBtBlue;
  result:= RGB(R,G,B);
end;

procedure TForm1.Button1Click(Sender: TObject);
var x,y,fx,fy,nbpoints,count:integer; d,sd,sd2,v,ppv:real;
begin

        Screen.Cursor:=crHourGlass;

        Grille.RowCount:=Image1.Picture.Bitmap.height*Image1.Picture.Bitmap.width;
        Image1.Picture.Bitmap.pixelformat := pf24bit;
        // précalcul des scanline
        For y := 0 to Image1.Picture.Bitmap.height-1 do Lignes[y] := Image1.Picture.Bitmap.scanline[y];

        count:=0;
        // pour chaque point de l'image

        for y:=0 to Image1.Picture.Height-1 do
                for x:=0 to Image1.Picture.Width-1 do begin

                    // (inutile si c'est un point de la frontière)
                    if (getpixel(x,y)=clWhite) then begin
                       Caption:=IntToStr(x)+' - '+IntToStr(y);
                       nbpoints:=0;  sd:=0; sd2:=0;
                       // pour chaque point de la frontière :
                       // (pour gagner du temps, on aurait pu chercher les points de la
                       //  frontière et les mettre dans un tableau avant cette boucle)
                       for fy:=0 to Image1.Picture.Height-1 do

                           for fx:=0 to Image1.Picture.Width-1 do
                               if (getpixel(fx,fy)=clBlack) then begin
                                    inc(nbpoints);
                                    d:=distance(Point(x,y),Point(fx,fy));
                                    sd:=sd+d; // somme des distances
                                    sd2:=sd2+d*d; // somme des carrés des distances
                               end;     
                       v:=sd2/nbpoints-(sd/nbpoints)*(sd/nbpoints); // variance
                       // on stocke le point et sa variance dans le tableau
                       Grille.Cells[0,count]:=IntToStr(x);
                       Grille.Cells[1,count]:=IntToStr(y);
                       Grille.Cells[2,count]:=FloatToStr(v);
                       inc(count);
                    end;
                end;
        
        // on cherche la plus petite variance.
        ppv:=StrToFloat(Grille.Cells[2,0]);   
        for y:=0 to count-1 do
             if StrToFloat(Grille.Cells[2,y])<ppv then ppv:=StrToFloat(Grille.Cells[2,y]);
        
        // on énumère tous les points qui ont la plus petite variance
        Memo1.Clear;       
        for y:=0 to count-1 do
                if Grille.Cells[2,y]=FloatToStr(ppv) then begin
                        memo1.Lines.Add(Grille.Cells[0,y]+#9+Grille.Cells[1,y]);
                        Image1.Picture.Bitmap.Canvas.Pixels[StrToInt(Grille.Cells[0,y]),StrToInt(Grille.Cells[1,y])]:=clRed;
                end;
        Screen.Cursor:=crDefault;
end;

end.