private void checkCrop( final Envelope envelope, final ImagePlus srcPixelProvider, final ImagePlus dstPixelProvider) throws IOException { // check metadata final RasterMetadata dstMetadata = geoRasterDst.getMetadata(); final float pixelSize_X = dstMetadata.getPixelSize_X(); final float pixelSize_Y = dstMetadata.getPixelSize_Y(); final float halfPixelSize_X = dstMetadata.getPixelSize_X() / 2; final float halfPixelSize_Y = Math.abs(dstMetadata.getPixelSize_Y()) / 2; final int ncols = dstMetadata.getNCols(); final int nrows = dstMetadata.getNRows(); final double xulcorner = dstMetadata.getXulcorner(); final double yulcorner = dstMetadata.getYulcorner(); assertTrue(pixelSize_X == geoRasterDst.getMetadata().getPixelSize_X()); assertTrue(ncols * pixelSize_X == dstMetadata.getEnvelope().getWidth()); assertTrue( xulcorner - (pixelSize_X / 2) + ncols * pixelSize_X == dstMetadata.getEnvelope().getMaxX()); assertTrue(pixelSize_Y == geoRasterDst.getMetadata().getPixelSize_Y()); assertTrue(nrows * Math.abs(pixelSize_Y) == dstMetadata.getEnvelope().getHeight()); assertTrue( yulcorner - (pixelSize_Y / 2) - nrows * Math.abs(pixelSize_Y) == dstMetadata.getEnvelope().getMinY()); // check raster values for (double y = envelope.getMinY() + halfPixelSize_Y; y < envelope.getMaxY(); y = y + 1) { for (double x = envelope.getMinX() + halfPixelSize_X; x < envelope.getMaxX(); x = x + 1) { final Point2D srcPixel = geoRasterSrc.fromRealWorldToPixel(x, y); final Point2D dstPixel = geoRasterDst.fromRealWorldToPixel(x, y); final float srcPixelValue = srcPixelProvider .getProcessor() .getPixelValue((int) srcPixel.getX(), (int) srcPixel.getY()); final float dstPixelValue = dstPixelProvider .getProcessor() .getPixelValue((int) dstPixel.getX(), (int) dstPixel.getY()); if (Float.isNaN(srcPixelValue)) { assertTrue(Float.isNaN(dstPixelValue)); } else { assertTrue("pixel[" + x + ", " + y + "]", srcPixelValue == dstPixelValue); } } } }