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;
}
}
public void dragOver(DropTargetDragEvent e) {
if (IJ.debugMode) IJ.log("DragOver: " + e.getLocation());
Point loc = e.getLocation();
int buttonSize = Toolbar.getButtonSize();
int width = IJ.getInstance().getSize().width;
openAsVirtualStack = width - loc.x <= buttonSize;
if (openAsVirtualStack) IJ.showStatus("<<Open as Virtual Stack>>");
else IJ.showStatus("<<Drag and Drop>>");
}
void loadParticleResults(String filename, ResultsTable res) {
try {
String line;
FileReader fr = new FileReader(filename);
BufferedReader br = new BufferedReader(fr);
java.lang.String header =
" Intensity X (px) Y (px) X (nm) Y (nm) Z (nm) Left-Width(px) Right-Width (px) Up-Height (px) Down-Height (px) X Symmetry (%) Y Symmetry (%) Width minus Height (px) Frame Number";
java.lang.String firstline = br.readLine();
if (!firstline.contains("X (px) Y (px) X (nm) Y (nm) Z (nm)")) {
IJ.error("File does not seam to be a Particles Table file");
IJ.log("Found header: " + firstline);
IJ.log("Expecting: " + header);
return;
}
res.reset();
int counter = 1;
java.util.concurrent.locks.Lock lock = new java.util.concurrent.locks.ReentrantLock();
ThreadedLoader tloader = new ThreadedLoader();
// java.lang.String txt = fr.read();
while ((line = br.readLine()) != null) {
tloader = new ThreadedLoader();
tloader.mysetup(res, lock, line);
tloader.start();
IJ.showStatus("Loading particle " + counter + "... sit back and relax.");
counter++;
}
try {
tloader.join();
} catch (Exception e) {
IJ.error("" + e);
}
if (res.getCounter() < 5000000) {
IJ.showStatus("Creating particle table, this should take a few seconds...");
res.show("Results");
} else
IJ.showMessage(
"Warning",
"Results table has too many particles, they will not be shown but the data still exists within it\nyou can still use all the plugin functionality or save table changes though the 'Save Particle Table' command.");
fr.close();
IJ.showStatus("Done loading table...");
} 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;
}
}
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;
}
}
public static byte[] download(String urlString, String name) {
int maxLength = 52428800; // 50MB
URL url = null;
boolean unknownLength = false;
byte[] data = null;
;
int n = 0;
try {
url = new URL(urlString);
if (IJ.debugMode) IJ.log("PluginInstaller: " + urlString + " " + url);
if (url == null) return null;
URLConnection uc = url.openConnection();
int len = uc.getContentLength();
unknownLength = len < 0;
if (unknownLength) len = maxLength;
if (name != null) IJ.showStatus("Downloading " + url.getFile());
InputStream in = uc.getInputStream();
data = new byte[len];
int lenk = len / 1024;
while (n < len) {
int count = in.read(data, n, len - n);
if (count < 0) break;
n += count;
if (name != null)
IJ.showStatus("Downloading " + name + " (" + (n / 1024) + "/" + lenk + "k)");
IJ.showProgress(n, len);
}
in.close();
} catch (Exception e) {
String msg = "" + e;
if (!msg.contains("://")) msg += "\n " + urlString;
IJ.error("Plugin Installer", msg);
return null;
} finally {
IJ.showProgress(1.0);
}
if (name != null) IJ.showStatus("");
if (unknownLength) {
byte[] data2 = data;
data = new byte[n];
for (int i = 0; i < n; i++) data[i] = data2[i];
}
return data;
}
private void fitglobal() {
int nparams = 11;
int nsel = 0;
for (int i = 0; i < ncurves; i++) {
if (include[i]) {
nsel++;
}
}
double[][] params = new double[nsel][nparams];
String[][] tempformulas = new String[nsel][nparams];
double[][][] constraints = new double[2][nsel][nparams];
int[][] vflmatrix = new int[nsel][nparams];
int counter = 0;
for (int i = 0; i < ncurves; i++) {
if (include[i]) {
for (int j = 0; j < nparams; j++) {
params[counter][j] = globalparams[i][j];
tempformulas[counter][j] = globalformulas[i][j];
constraints[0][counter][j] = globalconstraints[0][i][j];
constraints[1][counter][j] = globalconstraints[1][i][j];
vflmatrix[counter][j] = globalvflmatrix[i][j];
}
counter++;
}
for (int j = 0; j < nparams; j++) {
undoparams[i][j] = globalparams[i][j];
undoformulas[i][j] = globalformulas[i][j];
undovflmatrix[i][j] = globalvflmatrix[i][j];
}
for (int j = 0; j < xpts; j++) {
for (int k = 0; k < ypts; k++) {
undofit[i][j][k] = fit[i][j][k];
}
}
undoc2[i] = c2[i];
}
undoglobalc2 = globalc2;
if (showglobalfitdialog(params, tempformulas, vflmatrix)) {
counter = 0;
for (int i = 0; i < ncurves; i++) {
if (include[i]) {
for (int j = 0; j < nparams; j++) {
globalparams[i][j] = params[counter][j];
globalformulas[i][j] = tempformulas[counter][j];
globalvflmatrix[i][j] = vflmatrix[counter][j];
}
counter++;
}
}
double[] stats = new double[2];
float[][] tempdata = new float[nsel][xpts * ypts];
float[][] tempweights = new float[nsel][xpts * ypts];
counter = 0;
for (int i = 0; i < ncurves; i++) {
if (include[i]) {
for (int j = 0; j < xpts; j++) {
for (int k = 0; k < ypts; k++) {
tempdata[counter][j + k * xpts] = (float) ((double) pch[i][j][k] / (double) nmeas[i]);
tempweights[counter][j + k * xpts] = weights[i][j][k];
}
}
counter++;
}
}
int tempmaxiter = globalfitclass.maxiter;
if (checkc2) {
globalfitclass.changemaxiter(0);
}
double[] tempc2vals = new double[nsel];
IJ.showStatus("Fitting Globally");
float[][] tempfit =
globalfitclass.fitdata(
params,
vflmatrix,
tempformulas,
paramsnames,
constraints,
tempdata,
tempweights,
stats,
tempc2vals,
false);
IJ.showStatus("Fit Complete");
globalfitclass.changemaxiter(tempmaxiter);
globalc2 = stats[1];
globalc2label.setText("Global chi^2 = " + (float) globalc2);
counter = 0;
for (int i = 0; i < ncurves; i++) {
if (include[i]) {
for (int j = 0; j < xpts; j++) {
for (int k = 0; k < ypts; k++) {
fit[i][j][k] = tempfit[counter][j + xpts * k] * (float) nmeas[i];
}
}
if (i == dispcurve) {
pwfit.updateSeries(fit[dispcurve], 1, true);
}
for (int j = 0; j < nparams; j++) {
globalparams[i][j] = params[counter][j];
}
c2[i] = tempc2vals[counter];
c2array[i].setText("" + (float) c2[i]);
counter++;
}
}
float[] temp = pwfit.getLimits();
temp[4] = 1.0f;
pwfit.setLimits(temp);
}
}
public void showresults(String results) {
IJ.showStatus(results);
IJ.log(results);
}
public void mousePressed(MouseEvent e) {
Undo.reset();
if (!Prefs.noClickToGC) System.gc();
IJ.showStatus(version() + IJ.freeMemory());
if (IJ.debugMode) IJ.log("Windows: " + WindowManager.getWindowCount());
}
/**
* If 'applet' is not null, creates a new ImageJ frame that runs as an applet. If 'mode' is
* ImageJ.EMBEDDED and 'applet is null, creates an embedded (non-standalone) version of ImageJ.
*/
public ImageJ(java.applet.Applet applet, int mode) {
super("ImageJ");
if ((mode & DEBUG) != 0) IJ.setDebugMode(true);
mode = mode & 255;
if (IJ.debugMode) IJ.log("ImageJ starting in debug mode: " + mode);
embedded = applet == null && (mode == EMBEDDED || mode == NO_SHOW);
this.applet = applet;
String err1 = Prefs.load(this, applet);
setBackground(backgroundColor);
Menus m = new Menus(this, applet);
String err2 = m.addMenuBar();
m.installPopupMenu(this);
setLayout(new BorderLayout());
// Tool bar
toolbar = new Toolbar();
toolbar.addKeyListener(this);
add("Center", toolbar);
// Status bar
statusBar = new Panel();
statusBar.setLayout(new BorderLayout());
statusBar.setForeground(Color.black);
statusBar.setBackground(backgroundColor);
statusLine = new JLabel();
statusLine.setFont(new Font("SansSerif", Font.PLAIN, 13));
statusLine.addKeyListener(this);
statusLine.addMouseListener(this);
statusBar.add("Center", statusLine);
progressBar = new ProgressBar(120, 20);
progressBar.addKeyListener(this);
progressBar.addMouseListener(this);
statusBar.add("East", progressBar);
add("South", statusBar);
IJ.init(this, applet);
addKeyListener(this);
addWindowListener(this);
setFocusTraversalKeysEnabled(false);
m.installStartupMacroSet(); // add custom tools
runStartupMacro();
Point loc = getPreferredLocation();
Dimension tbSize = toolbar.getPreferredSize();
setCursor(Cursor.getDefaultCursor()); // work-around for JDK 1.1.8 bug
if (mode != NO_SHOW) {
if (IJ.isWindows())
try {
setIcon();
} catch (Exception e) {
}
setLocation(loc.x, loc.y);
setResizable(!IJ.isMacOSX());
pack();
setVisible(true);
}
if (err1 != null) IJ.error(err1);
if (err2 != null) {
IJ.error(err2);
IJ.runPlugIn("ij.plugin.ClassChecker", "");
}
if (IJ.isMacintosh() && applet == null) {
Object qh = null;
qh = IJ.runPlugIn("MacAdapter", "");
if (qh == null) IJ.runPlugIn("QuitHandler", "");
}
if (applet == null) IJ.runPlugIn("ij.plugin.DragAndDrop", "");
String str = m.getMacroCount() == 1 ? " macro" : " macros";
IJ.showStatus(version() + m.getPluginCount() + " commands; " + m.getMacroCount() + str);
configureProxy();
if (applet == null) loadCursors();
}
/**
* If 'applet' is not null, creates a new ImageJ frame that runs as an applet. If 'mode' is
* ImageJ.EMBEDDED and 'applet is null, creates an embedded (non-standalone) version of ImageJ.
*/
public ImageJ(java.applet.Applet applet, int mode) {
super("ImageJ");
embedded = applet == null && (mode == EMBEDDED || mode == NO_SHOW);
this.applet = applet;
String err1 = Prefs.load(this, applet);
if (IJ.isLinux()) {
backgroundColor = new Color(240, 240, 240);
setBackground(backgroundColor);
}
Menus m = new Menus(this, applet);
String err2 = m.addMenuBar();
m.installPopupMenu(this);
setLayout(new GridLayout(2, 1));
// Tool bar
toolbar = new Toolbar();
toolbar.addKeyListener(this);
add(toolbar);
// Status bar
statusBar = new Panel();
statusBar.setLayout(new BorderLayout());
statusBar.setForeground(Color.black);
statusBar.setBackground(backgroundColor);
statusLine = new Label();
statusLine.setFont(SansSerif12);
statusLine.addKeyListener(this);
statusLine.addMouseListener(this);
statusBar.add("Center", statusLine);
progressBar = new ProgressBar(120, 20);
progressBar.addKeyListener(this);
progressBar.addMouseListener(this);
statusBar.add("East", progressBar);
statusBar.setSize(toolbar.getPreferredSize());
add(statusBar);
IJ.init(this, applet);
addKeyListener(this);
addWindowListener(this);
setFocusTraversalKeysEnabled(false);
Point loc = getPreferredLocation();
Dimension tbSize = toolbar.getPreferredSize();
int ijWidth = tbSize.width + 10;
int ijHeight = 100;
setCursor(Cursor.getDefaultCursor()); // work-around for JDK 1.1.8 bug
if (mode != NO_SHOW) {
if (IJ.isWindows())
try {
setIcon();
} catch (Exception e) {
}
setBounds(loc.x, loc.y, ijWidth, ijHeight); // needed for pack to work
setLocation(loc.x, loc.y);
pack();
setResizable(!(IJ.isMacintosh() || IJ.isWindows())); // make resizable on Linux
show();
}
if (err1 != null) IJ.error(err1);
if (err2 != null) {
IJ.error(err2);
IJ.runPlugIn("ij.plugin.ClassChecker", "");
}
m.installStartupMacroSet();
if (IJ.isMacintosh() && applet == null) {
Object qh = null;
qh = IJ.runPlugIn("MacAdapter", "");
if (qh == null) IJ.runPlugIn("QuitHandler", "");
}
if (applet == null) IJ.runPlugIn("ij.plugin.DragAndDrop", "");
String str = m.getMacroCount() == 1 ? " macro" : " macros";
IJ.showStatus(version() + m.getPluginCount() + " commands; " + m.getMacroCount() + str);
// if (applet==null && !embedded && Prefs.runSocketListener)
// new SocketListener();
configureProxy();
if (applet == null) loadCursors();
}
public void dragExit(DropTargetEvent e) {
IJ.showStatus("");
}
public void dragEnter(DropTargetDragEvent e) {
IJ.showStatus("<<Drag and Drop>>");
if (IJ.debugMode) IJ.log("DragEnter: " + e.getLocation());
e.acceptDrag(DnDConstants.ACTION_COPY);
openAsVirtualStack = false;
}