run("Blobs (25K)");
  setAutoThreshold("Default");
  run("Set Measurements...", "area mean feret's decimal=4");
  run("Analyze Particles...", "  show=[Overlay Outlines] display exclude clear");
  resetThreshold();
  Overlay.drawLabels(false);
  for (i=0; i<nResults; i++) {
     x1 = getResult("FeretX",i);
     y1 = getResult("FeretY",i);
     length = getResult("Feret",i);
     degrees = getResult("FeretAngle",i);
     if (degrees>90)
        degrees -= 180; 
     angle = degrees*PI/180;
     x2 = x1 + cos(angle)*length;
     y2 = y1 - sin(angle)*length;
     makeLine(x1, y1, x2, y2);
     Overlay.addSelection("red");
  }