16位数据的双线性插值的Java错误

我有一个使用双线性插值16位数据的问题.我有两个图像,orig Image和display Image.我想使用AffineTransformOp通过AffineTransform将origImage过滤到displayImage,这是显示区域的大小. origImage的类型为BufferedImage.TYPE_USHORT_GRAY,并且具有sun.awt.image.ShortInterleavedRaster类型的栅格.这是我现在的代码

displayImage = new BufferedImage(getWidth(),getHeight(),origImage.getType());
try {
    op = new AffineTransformOp(atx,AffineTransformOp.TYPE_BILINEAR);
    op.filter(origImage,displayImage);
}
catch (Exception e) {
    e.printStackTrace();
}

为了显示错误,我创建了2个渐变图像.一个具有15位范围(最大32767)和16位范围(最大65535)中的值.以下是两张图片

15位图像

16位图像

这两个图像是以相同的时尚创建的,应该看起来是一样的,但注意到16位图像中间的行.起初我以为这是一个溢出的问题,但是它是奇怪的,它表现在梯度的中心,而不是在像素值越高的尽头.另外,如果这是一个溢出问题,我怀疑15位图像也会受到影响.

对此的任何帮助将不胜感激.

我只是想知道为什么没有人在回答,我提供了足够的信息吗？需要更多信息吗？

以下是我用于生成AffineTransform的代码.所有引用的变量都是基于用户输入(鼠标移动)计算的,应该是正确的(它已经被包括我在内的很多人所测试).希望这可以帮助错误.

AffineTransform panTranslate = new AffineTransform();
panTranslate.translate(imagePanOffset.x,imagePanOffset.y);

AffineTransform rotateCenterTranslate = new AffineTransform();
rotateCenterTranslate.translate(imageRotateCTR.x,imageRotateCTR.y);
AffineTransform rotateTransform = new AffineTransform();
rotateTransform.rotate(Math.toRadians(rotateValue));
AffineTransform rotateAntiCenterTranslate = new AffineTransform();
rotateAntiCenterTranslate.translate(-imageRotateCTR.x,-imageRotateCTR.y);

AffineTransform translateTransform = new AffineTransform();
translateTransform.translate(imageMagOffset.x,imageMagOffset.y);

AffineTransform flipMatrixTransform = new AffineTransform();

switch (flipState) {
    case ENV.FLIP_NORMAL: // NORMAL
        break;

    case ENV.FLIP_TOP_BOTTOM: // FLIP
        flipMatrixTransform.scale(1.0,-1.0);
        flipMatrixTransform.translate(0.0,-h);
        break;

    case ENV.FLIP_LEFT_RIGHT: // MIRROR
        flipMatrixTransform.scale(-1.0,1.0);
        flipMatrixTransform.translate(-w,0.0);
        break;

    case ENV.FLIP_TOP_BOTTOM_LEFT_RIGHT: // FLIP+MIRROR
        flipMatrixTransform.scale(-1.0,-1.0);
        flipMatrixTransform.translate(-w,-h);
        break;
}

scaleTransform = new AffineTransform();
scaleTransform.scale(magFactor,magFactor);

AffineTransform atx = new AffineTransform();
atx.concatenate(panTranslate);
atx.concatenate(rotateCenterTranslate);
atx.concatenate(rotateTransform);
atx.concatenate(rotateAntiCenterTranslate);
atx.concatenate(translateTransform);
atx.concatenate(flipMatrixTransform);
atx.concatenate(scaleTransform);

我仍然不知道这里发生了什么.我非常感谢可以提供的任何帮助.我还附上了一个我遇到的真实图像中发生的错误的例子,以供更多参考.

这是在手的X射线发生的错误

这是一个缩放版本,专注于拇指和第一个手指之间的区域.

请再次注意,在非常白色的区域,但是在动态范围中间的值,就像渐变图像中的错误一样.

我发现了更多的信息.我正在调整一些变换,发现如果我只是过滤一个单位矩阵,那么这个bug就不会发生.如果我翻译一个整数,它也不会发生.如果我用一个非整数的数字翻译,它会发生.如果我放大除1以外的任何数量(整数或不整数),也会发生这种情况.希望这有帮助.

经过更多的实验,该bug绝对表现在最大强度的一半之间的边界像素(65535/2 = 32767.5).它也只发生在此值.我希望这可能有助于诊断！

根据AlBlue的要求,这里是完全独立于我的应用程序的代码,可以产生错误.请注意,在原始帖子中,我包括使用以下代码生成的图像渐变,但是我放大了其中一个渐变以更好地显示效果.您应该在0.5个翻译图像上看到效果四次,而不是在其他两个图像中.还要注意,除了1以外的任何数量的缩放,这个错误出现.只需用AffineTransform.getScaleInstance(0.9,0.9)替换AffineTransform.getTranslateInstance()即可查看错误.

private static class MyJPanel extends JPanel {
    BufferedImage displayImage = null;
    public MyJPanel(double translateValue) {
        super();
        BufferedImage bi = new BufferedImage(1024,1024,BufferedImage.TYPE_USHORT_GRAY);

        int dataRange = (int)Math.pow(2,16);
        double step = dataRange/(bi.getRaster().getDataBuffer().getSize()/4.0);
        double value = 0;
        for (int i=0; i<bi.getRaster().getDataBuffer().getSize(); i++) {
            bi.getRaster().getDataBuffer().setElem(i,(int)value);
            if (value >= dataRange)
                value = 0;
            else
                value += step;
        }
        displayImage = new BufferedImage(bi.getWidth(),bi.getHeight(),bi.getType());
        AffineTransform tx = AffineTransform.getTranslateInstance(translateValue,translateValue);
        AffineTransformOp op = new AffineTransformOp(tx,AffineTransformOp.TYPE_BILINEAR);
        op.filter(bi,displayImage);
    }

    public void paint(Graphics g) {
        super.paint(g);
        g.drawImage(displayImage,this);
    }
}

private static void showDisplayError() {
    JDialog dialog1 = new JDialog();
    dialog1.setTitle("No Translation");
    MyJPanel panel1 = new MyJPanel(0);
    dialog1.getContentPane().add(panel1);
    dialog1.setSize(1024,1024);
    dialog1.setVisible(true);

    JDialog dialog2 = new JDialog();
    dialog2.setTitle("Translation of 0.5");
    MyJPanel panel2 = new MyJPanel(0.5);
    dialog2.getContentPane().add(panel2);
    dialog2.setSize(1024,1024);
    dialog2.setVisible(true);

    JDialog dialog3 = new JDialog();
    dialog3.setTitle("Translation of 1.0");
    MyJPanel panel3 = new MyJPanel(1.0);
    dialog3.getContentPane().add(panel3);
    dialog3.setSize(1024,1024);
    dialog3.setVisible(true);
}

作为另一个更新,我刚刚在Fedora 10上尝试过,看到bug仍然存在.