Changeset 7375

Timestamp:

12/09/10 13:15:35 (9 years ago)

Author:

melissa

Message:

Added more sophisticated image name parsing to OpenlabReader. Closes #615, closes #616.

File:

• trunk/components/bio-formats/src/loci/formats/in/OpenlabReader.java (modified) (10 diffs)

trunk/components/bio-formats/src/loci/formats/in/OpenlabReader.java

@@ -108 +108 @@
 
   private String gain, detectorOffset, xPos, yPos, zPos;
+  private boolean specialPlateNames = false;
 
   // -- Constructor --
@@ -142 +143 @@
     lastPlane = no;
 
-    if (no >= planeOffsets[series].length) return buf;
-    int index = planeOffsets[series][no];
-
-    long first = planes[index].planeOffset;
-    long last = no == planeOffsets[series].length - 1 ||
-      planeOffsets[series][no + 1] >= planes.length ? in.length() :
-      planes[planeOffsets[series][no + 1]].planeOffset;
+    PlaneInfo planeInfo = null;
+    if (specialPlateNames) {
+      planeInfo = getPlane(getZCTCoords(no));
+    }
+    else if (no < planeOffsets[series].length) {
+      int index = planeOffsets[series][no];
+      planeInfo = planes[index];
+    }
+    if (planeInfo == null) return buf;
+
+    long first = planeInfo.planeOffset;
+    long last = first + FormatTools.getPlaneSize(this) * 2;
     in.seek(first);
 
     int bpp = FormatTools.getBytesPerPixel(getPixelType());
 
-    if (!planes[index].pict) {
+    if (!planeInfo.pict) {
       if (version == 2) {
         readPlane(in, x, y, w, h, buf);
@@ -188 +194 @@
         b = null;
       }
-      if (planes[index].volumeType == MAC_256_GREYS ||
-        planes[index].volumeType == MAC_256_COLORS)
+      if (planeInfo.volumeType == MAC_256_GREYS ||
+        planeInfo.volumeType == MAC_256_COLORS)
       {
         for (int i=0; i<buf.length; i++) {
@@ -229 +235 @@
         r.close();
         LOGGER.debug("", exc);
-        in.seek(planes[index].planeOffset - 298);
+        in.seek(planeInfo.planeOffset - 298);
 
         if (in.readByte() == 1) in.skipBytes(128);
@@ -298 +304 @@
       gain = detectorOffset = null;
       xPos = yPos = zPos = null;
+      specialPlateNames = false;
     }
   }
@@ -684 +691 @@
     core[s].dimensionOrder = "XY";
     for (PlaneInfo plane : planes) {
-      if (plane == null) continue;
-      if (plane.series == s) {
-        String name = plane.planeName;
+      if (plane == null || plane.series != s) continue;
+
+      String name = plane.planeName;
+
+      // check for a specific name format:
+      // <channel name><optional timepoint>_<plate>_<well>_<Z section>
+
+      String[] tokens = name.split("_");
+      if (tokens.length == 4) {
+        specialPlateNames = true;
+
+        if (!uniqueZ.contains(tokens[3])) {
+          uniqueZ.add(tokens[3]);
+        }
+        plane.channelName = tokens[0];
+        int endIndex = 0;
+        while (endIndex < plane.channelName.length() &&
+          !Character.isDigit(plane.channelName.charAt(endIndex)))
+        {
+          endIndex++;
+        }
+        String timepoint = plane.channelName.substring(endIndex);
+        if (timepoint.equals("")) timepoint = "1";
+        plane.channelName = plane.channelName.substring(0, endIndex);
+
+        if (!uniqueC.contains(plane.channelName)) {
+          uniqueC.add(plane.channelName);
+        }
+        if (!uniqueT.contains(timepoint)) {
+          uniqueT.add(timepoint);
+        }
+
+        core[s].dimensionOrder = "XYCTZ";
+        plane.wavelength = uniqueC.indexOf(plane.channelName);
+        plane.timepoint = uniqueT.indexOf(timepoint);
+        plane.zPosition = uniqueZ.indexOf(tokens[3]);
+      }
+      else {
         for (String axis : axes) {
           Vector<String> unique = null;
@@ -717 +759 @@
     }
 
+    if (specialPlateNames) {
+      core[s].sizeC *= uniqueC.size();
+      core[s].sizeT = uniqueT.size();
+      core[s].sizeZ = uniqueZ.size();
+      core[s].imageCount = core[s].sizeC * core[s].sizeZ * core[s].sizeT;
+      return;
+    }
+
     if (core[s].rgb && uniqueC.size() <= 1) {
       core[s].dimensionOrder = core[s].dimensionOrder.replaceAll("C", "");
@@ -750 +800 @@
   }
 
+  private PlaneInfo getPlane(int[] zct) {
+    for (PlaneInfo plane : planes) {
+      if (plane != null && plane.zPosition == zct[0] &&
+        plane.wavelength == zct[1] && plane.timepoint == zct[2] &&
+        plane.series == getSeries())
+      {
+        return plane;
+      }
+    }
+    return null;
+  }
+
   // -- Helper classes --
 
@@ -757 +819 @@
     protected int zPosition;
     protected int wavelength;
+    protected int timepoint;
     protected String planeName;
     protected long timestamp;
@@ -765 +828 @@
     protected int height;
     protected int series = -1;
+    protected String channelName;
   }