我感兴趣的主要领域是生物医学图像处理。我有一套生物医学图像,我需要找到它们之间的相位。
为了计算图像之间的相位,我使用相位相关,但我不确定我是如何做到的。我正在使用“Java Advanced Imaging”库,我的代码是:
BufferedImage first = grayscale(first);
BufferedImage second = grayscale(second);
RenderedOp dftFirst = DFTDescriptor.create(PlanarImage.wrapRenderedImage(first), DFTDescriptor.SCALING_NONE, DFTDescriptor.REAL_TO_COMPLEX, null);
RenderedOp dftSecond = DFTDescriptor.create(PlanarImage.wrapRenderedImage(second), DFTDescriptor.SCALING_NONE, DFTDescriptor.REAL_TO_COMPLEX, null);
RenderedOp conj = ConjugateDescriptor.create(dftSecond, null);
RenderedOp conv = MultiplyComplexDescriptor.create(dftFirst, conj, null);
RenderedOp abs = AbsoluteDescriptor.create(conv, null);
RenderedOp R = DivideComplexDescriptor.create(conv, abs, null);
RenderedOp idft = IDFTDescriptor.create(R, IDFTDescriptor.SCALING_DIMENSIONS, IDFTDescriptor.COMPLEX_TO_COMPLEX, null);
RenderedOp magnitude = MagnitudeDescriptor.create(idft, null);
RenderedOp shift = PeriodicShiftDescriptor.create(magnitude, new Integer(magnitude.getWidth() / 2), new Integer(magnitude.getHeight() / 2), null);
RenderedOp opExtrema = ExtremaDescriptor.create(shift, new ROIShape(shift.getBounds()), 1, 1, Boolean.TRUE, 1, null);
int[] maxLocation = (int[])((List[])opExtrema.getProperty("maxLocations"))[0].get(0);
double[][] extrema = (double[][])opExtrema.getProperty("extrema");
double c = 255.0 / (extrema[1][0] - extrema[0][0]);
double o = -extrema[0][0] * c;
System.out.println("Maximum value " + extrema[1][0] + " found at (" + maxLocation[0] + "," + maxLocation[1] + ")");
Test images and Image result of Phase Correlation
我的问题是:
谢谢
第2点解决了
int deltaX = maxLocation[0] - (conv.getWidth() / 2);
int deltaY = maxLocation[1] - (conv.getHeight() / 2);