2

我的主要兴趣领域是生物医学图像处理。我有一组生物医学图像，我需要在它们之间找到相位。

对于计算图像之间的相位，我使用相位相关，但我不确定我是如何做到的。我正在使用“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] + ")");

相位相关的测试图像和图像结果

我的问题是：

  1. JAI中相位相关的步骤是否正确？
  2. 如何从 JAI 中的图像结果中获取相位大小？

谢谢

更新

第 2 点。解决了

int deltaX = maxLocation[0] - (conv.getWidth() / 2);
int deltaY = maxLocation[1] - (conv.getHeight() / 2);
4

0 回答 0