Ускорить программу, которая выполняет преобразование Фурье - Delphi - Ответ 15667645

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, ExtCtrls, Jpeg;
 
  type
 Tinf = record
    Im  :real;
    Re  :real;
    modul:real;
  end;
 ImRemage = array of array of Tinf;
 
type
  TForm1 = class(TForm)
    GroupBox1: TGroupBox;
    Image1: TImage;
    GroupBox2: TGroupBox;
    GroupBox3: TGroupBox;
    Image2: TImage;
    Image3: TImage;
    Button1: TButton;
    Button2: TButton;
    Button3: TButton;
    Button4: TButton;
    OpenDialog1: TOpenDialog;
    SaveDialog1: TSaveDialog;
    procedure Button1Click(Sender: TObject);
    procedure Button2Click(Sender: TObject);
    procedure Button3Click(Sender: TObject);
    procedure Button4Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
var
  Form1: TForm1;
  bmp1,bmp2         : TBitMap ;
  m,n               : integer;
  IRmage,IRmage2    : ImRemage;
  jpg               : TJpegImage;
 
implementation
 
{$R *.dfm}
 
procedure TForm1.Button1Click(Sender: TObject);
begin
 if opendialog1.Execute() then
  begin
     image1.Picture.LoadFromFile(opendialog1.FileName);
     bmp1:=TBitmap.Create;
     bmp2:=TBitmap.Create;
 
     Jpg := TJPEGImage.Create;
     Jpg.LoadFromFile(opendialog1.FileName);
     Bmp1.Assign(jpg);
     Bmp2.Assign(jpg);
     setlength(irmage,jpg.Height, jpg.Width);
       setlength(irmage2,jpg.Height, jpg.Width);
  end;
end;
 
procedure TForm1.Button2Click(Sender: TObject);
var
q,s,i,j           : integer;
v,u,x,y,f         : integer;
g                 : boolean;
scan              : pByteArray;
img               : Timage;
md                : real;
re,im             : real;
Pixel             : TColor;
begin
 
  M := bmp1.Height;
  N := bmp1.Width;
 
    begin
    for v := 0 to M - 1 do begin
       for u := 0 to N - 1 do begin
          md := 0;  Re := 0;  Im := 0;
          irmage[v,u].Re := 0;
          irmage[v,u].Im := 0;
          for y := 0 to M - 1 do begin
            scan := bmp1.ScanLine[y];
            for x := 0 to N - 1 do begin
              f:= 0;
              f := {round}(scan[x]);
              if odd(x+y) then f := f*(-1);
              Re := Re + cos((-1)*2*pi*(u*x/M+v*y/N))*f;
              Im := Im + sin((-1)*2*pi*(u*x/M+v*y/N))*f;
            end;
          end;
          irmage[v,u].Re := Re/(M*N);
          irmage[v,u].Im := Im/(M*N);
          irmage2[v,u].Re := Re/(M*N);
          irmage2[v,u].Im := Im/(M*N);
          md := sqrt( ((Re*Re)/(m*n)) + ((Im*Im)/(m*n)) );
          irmage[v,u].Modul := md;
          i := round (20 * ln(1+md));
          if i >= 255 then i := 255;
scan:=bmp1.ScanLine[x];         //<-----Этот момент!
          scan[y]:=f;
                                        
       end;
    end;
    image2.Picture.Bitmap.Assign(bmp2);
    bmp1.Free;
    bmp2.Free;
    end;
end;