protected static double[] createArray(final Interval interval, final double sigma) {
final double[] sigmas = new double[interval.numDimensions()];
for (int d = 0; d < interval.numDimensions(); ++d) sigmas[d] = sigma;
return sigmas;
}
private long numPixels(final Interval i) {
assert i != null;
long[] dims = new long[i.numDimensions()];
i.dimensions(dims);
long acc = 1;
for (long l : dims) {
acc *= l;
}
return acc;
}
/**
* Computes a Gaussian convolution with the precision of the type provided on an infinite {@link
* RandomAccessible}
*
* @param sigma - the sigma for the convolution
* @param img - the img {@link RandomAccessible} (infinite -> Views.extend( ... ) )
* @param interval - the interval which should be convolved
* @param output - the output {@link RandomAccessible} (img and output can be the same)
* @param origin - the origin in the output where the result should be placed
* @param imgFactory - the {@link ImgFactory} for T which is needed for temporary imgs
*/
@SuppressWarnings("rawtypes")
public static <T extends NumericType<T>> void inNumericType(
final double[] sigma,
final RandomAccessible<T> img,
final Interval interval,
final RandomAccessible<T> output,
final Localizable origin,
final ImgFactory<T> imgFactory) {
// find out if it is a NativeType, therefore we get the first value
// defined by the
// interval in the img - we also need it in this case for the
// convolution
final long[] tmpCoordinate = new long[img.numDimensions()];
interval.min(tmpCoordinate);
final RandomAccess<T> tmp = img.randomAccess();
tmp.setPosition(tmpCoordinate);
if (NativeType.class.isInstance(tmp.get())) {
// we have to call it in an extra method because we cannot cast to
// (NumericType<T> & NativeType<T>)
computeInNativeType(
sigma,
(RandomAccessible) img,
interval,
(RandomAccessible) output,
origin,
(ImgFactory) imgFactory,
tmp.get());
} else {
final GaussGeneral<T> gauss =
new GaussGeneral<T>(sigma, img, interval, output, origin, imgFactory, tmp.get());
gauss.call();
}
}
/**
* Computes a Gaussian convolution with double precision on an infinite {@link RandomAccessible}
*
* @param sigma - the sigma for the convolution
* @param img - the img {@link RandomAccessible} (infinite -> Views.extend( ... ) )
* @param interval - the interval which should be convolved
* @param output - the output {@link RandomAccessible} (img and output can be the same)
* @param origin - the origin in the output where the result should be placed
* @param imgFactory - the {@link ImgFactory} for {@link DoubleType} which is needed for temporary
* imgs
*/
public static <T extends RealType<T>> void inDouble(
final double[] sigma,
final RandomAccessible<T> img,
final Interval interval,
final RandomAccessible<T> output,
final Localizable origin,
final ImgFactory<DoubleType> imgFactory) {
// find out if it is a DoubleType, therefore we get the first value
// defined by the
// interval in the img
final long[] tmpCoordinate = new long[img.numDimensions()];
interval.min(tmpCoordinate);
final RandomAccess<T> tmp = img.randomAccess();
tmp.setPosition(tmpCoordinate);
if (DoubleType.class.isInstance(tmp.get())) {
@SuppressWarnings({"rawtypes", "unchecked"})
final RandomAccessible<DoubleType> rIn = (RandomAccessible) img;
@SuppressWarnings({"rawtypes", "unchecked"})
final RandomAccessible<DoubleType> rOut = (RandomAccessible) output;
new GaussDouble(sigma, rIn, interval, rOut, origin, imgFactory).call();
} else {
final RandomAccessible<DoubleType> rIn =
new WriteConvertedRandomAccessible<T, DoubleType>(
img, new RealDoubleSamplerConverter<T>());
final RandomAccessible<DoubleType> rOut =
new WriteConvertedRandomAccessible<T, DoubleType>(
output, new RealDoubleSamplerConverter<T>());
new GaussDouble(sigma, rIn, interval, rOut, origin, imgFactory).call();
}
}
/**
* Test whether the {@code containing} interval completely contains the {@code contained}
* interval.
*
* <p>TODO: move to {@link Intervals}
*
* @param containing
* @param contained
* @return
*/
public static final boolean contains(final Interval containing, final Interval contained) {
assert containing.numDimensions() == contained.numDimensions();
final int n = containing.numDimensions();
for (int d = 0; d < n; ++d) {
if (containing.min(d) > contained.min(d) || containing.max(d) < contained.max(d))
return false;
}
return true;
}
public static String printInterval(final Interval interval) {
String out = "(Interval empty)";
if (interval == null || interval.numDimensions() == 0) return out;
out = "[" + interval.min(0);
for (int i = 1; i < interval.numDimensions(); i++) out += ", " + interval.min(i);
out += "] -> [" + interval.max(0);
for (int i = 1; i < interval.numDimensions(); i++) out += ", " + interval.max(i);
out += "], dimensions (" + interval.dimension(0);
for (int i = 1; i < interval.numDimensions(); i++) out += ", " + interval.dimension(i);
out += ")";
return out;
}
/**
* Walk through an Interval on an Img using a LocalizingCursor, localizing on every step.
*
* @param img
* @param interval
*/
protected static void localizingWalkThrough(final Img<IntType> img, final Interval interval) {
final Cursor<IntType> c = Views.interval(img, interval).localizingCursor();
final long[] pos = new long[interval.numDimensions()];
int i = 0;
while (c.hasNext()) {
c.fwd();
i += c.get().get();
c.localize(pos);
}
j = (int) pos[0] + i;
}
/**
* Create a long[] with the min coordinates of an {@link Interval}.
*
* <p>Keep in mind that creating arrays wildly is not good practice and consider using the
* interval directly.
*
* @param interval
* @return dimensions of the interval as a new long[]
*/
public static final long[] intervalMin(final Interval interval) {
final long[] min = new long[interval.numDimensions()];
interval.min(min);
return min;
}
/**
* Create a long[] with the max coordinates of an {@link Interval}.
*
* <p>Keep in mind that creating arrays wildly is not good practice and consider using the
* interval directly.
*
* @param interval
* @return dimensions of the interval as a new long[]
*/
public static final long[] intervalMax(final Interval interval) {
final long[] max = new long[interval.numDimensions()];
interval.max(max);
return max;
}
/**
* Create a long[] with the dimensions of an {@link Interval}.
*
* <p>Keep in mind that creating arrays wildly is not good practice and consider using the
* interval directly.
*
* @param interval
* @return dimensions of the interval as a new long[]
*/
public static final long[] intervalDimensions(final Interval interval) {
final long[] dimensions = new long[interval.numDimensions()];
interval.dimensions(dimensions);
return dimensions;
}
private long[] resultDims(final Interval src) {
long[] dims = new long[src.numDimensions()];
src.dimensions(dims);
return dims;
}
@Override
public int numDimensions() {
return interval.numDimensions();
}
/**
* Get a "good" initial viewer transform. The viewer transform is chosen such that for the first
* source,
*
* <ul>
* <li>the XY plane is aligned with the screen plane,
* <li>the <em>z = dim_z / 2</em> slice is shown,
* <li>centered and scaled such that the full <em>dim_x</em> by <em>dim_y</em> is visible.
* </ul>
*
* @param viewerWidth width of the viewer display
* @param viewerHeight height of the viewer display
* @param state the {@link ViewerState} containing at least one source.
* @return proposed initial viewer transform.
*/
public static AffineTransform3D initTransform(
final int viewerWidth,
final int viewerHeight,
final boolean zoomedIn,
final ViewerState state) {
final int cX = viewerWidth / 2;
final int cY = viewerHeight / 2;
final Source<?> source = state.getSources().get(state.getCurrentSource()).getSpimSource();
final int timepoint = state.getCurrentTimepoint();
if (!source.isPresent(timepoint)) return new AffineTransform3D();
final AffineTransform3D sourceTransform = new AffineTransform3D();
source.getSourceTransform(timepoint, 0, sourceTransform);
final Interval sourceInterval = source.getSource(timepoint, 0);
final double sX0 = sourceInterval.min(0);
final double sX1 = sourceInterval.max(0);
final double sY0 = sourceInterval.min(1);
final double sY1 = sourceInterval.max(1);
final double sZ0 = sourceInterval.min(2);
final double sZ1 = sourceInterval.max(2);
final double sX = (sX0 + sX1 + 1) / 2;
final double sY = (sY0 + sY1 + 1) / 2;
final double sZ = (sZ0 + sZ1 + 1) / 2;
final double[][] m = new double[3][4];
// rotation
final double[] qSource = new double[4];
final double[] qViewer = new double[4];
Affine3DHelpers.extractApproximateRotationAffine(sourceTransform, qSource, 2);
LinAlgHelpers.quaternionInvert(qSource, qViewer);
LinAlgHelpers.quaternionToR(qViewer, m);
// translation
final double[] centerSource = new double[] {sX, sY, sZ};
final double[] centerGlobal = new double[3];
final double[] translation = new double[3];
sourceTransform.apply(centerSource, centerGlobal);
LinAlgHelpers.quaternionApply(qViewer, centerGlobal, translation);
LinAlgHelpers.scale(translation, -1, translation);
LinAlgHelpers.setCol(3, translation, m);
final AffineTransform3D viewerTransform = new AffineTransform3D();
viewerTransform.set(m);
// scale
final double[] pSource = new double[] {sX1 + 0.5, sY1 + 0.5, sZ};
final double[] pGlobal = new double[3];
final double[] pScreen = new double[3];
sourceTransform.apply(pSource, pGlobal);
viewerTransform.apply(pGlobal, pScreen);
final double scaleX = cX / pScreen[0];
final double scaleY = cY / pScreen[1];
final double scale;
if (zoomedIn) scale = Math.max(scaleX, scaleY);
else scale = Math.min(scaleX, scaleY);
viewerTransform.scale(scale);
// window center offset
viewerTransform.set(viewerTransform.get(0, 3) + cX, 0, 3);
viewerTransform.set(viewerTransform.get(1, 3) + cY, 1, 3);
return viewerTransform;
}
