void loadTransformation(String filename, ResultsTable res) {
try {
String line;
FileReader fr = new FileReader(filename);
BufferedReader br = new BufferedReader(fr);
if (!br.readLine().equals(" Z-Step Raw Width minus Heigh Calibration Width minus Height")) {
IJ.error("File does not seam to be an Astigmatism calibration file");
return;
}
// java.lang.String [] elements = new java.lang.String [3];
java.lang.String[] elements;
int counter = 1;
res.reset();
while ((line = br.readLine()) != null) {
IJ.showStatus("Loading element " + counter + "... sit back and relax.");
counter++;
line.trim();
elements = line.split("\t");
res.incrementCounter();
res.addValue("Z-Step", Double.parseDouble(elements[1]));
res.addValue("Raw Width minus Heigh", Double.parseDouble(elements[2]));
res.addValue("Calibration Width minus Height", Double.parseDouble(elements[3]));
}
fr.close();
} catch (FileNotFoundException e) {
IJ.error("File not found exception" + e);
return;
} catch (IOException e) {
IJ.error("IOException exception" + e);
return;
} catch (NumberFormatException e) {
IJ.error("Number format exception" + e);
return;
}
}
/** Called by ImageJ when the user selects Quit. */
public void quit() {
quitMacro = IJ.macroRunning();
Thread thread = new Thread(this, "Quit");
thread.setPriority(Thread.NORM_PRIORITY);
thread.start();
IJ.wait(10);
}
void createEllipse(ImagePlus imp) {
IJ.showStatus("Fitting ellipse");
Roi roi = imp.getRoi();
if (roi == null) {
noRoi("Fit Ellipse");
return;
}
if (roi.isLine()) {
IJ.error("Fit Ellipse", "\"Fit Ellipse\" does not work with line selections");
return;
}
ImageProcessor ip = imp.getProcessor();
ip.setRoi(roi);
int options = Measurements.CENTROID + Measurements.ELLIPSE;
ImageStatistics stats = ImageStatistics.getStatistics(ip, options, null);
double dx = stats.major * Math.cos(stats.angle / 180.0 * Math.PI) / 2.0;
double dy = -stats.major * Math.sin(stats.angle / 180.0 * Math.PI) / 2.0;
double x1 = stats.xCentroid - dx;
double x2 = stats.xCentroid + dx;
double y1 = stats.yCentroid - dy;
double y2 = stats.yCentroid + dy;
double aspectRatio = stats.minor / stats.major;
imp.killRoi();
imp.setRoi(new EllipseRoi(x1, y1, x2, y2, aspectRatio));
}
/** Handle menu events. */
public void actionPerformed(ActionEvent e) {
if ((e.getSource() instanceof MenuItem)) {
MenuItem item = (MenuItem) e.getSource();
String cmd = e.getActionCommand();
commandName = cmd;
ImagePlus imp = null;
if (item.getParent() == Menus.getOpenRecentMenu()) {
new RecentOpener(cmd); // open image in separate thread
return;
} else if (item.getParent() == Menus.getPopupMenu()) {
Object parent = Menus.getPopupMenu().getParent();
if (parent instanceof ImageCanvas) imp = ((ImageCanvas) parent).getImage();
}
int flags = e.getModifiers();
hotkey = false;
actionPerformedTime = System.currentTimeMillis();
long ellapsedTime = actionPerformedTime - keyPressedTime;
if (cmd != null && (ellapsedTime >= 200L || !cmd.equals(lastKeyCommand))) {
if ((flags & Event.ALT_MASK) != 0) IJ.setKeyDown(KeyEvent.VK_ALT);
if ((flags & Event.SHIFT_MASK) != 0) IJ.setKeyDown(KeyEvent.VK_SHIFT);
new Executer(cmd, imp);
}
lastKeyCommand = null;
if (IJ.debugMode) IJ.log("actionPerformed: time=" + ellapsedTime + ", " + e);
}
}
void interpolate() {
Roi roi = imp.getRoi();
if (roi == null) {
noRoi("Interpolate");
return;
}
if (roi.getType() == Roi.POINT) return;
if (IJ.isMacro() && Macro.getOptions() == null) Macro.setOptions("interval=1");
GenericDialog gd = new GenericDialog("Interpolate");
gd.addNumericField("Interval:", 1.0, 1, 4, "pixel");
gd.addCheckbox("Smooth", IJ.isMacro() ? false : smooth);
gd.showDialog();
if (gd.wasCanceled()) return;
double interval = gd.getNextNumber();
smooth = gd.getNextBoolean();
Undo.setup(Undo.ROI, imp);
FloatPolygon poly = roi.getInterpolatedPolygon(interval, smooth);
int t = roi.getType();
int type = roi.isLine() ? Roi.FREELINE : Roi.FREEROI;
if (t == Roi.POLYGON && interval > 1.0) type = Roi.POLYGON;
if ((t == Roi.RECTANGLE || t == Roi.OVAL || t == Roi.FREEROI) && interval >= 5.0)
type = Roi.POLYGON;
if ((t == Roi.LINE || t == Roi.FREELINE) && interval >= 5.0) type = Roi.POLYLINE;
if (t == Roi.POLYLINE && interval >= 1.0) type = Roi.POLYLINE;
ImageCanvas ic = imp.getCanvas();
if (poly.npoints <= 150 && ic != null && ic.getMagnification() >= 12.0)
type = roi.isLine() ? Roi.POLYLINE : Roi.POLYGON;
Roi p = new PolygonRoi(poly, type);
if (roi.getStroke() != null) p.setStrokeWidth(roi.getStrokeWidth());
p.setStrokeColor(roi.getStrokeColor());
p.setName(roi.getName());
transferProperties(roi, p);
imp.setRoi(p);
}
void runMacro(String arg) {
Roi roi = imp.getRoi();
if (IJ.macroRunning()) {
String options = Macro.getOptions();
if (options != null
&& (options.indexOf("grid=") != -1 || options.indexOf("interpolat") != -1)) {
IJ.run("Rotate... ", options); // run Image>Transform>Rotate
return;
}
}
if (roi == null) {
noRoi("Rotate>Selection");
return;
}
roi = (Roi) roi.clone();
if (arg.equals("rotate")) {
double d = Tools.parseDouble(angle);
if (Double.isNaN(d)) angle = "15";
String value = IJ.runMacroFile("ij.jar:RotateSelection", angle);
if (value != null) angle = value;
} else if (arg.equals("enlarge")) {
String value = IJ.runMacroFile("ij.jar:EnlargeSelection", enlarge);
if (value != null) enlarge = value;
Roi.previousRoi = roi;
}
}
public boolean beadCalibration3d() {
imp = IJ.getImage();
if (imp == null) {
IJ.noImage();
return false;
} else if (imp.getStackSize() == 1) {
IJ.error("Stack required");
return false;
} else if (imp.getType() != ImagePlus.GRAY8 && imp.getType() != ImagePlus.GRAY16) {
// In order to support 32bit images, pict[] must be changed to float[], and getPixel(x, y);
// requires a Float.intBitsToFloat() conversion
IJ.error("8 or 16 bit greyscale image required");
return false;
}
width = imp.getWidth();
height = imp.getHeight();
nslices = imp.getStackSize();
imtitle = imp.getTitle();
models[0] = "*None*";
models[1] = "line";
models[2] = "2nd degree polynomial";
models[3] = "3rd degree polynomial";
models[4] = "4th degree polynomial";
GenericDialog gd = new GenericDialog("3D PALM calibration");
gd.addNumericField("Maximum FWHM (in px)", prefs.get("QuickPALM.3Dcal_fwhm", 20), 0);
gd.addNumericField(
"Particle local threshold (% maximum intensity)", prefs.get("QuickPALM.pthrsh", 20), 0);
gd.addNumericField("Z-spacing (nm)", prefs.get("QuickPALM.z-step", 10), 2);
gd.addNumericField("Calibration Z-smoothing (radius)", prefs.get("QuickPALM.window", 1), 0);
gd.addChoice("Model", models, prefs.get("QuickPALM.model", models[3]));
gd.addCheckbox(
"Show divergence of bead positions against model",
prefs.get("QuickPALM.3Dcal_showDivergence", false));
gd.addCheckbox("Show extra particle info", prefs.get("QuickPALM.3Dcal_showExtraInfo", false));
gd.addMessage("\n\nDon't forget to save the table in the end...");
gd.showDialog();
if (gd.wasCanceled()) return false;
fwhm = gd.getNextNumber();
prefs.set("QuickPALM.QuickPALM.3Dcal_fwhm", fwhm);
pthrsh = gd.getNextNumber() / 100;
prefs.set("QuickPALM.pthrsh", pthrsh * 100);
cal_z = gd.getNextNumber();
prefs.set("QuickPALM.z-step", cal_z);
window = (int) gd.getNextNumber();
prefs.set("QuickPALM.window", window);
model = gd.getNextChoice();
prefs.set("QuickPALM.model", model);
part_divergence = gd.getNextBoolean();
prefs.set("QuickPALM.3Dcal_showDivergence", part_divergence);
part_extrainfo = gd.getNextBoolean();
prefs.set("QuickPALM.3Dcal_showExtraInfo", part_extrainfo);
return true;
}
/** Performs actual projection using specified method. */
public void doProjection() {
if (imp == null) return;
sliceCount = 0;
if (method < AVG_METHOD || method > MEDIAN_METHOD) method = AVG_METHOD;
for (int slice = startSlice; slice <= stopSlice; slice += increment) sliceCount++;
if (method == MEDIAN_METHOD) {
projImage = doMedianProjection();
return;
}
// Create new float processor for projected pixels.
FloatProcessor fp = new FloatProcessor(imp.getWidth(), imp.getHeight());
ImageStack stack = imp.getStack();
RayFunction rayFunc = getRayFunction(method, fp);
if (IJ.debugMode == true) {
IJ.log("\nProjecting stack from: " + startSlice + " to: " + stopSlice);
}
// Determine type of input image. Explicit determination of
// processor type is required for subsequent pixel
// manipulation. This approach is more efficient than the
// more general use of ImageProcessor's getPixelValue and
// putPixel methods.
int ptype;
if (stack.getProcessor(1) instanceof ByteProcessor) ptype = BYTE_TYPE;
else if (stack.getProcessor(1) instanceof ShortProcessor) ptype = SHORT_TYPE;
else if (stack.getProcessor(1) instanceof FloatProcessor) ptype = FLOAT_TYPE;
else {
IJ.error("Z Project", "Non-RGB stack required");
return;
}
// Do the projection.
for (int n = startSlice; n <= stopSlice; n += increment) {
IJ.showStatus("ZProjection " + color + ": " + n + "/" + stopSlice);
IJ.showProgress(n - startSlice, stopSlice - startSlice);
projectSlice(stack.getPixels(n), rayFunc, ptype);
}
// Finish up projection.
if (method == SUM_METHOD) {
fp.resetMinAndMax();
projImage = new ImagePlus(makeTitle(), fp);
} else if (method == SD_METHOD) {
rayFunc.postProcess();
fp.resetMinAndMax();
projImage = new ImagePlus(makeTitle(), fp);
} else {
rayFunc.postProcess();
projImage = makeOutputImage(imp, fp, ptype);
}
if (projImage == null) IJ.error("Z Project", "Error computing projection.");
}
boolean setProperties(String title, Roi roi) {
Frame f = WindowManager.getFrontWindow();
if (f != null && f.getTitle().indexOf("3D Viewer") != -1) return false;
if (roi == null) {
IJ.error("This command requires a selection.");
return false;
}
RoiProperties rp = new RoiProperties(title, roi);
boolean ok = rp.showDialog();
if (IJ.debugMode) IJ.log(roi.getDebugInfo());
return ok;
}
private void rotate(ImagePlus imp) {
Roi roi = imp.getRoi();
if (IJ.macroRunning()) {
String options = Macro.getOptions();
if (options != null
&& (options.indexOf("grid=") != -1 || options.indexOf("interpolat") != -1)) {
IJ.run("Rotate... ", options); // run Image>Transform>Rotate
return;
}
}
(new RoiRotator()).run("");
}
/** Quit using a separate thread, hopefully avoiding thread deadlocks. */
public void run() {
quitting = true;
boolean changes = false;
int[] wList = WindowManager.getIDList();
if (wList != null) {
for (int i = 0; i < wList.length; i++) {
ImagePlus imp = WindowManager.getImage(wList[i]);
if (imp != null && imp.changes == true) {
changes = true;
break;
}
}
}
Frame[] frames = WindowManager.getNonImageWindows();
if (frames != null) {
for (int i = 0; i < frames.length; i++) {
if (frames[i] != null && (frames[i] instanceof Editor)) {
if (((Editor) frames[i]).fileChanged()) {
changes = true;
break;
}
}
}
}
if (windowClosed
&& !changes
&& Menus.window.getItemCount() > Menus.WINDOW_MENU_ITEMS
&& !(IJ.macroRunning() && WindowManager.getImageCount() == 0)) {
GenericDialog gd = new GenericDialog("ImageJ", this);
gd.addMessage("Are you sure you want to quit ImageJ?");
gd.showDialog();
quitting = !gd.wasCanceled();
windowClosed = false;
}
if (!quitting) return;
if (!WindowManager.closeAllWindows()) {
quitting = false;
return;
}
// IJ.log("savePreferences");
if (applet == null) {
saveWindowLocations();
Prefs.savePreferences();
}
IJ.cleanup();
// setVisible(false);
// IJ.log("dispose");
dispose();
if (exitWhenQuitting) System.exit(0);
}
void createMaskFromThreshold(ImagePlus imp) {
ImageProcessor ip = imp.getProcessor();
if (ip.getMinThreshold() == ImageProcessor.NO_THRESHOLD) {
IJ.error("Create Mask", "Area selection or thresholded image required");
return;
}
double t1 = ip.getMinThreshold();
double t2 = ip.getMaxThreshold();
IJ.run("Duplicate...", "title=mask");
ImagePlus imp2 = WindowManager.getCurrentImage();
ImageProcessor ip2 = imp2.getProcessor();
ip2.setThreshold(t1, t2, ImageProcessor.NO_LUT_UPDATE);
IJ.run("Convert to Mask");
}
void fitSpline() {
Roi roi = imp.getRoi();
if (roi == null) {
noRoi("Spline");
return;
}
int type = roi.getType();
boolean segmentedSelection = type == Roi.POLYGON || type == Roi.POLYLINE;
if (!(segmentedSelection
|| type == Roi.FREEROI
|| type == Roi.TRACED_ROI
|| type == Roi.FREELINE)) {
IJ.error("Spline", "Polygon or polyline selection required");
return;
}
if (roi instanceof EllipseRoi) return;
PolygonRoi p = (PolygonRoi) roi;
if (!segmentedSelection) {
if (p.subPixelResolution()) p = trimFloatPolygon(p, p.getUncalibratedLength());
else p = trimPolygon(p, p.getUncalibratedLength());
}
String options = Macro.getOptions();
if (options != null && options.indexOf("straighten") != -1) p.fitSplineForStraightening();
else if (options != null && options.indexOf("remove") != -1) p.removeSplineFit();
else p.fitSpline();
imp.draw();
LineWidthAdjuster.update();
}
public void run(String arg) {
imp = WindowManager.getCurrentImage();
if (arg.equals("add")) {
addToRoiManager(imp);
return;
}
if (imp == null) {
IJ.noImage();
return;
}
if (arg.equals("all")) imp.setRoi(0, 0, imp.getWidth(), imp.getHeight());
else if (arg.equals("none")) imp.killRoi();
else if (arg.equals("restore")) imp.restoreRoi();
else if (arg.equals("spline")) fitSpline();
else if (arg.equals("circle")) fitCircle(imp);
else if (arg.equals("ellipse")) createEllipse(imp);
else if (arg.equals("hull")) convexHull(imp);
else if (arg.equals("mask")) createMask(imp);
else if (arg.equals("from")) createSelectionFromMask(imp);
else if (arg.equals("inverse")) invert(imp);
else if (arg.equals("toarea")) lineToArea(imp);
else if (arg.equals("toline")) areaToLine(imp);
else if (arg.equals("properties")) {
setProperties("Properties ", imp.getRoi());
imp.draw();
} else if (arg.equals("band")) makeBand(imp);
else if (arg.equals("tobox")) toBoundingBox(imp);
else runMacro(arg);
}
void lineToArea(ImagePlus imp) {
Roi roi = imp.getRoi();
if (roi == null || !roi.isLine()) {
IJ.error("Line to Area", "Line selection required");
return;
}
Undo.setup(Undo.ROI, imp);
Roi roi2 = null;
if (roi.getType() == Roi.LINE) {
double width = roi.getStrokeWidth();
if (width <= 1.0) roi.setStrokeWidth(1.0000001);
FloatPolygon p = roi.getFloatPolygon();
roi.setStrokeWidth(width);
roi2 = new PolygonRoi(p, Roi.POLYGON);
roi2.setDrawOffset(roi.getDrawOffset());
} else {
ImageProcessor ip2 = new ByteProcessor(imp.getWidth(), imp.getHeight());
ip2.setColor(255);
roi.drawPixels(ip2);
// new ImagePlus("ip2", ip2.duplicate()).show();
ip2.setThreshold(255, 255, ImageProcessor.NO_LUT_UPDATE);
ThresholdToSelection tts = new ThresholdToSelection();
roi2 = tts.convert(ip2);
}
transferProperties(roi, roi2);
roi2.setStrokeWidth(0);
Color c = roi2.getStrokeColor();
if (c != null) // remove any transparency
roi2.setStrokeColor(new Color(c.getRed(), c.getGreen(), c.getBlue()));
imp.setRoi(roi2);
Roi.previousRoi = (Roi) roi.clone();
}
public boolean checkBeads() {
rmanager = RoiManager.getInstance();
if (IJ.versionLessThan("1.26i")) return false;
else if (rmanager == null) {
IJ.error("Add bead ROIs to the RoiManager first (select region then press [t]).");
return false;
}
nrois = rmanager.getCount();
if (nrois == 0) {
IJ.error("Add bead ROIs to the RoiManager first (select region then press [t]).");
return false;
}
rois = rmanager.getRoisAsArray();
return true;
}
public String getInfo() {
return version()
+ System.getProperty("os.name")
+ " "
+ System.getProperty("os.version")
+ "; "
+ IJ.freeMemory();
}
/**
* Execute the plugin functionality: duplicate and scale the given image.
*
* @return an Object[] array with the name and the scaled ImagePlus. Does NOT show the new, image;
* just returns it.
*/
public Object[] exec(
ImagePlus imp, String myMethod, int radius, double par1, double par2, boolean doIwhite) {
// 0 - Check validity of parameters
if (null == imp) return null;
ImageProcessor ip = imp.getProcessor();
int xe = ip.getWidth();
int ye = ip.getHeight();
// int [] data = (ip.getHistogram());
IJ.showStatus("Thresholding...");
long startTime = System.currentTimeMillis();
// 1 Do it
if (imp.getStackSize() == 1) {
ip.snapshot();
Undo.setup(Undo.FILTER, imp);
}
// Apply the selected algorithm
if (myMethod.equals("Bernsen")) {
Bernsen(imp, radius, par1, par2, doIwhite);
} else if (myMethod.equals("Contrast")) {
Contrast(imp, radius, par1, par2, doIwhite);
} else if (myMethod.equals("Mean")) {
Mean(imp, radius, par1, par2, doIwhite);
} else if (myMethod.equals("Median")) {
Median(imp, radius, par1, par2, doIwhite);
} else if (myMethod.equals("MidGrey")) {
MidGrey(imp, radius, par1, par2, doIwhite);
} else if (myMethod.equals("Niblack")) {
Niblack(imp, radius, par1, par2, doIwhite);
} else if (myMethod.equals("Otsu")) {
Otsu(imp, radius, par1, par2, doIwhite);
} else if (myMethod.equals("Phansalkar")) {
Phansalkar(imp, radius, par1, par2, doIwhite);
} else if (myMethod.equals("Sauvola")) {
Sauvola(imp, radius, par1, par2, doIwhite);
}
// IJ.showProgress((double)(255-i)/255);
imp.updateAndDraw();
imp.getProcessor().setThreshold(255, 255, ImageProcessor.NO_LUT_UPDATE);
// 2 - Return the threshold and the image
IJ.showStatus("\nDone " + (System.currentTimeMillis() - startTime) / 1000.0);
return new Object[] {imp};
}
private String version() {
return "ImageJ "
+ VERSION
+ BUILD
+ "; "
+ "Java "
+ System.getProperty("java.version")
+ (IJ.is64Bit() ? " [64-bit]; " : " [32-bit]; ");
}
/**
* Builds dialog to query users for projection parameters.
*
* @param start starting slice to display
* @param stop last slice
*/
protected GenericDialog buildControlDialog(int start, int stop) {
GenericDialog gd = new GenericDialog("ZProjection", IJ.getInstance());
gd.addNumericField("Start slice:", startSlice, 0 /*digits*/);
gd.addNumericField("Stop slice:", stopSlice, 0 /*digits*/);
gd.addChoice("Projection type", METHODS, METHODS[method]);
if (isHyperstack && imp.getNFrames() > 1 && imp.getNSlices() > 1)
gd.addCheckbox("All time frames", allTimeFrames);
return gd;
}
// Input/Output options
void io() {
GenericDialog gd = new GenericDialog("I/O Options");
gd.addNumericField("JPEG quality (0-100):", FileSaver.getJpegQuality(), 0, 3, "");
gd.addNumericField("GIF and PNG transparent index:", Prefs.getTransparentIndex(), 0, 3, "");
gd.addStringField(
"File extension for tables (.txt, .xls or .csv):", Prefs.get("options.ext", ".csv"), 4);
gd.addCheckbox("Use JFileChooser to open/save", Prefs.useJFileChooser);
if (!IJ.isMacOSX())
gd.addCheckbox("Use_file chooser to import sequences", Prefs.useFileChooser);
gd.addCheckbox("Save TIFF and raw in Intel byte order", Prefs.intelByteOrder);
gd.addCheckbox("Skip dialog when opening .raw files", Prefs.skipRawDialog);
gd.setInsets(15, 20, 0);
gd.addMessage("Results Table Options");
gd.setInsets(3, 40, 0);
gd.addCheckbox("Copy_column headers", Prefs.copyColumnHeaders);
gd.setInsets(0, 40, 0);
gd.addCheckbox("Copy_row numbers", !Prefs.noRowNumbers);
gd.setInsets(0, 40, 0);
gd.addCheckbox("Save_column headers", !Prefs.dontSaveHeaders);
gd.setInsets(0, 40, 0);
gd.addCheckbox("Save_row numbers", !Prefs.dontSaveRowNumbers);
gd.showDialog();
if (gd.wasCanceled()) return;
int quality = (int) gd.getNextNumber();
if (quality < 0) quality = 0;
if (quality > 100) quality = 100;
FileSaver.setJpegQuality(quality);
int transparentIndex = (int) gd.getNextNumber();
Prefs.setTransparentIndex(transparentIndex);
String extension = gd.getNextString();
if (!extension.startsWith(".")) extension = "." + extension;
Prefs.set("options.ext", extension);
Prefs.useJFileChooser = gd.getNextBoolean();
if (!IJ.isMacOSX()) Prefs.useFileChooser = gd.getNextBoolean();
Prefs.intelByteOrder = gd.getNextBoolean();
Prefs.skipRawDialog = gd.getNextBoolean();
Prefs.copyColumnHeaders = gd.getNextBoolean();
Prefs.noRowNumbers = !gd.getNextBoolean();
Prefs.dontSaveHeaders = !gd.getNextBoolean();
Prefs.dontSaveRowNumbers = !gd.getNextBoolean();
return;
}
void Bernsen(ImagePlus imp, int radius, double par1, double par2, boolean doIwhite) {
// Bernsen recommends WIN_SIZE = 31 and CONTRAST_THRESHOLD = 15.
// 1) Bernsen J. (1986) "Dynamic Thresholding of Grey-Level Images"
// Proc. of the 8th Int. Conf. on Pattern Recognition, pp. 1251-1255
// 2) Sezgin M. and Sankur B. (2004) "Survey over Image Thresholding
// Techniques and Quantitative Performance Evaluation" Journal of
// Electronic Imaging, 13(1): 146-165
// http://citeseer.ist.psu.edu/sezgin04survey.html
// Ported to ImageJ plugin from E Celebi's fourier_0.8 routines
// This version uses a circular local window, instead of a rectagular one
ImagePlus Maximp, Minimp;
ImageProcessor ip = imp.getProcessor(), ipMax, ipMin;
int contrast_threshold = 15;
int local_contrast;
int mid_gray;
byte object;
byte backg;
int temp;
if (par1 != 0) {
IJ.log("Bernsen: changed contrast_threshold from :" + contrast_threshold + " to:" + par1);
contrast_threshold = (int) par1;
}
if (doIwhite) {
object = (byte) 0xff;
backg = (byte) 0;
} else {
object = (byte) 0;
backg = (byte) 0xff;
}
Maximp = duplicateImage(ip);
ipMax = Maximp.getProcessor();
RankFilters rf = new RankFilters();
rf.rank(ipMax, radius, rf.MAX); // Maximum
// Maximp.show();
Minimp = duplicateImage(ip);
ipMin = Minimp.getProcessor();
rf.rank(ipMin, radius, rf.MIN); // Minimum
// Minimp.show();
byte[] pixels = (byte[]) ip.getPixels();
byte[] max = (byte[]) ipMax.getPixels();
byte[] min = (byte[]) ipMin.getPixels();
for (int i = 0; i < pixels.length; i++) {
local_contrast = (int) ((max[i] & 0xff) - (min[i] & 0xff));
mid_gray = (int) ((min[i] & 0xff) + (max[i] & 0xff)) / 2;
temp = (int) (pixels[i] & 0x0000ff);
if (local_contrast < contrast_threshold)
pixels[i] = (mid_gray >= 128) ? object : backg; // Low contrast region
else pixels[i] = (temp >= mid_gray) ? object : backg;
}
// imp.updateAndDraw();
return;
}
void createSelectionFromMask(ImagePlus imp) {
ImageProcessor ip = imp.getProcessor();
if (ip.getMinThreshold() != ImageProcessor.NO_THRESHOLD) {
IJ.runPlugIn("ij.plugin.filter.ThresholdToSelection", "");
return;
}
if (!ip.isBinary()) {
IJ.error(
"Create Selection",
"This command creates a composite selection from\n"
+ "a mask (8-bit binary image with white background)\n"
+ "or from an image that has been thresholded using\n"
+ "the Image>Adjust>Threshold tool. The current\n"
+ "image is not a mask and has not been thresholded.");
return;
}
int threshold = ip.isInvertedLut() ? 255 : 0;
ip.setThreshold(threshold, threshold, ImageProcessor.NO_LUT_UPDATE);
IJ.runPlugIn("ij.plugin.filter.ThresholdToSelection", "");
}
void areaToLine(ImagePlus imp) {
Roi roi = imp.getRoi();
if (roi == null || !roi.isArea()) {
IJ.error("Area to Line", "Area selection required");
return;
}
Polygon p = roi.getPolygon();
if (p == null) return;
int type1 = roi.getType();
if (type1 == Roi.COMPOSITE) {
IJ.error("Area to Line", "Composite selections cannot be converted to lines.");
return;
}
int type2 = Roi.POLYLINE;
if (type1 == Roi.OVAL
|| type1 == Roi.FREEROI
|| type1 == Roi.TRACED_ROI
|| ((roi instanceof PolygonRoi) && ((PolygonRoi) roi).isSplineFit())) type2 = Roi.FREELINE;
Roi roi2 = new PolygonRoi(p.xpoints, p.ypoints, p.npoints, type2);
imp.setRoi(roi2);
}
void invert(ImagePlus imp) {
Roi roi = imp.getRoi();
if (roi == null || !roi.isArea()) {
IJ.error("Inverse", "Area selection required");
return;
}
ShapeRoi s1, s2;
if (roi instanceof ShapeRoi) s1 = (ShapeRoi) roi;
else s1 = new ShapeRoi(roi);
s2 = new ShapeRoi(new Roi(0, 0, imp.getWidth(), imp.getHeight()));
imp.setRoi(s1.xor(s2));
}
public void run(String arg) {
if (IJ.versionLessThan("1.49d")) return;
if (!showDialog()) return;
SaveDialog sd = new SaveDialog("Save as Bricks...", "", "");
basename = sd.getFileName();
directory = sd.getDirectory();
if (basename == null || directory == null) return;
build_bricks();
}
// Conversion Options
void conversions() {
double[] weights = ColorProcessor.getWeightingFactors();
boolean weighted = !(weights[0] == 1d / 3d && weights[1] == 1d / 3d && weights[2] == 1d / 3d);
// boolean weighted = !(Math.abs(weights[0]-1d/3d)<0.0001 && Math.abs(weights[1]-1d/3d)<0.0001
// && Math.abs(weights[2]-1d/3d)<0.0001);
GenericDialog gd = new GenericDialog("Conversion Options");
gd.addCheckbox("Scale when converting", ImageConverter.getDoScaling());
String prompt = "Weighted RGB conversions";
if (weighted)
prompt +=
" (" + IJ.d2s(weights[0]) + "," + IJ.d2s(weights[1]) + "," + IJ.d2s(weights[2]) + ")";
gd.addCheckbox(prompt, weighted);
gd.showDialog();
if (gd.wasCanceled()) return;
ImageConverter.setDoScaling(gd.getNextBoolean());
Prefs.weightedColor = gd.getNextBoolean();
if (!Prefs.weightedColor) ColorProcessor.setWeightingFactors(1d / 3d, 1d / 3d, 1d / 3d);
else if (Prefs.weightedColor && !weighted)
ColorProcessor.setWeightingFactors(0.299, 0.587, 0.114);
return;
}
void abortPluginOrMacro(ImagePlus imp) {
if (imp != null) {
ImageWindow win = imp.getWindow();
if (win != null) {
win.running = false;
win.running2 = false;
}
}
Macro.abort();
Interpreter.abort();
if (Interpreter.getInstance() != null) IJ.beep();
}
public void run(ImageProcessor ip) {
dimz = stack.getSize();
dimy = stack.getWidth();
dimx = stack.getHeight();
SaveDialog sd = new SaveDialog("Save Measurements as Text...", "res", ".dat");
String name = sd.getFileName();
if (name == null) return;
String directory = sd.getDirectory();
nb = calculnb(stack); // -1;
IJ.showStatus("Measure parameters for the " + nb + " objects ...");
if (nb < 1) {
IJ.showMessage("volume must be labeled");
} else {
double[] volume_m = new double[nb];
int[] volume_p = new int[nb];
double[] surface = new double[nb];
double[] surfacenb = new double[nb];
double[][][] I = new double[3][3][nb];
double[][] J = new double[3][nb];
double[][][] dir = new double[3][3][nb];
double[] xg = new double[nb];
double[] yg = new double[nb];
double[] zg = new double[nb];
byte[] bord = new byte[nb];
// double[] a = new double[nb];
// double[] b = new double[nb];
// double[] c = new double[nb];
// double[] Fab = new double[nb];
// double[] Fac = new double[nb];
// double[] Fbc = new double[nb];
// double[] sp = new double[nb];
double[][] lmin = new double[nb][3];
double[][] lmax = new double[nb][3];
IJ.showStatus("Measure surfaces ...");
calculmarchsurfstack(stack, nb, surface, volume_m);
calculmarchsurfstacknb(stack, nb, surfacenb);
// calculvolumestack(stack,nb,volume_p);
IJ.showStatus("Measure volumes and inertia ...");
calculcgstack(stack, nb, volume_p, xg, yg, zg);
calculinertiestack(stack, nb, xg, yg, zg, I);
inertie(nb, I, J, dir);
IJ.showStatus("Measure bounding boxes ...");
boitestack(stack, nb, xg, yg, zg, dir, lmin, lmax);
borderstack(stack, nb, bord);
IJ.showStatus("Save results ...");
sauvegarde(
volume_p, volume_m, surface, surfacenb, xg, yg, zg, J, dir, nb, bord, lmin, lmax,
directory, name);
volume_m = null;
volume_p = null;
surface = null;
xg = null;
yg = null;
zg = null;
}
}
void lineWidth() {
int width = (int) IJ.getNumber("Line Width:", Line.getWidth());
if (width == IJ.CANCELED) return;
Line.setWidth(width);
LineWidthAdjuster.update();
ImagePlus imp = WindowManager.getCurrentImage();
if (imp != null && imp.isProcessor()) {
ImageProcessor ip = imp.getProcessor();
ip.setLineWidth(Line.getWidth());
Roi roi = imp.getRoi();
if (roi != null && roi.isLine()) imp.draw();
}
}
