public void computeData() {
computing = true;
cancelCompute = false;
try {
maxData = new double[sequence.getSizeC()];
minData = new double[sequence.getSizeC()];
drawPolygon = new Polygon[sequence.getSizeC()];
if (paintMode == PaintMode.line) {
Line2D line = ((ROI2DLine) guideRoi).getLine();
dataCount = (int) line.getP1().distance(line.getP2());
} else {
dataCount = sequence.getSizeZ();
}
dataArr = new ArrayList<double[]>();
for (int component = 0; component < sequence.getSizeC(); component++) {
double[] data = new double[dataCount];
dataArr.add(data);
}
if (paintMode == PaintMode.line) {
Line2D line = ((ROI2DLine) guideRoi).getLine();
dataCount = (int) line.getP1().distance(line.getP2());
ShapeUtil.consumeShapeFromPath(
((ROI2DShape) guideRoi).getPathIterator(null),
new ShapeConsumer() {
@Override
public boolean consume(Shape shape) {
if (shape instanceof Line2D) {
Line2D line = (Line2D) shape;
Point2D Lp;
Point2D Rp;
if (line.getX2() > line.getX1()) {
Lp = line.getP1();
Rp = line.getP2();
} else {
Lp = line.getP2();
Rp = line.getP1();
}
for (int component = 0; component < sequence.getSizeC(); component++) {
// create histo data
int distance = dataCount;
double vx = (Rp.getX() - Lp.getX()) / distance;
double vy = (Rp.getY() - Lp.getY()) / distance;
double[] data = dataArr.get(component);
double x = Lp.getX();
double y = Lp.getY();
IcyBufferedImage image = canvas.getCurrentImage();
if (image.isInside((int) x, (int) y)) {
maxData[component] =
Array1DUtil.getValue(
image.getDataXY(component),
image.getOffset((int) x, (int) y),
image.isSignedDataType());
} else {
maxData[component] = 0;
}
minData[component] = maxData[component];
for (int i = 0; i < dataCount; i++) {
if (cancelCompute) break;
if (image.isInside((int) x, (int) y)) {
data[i] =
Array1DUtil.getValue(
image.getDataXY(component),
image.getOffset((int) x, (int) y),
image.isSignedDataType());
} else {
data[i] = 0;
}
if (data[i] > maxData[component]) maxData[component] = data[i];
if (data[i] < minData[component]) minData[component] = data[i];
x += vx;
y += vy;
}
Polygon polygon = new Polygon();
polygon.addPoint(0, 0);
for (int i = 0; i < dataCount; i++) {
polygon.addPoint(i, (int) (data[i] - minData[component]));
}
polygon.addPoint(dataCount, 0);
drawPolygon[component] = polygon;
}
}
return true; // continue
}
});
} else {
for (int component = 0; component < sequence.getSizeC(); component++) {
double[] data = dataArr.get(component);
if (paintMode == PaintMode.point) {
Point p = guideRoi.getPosition();
if (p.x < sequence.getSizeX() && p.y < sequence.getSizeY()) {
maxData[component] = sequence.getData(0, 0, component, p.y, p.x);
minData[component] = maxData[component];
for (int i = 0; i < dataCount; i++) {
if (cancelCompute) break;
data[i] = sequence.getData(0, i, component, p.y, p.x);
if (data[i] > maxData[component]) maxData[component] = data[i];
if (data[i] < minData[component]) minData[component] = data[i];
}
}
} else {
maxData[component] = ROIUtil.getMeanIntensity(sequence, guideRoi, 0, -1, component);
;
minData[component] = maxData[component];
for (int i = 0; i < dataCount; i++) {
if (cancelCompute) break;
data[i] = ROIUtil.getMeanIntensity(sequence, guideRoi, i, -1, component);
if (data[i] > maxData[component]) maxData[component] = data[i];
if (data[i] < minData[component]) minData[component] = data[i];
}
}
Polygon polygon = new Polygon();
polygon.addPoint(0, 0);
for (int i = 0; i < dataCount; i++)
// pity polygon does not support this with double...
polygon.addPoint(i, (int) (data[i] - minData[component]));
polygon.addPoint(dataCount, 0);
drawPolygon[component] = polygon;
}
}
} catch (Exception e) {
System.out.print(e);
}
computing = false;
}
void drawHisto(ROI2D roi, Graphics2D g, Sequence sequence, final IcyCanvas canvas) {
if (!roiPairDict.containsKey(roi)) return;
IntensityPaint ip = roiPairDict.get(roi);
String currentValueX = "";
String currentValueV = "";
String maxValue = "";
String minValue = "";
for (int component = 0; component < sequence.getSizeC(); component++) {
AffineTransform originalTransform = g.getTransform();
g.setColor(new Color(236, 10, 170));
if (sequence.getSizeC() != 1) {
if (component == 0) g.setColor(Color.red);
if (component == 1) g.setColor(Color.green);
if (component == 2) g.setColor(Color.blue);
}
Rectangle2D rectBox = ((ROI2DRectangle) ip.displayRectangle).getRectangle();
Rectangle2D polyBox = ip.drawPolygon[component].getBounds2D();
try {
if (ip.paintMode == PaintMode.line) {
Line2D line = ((ROI2DLine) roi).getLine();
Point2D Lp;
Point2D Rp;
if (line.getX2() > line.getX1()) {
Lp = line.getP1();
Rp = line.getP2();
} else {
Lp = line.getP2();
Rp = line.getP1();
}
int pos;
if (Math.min(line.getX1(), line.getX2()) >= cursorPos.x) pos = 0;
else if (Math.max(line.getX1(), line.getX2()) <= cursorPos.x) pos = ip.dataCount;
else {
pos = (int) ((cursorPos.x - Lp.getX()) / (Rp.getX() - Lp.getX()) * ip.dataCount);
try {
currentValueX = String.format("X:%.1f", cursorPos.x);
currentValueV += String.format("%.1f ", ip.dataArr.get(component)[pos]);
} catch (Exception e2) {
}
}
ip.cursor1.setLine(pos, 0, pos, polyBox.getHeight());
} else {
int pos = (int) cursorPos.z;
ip.cursor1.setLine(pos, 0, pos, polyBox.getHeight());
try {
currentValueX = String.format("Z:%.1f", cursorPos.z);
currentValueV += String.format("%.1f ", ip.dataArr.get(component)[pos]);
} catch (Exception e2) {
}
}
double sx = rectBox.getWidth() / polyBox.getWidth();
double sy = rectBox.getHeight() / polyBox.getHeight();
if (sx < 100 && sy < 100) {
g.translate(rectBox.getMinX(), rectBox.getMaxY());
g.scale(sx, -sy);
g.draw(ip.drawPolygon[component]);
g.setColor(new Color(100, 100, 170));
g.draw(ip.cursor1);
g.setColor(new Color(236, 10, 170));
} else {
char[] c = "Exceeding display limit!".toCharArray();
g.drawChars(c, 0, c.length, (int) rectBox.getCenterX() - 10, (int) rectBox.getCenterY());
}
} finally {
g.setTransform(originalTransform);
// min,max
double xStart, xEnd;
if (ip.paintMode == PaintMode.line) {
Line2D line = ((ROI2DLine) roi).getLine();
Point2D Lp;
Point2D Rp;
if (line.getX2() > line.getX1()) {
Lp = line.getP1();
Rp = line.getP2();
} else {
Lp = line.getP2();
Rp = line.getP1();
}
xStart = Lp.getX();
xEnd = Rp.getX();
int pos;
double yp;
if (Math.min(line.getX1(), line.getX2()) >= cursorPos.x) {
pos = (int) Lp.getX();
yp = Lp.getY();
} else if (Math.max(line.getX1(), line.getX2()) <= cursorPos.x) {
pos = (int) Rp.getX();
yp = Rp.getY();
} else {
pos = (int) cursorPos.x;
yp =
(cursorPos.x - Lp.getX())
/ (line.getX2() - line.getX1())
* (line.getY2() - line.getY1())
+ Lp.getY();
}
ip.cursor2.setLine(pos, yp + 10, pos, yp - 10);
g.draw(ip.cursor2);
} else {
xStart = 0;
xEnd = ip.dataCount;
}
maxValue += String.format("%.1f ", ip.maxData[component]);
minValue += String.format("%.1f ", ip.minData[component]);
if (component == sequence.getSizeC() - 1) {
char[] c = String.format("%.1f", xStart).toCharArray();
// x1
g.drawChars(c, 0, c.length, (int) rectBox.getMinX(), (int) rectBox.getMaxY() + 30);
c = String.format("%.1f", xEnd).toCharArray();
// x2
g.drawChars(c, 0, c.length, (int) rectBox.getMaxX(), (int) rectBox.getMaxY() + 30);
c = maxValue.toCharArray();
g.drawChars(c, 0, c.length, (int) rectBox.getMaxX() + 10, (int) rectBox.getMinY() + 10);
c = minValue.toCharArray();
g.drawChars(c, 0, c.length, (int) rectBox.getMaxX() + 10, (int) rectBox.getMaxY() - 5);
c = currentValueX.toCharArray();
g.drawChars(
c,
0,
c.length,
(int) (rectBox.getMinX() + (ip.cursor1.x1 / ip.dataCount) * rectBox.getWidth()) - 20,
(int) rectBox.getMaxY() + 15);
c = currentValueV.toCharArray();
g.drawChars(
c,
0,
c.length,
(int) (rectBox.getMinX() + (ip.cursor1.x1 / ip.dataCount) * rectBox.getWidth()) - 20,
(int) rectBox.getMinY() - 5);
}
}
}
}
/**
* Save the specified sequence in the specified file.<br>
* When the sequence contains severals image the multipleFile flag is used to indicate<br>
* if images are saved as separate files (file then specify a directory) or not.<br>
* zMin - zMax and tMin - tMax define the Z and T images range to save.<br>
*
* @param formatWriter writer used to save sequence (define the image format)
* @param sequence sequence to save
* @param filePath file name where we want to save sequence
* @param zMin start Z position to save
* @param zMax end Z position to save
* @param tMin start T position to save
* @param tMax end T position to save
* @param fps frame rate for AVI sequence save
* @param saveFrame progress frame for save operation (can be null)
* @throws ServiceException
* @throws IOException
* @throws FormatException
*/
private static void save(
IFormatWriter formatWriter,
Sequence sequence,
String filePath,
int zMin,
int zMax,
int tMin,
int tMax,
int fps,
FileFrame saveFrame)
throws ServiceException, FormatException, IOException {
final File file = new File(filePath);
final IFormatWriter writer;
if (formatWriter == null) writer = getWriter(file, ImageFileFormat.TIFF);
else writer = formatWriter;
// TODO: temporary fix for the "incorrect close operation" bug in Bio-Formats
// with OME TIF writer, remove it when fixed.
// {
// try
// {
// writer = formatWriter.getClass().newInstance();
// }
// catch (Exception e)
// {
// throw new ServiceException("Can't create new writer instance: " + e);
// }
// }
if (writer == null)
throw new UnknownFormatException(
"Can't find a valid image writer for the specified file: " + filePath);
// first delete the file else LOCI won't save it correctly
if (file.exists()) file.delete();
// ensure parent directory exist
FileUtil.ensureParentDirExist(file);
final int sizeC = sequence.getSizeC();
final boolean separateChannel = getSeparateChannelFlag(writer, sequence.getColorModel());
// set settings
writer.setFramesPerSecond(fps);
// generate metadata
writer.setMetadataRetrieve(
MetaDataUtil.generateMetaData(
sequence, (zMax - zMin) + 1, (tMax - tMin) + 1, separateChannel));
// no interleave (XP default viewer want interleaved channel to correctly read image)
writer.setInterleaved(false);
// set id
writer.setId(filePath);
// init
writer.setSeries(0);
// usually give better save performance
writer.setWriteSequentially(true);
// get endianess
final boolean littleEndian =
!writer.getMetadataRetrieve().getPixelsBinDataBigEndian(0, 0).booleanValue();
byte[] data = null;
try {
int imageIndex = 0;
// XYCZT order is important here (see metadata)
for (int t = tMin; t <= tMax; t++) {
for (int z = zMin; z <= zMax; z++) {
if ((saveFrame != null) && saveFrame.isCancelRequested()) return;
final IcyBufferedImage image = sequence.getImage(t, z);
// separated channel data
if (separateChannel) {
for (int c = 0; c < sizeC; c++) {
if (image != null) {
// avoid multiple allocation
data = image.getRawData(c, data, 0, littleEndian);
writer.saveBytes(imageIndex, data);
}
imageIndex++;
}
} else {
if (image != null) {
// avoid multiple allocation
data = image.getRawData(data, 0, littleEndian);
writer.saveBytes(imageIndex, data);
}
imageIndex++;
}
if (saveFrame != null) saveFrame.incPosition();
}
}
} finally {
// always close writer after a file has been saved
writer.close();
}
}
