/** Displays a modal options dialog. */
public boolean showDialog() {
Calibration cal = imp != null ? imp.getCalibration() : (new Calibration());
double unitSquared = cal.pixelWidth * cal.pixelHeight;
if (pixelUnits) unitSquared = 1.0;
if (Macro.getOptions() != null) {
boolean oldMacro = updateMacroOptions();
if (oldMacro) unitSquared = 1.0;
staticMinSize = 0.0;
staticMaxSize = DEFAULT_MAX_SIZE;
staticMinCircularity = 0.0;
staticMaxCircularity = 1.0;
staticShowChoice = NOTHING;
}
GenericDialog gd = new GenericDialog("Analyze Particles");
minSize = staticMinSize;
maxSize = staticMaxSize;
minCircularity = staticMinCircularity;
maxCircularity = staticMaxCircularity;
showChoice = staticShowChoice;
if (maxSize == 999999) maxSize = DEFAULT_MAX_SIZE;
options = staticOptions;
String unit = cal.getUnit();
boolean scaled = cal.scaled();
if (unit.equals("inch")) {
unit = "pixel";
unitSquared = 1.0;
scaled = false;
pixelUnits = true;
}
String units = unit + "^2";
int places = 0;
double cmin = minSize * unitSquared;
if ((int) cmin != cmin) places = 2;
double cmax = maxSize * unitSquared;
if ((int) cmax != cmax && cmax != DEFAULT_MAX_SIZE) places = 2;
String minStr = ResultsTable.d2s(cmin, places);
if (minStr.indexOf("-") != -1) {
for (int i = places; i <= 6; i++) {
minStr = ResultsTable.d2s(cmin, i);
if (minStr.indexOf("-") == -1) break;
}
}
String maxStr = ResultsTable.d2s(cmax, places);
if (maxStr.indexOf("-") != -1) {
for (int i = places; i <= 6; i++) {
maxStr = ResultsTable.d2s(cmax, i);
if (maxStr.indexOf("-") == -1) break;
}
}
if (scaled) gd.setInsets(5, 0, 0);
gd.addStringField("Size (" + units + "):", minStr + "-" + maxStr, 12);
if (scaled) {
gd.setInsets(0, 40, 5);
gd.addCheckbox("Pixel units", pixelUnits);
}
gd.addStringField("Circularity:", IJ.d2s(minCircularity) + "-" + IJ.d2s(maxCircularity), 12);
gd.addChoice("Show:", showStrings, showStrings[showChoice]);
String[] labels = new String[8];
boolean[] states = new boolean[8];
labels[0] = "Display results";
states[0] = (options & SHOW_RESULTS) != 0;
labels[1] = "Exclude on edges";
states[1] = (options & EXCLUDE_EDGE_PARTICLES) != 0;
labels[2] = "Clear results";
states[2] = (options & CLEAR_WORKSHEET) != 0;
labels[3] = "Include holes";
states[3] = (options & INCLUDE_HOLES) != 0;
labels[4] = "Summarize";
states[4] = (options & DISPLAY_SUMMARY) != 0;
labels[5] = "Record starts";
states[5] = (options & RECORD_STARTS) != 0;
labels[6] = "Add to Manager";
states[6] = (options & ADD_TO_MANAGER) != 0;
labels[7] = "In_situ Show";
states[7] = (options & IN_SITU_SHOW) != 0;
gd.addCheckboxGroup(4, 2, labels, states);
gd.addHelp(IJ.URL + "/docs/menus/analyze.html#ap");
gd.showDialog();
if (gd.wasCanceled()) return false;
String size = gd.getNextString(); // min-max size
if (scaled) pixelUnits = gd.getNextBoolean();
if (pixelUnits) unitSquared = 1.0;
else unitSquared = cal.pixelWidth * cal.pixelHeight;
String[] minAndMax = Tools.split(size, " -");
double mins = gd.parseDouble(minAndMax[0]);
double maxs = minAndMax.length == 2 ? gd.parseDouble(minAndMax[1]) : Double.NaN;
minSize = Double.isNaN(mins) ? DEFAULT_MIN_SIZE : mins / unitSquared;
maxSize = Double.isNaN(maxs) ? DEFAULT_MAX_SIZE : maxs / unitSquared;
if (minSize < DEFAULT_MIN_SIZE) minSize = DEFAULT_MIN_SIZE;
if (maxSize < minSize) maxSize = DEFAULT_MAX_SIZE;
staticMinSize = minSize;
staticMaxSize = maxSize;
minAndMax = Tools.split(gd.getNextString(), " -"); // min-max circularity
double minc = gd.parseDouble(minAndMax[0]);
double maxc = minAndMax.length == 2 ? gd.parseDouble(minAndMax[1]) : Double.NaN;
minCircularity = Double.isNaN(minc) ? 0.0 : minc;
maxCircularity = Double.isNaN(maxc) ? 1.0 : maxc;
if (minCircularity < 0.0 || minCircularity > 1.0) minCircularity = 0.0;
if (maxCircularity < minCircularity || maxCircularity > 1.0) maxCircularity = 1.0;
if (minCircularity == 1.0 && maxCircularity == 1.0) minCircularity = 0.0;
staticMinCircularity = minCircularity;
staticMaxCircularity = maxCircularity;
if (gd.invalidNumber()) {
IJ.error("Bins invalid.");
canceled = true;
return false;
}
showChoice = gd.getNextChoiceIndex();
staticShowChoice = showChoice;
if (gd.getNextBoolean()) options |= SHOW_RESULTS;
else options &= ~SHOW_RESULTS;
if (gd.getNextBoolean()) options |= EXCLUDE_EDGE_PARTICLES;
else options &= ~EXCLUDE_EDGE_PARTICLES;
if (gd.getNextBoolean()) options |= CLEAR_WORKSHEET;
else options &= ~CLEAR_WORKSHEET;
if (gd.getNextBoolean()) options |= INCLUDE_HOLES;
else options &= ~INCLUDE_HOLES;
if (gd.getNextBoolean()) options |= DISPLAY_SUMMARY;
else options &= ~DISPLAY_SUMMARY;
if (gd.getNextBoolean()) options |= RECORD_STARTS;
else options &= ~RECORD_STARTS;
if (gd.getNextBoolean()) options |= ADD_TO_MANAGER;
else options &= ~ADD_TO_MANAGER;
if (gd.getNextBoolean()) options |= IN_SITU_SHOW;
else options &= ~IN_SITU_SHOW;
staticOptions = options;
options |= SHOW_PROGRESS;
if ((options & DISPLAY_SUMMARY) != 0)
Analyzer.setMeasurements(Analyzer.getMeasurements() | AREA);
return true;
}
public boolean get_errors(double[] params, int[] fixes) {
GenericDialog gd = new GenericDialog("Error Options");
String[] methods = {"Support Plane", "Monte Carlo"};
gd.addChoice("Method", methods, methods[0]);
float conf = 0.67f;
gd.addNumericField("SP_Confidence Limit (%)", (int) (conf * 100.0f), 5, 10, null);
String[] labels = {"P1", "P2", "P3", "P4", "P5", "P6", "P7", "P8", "P9", "P10"};
gd.addChoice("SP_Parameter", labels, labels[0]);
double spacing = 0.01;
gd.addNumericField("SP_Chi^2_plot_spacing (% of value)?", spacing * 100.0, 2, 10, null);
int ntrials = 100;
gd.addNumericField("MC_#_Trials", ntrials, 0);
gd.showDialog();
if (gd.wasCanceled()) {
return false;
}
int methodindex = gd.getNextChoiceIndex();
conf = 0.01f * (float) gd.getNextNumber();
int paramindex = gd.getNextChoiceIndex();
spacing = 0.01 * gd.getNextNumber();
ntrials = (int) gd.getNextNumber();
if (methodindex == 0) {
support_plane_errors_v2 erclass = new support_plane_errors_v2(this, 0.0001, 50, false, 0.1);
int errindex = paramindex;
int nfit = 0;
for (int i = 0; i < labels.length; i++) {
if (fixes[i] == 0) {
nfit++;
}
}
int npts = tempdata.length;
int dofnum = npts - (nfit - 1) - 1;
int dofden = npts - nfit - 1;
double flim = (new jdist()).FLimit(dofnum, dofden, (double) conf);
IJ.log("FLimit = " + (float) flim);
if (flim == Double.NaN && flim < 1.0) {
IJ.showMessage("Invalid Limiting F Value");
return false;
}
double truespacing = Math.abs(params[errindex] * spacing);
double[][] c2plot =
erclass.geterrors(
params, fixes, constraints, tempdata, weights, flim, truespacing, errindex);
IJ.log("upper limit = " + c2plot[1][0] + " lower limit = " + c2plot[0][0]);
IJ.log(
"upper error = "
+ (c2plot[1][0] - params[errindex])
+ " lower error = "
+ (params[errindex] - c2plot[0][0]));
int templength = c2plot[0].length;
float[][] c2plotf = new float[2][templength - 1];
for (int i = 0; i < (templength - 1); i++) {
c2plotf[0][i] = (float) c2plot[0][i + 1];
c2plotf[1][i] = (float) c2plot[1][i + 1];
}
new PlotWindow4("c2 plot", labels[errindex], "Chi^2", c2plotf[0], c2plotf[1]).draw();
} else {
StringBuffer sb = new StringBuffer();
sb.append("Trial\t");
for (int i = 0; i < labels.length; i++) {
if (fixes[i] == 0) sb.append(labels[i] + "\t");
}
sb.append("chi^2");
tw = new TextWindow("Monte Carlo Results", sb.toString(), "", 400, 400);
redirect = true;
monte_carlo_errors_v2 erclass = new monte_carlo_errors_v2(this, 0.0001, 50, false, 0.1);
double[][] errors = erclass.geterrors(params, fixes, constraints, tempdata, weights, ntrials);
sb = new StringBuffer();
sb.append("StDev\t");
for (int i = 0; i < errors.length; i++) {
float[] ferr = new float[errors[0].length];
for (int j = 0; j < ferr.length; j++) ferr[j] = (float) errors[i][j];
float stdev = jstatistics.getstatistic("StDev", ferr, null);
sb.append("" + stdev);
if (i < (errors.length - 1)) sb.append("\t");
}
tw.append(sb.toString());
redirect = false;
}
return true;
}
public void run(String arg) {
Frame[] niframes = WindowManager.getNonImageWindows();
String[] titles = new String[niframes.length + 1];
for (int i = 0; i < niframes.length; i++) {
titles[i] = niframes[i].getTitle();
}
titles[niframes.length] = "Clipboard";
GenericDialog gd = new GenericDialog("Windows");
boolean importfile = false;
gd.addCheckbox("Import from file?", importfile);
gd.addChoice("Windows", titles, titles[0]);
boolean hasxvals = false;
gd.addCheckbox("X Vals Column?", hasxvals);
boolean multix = false;
gd.addCheckbox("Multi_X_Columns?", multix);
boolean skipendzeros = false;
gd.addCheckbox("Skip_end_zeros?", skipendzeros);
String[] delimiters = {"Tab", "Comma", "Space"};
gd.addChoice("Delimiter", delimiters, delimiters[0]);
gd.showDialog();
if (gd.wasCanceled()) {
return;
}
importfile = gd.getNextBoolean();
int index = gd.getNextChoiceIndex();
hasxvals = gd.getNextBoolean();
multix = gd.getNextBoolean();
skipendzeros = gd.getNextBoolean();
int delimindex = gd.getNextChoiceIndex();
if (multix) hasxvals = true;
String textdata = "";
if (importfile) {
OpenDialog od = new OpenDialog("Open File", "", ".txt");
String directory = od.getDirectory();
String name = od.getFileName();
if (name == null) {
return;
}
try {
File infile = new File(directory + name);
BufferedReader b = new BufferedReader(new FileReader(infile));
textdata = (new jdataio()).readstringfile(b);
b.close();
} catch (IOException e) {
return;
}
} else {
if (index == niframes.length) {
// here we get the data from the clipboard
Transferable t = Toolkit.getDefaultToolkit().getSystemClipboard().getContents(null);
try {
if (t != null && t.isDataFlavorSupported(DataFlavor.stringFlavor)) {
textdata = (String) t.getTransferData(DataFlavor.stringFlavor);
}
} catch (UnsupportedFlavorException e) {
} catch (IOException e) {
}
if (textdata.equals("")) {
IJ.error("Error copying from clipboard.");
return;
}
} else {
if (niframes[index] instanceof Editor) {
Editor tw = (Editor) niframes[index];
textdata = tw.getText();
} else {
if (niframes[index] instanceof TextWindow) {
TextWindow tw = (TextWindow) niframes[index];
textdata = tw.getTextPanel().getText();
} else {
IJ.showMessage("Not a valid text window");
return;
}
}
}
}
if (textdata == null) {
IJ.showMessage("Error in Obtaining String");
return;
}
if (textdata.indexOf("\r") >= 0) {
textdata = textdata.replace('\r', '\n');
}
char[] delims = {'\t', ',', ' '};
delimit_string ds = new delimit_string(delims[delimindex]);
String[] rows = ds.getrows(textdata);
int lines = rows.length;
int columns = ds.getnumcolumns(rows[0]);
int ycolumns = columns;
if (hasxvals) {
if (multix) {
ycolumns /= 2;
} else {
ycolumns--;
}
}
if (multix) {
float[][] ydata = new float[ycolumns][lines];
float[][] xdata = new float[ycolumns][lines];
for (int i = 0; i < lines; i++) {
float[] temp = ds.delim2float(rows[i], columns);
for (int j = 0; j < ycolumns; j++) {
ydata[j][i] = temp[2 * j + 1];
xdata[j][i] = temp[2 * j];
}
}
int[] npts = new int[ycolumns];
for (int i = 0; i < ycolumns; i++) {
npts[i] = lines;
}
if (skipendzeros) {
for (int i = 0; i < ycolumns; i++) {
int counter = lines - 1;
while ((xdata[i][counter] == 0.0f || Float.isNaN(xdata[i][counter])) && counter > 0) {
xdata[i][counter] = 0.0f;
ydata[i][counter] = 0.0f;
npts[i]--;
counter--;
}
}
}
(new PlotWindow4("Text Plot", "x", "y", xdata, ydata, npts)).draw();
} else {
float[][] tempydata = new float[ycolumns][lines];
float[] tempxdata = new float[lines];
float[][] xdata = null;
float[][] ydata = null;
int startcolumn = 0;
if (hasxvals) startcolumn = 1;
for (int i = 0; i < lines; i++) {
float[] temp = ds.delim2float(rows[i], columns);
if (hasxvals) {
tempxdata[i] = temp[0];
} else {
tempxdata[i] = (float) (i + 1);
}
for (int j = 0; j < ycolumns; j++) {
tempydata[j][i] = temp[j + startcolumn];
}
}
int[] npts = new int[ycolumns];
npts[0] = lines;
if (skipendzeros) {
int maxpts = 0;
for (int i = 0; i < ycolumns; i++) {
int counter = lines - 1;
npts[i] = lines;
while ((tempydata[i][counter] == 0.0f || Float.isNaN(tempydata[i][counter]))
&& counter > 0) {
npts[i]--;
counter--;
}
if (npts[i] > maxpts) maxpts = npts[i];
IJ.log("" + npts[i]);
}
ydata = new float[ycolumns][maxpts];
xdata = new float[ycolumns][maxpts];
for (int i = 0; i < ycolumns; i++) {
// npts[i]=npts[0];
System.arraycopy(tempxdata, 0, xdata[i], 0, npts[i]);
System.arraycopy(tempydata[i], 0, ydata[i], 0, npts[i]);
}
} else {
ydata = tempydata;
xdata = new float[ycolumns][];
for (int i = 0; i < ycolumns; i++) {
npts[i] = npts[0];
xdata[i] = tempxdata.clone();
}
}
(new PlotWindow4("Text Plot", "x", "y", xdata, ydata, npts)).draw();
}
}