//////////////////////////////////
  img2m;  //convert to data matrix
  string mat$=%H;
  m2w method:=xyz;  //convert to a worksheet
  string wkb$=%H;
  sum(col(3));  //Get stats
  minr=sum.imin;  //imin for darkest, and imax for brightest
  minx=col(1)[minr];
  miny=col(2)[minr];
  minz=col(3)[minr];
  meanz=sum.mean;
  sdz=sum.sd;
  maxz=sum.max;
  minz=sum.min;
  win -t plot image;  //make an empty image plot
  string gr$=%H;
  layer -i %(mat$);
  layer -a;
  speedmode;
  window -a %(wkb$);
  newsheet;
  col(1)[1]=minx;
  col(2)[1]=miny;
  plotxy iy:=(1,2) plot:=201 ogl:=[gr$]1!1;  //plot a marker in the image plot
  window -a %(gr$);
  range ss=[%(wkb$)]2!2;
  set ss -c 2;  //RED
  set ss -k 2;  //Circle
  set ss -kf 4; //Center Plus
  type -a Z value at point is $(minz) at ( $(minx) , $(miny) );
  type -a Region Statistics: Mean: $(meanz), SD: $(sdz), Max: $(maxz), Min: $(minz);
//////////////////////////////////