Ticket #3627: Test.java

import java.security.MessageDigest;
import java.security.NoSuchAlgorithmException;

import loci.common.DataTools;
import loci.formats.ChannelFiller;
import loci.formats.ChannelSeparator;
import loci.formats.FormatTools;
import loci.formats.IFormatReader;
import loci.formats.ImageReader;
import loci.formats.MinMaxCalculator;


public class Test
{
  private static final String TARGET = 
    "/home/callan/testimages/lsm/2chZT.lsm";

  public void parse(String id) throws Exception
  {
    IFormatReader reader = new ImageReader();
    reader = new ChannelFiller(reader);
    reader = new ChannelSeparator(reader);
    reader = new MinMaxCalculator(reader);
    reader.setId(id);

    int sizeX = reader.getSizeX();
    int sizeY = reader.getSizeY();
    int sizeZ = reader.getSizeZ();
    int sizeC = reader.getSizeC();
    int sizeT = reader.getSizeT();
    int bpp = FormatTools.getBytesPerPixel(reader.getPixelType());
    int planeSize = sizeX * sizeY * bpp;

    byte[] plane = new byte[planeSize];
    byte[] halfPlane = new byte[planeSize / 2];
    String planeDigest, halfPlaneDigest;
    for (int i = 0; i < reader.getImageCount(); i++) {
      reader.openBytes(i, plane);
      reader.openBytes(i, halfPlane, 0, 0, sizeX, sizeY / 2);
      planeDigest = md5(plane, 0, planeSize / 2);
      halfPlaneDigest = md5(halfPlane, 0, planeSize / 2);
      if (!planeDigest.equals(halfPlaneDigest)) {
        System.err.println(String.format("MD5:%d %s != %s",
            i, planeDigest, halfPlaneDigest));
      }
    }
  }

  /** Calculate the MD5 of a byte array. */
  public static String md5(byte[] b, int offset, int len) {
    try {
      MessageDigest md = MessageDigest.getInstance("MD5");
      md.reset();
      md.update(b, offset, len);
      byte[] digest = md.digest();
      return DataTools.bytesToHex(digest);
    }
    catch (NoSuchAlgorithmException e) { }
    return null;
  }

  public static void main(String[] args) throws Exception
  {
    Test t = new Test();
    t.parse(TARGET);
  }

}