1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
|
package ij.plugin.filter;
import ij.*;
import ij.process.*;
import ij.gui.*;
import ij.measure.*;
import ij.plugin.filter.Analyzer;
import ij.util.Tools;
import java.awt.Rectangle;
/** Implements the Image/Stack/Plot Z-axis Profile command. */
public class ZAxisProfiler implements PlugInFilter, Measurements {
ImagePlus imp;
public int setup(String arg, ImagePlus imp) {
this.imp = imp;
return DOES_ALL+NO_CHANGES;
}
public void run(ImageProcessor ip) {
if (imp.getStackSize()<2) {
IJ.error("ZAxisProfiler", "This command requires a stack.");
return;
}
Roi roi = imp.getRoi();
if (roi!=null && roi.isLine()) {
IJ.error("ZAxisProfiler", "This command does not work with line selections.");
return;
}
double minThreshold = ip.getMinThreshold();
double maxThreshold = ip.getMaxThreshold();
float[] y = getZAxisProfile(roi, minThreshold, maxThreshold);
if (y!=null) {
float[] x = new float[y.length];
for (int i=0; i<x.length; i++)
x[i] = i+1;
String title;
if (roi!=null) {
Rectangle r = imp.getRoi().getBounds();
title = imp.getTitle()+"-"+r.x+"-"+r.y;
} else
title = imp.getTitle()+"-0-0";
Plot plot = new Plot(title, "Slice", "Mean", x, y);
double ymin = ProfilePlot.getFixedMin();
double ymax= ProfilePlot.getFixedMax();
if (!(ymin==0.0 && ymax==0.0)) {
double[] a = Tools.getMinMax(x);
double xmin=a[0]; double xmax=a[1];
plot.setLimits(xmin, xmax, ymin, ymax);
}
plot.show();
}
}
float[] getZAxisProfile(Roi roi, double minThreshold, double maxThreshold) {
ImageStack stack = imp.getStack();
int size = stack.getSize();
float[] values = new float[size];
Calibration cal = imp.getCalibration();
Analyzer analyzer = new Analyzer(imp);
int measurements = analyzer.getMeasurements();
boolean showResults = measurements!=0 && measurements!=LIMIT;
boolean showingLabels = (measurements&LABELS)!=0 || (measurements&SLICE)!=0;
measurements |= MEAN;
if (showResults) {
if (!analyzer.resetCounter())
return null;
}
int current = imp.getCurrentSlice();
for (int i=1; i<=size; i++) {
if (showingLabels) imp.setSlice(i);
ImageProcessor ip = stack.getProcessor(i);
if (minThreshold!=ImageProcessor.NO_THRESHOLD)
ip.setThreshold(minThreshold,maxThreshold,ImageProcessor.NO_LUT_UPDATE);
ip.setRoi(roi);
ImageStatistics stats = ImageStatistics.getStatistics(ip, measurements, cal);
analyzer.saveResults(stats, roi);
if (showResults)
analyzer.displayResults();
values[i-1] = (float)stats.mean;
}
if (showingLabels) imp.setSlice(current);
return values;
}
}
|
