compare center of mass and gaussian fit.ijm

compare center of mass and gaussian fit.ijm

// see https://twitter.com/Maurice_Y_Lee/status/1487093624982499330

// show center of mass
makePoint(getValue('XM'),getValue('YM'),"medium green hybrid");
Overlay.addSelection('green');

// gaussian fit
makePoint(gaussianFit(true),gaussianFit(false),"medium red hybrid");
Overlay.addSelection('red');

run("Select None");

function gaussianFit(horizontal) {
  run("Select All");
  if (!horizontal) setKeyDown('alt');
  k=getProfile();
  Fit.doFit("Gaussian",Array.getSequence(k.length), k);
  return Fit.p(2);
}