var {Global variables, initially zero}
  pid1,pid2:integer;

procedure DoGammaTransform(gamma:real);
var
  i,v:integer;
  n,mode,min,max:integer
  mean:real;
BEGIN
  measure;
  GetResults(n,mean,mode,min,max);
  ShowMessage('min=',min:1,'\max=',max:1);
  for i:=1 to 254 DO begin
    if (i>min) and (i<max)
      then v:=exp(gamma*ln((i-min)/(max-min)))*255 {x^y=exp(y*ln(x)}
      else begin
        if i<=min then v:=0 else v:=255;
      end;
    RedLUT[i]:=255-v;
    GreenLUT[i]:=255-v;
    BlueLUT[i]:=255-v;
  end;
  UpdateLUT;
END.


procedure MarkPeaks(PeakRadius:real);
var
  gamma,Peakedness:real;
  i,width,height,threshold:integer;
begin
  SetforegroundColor(0);
  SetbackgroundColor(255);
  gamma:=1.6;
  {Invert;}
  SubtractBackground('2D Rolling Ball',3);
  GetPicSize(width,height);
  MakeRoi(0.05*width,0.05*height,0.9*width,0.9*height);
  EnhanceContrast;
  KillRoi;
  ApplyLUT;
  DoGammaTransform(gamma);
  ApplyLUT;
  Smooth;
  AutoThreshold;
  i:=0;
  Repeat
    i:=i+1;
  until RedLUT[i]<255;
  threshold:=round(0.5*i);
  SetThreshold(-1);
  peakedness:=0.4;  {0.05-0.95}
  UserCode('Mark Peaks', threshold, PeakRadius, Peakedness); 
end;


macro 'Mark Grains [M]';
var
  radius:real;
begin
  RequiresVersion(1.53);
  if nPics=999 then begin
    PutMessage('Exactly one window must be open.');
    exit;
  end;
  KillRoi;
  pid1:=PidNumber;
  Duplicate(concat(WindowTitle,'-2'));
  pid2:=PidNumber;
   radius:=3;
  {radius:=GetNumber('Peak Radius[1-50]:',radius);}
   Invert;
  MarkPeaks(radius);
  Invert;
  ChangeValues(254, 254, 253);
  ChangeValues(255, 255, 254);
  RedLUT[254]:=255;
  GreenLUT[254]:=255;
  BlueLUT[254]:=0;
  UpdateLUT;
  ResetCounter;
end;

macro 'Count Marked Grains in ROI [C]';
var
  t1,t2:integer;
begin
  if pid1=0 then begin
    PutMessage('You must first use the Count Grains macro.');
    exit;
  end;
  SetOptions('Area, Mean, X-Y Center, User1');
  SetUser1Label('Count');
  if PidNumber<>pid1 then begin
    SelectPic(pid1);
    RestoreRoi;
  end;
  GetThresholds(t1,t2);
  Measure;
  SelectPic(pid2);
  RestoreRoi;
  Measure;
  SetCounter(rCount-1);
  rUser1[rCount]:=histogram[254];
  UpdateResults;
  SelectPic(pid1);
end;