/** Called from io/Opener.java. */
public void run(String path) {
if (path.equals("")) return;
File theFile = new File(path);
String directory = theFile.getParent();
String fileName = theFile.getName();
if (directory == null) directory = "";
// Try and recognise file type and load the file if recognised
ImagePlus imp = openImage(directory, fileName, path);
if (imp == null) {
IJ.showStatus("");
return; // failed to load file or plugin has opened and displayed it
}
ImageStack stack = imp.getStack();
// get the title from the stack (falling back to the fileName)
String title = imp.getTitle().equals("") ? fileName : imp.getTitle();
// set the stack of this HandleExtraFileTypes object
// to that attached to the ImagePlus object returned by openImage()
setStack(title, stack);
// copy over the calibration info since it doesn't come with the ImageProcessor
setCalibration(imp.getCalibration());
// also copy the Show Info field over if it exists
if (imp.getProperty("Info") != null) setProperty("Info", imp.getProperty("Info"));
// copy the FileInfo
setFileInfo(imp.getOriginalFileInfo());
// copy dimensions
if (IJ.getVersion().compareTo("1.38s") >= 0)
setDimensions(imp.getNChannels(), imp.getNSlices(), imp.getNFrames());
if (IJ.getVersion().compareTo("1.41o") >= 0) setOpenAsHyperStack(imp.getOpenAsHyperStack());
}
/** Set display range for each channel to min-max. */
private static void displayMinToMax(ImagePlus imp) {
if (imp.isComposite()) {
for (int c = 1; c <= imp.getNChannels(); c++) {
double min = I1l.getStatsForChannel(imp, c).min;
double max = I1l.getStatsForChannel(imp, c).max;
imp.setC(c);
IJ.setMinAndMax(imp, min, max);
}
imp.setC(1);
} else {
double min = imp.getStatistics().min;
double max = imp.getStatistics().max;
IJ.setMinAndMax(imp, min, max);
}
}
/** Mask central offset / low freq spike to better use display range. */
private void maskCentralRegion(ImagePlus imp) {
int w = imp.getWidth();
double d = this.maskDiameter * w;
if (w % 2 == 0) {
d += 1; // tweak to centre the mask about zero freq stripes
}
// N.B. we assume the image is square! (FFT result)
OvalRoi mask = new OvalRoi(w / 2 - d / 2 + 1, w / 2 - d / 2 + 1, d, d);
imp.setRoi(mask);
for (int c = 1; c <= imp.getNChannels(); c++) {
imp.setC(c);
double min = I1l.getStatsForChannel(imp, c).min;
IJ.run(imp, "Set...", "value=" + min + " slice");
}
imp.setC(1);
imp.deleteRoi();
}
public void run(String arg) {
ImagePlus imp = WindowManager.getCurrentImage();
Calibration cal = imp.getCalibration();
GenericDialog gd = new GenericDialog("Options");
int subsize = 32;
gd.addNumericField("Subregion Size (pixels)?", subsize, 0);
int stepsize = 16;
gd.addNumericField("Step Size?", stepsize, 0);
int shift = 3;
gd.addNumericField("STICS temporal Shift?", shift, 0);
float xoffset = 0.0f;
gd.addNumericField("X_Offset", xoffset, 5, 15, null);
float yoffset = 0.0f;
gd.addNumericField("Y_Offset", yoffset, 5, 15, null);
float multiplier = 8.0f;
gd.addNumericField("Velocity Multiplier", multiplier, 5, 15, null);
float ftime = 1.0f;
gd.addNumericField("Frame_Time(min)", ftime, 5, 15, null);
float scaling = (float) cal.pixelWidth;
gd.addNumericField("Pixel_Size(um)", scaling, 5, 15, null);
boolean norm = true;
gd.addCheckbox("Normalize_Vector_lengths?", norm);
boolean centered = true;
gd.addCheckbox("Center_Vectors?", centered);
float magthresh = 0.0f;
gd.addNumericField("Magnitude_Threshhold?", magthresh, 5, 15, null);
int rlength = 10;
gd.addNumericField("Running_avg_length", rlength, 0);
int inc = 5;
gd.addNumericField("Start_frame_increment", inc, 0);
gd.showDialog();
if (gd.wasCanceled()) {
return;
}
subsize = (int) gd.getNextNumber();
stepsize = (int) gd.getNextNumber();
shift = (int) gd.getNextNumber();
xoffset = (float) gd.getNextNumber();
yoffset = (float) gd.getNextNumber();
multiplier = (float) gd.getNextNumber();
ftime = (float) gd.getNextNumber();
scaling = (float) gd.getNextNumber();
norm = gd.getNextBoolean();
centered = gd.getNextBoolean();
magthresh = (float) gd.getNextNumber();
rlength = (int) gd.getNextNumber();
inc = (int) gd.getNextNumber();
int width = imp.getWidth();
int xregions = 1 + (int) (((float) width - (float) subsize) / (float) stepsize);
int newwidth = xregions * subsize;
int height = imp.getHeight();
int yregions = 1 + (int) (((float) height - (float) subsize) / (float) stepsize);
int newheight = yregions * subsize;
ImageStack stack = imp.getStack();
int slices = imp.getNSlices();
int channels = imp.getNChannels();
int frames = imp.getNFrames();
if (frames == 1) {
frames = slices;
slices = 1;
}
Roi roi = imp.getRoi();
if (roi == null) {
roi = new Roi(0, 0, width, height);
}
STICS_map map = new STICS_map(subsize, stepsize);
Object[] tseries = jutils.get3DTSeries(stack, 0, 0, frames, slices, channels);
map.update_STICS_map(tseries, width, height, 0, rlength, roi.getPolygon(), shift);
FloatProcessor fp =
map.get_map(scaling, ftime, stepsize, centered, norm, multiplier, stepsize, magthresh);
ImageStack vector_stack = new ImageStack(fp.getWidth(), fp.getHeight());
vector_stack.addSlice("", fp);
float[][] vel = map.get_scaled_velocities(scaling, ftime, stepsize);
ImageStack velstack = new ImageStack(map.xregions, map.yregions);
velstack.addSlice("", vel[0]);
velstack.addSlice("", vel[1]);
int velframes = 2;
IJ.showStatus("frame " + 0 + " calculated");
for (int i = inc; i < (frames - rlength); i += inc) {
map.update_STICS_map(tseries, width, height, i, rlength, roi.getPolygon(), shift);
FloatProcessor fp2 =
map.get_map(scaling, ftime, stepsize, centered, norm, multiplier, stepsize, magthresh);
vector_stack.addSlice("", fp2);
vel = map.get_scaled_velocities(scaling, ftime, stepsize);
velstack.addSlice("", vel[0]);
velstack.addSlice("", vel[1]);
velframes += 2;
IJ.showStatus("frame " + i + " calculated");
}
(new ImagePlus("STICS Vectors", vector_stack)).show();
ImagePlus imp3 = new ImagePlus("Velocities", velstack);
imp3.setOpenAsHyperStack(true);
imp3.setDimensions(2, 1, velframes / 2);
new CompositeImage(imp3, CompositeImage.COLOR).show();
}
public void itemStateChanged(ItemEvent e) {
ImagePlus imp = WindowManager.getCurrentImage();
if (imp == null) return;
if (!imp.isComposite()) {
int channels = imp.getNChannels();
if (channels == 1 && imp.getStackSize() <= 4) channels = imp.getStackSize();
if (imp.getBitDepth() == 24 || (channels > 1 && channels < CompositeImage.MAX_CHANNELS)) {
GenericDialog gd = new GenericDialog(imp.getTitle());
gd.addMessage("Convert to multi-channel composite image?");
gd.showDialog();
if (gd.wasCanceled()) return;
else IJ.doCommand("Make Composite");
} else {
IJ.error(
"Channels",
"A composite image is required (e.g., "
+ moreLabel
+ " Open HeLa Cells),\nor create one using "
+ moreLabel
+ " Make Composite.");
return;
}
}
if (!imp.isComposite()) return;
CompositeImage ci = (CompositeImage) imp;
Object source = e.getSource();
if (source == choice) {
int index = ((Choice) source).getSelectedIndex();
switch (index) {
case 0:
ci.setMode(IJ.COMPOSITE);
break;
case 1:
ci.setMode(IJ.COLOR);
break;
case 2:
ci.setMode(IJ.GRAYSCALE);
break;
}
ci.updateAndDraw();
if (Recorder.record) {
String mode = null;
switch (index) {
case 0:
mode = "composite";
break;
case 1:
mode = "color";
break;
case 2:
mode = "grayscale";
break;
}
Recorder.record("Stack.setDisplayMode", mode);
}
} else if (source instanceof Checkbox) {
for (int i = 0; i < checkbox.length; i++) {
Checkbox cb = (Checkbox) source;
if (cb == checkbox[i]) {
if (ci.getMode() == IJ.COMPOSITE) {
boolean[] active = ci.getActiveChannels();
active[i] = cb.getState();
if (Recorder.record) {
String str = "";
for (int c = 0; c < ci.getNChannels(); c++) str += active[c] ? "1" : "0";
Recorder.record("Stack.setActiveChannels", str);
Recorder.record("//Stack.toggleChannel", imp.getChannel());
}
} else {
imp.setPosition(i + 1, imp.getSlice(), imp.getFrame());
if (Recorder.record) Recorder.record("Stack.setChannel", i + 1);
}
ci.updateAndDraw();
return;
}
}
}
}