program TestColorTransform;
{$APPTYPE CONSOLE}
uses
  Windows, SysUtils, math;

const
  MaxLen=100000000;

type
  TTransTab= array[byte] of byte;
  TDataTab = array[0..MaxLen-1] of byte;

  procedure CreateTransformTab(const Brighness, contrast, gamma: single;
                                 out PixTab, NoizeTab: TTransTab);
  var
    i: byte;
    Val: single;
  begin
    for i:=0 to maxByte do begin
      val:= MaxByte*(Power(i/MaxByte, gamma)+Brighness)*contrast;
      if Val<0
        then PixTab[i]:=0
      else if Val>MaxByte
        then PixTab[i]:=MaxByte
      else PixTab[i]:=Trunc(Val);
    end;

    for i:=0 to maxByte-1 do
      NoizeTab[i]:=pixTab[i+1]-pixTab[i]+1;
    NoizeTab[maxByte]:=1;
  end;

  procedure InitData(out DataTab: TDataTab);
  var
    i: integer;
  begin
    for i:=0 to High(DataTab) do DataTab[i]:=random(256);
  end;


  procedure ProcDataTabNoNoize(var DataTab: TDataTab; const TransTab: TTransTab);
  var
    i: integer;
  begin
    for i:=0 to High(DataTab) do begin
      DataTab[i]:=TransTab[DataTab[i]];
    end;
  end;

  procedure ProcDataTabRegNoize(var DataTab: TDataTab; const TransTab, NoizeTab: TTransTab);
  var
    i: integer;
  begin
    for i:=0 to High(DataTab) do begin
      DataTab[i]:=TransTab[DataTab[i]]+I mod NoizeTab[DataTab[i]];
    end;
  end;

  procedure ProcDataTabRandNoize(var DataTab: TDataTab; const TransTab, NoizeTab: TTransTab);
  var
    i: integer;
   begin
    for i:=0 to High(DataTab) do begin
      DataTab[i]:=TransTab[DataTab[i]]+ Random(NoizeTab[DataTab[i]]);
    end;
  end;



//собственно тест
var
  PixTab, NoizeTab: TTransTab;
  TC: cardinal;
  Data: TDataTab;

begin
  CreateTransformTab(-0.05, 1.1, 0.5, PixTab, NoizeTab);

  InitData(Data);
  TC:=gettickcount;
  ProcDataTabNoNoize(Data, PixTab);
  Writeln('No noize: ', Gettickcount-TC);

  InitData(Data);
  TC:=gettickcount;
  ProcDataTabRegNoize(Data, PixTab, NoizeTab);
  Writeln('Regular noize: ', Gettickcount-TC);

  InitData(Data);
  TC:=gettickcount;
  ProcDataTabRandNoize(Data, PixTab, NoizeTab);
  Writeln('Random noize: ', Gettickcount-TC);
  Readln;
end.

const
  BaseShift=8;
type
  TRandTab = array[0..MaxLen-1] of byte; //значения в диапазоне 0.. 1 shl (BaseShift-1)

  procedure ProcDataTabNoize_LK(var DataTab: TDataTab; const TransTab, NoizeTab: TTransTab; const RandTab: TRandTab);
  var
    i: integer;
   begin
    for i:=0 to High(DataTab) do begin
      DataTab[i]:=TransTab[DataTab[i]]+  randtab[i]*NoizeTab[DataTab[i]] shr BaseShift;
    end;
  end;