bool
RawSlabsPlugin::setFile(QStringList files)
{
m_fileName = files;
{
m_slices.clear();
QFileInfo fi(files[0]);
QString hdrflnm;
hdrflnm = QFileDialog::getOpenFileName(0,
"Load text header file",
fi.absolutePath(),
"Files (*.*)");
int nfX0, nfY0, nfZ0;
if (hdrflnm.isEmpty())
{
uchar nvt0;
QFile fd(m_fileName[0]);
fd.open(QFile::ReadOnly);
fd.read((char*)&nvt0, sizeof(unsigned char));
fd.read((char*)&nfX0, sizeof(int));
fd.read((char*)&nfY0, sizeof(int));
fd.read((char*)&nfZ0, sizeof(int));
fd.close();
m_voxelType = _UChar;
if (nvt0 == 0) m_voxelType = _UChar;
if (nvt0 == 1) m_voxelType = _Char;
if (nvt0 == 2) m_voxelType = _UShort;
if (nvt0 == 3) m_voxelType = _Short;
if (nvt0 == 4) m_voxelType = _Int;
if (nvt0 == 8) m_voxelType = _Float;
m_headerBytes = m_skipBytes = 13;
int nslices = nfX0;
m_slices.append(nslices);
for (int nf=1; nf<m_fileName.count(); nf++)
{
uchar nvt;
int nfX, nfY, nfZ;
QFile fd(m_fileName[nf]);
fd.open(QFile::ReadOnly);
fd.read((char*)&nvt, sizeof(unsigned char));
fd.read((char*)&nfX, sizeof(int));
fd.read((char*)&nfY, sizeof(int));
fd.read((char*)&nfZ, sizeof(int));
fd.close();
if (nvt != nvt0 || nfY != nfY0 || nfZ != nfZ0)
{
QMessageBox::information(0, "", "Raw File format does not match");
return false;
}
nslices += nfX;
m_slices.append(nslices);
}
}
else
{
// load header information from a text file
int nvt0, nvt;
int nfX, nfY, nfZ;
QFile fd(hdrflnm);
fd.open(QFile::ReadOnly | QFile::Text);
QTextStream in(&fd);
QString line = (in.readLine()).simplified();
QStringList words = line.split(" ", QString::SkipEmptyParts);
if (words.count() >= 2)
{
int nvt0 = words[0].toInt();
m_voxelType = _UChar;
if (nvt0 == 0) m_voxelType = _UChar;
if (nvt0 == 1) m_voxelType = _Char;
if (nvt0 == 2) m_voxelType = _UShort;
if (nvt0 == 3) m_voxelType = _Short;
if (nvt0 == 4) m_voxelType = _Int;
if (nvt0 == 8) m_voxelType = _Float;
m_headerBytes = m_skipBytes = words[1].toInt();
}
else
{
QMessageBox::information(0, "",
QString("Expecting voxeltype and headerbytes\nGot this %1").arg(line));
return false;
}
int nslices = 0;
while (!in.atEnd())
{
line = (in.readLine()).simplified();
words = line.split(" ", QString::SkipEmptyParts);
if (words.count() >= 3)
{
nfX = words[0].toInt();
nfY = words[1].toInt();
nfZ = words[2].toInt();
if (nslices == 0)
{
nfY0 = nfY;
nfZ0 = nfZ;
}
else
{
if (nfY0 != nfY || nfZ0 != nfZ)
{
QMessageBox::information(0, "",
QString("Slice size not same :: %1 %2 : %3 %4").\
arg(nfY0).arg(nfZ0).arg(nfY).arg(nfZ));
return false;
}
}
nslices += nfX;
m_slices.append(nslices);
}
else
{
QMessageBox::information(0, "",
QString("Expecting dimensions\nGot this %1").arg(line));
return false;
}
}
}
m_depth = m_slices[m_slices.count()-1];
m_width = nfY0;
m_height = nfZ0;
}
//------------------------------
m_bytesPerVoxel = 1;
if (m_voxelType == _UChar) m_bytesPerVoxel = 1;
else if (m_voxelType == _Char) m_bytesPerVoxel = 1;
else if (m_voxelType == _UShort) m_bytesPerVoxel = 2;
else if (m_voxelType == _Short) m_bytesPerVoxel = 2;
else if (m_voxelType == _Int) m_bytesPerVoxel = 4;
else if (m_voxelType == _Float) m_bytesPerVoxel = 4;
if (m_voxelType == _UChar ||
m_voxelType == _Char ||
m_voxelType == _UShort ||
m_voxelType == _Short)
{
findMinMaxandGenerateHistogram();
}
else
{
findMinMax();
generateHistogram();
}
return true;
}
bool
RawPlugin::setFile(QStringList files)
{
m_fileName = files;
int nX, nY, nZ;
{
// --- load various parameters from the raw file ---
LoadRawDialog loadRawDialog(0,
(char *)m_fileName[0].toAscii().data());
if (!m_4dvol)
{
loadRawDialog.exec();
if (loadRawDialog.result() == QDialog::Rejected)
return false;
}
m_voxelType = loadRawDialog.voxelType();
m_skipBytes = loadRawDialog.skipHeaderBytes();
loadRawDialog.gridSize(nX, nY, nZ);
m_depth = nX;
m_width = nY;
m_height = nZ;
}
//-----------------------------------
QFile fin(m_fileName[0]);
fin.open(QFile::ReadOnly);
//-- recheck the information (for backward compatibility) ----
if (m_skipBytes == 0)
{
int bpv = 1;
if (m_voxelType == _UChar) bpv = 1;
else if (m_voxelType == _Char) bpv = 1;
else if (m_voxelType == _UShort) bpv = 2;
else if (m_voxelType == _Short) bpv = 2;
else if (m_voxelType == _Int) bpv = 4;
else if (m_voxelType == _Float) bpv = 4;
if (fin.size() == 13+nX*nY*nZ*bpv)
m_skipBytes = 13;
else if (fin.size() == 12+nX*nY*nZ*bpv)
m_skipBytes = 12;
else
m_skipBytes = 0;
}
m_headerBytes = m_skipBytes;
fin.close();
//------------------------------
m_bytesPerVoxel = 1;
if (m_voxelType == _UChar) m_bytesPerVoxel = 1;
else if (m_voxelType == _Char) m_bytesPerVoxel = 1;
else if (m_voxelType == _UShort) m_bytesPerVoxel = 2;
else if (m_voxelType == _Short) m_bytesPerVoxel = 2;
else if (m_voxelType == _Int) m_bytesPerVoxel = 4;
else if (m_voxelType == _Float) m_bytesPerVoxel = 4;
if (m_4dvol) // do not perform further calculations.
return true;
if (m_voxelType == _UChar ||
m_voxelType == _Char ||
m_voxelType == _UShort ||
m_voxelType == _Short)
{
findMinMaxandGenerateHistogram();
}
else
{
findMinMax();
generateHistogram();
}
return true;
}
bool
GrdPlugin::setFile(QStringList files)
{
m_fileName = files;
m_imageList.clear();
QFileInfo f(m_fileName[0]);
if (f.isDir())
{
// list all image files in the directory
QStringList imageNameFilter;
imageNameFilter << "*";
QStringList imgfiles= QDir(m_fileName[0]).entryList(imageNameFilter,
QDir::NoSymLinks|
QDir::NoDotAndDotDot|
QDir::Readable|
QDir::Files);
m_imageList.clear();
for(uint i=0; i<imgfiles.size(); i++)
{
QFileInfo fileInfo(m_fileName[0], imgfiles[i]);
QString imgfl = fileInfo.absoluteFilePath();
m_imageList.append(imgfl);
}
}
else
m_imageList = files;
// --- load various parameters from the raw file ---
LoadRawDialog loadRawDialog(0,
(char *)m_imageList[0].toAscii().data());
//(char *)m_fileName[0].toAscii().data());
loadRawDialog.exec();
if (loadRawDialog.result() == QDialog::Rejected)
return false;
m_voxelType = loadRawDialog.voxelType();
m_headerBytes = loadRawDialog.skipHeaderBytes();
int nX, nY, nZ;
loadRawDialog.gridSize(nX, nY, nZ);
m_depth = m_imageList.size();
m_width = nX;
m_height = nY;
m_bytesPerVoxel = 1;
if (m_voxelType == _UChar) m_bytesPerVoxel = 1;
else if (m_voxelType == _Char) m_bytesPerVoxel = 1;
else if (m_voxelType == _UShort) m_bytesPerVoxel = 2;
else if (m_voxelType == _Short) m_bytesPerVoxel = 2;
else if (m_voxelType == _Int) m_bytesPerVoxel = 4;
else if (m_voxelType == _Float) m_bytesPerVoxel = 4;
if (m_voxelType == _UChar ||
m_voxelType == _Char ||
m_voxelType == _UShort ||
m_voxelType == _Short)
{
findMinMaxandGenerateHistogram();
}
else
{
findMinMax();
generateHistogram();
}
return true;
}
bool
NrrdPlugin::setFile(QStringList files)
{
m_fileName = files;
typedef itk::Image<unsigned char, 3> ImageType;
typedef itk::ImageFileReader<ImageType> ReaderType;
ReaderType::Pointer reader = ReaderType::New();
reader->SetFileName(m_fileName[0].toAscii().data());
typedef itk::NrrdImageIO NrrdIOType;
NrrdIOType::Pointer nrrdIO = NrrdIOType::New();
reader->SetImageIO(nrrdIO);
reader->Update();
itk::ImageIOBase::Pointer imageIO = reader->GetImageIO();
m_height = imageIO->GetDimensions(0);
m_width = imageIO->GetDimensions(1);
m_depth = imageIO->GetDimensions(2);
m_voxelSizeX = imageIO->GetSpacing(0);
m_voxelSizeY = imageIO->GetSpacing(1);
m_voxelSizeZ = imageIO->GetSpacing(2);
int et = imageIO->GetComponentType();
if (et == itk::ImageIOBase::UCHAR) m_voxelType = _UChar;
if (et == itk::ImageIOBase::CHAR) m_voxelType = _Char;
if (et == itk::ImageIOBase::USHORT) m_voxelType = _UShort;
if (et == itk::ImageIOBase::SHORT) m_voxelType = _Short;
if (et == itk::ImageIOBase::INT) m_voxelType = _Int;
if (et == itk::ImageIOBase::FLOAT) m_voxelType = _Float;
m_skipBytes = m_headerBytes = 0;
m_bytesPerVoxel = 1;
if (m_voxelType == _UChar) m_bytesPerVoxel = 1;
else if (m_voxelType == _Char) m_bytesPerVoxel = 1;
else if (m_voxelType == _UShort) m_bytesPerVoxel = 2;
else if (m_voxelType == _Short) m_bytesPerVoxel = 2;
else if (m_voxelType == _Int) m_bytesPerVoxel = 4;
else if (m_voxelType == _Float) m_bytesPerVoxel = 4;
if (m_4dvol) // do not perform further calculations.
return true;
if (m_voxelType == _UChar ||
m_voxelType == _Char ||
m_voxelType == _UShort ||
m_voxelType == _Short)
{
findMinMaxandGenerateHistogram();
}
else
{
findMinMax();
generateHistogram();
}
return true;
}
bool
RemapHDF4::setFile(QList<QString> fl)
{
m_fileName = fl;
// list all image files in the directory
QStringList imageNameFilter;
imageNameFilter << "*.hdf";
QStringList files= QDir(m_fileName[0]).entryList(imageNameFilter,
QDir::NoSymLinks|
QDir::NoDotAndDotDot|
QDir::Readable|
QDir::Files);
m_imageList.clear();
for(uint i=0; i<files.size(); i++)
{
QFileInfo fileInfo(m_fileName[0], files[i]);
QString imgfl = fileInfo.absoluteFilePath();
m_imageList.append(imgfl);
}
m_depth = m_imageList.size();
/* Open the file and initiate the SD interface. */
int32 sd_id = SDstart(strdup(m_imageList[0].toAscii().data()),
DFACC_READ);
if (sd_id < 0) {
QMessageBox::information(0,
"Error",
QString("Failed to open %1").arg(m_imageList[0]));
return false;
}
/* Determine the contents of the file. */
int32 dim_sizes[MAX_VAR_DIMS];
int32 rank, num_type, attributes, istat;
char name[64];
int32 n_datasets, n_file_attrs;
istat = SDfileinfo(sd_id, &n_datasets, &n_file_attrs);
/* Access the name of every data set in the file. */
QStringList varNames;
for (int32 index = 0; index < n_datasets; index++)
{
int32 sds_id = SDselect(sd_id, index);
istat = SDgetinfo(sds_id, name, &rank, dim_sizes, \
&num_type, &attributes);
istat = SDendaccess(sds_id);
if (rank == 2)
varNames.append(name);
}
QString var;
if (varNames.size() == 0) {
QMessageBox::information(0,
"Error",
QString("No variables in file with rank of 2"));
return false;
} else if (varNames.size() == 1)
{
var = varNames[0];
}
else
{
var = QInputDialog::getItem(0,
"Select Variable to Extract",
"Variable Names",
varNames,
0,
false);
}
m_Index = 0;
for (int32 index = 0; index < n_datasets; index++)
{
int32 sds_id = SDselect(sd_id, index);
istat = SDgetinfo(sds_id, name, &rank, dim_sizes, \
&num_type, &attributes);
istat = SDendaccess(sds_id);
if (var == QString(name))
{
m_Index = index;
break;
}
}
{
int32 sds_id = SDselect(sd_id, m_Index);
istat = SDgetinfo(sds_id,
name,
&rank,
dim_sizes,
&num_type,
&attributes);
istat = SDendaccess(sds_id);
}
/* Terminate access to the SD interface and close the file. */
istat = SDend(sd_id);
if (num_type == DFNT_CHAR8)
m_voxelType = _Char;
else if (num_type == DFNT_UCHAR8)
m_voxelType = _UChar;
else if (num_type == DFNT_INT8)
m_voxelType = _Char;
else if (num_type == DFNT_UINT8)
m_voxelType = _UChar;
else if (num_type == DFNT_INT16)
m_voxelType = _Short;
else if (num_type == DFNT_UINT16)
m_voxelType = _UShort;
else if (num_type == DFNT_INT32)
m_voxelType = _Int;
else if (num_type == DFNT_FLOAT32)
m_voxelType = _Float;
else
{
QMessageBox::information(0, "Error",
QString("Cannot handle datatype %1").arg(num_type));
return false;
}
m_width = dim_sizes[0];
m_height = dim_sizes[1];
m_bytesPerVoxel = 1;
if (m_voxelType == _UChar) m_bytesPerVoxel = 1;
else if (m_voxelType == _Char) m_bytesPerVoxel = 1;
else if (m_voxelType == _UShort) m_bytesPerVoxel = 2;
else if (m_voxelType == _Short) m_bytesPerVoxel = 2;
else if (m_voxelType == _Int) m_bytesPerVoxel = 4;
else if (m_voxelType == _Float) m_bytesPerVoxel = 4;
m_headerBytes = 0;
if (m_voxelType == _UChar ||
m_voxelType == _Char ||
m_voxelType == _UShort ||
m_voxelType == _Short)
{
findMinMaxandGenerateHistogram();
}
else
{
findMinMax();
generateHistogram();
}
m_rawMap.append(m_rawMin);
m_rawMap.append(m_rawMax);
m_pvlMap.append(0);
m_pvlMap.append(255);
return true;
}
