use of loci.common.Location in project bioformats by openmicroscopy.
the class DeltavisionReader method initFile.
// -- Internal FormatReader API methods --
/* @see loci.formats.FormatReader#initFile(String) */
@Override
protected void initFile(String id) throws FormatException, IOException {
if (!checkSuffix(id, "dv")) {
if (checkSuffix(id, "dv.log") || checkSuffix(id, "r3d.log")) {
id = id.substring(0, id.lastIndexOf("."));
} else if (id.endsWith("_log.txt")) {
id = id.substring(0, id.lastIndexOf("_")) + ".dv";
}
Location file = new Location(id).getAbsoluteFile();
if (!file.exists()) {
Location dir = file.getParentFile();
String[] list = dir.list(true);
String name = file.getName();
name = name.substring(0, name.lastIndexOf("."));
for (String f : list) {
if (checkSuffix(f, "dv") && f.startsWith(name)) {
id = new Location(dir, f).getAbsolutePath();
break;
}
}
}
}
super.initFile(id);
findLogFiles();
in = new RandomAccessInputStream(currentId);
initPixels();
MetadataLevel metadataLevel = metadataOptions.getMetadataLevel();
if (metadataLevel != MetadataLevel.MINIMUM) {
initExtraMetadata();
}
}
use of loci.common.Location in project bioformats by openmicroscopy.
the class FV1000Reader method initFile.
// -- Internal FormatReader API methods --
/* @see loci.formats.FormatReader#initFile(String) */
@Override
protected void initFile(String id) throws FormatException, IOException {
super.initFile(id);
parser.setCommentDelimiter(null);
isOIB = checkSuffix(id, OIB_SUFFIX);
if (isOIB) {
initPOIService();
}
// mappedOIF is used to distinguish between datasets that are being read
// directly (e.g. using ImageJ or showinf), and datasets that are being
// imported through omebf. In the latter case, the necessary directory
// structure is not preserved (only relative file names are stored in
// OMEIS), so we will need to use slightly different logic to build the
// list of associated files.
boolean mappedOIF = !isOIB && !new File(id).getAbsoluteFile().exists();
wavelengths = new ArrayList<Double>();
illuminations = new ArrayList<String>();
channels = new ArrayList<ChannelData>();
planes = new ArrayList<PlaneData>();
String oifName = null;
if (isOIB) {
oifName = mapOIBFiles();
} else {
// make sure we have the OIF file, not a TIFF
if (!checkSuffix(id, OIF_SUFFIX)) {
currentId = findOIFFile(id);
initFile(currentId);
}
oifName = currentId;
}
String oifPath = new Location(oifName).getAbsoluteFile().getAbsolutePath();
if (mappedOIF) {
oifPath = oifName.substring(0, oifName.lastIndexOf(File.separator) + 1);
}
String path = isOIB ? "" : oifPath.substring(0, oifPath.lastIndexOf(File.separator) + 1);
try {
RandomAccessInputStream s = getFile(oifName);
s.close();
} catch (IOException e) {
oifName = oifName.replaceAll(".oif", ".OIF");
}
// parse key/value pairs from the OIF file
code = new String[NUM_DIMENSIONS];
size = new String[NUM_DIMENSIONS];
pixelSize = new Double[NUM_DIMENSIONS];
previewNames = new ArrayList<String>();
SortedMap<Integer, String> previewFileNames = new TreeMap<Integer, String>();
boolean laserEnabled = true;
IniList f = getIniFile(oifName);
IniTable saveInfo = f.getTable("ProfileSaveInfo");
String[] saveKeys = saveInfo.keySet().toArray(new String[saveInfo.size()]);
for (String key : saveKeys) {
String value = saveInfo.get(key).toString();
value = sanitizeValue(value).trim();
if (key.startsWith("IniFileName") && key.indexOf("Thumb") == -1 && !isPreviewName(value)) {
filenames.put(new Integer(key.substring(11)), value);
} else if (key.startsWith("RoiFileName") && key.indexOf("Thumb") == -1 && !isPreviewName(value)) {
try {
roiFilenames.put(new Integer(key.substring(11)), value);
} catch (NumberFormatException e) {
}
} else if (key.equals("PtyFileNameS"))
ptyStart = value;
else if (key.equals("PtyFileNameE"))
ptyEnd = value;
else if (key.equals("PtyFileNameT2"))
ptyPattern = value;
else if (key.indexOf("Thumb") != -1) {
if (thumbId == null)
thumbId = value.trim();
} else if (key.startsWith("LutFileName")) {
lutNames.add(path + value);
} else if (isPreviewName(value)) {
try {
RandomAccessInputStream s = getFile(path + value.trim());
if (s != null) {
s.close();
Integer previewIndex = getPreviewNameIndex(key);
if (previewIndex != null) {
previewFileNames.put(previewIndex, path + value.trim());
} else {
previewNames.add(path + value.trim());
}
}
} catch (FormatException e) {
LOGGER.debug("Preview file not found", e);
} catch (IOException e) {
LOGGER.debug("Preview file not found", e);
}
}
}
// Store sorted list of preview names
previewNames.addAll(previewFileNames.values());
if (filenames.isEmpty())
addPtyFiles();
for (int i = 0; i < NUM_DIMENSIONS; i++) {
IniTable commonParams = f.getTable("Axis " + i + " Parameters Common");
code[i] = commonParams.get("AxisCode");
size[i] = commonParams.get("MaxSize");
double end = Double.parseDouble(commonParams.get("EndPosition"));
double start = Double.parseDouble(commonParams.get("StartPosition"));
pixelSize[i] = end - start;
}
IniTable referenceParams = f.getTable("Reference Image Parameter");
imageDepth = Integer.parseInt(referenceParams.get("ImageDepth"));
pixelSizeX = referenceParams.get("WidthConvertValue");
pixelSizeY = referenceParams.get("HeightConvertValue");
String ripValidBitCounts = referenceParams.get("ValidBitCounts");
if (ripValidBitCounts != null) {
validBits = Integer.parseInt(ripValidBitCounts);
}
int index = 0;
IniTable laser = f.getTable("Laser " + index + " Parameters");
while (laser != null) {
laserEnabled = laser.get("Laser Enable").equals("1");
if (laserEnabled) {
wavelengths.add(new Double(laser.get("LaserWavelength")));
}
creationDate = laser.get("ImageCaputreDate");
if (creationDate == null) {
creationDate = laser.get("ImageCaptureDate");
}
index++;
laser = f.getTable("Laser " + index + " Parameters");
}
if (getMetadataOptions().getMetadataLevel() != MetadataLevel.MINIMUM) {
index = 1;
IniTable guiChannel = f.getTable("GUI Channel " + index + " Parameters");
while (guiChannel != null) {
ChannelData channel = new ChannelData();
channel.gain = DataTools.parseDouble(guiChannel.get("AnalogPMTGain"));
channel.voltage = DataTools.parseDouble(guiChannel.get("AnalogPMTVoltage"));
channel.barrierFilter = channel.getFilter(guiChannel.get("BF Name"));
channel.active = Integer.parseInt(guiChannel.get("CH Activate")) != 0;
channel.name = guiChannel.get("CH Name");
channel.dyeName = guiChannel.get("DyeName");
channel.emissionFilter = channel.getFilter(guiChannel.get("EmissionDM Name"));
channel.emWave = DataTools.parseDouble(guiChannel.get("EmissionWavelength"));
channel.excitationFilter = channel.getFilter(guiChannel.get("ExcitationDM Name"));
channel.exWave = DataTools.parseDouble(guiChannel.get("ExcitationWavelength"));
channels.add(channel);
index++;
guiChannel = f.getTable("GUI Channel " + index + " Parameters");
}
index = 1;
IniTable channel = f.getTable("Channel " + index + " Parameters");
while (channel != null) {
String illumination = channel.get("LightType");
if (illumination != null)
illumination = illumination.toLowerCase();
if (illumination == null) {
// Ignored
} else if (illumination.indexOf("fluorescence") != -1) {
illumination = "Epifluorescence";
} else if (illumination.indexOf("transmitted") != -1) {
illumination = "Transmitted";
} else
illumination = null;
illuminations.add(illumination);
index++;
channel = f.getTable("Channel " + index + " Parameters");
}
HashMap<String, String> iniMap = f.flattenIntoHashMap();
metadata.putAll(iniMap);
}
LOGGER.info("Initializing helper readers");
if (previewNames.size() > 0) {
final List<String> v = new ArrayList<String>();
for (int i = 0; i < previewNames.size(); i++) {
String ss = previewNames.get(i);
ss = replaceExtension(ss, "pty", "tif");
if (ss.endsWith(".tif"))
v.add(ss);
}
previewNames = v;
if (previewNames.size() > 0) {
core.clear();
core.add(new CoreMetadata());
core.add(new CoreMetadata());
IFDList ifds = null;
CoreMetadata ms1 = core.get(1);
for (String previewName : previewNames) {
RandomAccessInputStream preview = getFile(previewName);
TiffParser tp = new TiffParser(preview);
ifds = tp.getIFDs();
preview.close();
tp = null;
ms1.imageCount += ifds.size();
}
ms1.sizeX = (int) ifds.get(0).getImageWidth();
ms1.sizeY = (int) ifds.get(0).getImageLength();
ms1.sizeZ = 1;
ms1.sizeT = 1;
ms1.sizeC = ms1.imageCount;
ms1.rgb = false;
int bits = ifds.get(0).getBitsPerSample()[0];
while ((bits % 8) != 0) bits++;
bits /= 8;
ms1.pixelType = FormatTools.pixelTypeFromBytes(bits, false, false);
ms1.dimensionOrder = "XYCZT";
ms1.indexed = false;
}
}
CoreMetadata ms0 = core.get(0);
ms0.imageCount = filenames.size();
tiffs = new ArrayList<String>(getImageCount());
thumbReader = new BMPReader();
if (thumbId != null) {
thumbId = replaceExtension(thumbId, "pty", "bmp");
thumbId = sanitizeFile(thumbId, path);
}
LOGGER.info("Reading additional metadata");
// open each INI file (.pty extension) and build list of TIFF files
String tiffPath = null;
ms0.dimensionOrder = "XY";
final Map<String, String> values = new HashMap<String, String>();
final List<String> baseKeys = new ArrayList<String>();
for (int i = 0, ii = 0; ii < getImageCount(); i++, ii++) {
String file = filenames.get(i);
while (file == null) file = filenames.get(++i);
file = sanitizeFile(file, path);
if (file.indexOf(File.separator) != -1) {
tiffPath = file.substring(0, file.lastIndexOf(File.separator));
} else
tiffPath = file;
Location ptyFile = new Location(file);
if (!isOIB && !ptyFile.exists()) {
LOGGER.warn("Could not find .pty file ({}); guessing at the " + "corresponding TIFF file.", file);
String tiff = replaceExtension(file, "pty", "tif");
Location tiffFile = new Location(tiff);
if (tiffFile.exists()) {
tiffs.add(ii, tiffFile.getAbsolutePath());
continue;
} else {
if (!tiffFile.getParentFile().exists()) {
String realOIFName = new Location(currentId).getName();
String basePath = tiffFile.getParentFile().getParent();
if (mappedOIF) {
tiffPath = basePath + File.separator + realOIFName + ".files";
ptyFile = new Location(tiffPath, ptyFile.getName());
file = ptyFile.getAbsolutePath();
} else {
Location newFile = new Location(basePath, realOIFName + ".files");
ptyFile = new Location(newFile, ptyFile.getName());
file = ptyFile.getAbsolutePath();
tiffPath = newFile.getAbsolutePath();
}
}
}
} else if (!isOIB) {
file = ptyFile.getAbsolutePath();
}
IniList pty = getIniFile(file);
IniTable fileInfo = pty.getTable("File Info");
file = sanitizeValue(fileInfo.get("DataName"));
if (!isPreviewName(file)) {
while (file.indexOf("GST") != -1) {
file = removeGST(file);
}
if (isOIB) {
file = tiffPath + File.separator + file;
} else
file = new Location(tiffPath, file).getAbsolutePath();
file = replaceExtension(file, "pty", "tif");
tiffs.add(ii, file);
}
PlaneData plane = new PlaneData();
for (int dim = 0; dim < NUM_DIMENSIONS; dim++) {
IniTable axis = pty.getTable("Axis " + dim + " Parameters");
if (axis == null)
break;
boolean addAxis = Integer.parseInt(axis.get("Number")) > 1;
if (dim == 2) {
if (addAxis && getDimensionOrder().indexOf('C') == -1) {
ms0.dimensionOrder += 'C';
}
} else if (dim == 3) {
if (addAxis && getDimensionOrder().indexOf('Z') == -1) {
ms0.dimensionOrder += 'Z';
}
final Double number = Double.valueOf(axis.get("AbsPositionValue"));
plane.positionZ = new Length(number, UNITS.REFERENCEFRAME);
} else if (dim == 4) {
if (addAxis && getDimensionOrder().indexOf('T') == -1) {
ms0.dimensionOrder += 'T';
}
// divide by 1000, as the position is in milliseconds
// and DeltaT is in seconds
plane.deltaT = Double.parseDouble(axis.get("AbsPositionValue")) / 1000;
} else if (dim == 7) {
try {
String xPos = axis.get("AbsPositionValueX");
if (xPos != null) {
final Double number = Double.valueOf(xPos);
plane.positionX = new Length(number, UNITS.REFERENCEFRAME);
}
} catch (NumberFormatException e) {
}
try {
String yPos = axis.get("AbsPositionValueY");
if (yPos != null) {
final Double number = Double.valueOf(yPos);
plane.positionY = new Length(number, UNITS.REFERENCEFRAME);
}
} catch (NumberFormatException e) {
}
}
}
ms0.bitsPerPixel = validBits;
planes.add(plane);
IniTable acquisition = pty.getTable("Acquisition Parameters Common");
if (acquisition != null) {
magnification = acquisition.get("Magnification");
lensNA = acquisition.get("ObjectiveLens NAValue");
objectiveName = acquisition.get("ObjectiveLens Name");
workingDistance = acquisition.get("ObjectiveLens WDValue");
pinholeSize = acquisition.get("PinholeDiameter");
String validBitCounts = acquisition.get("ValidBitCounts");
if (validBitCounts != null) {
ms0.bitsPerPixel = Integer.parseInt(validBitCounts);
}
}
if (getMetadataOptions().getMetadataLevel() != MetadataLevel.MINIMUM) {
for (IniTable table : pty) {
String[] keys = table.keySet().toArray(new String[table.size()]);
for (String key : keys) {
values.put("Image " + ii + " : " + key, table.get(key));
if (!baseKeys.contains(key))
baseKeys.add(key);
}
}
}
}
for (String key : baseKeys) {
if (key.equals("DataName") || key.indexOf("FileName") >= 0)
break;
boolean equal = true;
String first = values.get("Image 0 : " + key);
for (int i = 1; i < getImageCount(); i++) {
if (!first.equals(values.get("Image " + i + " : " + key))) {
equal = false;
break;
}
}
if (equal) {
addGlobalMeta(key, first);
} else {
for (int i = 0; i < getImageCount(); i++) {
String k = "Image " + i + " : " + key;
addGlobalMeta(k, values.get(k));
}
}
}
if (tiffs.size() != getImageCount()) {
ms0.imageCount = tiffs.size();
}
usedFiles = new ArrayList<String>();
if (tiffPath != null) {
usedFiles.add(isOIB ? id : oifName);
if (!isOIB) {
Location dir = new Location(tiffPath);
if (!mappedOIF && !dir.exists()) {
throw new FormatException("Required directory " + tiffPath + " was not found.");
}
String[] list = mappedOIF ? Location.getIdMap().keySet().toArray(new String[0]) : dir.list(true);
for (int i = 0; i < list.length; i++) {
if (mappedOIF)
usedFiles.add(list[i]);
else {
String p = new Location(tiffPath, list[i]).getAbsolutePath();
String check = p.toLowerCase();
if (!check.endsWith(".tif") && !check.endsWith(".pty") && !check.endsWith(".roi") && !check.endsWith(".lut") && !check.endsWith(".bmp")) {
continue;
}
usedFiles.add(p);
}
}
}
}
LOGGER.info("Populating metadata");
// calculate axis sizes
int realChannels = 0;
for (int i = 0; i < NUM_DIMENSIONS; i++) {
int ss = Integer.parseInt(size[i]);
if (pixelSize[i] == null)
pixelSize[i] = 1.0;
if (code[i].equals("X"))
ms0.sizeX = ss;
else if (code[i].equals("Y") && ss > 1)
ms0.sizeY = ss;
else if (code[i].equals("Z")) {
if (getSizeY() == 0) {
ms0.sizeY = ss;
} else {
ms0.sizeZ = ss;
// Z size stored in nm
pixelSizeZ = Math.abs((pixelSize[i].doubleValue() / (getSizeZ() - 1)) / 1000);
}
} else if (code[i].equals("T")) {
if (getSizeY() == 0) {
ms0.sizeY = ss;
} else {
ms0.sizeT = ss;
pixelSizeT = Math.abs((pixelSize[i].doubleValue() / (getSizeT() - 1)) / 1000);
}
} else if (ss > 0) {
if (getSizeC() == 0)
ms0.sizeC = ss;
else
ms0.sizeC *= ss;
if (code[i].equals("C"))
realChannels = ss;
}
}
if (getSizeZ() == 0)
ms0.sizeZ = 1;
if (getSizeC() == 0)
ms0.sizeC = 1;
if (getSizeT() == 0)
ms0.sizeT = 1;
if (getImageCount() == getSizeC() && getSizeY() == 1) {
ms0.imageCount *= getSizeZ() * getSizeT();
} else if (getImageCount() == getSizeC()) {
ms0.sizeZ = 1;
ms0.sizeT = 1;
}
if (getSizeZ() * getSizeT() * getSizeC() != getImageCount()) {
int diff = (getSizeZ() * getSizeC() * getSizeT()) - getImageCount();
if (diff == previewNames.size() || diff < 0) {
diff /= getSizeC();
if (getSizeT() > 1 && getSizeZ() == 1)
ms0.sizeT -= diff;
else if (getSizeZ() > 1 && getSizeT() == 1)
ms0.sizeZ -= diff;
} else
ms0.imageCount += diff;
}
if (getSizeC() > 1 && getSizeZ() == 1 && getSizeT() == 1) {
if (getDimensionOrder().indexOf('C') == -1)
ms0.dimensionOrder += 'C';
}
if (getDimensionOrder().indexOf('Z') == -1)
ms0.dimensionOrder += 'Z';
if (getDimensionOrder().indexOf('C') == -1)
ms0.dimensionOrder += 'C';
if (getDimensionOrder().indexOf('T') == -1)
ms0.dimensionOrder += 'T';
ms0.pixelType = FormatTools.pixelTypeFromBytes(imageDepth, false, false);
try {
RandomAccessInputStream thumb = getFile(thumbId);
byte[] b = new byte[(int) thumb.length()];
thumb.read(b);
thumb.close();
Location.mapFile("thumbnail.bmp", new ByteArrayHandle(b));
thumbReader.setId("thumbnail.bmp");
for (int i = 0; i < getSeriesCount(); i++) {
core.get(i).thumbSizeX = thumbReader.getSizeX();
core.get(i).thumbSizeY = thumbReader.getSizeY();
}
thumbReader.close();
Location.mapFile("thumbnail.bmp", null);
} catch (IOException e) {
LOGGER.debug("Could not read thumbnail", e);
} catch (FormatException e) {
LOGGER.debug("Could not read thumbnail", e);
}
// initialize lookup table
lut = new short[getSizeC()][3][65536];
byte[] buffer = new byte[65536 * 4];
int count = (int) Math.min(getSizeC(), lutNames.size());
for (int c = 0; c < count; c++) {
Exception exc = null;
try {
RandomAccessInputStream stream = getFile(lutNames.get(c));
stream.seek(stream.length() - 65536 * 4);
stream.read(buffer);
stream.close();
for (int q = 0; q < buffer.length; q += 4) {
lut[c][0][q / 4] = (short) ((buffer[q + 2] & 0xff) * 257);
lut[c][1][q / 4] = (short) ((buffer[q + 1] & 0xff) * 257);
lut[c][2][q / 4] = (short) ((buffer[q] & 0xff) * 257);
}
} catch (IOException e) {
exc = e;
} catch (FormatException e) {
exc = e;
}
if (exc != null) {
LOGGER.debug("Could not read LUT", exc);
lut = null;
break;
}
}
for (int i = 0; i < getSeriesCount(); i++) {
CoreMetadata ms = core.get(i);
ms.rgb = false;
ms.littleEndian = true;
ms.interleaved = false;
ms.metadataComplete = true;
ms.indexed = lut != null;
ms.falseColor = true;
int nFiles = i == 0 ? tiffs.size() : previewNames.size();
for (int file = 0; file < nFiles; file++) {
RandomAccessInputStream plane = getFile(i == 0 ? tiffs.get(file) : previewNames.get(file));
if (plane == null) {
ifds.add(null);
continue;
}
try {
TiffParser tp = new TiffParser(plane);
IFDList ifd = tp.getIFDs();
ifds.add(ifd);
} finally {
plane.close();
}
}
}
// populate MetadataStore
MetadataStore store = makeFilterMetadata();
MetadataTools.populatePixels(store, this, true);
if (creationDate != null) {
creationDate = creationDate.replaceAll("'", "");
creationDate = DateTools.formatDate(creationDate, DATE_FORMAT);
}
for (int i = 0; i < getSeriesCount(); i++) {
// populate Image data
store.setImageName("Series " + (i + 1), i);
if (creationDate != null)
store.setImageAcquisitionDate(new Timestamp(creationDate), i);
}
if (getMetadataOptions().getMetadataLevel() != MetadataLevel.MINIMUM) {
populateMetadataStore(store, path);
}
}
use of loci.common.Location in project bioformats by openmicroscopy.
the class FV1000Reader method findOIFFile.
// -- Helper methods --
private String findOIFFile(String baseFile) throws FormatException {
Location current = new Location(baseFile).getAbsoluteFile();
String parent = current.getParent();
Location tmp = new Location(parent).getParentFile();
parent = tmp.getAbsolutePath();
baseFile = current.getName();
if (baseFile == null || baseFile.indexOf('_') == -1)
return null;
baseFile = baseFile.substring(0, baseFile.lastIndexOf("_"));
if (checkSuffix(current.getName(), new String[] { "roi", "lut" })) {
if (!new Location(tmp, baseFile + ".oif").exists() && !new Location(tmp, baseFile + ".OIF").exists() && baseFile.indexOf('_') >= 0) {
// some metadata files have an extra underscore
baseFile = baseFile.substring(0, baseFile.lastIndexOf("_"));
}
}
baseFile += ".oif";
tmp = new Location(tmp, baseFile);
String oifFile = tmp.getAbsolutePath();
if (!tmp.exists()) {
oifFile = oifFile.substring(0, oifFile.lastIndexOf(".")) + ".OIF";
tmp = new Location(oifFile);
if (!tmp.exists()) {
baseFile = current.getParent();
baseFile = baseFile.substring(0, baseFile.lastIndexOf("."));
baseFile = baseFile.substring(0, baseFile.lastIndexOf("."));
tmp = new Location(baseFile + ".oif");
oifFile = tmp.getAbsolutePath();
if (!tmp.exists()) {
tmp = new Location(tmp.getParent(), tmp.getName().toUpperCase());
oifFile = tmp.getAbsolutePath();
if (!tmp.exists()) {
// check in parent directory
if (parent.endsWith(File.separator)) {
parent = parent.substring(0, parent.length() - 1);
}
String dir = parent.substring(parent.lastIndexOf(File.separator));
dir = dir.substring(0, dir.lastIndexOf("."));
tmp = new Location(parent);
oifFile = new Location(tmp, dir).getAbsolutePath();
if (!new Location(oifFile).exists()) {
throw new FormatException("OIF file not found");
}
}
}
}
}
return oifFile;
}
use of loci.common.Location in project bioformats by openmicroscopy.
the class BioRadReader method fileGroupOption.
/* @see loci.formats.IFormatReader#fileGroupOption(String) */
@Override
public int fileGroupOption(String id) throws FormatException, IOException {
Location thisFile = new Location(id).getAbsoluteFile();
Location parent = thisFile.getParentFile();
String[] list = parent.list(true);
for (String f : list) {
if (checkSuffix(f, "raw") || checkSuffix(f, "xml")) {
return FormatTools.MUST_GROUP;
}
}
return FormatTools.CAN_GROUP;
}
use of loci.common.Location in project bioformats by openmicroscopy.
the class AnalyzeReader method initFile.
// -- Internal FormatReader API methods --
/* @see loci.formats.FormatReader#initFile(String) */
@Override
protected void initFile(String id) throws FormatException, IOException {
// the dataset has two files - we want the one ending in '.hdr'
if (id.endsWith(".img")) {
LOGGER.info("Looking for header file");
String header = id.substring(0, id.lastIndexOf(".")) + ".hdr";
if (new Location(header).exists()) {
setId(header);
return;
} else
throw new FormatException("Header file not found.");
}
super.initFile(id);
in = new RandomAccessInputStream(id);
pixelsFilename = id.substring(0, id.lastIndexOf(".")) + ".img";
pixelFile = new RandomAccessInputStream(pixelsFilename);
LOGGER.info("Reading header");
int fileSize = in.readInt();
boolean little = fileSize != in.length();
in.order(little);
pixelFile.order(little);
in.skipBytes(10);
String imageName = in.readString(18);
in.skipBytes(8);
int ndims = in.readShort();
int x = in.readShort();
int y = in.readShort();
int z = in.readShort();
int t = in.readShort();
in.skipBytes(20);
int dataType = in.readShort();
int nBitsPerPixel = in.readShort();
String description = null;
double voxelWidth = 0d, voxelHeight = 0d, sliceThickness = 0d, deltaT = 0d;
if (getMetadataOptions().getMetadataLevel() != MetadataLevel.MINIMUM) {
in.skipBytes(6);
voxelWidth = in.readFloat();
voxelHeight = in.readFloat();
sliceThickness = in.readFloat();
deltaT = in.readFloat();
in.skipBytes(12);
pixelOffset = (int) in.readFloat();
in.skipBytes(12);
float calibratedMax = in.readFloat();
float calibratedMin = in.readFloat();
float compressed = in.readFloat();
float verified = in.readFloat();
float pixelMax = in.readFloat();
float pixelMin = in.readFloat();
description = in.readString(80);
String auxFile = in.readString(24);
char orient = (char) in.readByte();
String originator = in.readString(10);
String generated = in.readString(10);
String scannum = in.readString(10);
String patientID = in.readString(10);
String expDate = in.readString(10);
String expTime = in.readString(10);
in.skipBytes(3);
int views = in.readInt();
int volsAdded = in.readInt();
int startField = in.readInt();
int fieldSkip = in.readInt();
int omax = in.readInt();
int omin = in.readInt();
int smax = in.readInt();
int smin = in.readInt();
addGlobalMeta("Database name", imageName);
addGlobalMeta("Number of dimensions", ndims);
addGlobalMeta("Data type", dataType);
addGlobalMeta("Number of bits per pixel", nBitsPerPixel);
addGlobalMeta("Voxel width", voxelWidth);
addGlobalMeta("Voxel height", voxelHeight);
addGlobalMeta("Slice thickness", sliceThickness);
addGlobalMeta("Exposure time", deltaT);
addGlobalMeta("Pixel offset", pixelOffset);
addGlobalMeta("Calibrated maximum", calibratedMax);
addGlobalMeta("Calibrated minimum", calibratedMin);
addGlobalMeta("Compressed", compressed);
addGlobalMeta("Verified", verified);
addGlobalMeta("Pixel maximum", pixelMax);
addGlobalMeta("Pixel minimum", pixelMin);
addGlobalMeta("Description", description);
addGlobalMeta("Auxiliary file", auxFile);
addGlobalMeta("Orientation", orient);
addGlobalMeta("Originator", originator);
addGlobalMeta("Generated", generated);
addGlobalMeta("Scan Number", scannum);
addGlobalMeta("Patient ID", patientID);
addGlobalMeta("Acquisition Date", expDate);
addGlobalMeta("Acquisition Time", expTime);
} else {
in.skipBytes(34);
pixelOffset = (int) in.readFloat();
}
LOGGER.info("Populating core metadata");
CoreMetadata m = core.get(0);
m.littleEndian = little;
m.sizeX = x;
m.sizeY = y;
m.sizeZ = z;
m.sizeT = t;
m.sizeC = 1;
if (getSizeZ() == 0)
m.sizeZ = 1;
if (getSizeT() == 0)
m.sizeT = 1;
m.imageCount = getSizeZ() * getSizeT();
m.rgb = false;
m.interleaved = false;
m.indexed = false;
m.dimensionOrder = "XYZTC";
switch(dataType) {
case 1:
case 2:
m.pixelType = FormatTools.UINT8;
break;
case 4:
m.pixelType = FormatTools.INT16;
break;
case 8:
m.pixelType = FormatTools.INT32;
break;
case 16:
m.pixelType = FormatTools.FLOAT;
break;
case 64:
m.pixelType = FormatTools.DOUBLE;
break;
case 128:
m.pixelType = FormatTools.UINT8;
m.sizeC = 3;
m.rgb = true;
m.interleaved = true;
m.dimensionOrder = "XYCZT";
default:
throw new FormatException("Unsupported data type: " + dataType);
}
LOGGER.info("Populating MetadataStore");
MetadataStore store = makeFilterMetadata();
MetadataTools.populatePixels(store, this);
store.setImageName(imageName, 0);
if (getMetadataOptions().getMetadataLevel() != MetadataLevel.MINIMUM) {
store.setImageDescription(description, 0);
Length sizeX = FormatTools.getPhysicalSizeX(voxelWidth, UNITS.MILLIMETER);
Length sizeY = FormatTools.getPhysicalSizeY(voxelHeight, UNITS.MILLIMETER);
Length sizeZ = FormatTools.getPhysicalSizeZ(sliceThickness, UNITS.MILLIMETER);
if (sizeX != null) {
store.setPixelsPhysicalSizeX(sizeX, 0);
}
if (sizeY != null) {
store.setPixelsPhysicalSizeY(sizeY, 0);
}
if (sizeZ != null) {
store.setPixelsPhysicalSizeZ(sizeZ, 0);
}
store.setPixelsTimeIncrement(new Time(deltaT, UNITS.MILLISECOND), 0);
}
}
Aggregations