bool DCMTKImageIO::CanReadFile(const char* filename)
{
DcmFileFormat dicomFile;
OFCondition condition = dicomFile.loadFile( filename );
if ( !condition.good() ) {
return false;
}
E_TransferSyntax xfer = dicomFile.getDataset()->getOriginalXfer();
if( xfer == EXS_JPEG2000LosslessOnly ||
xfer == EXS_JPEG2000 ||
xfer == EXS_JPEG2000MulticomponentLosslessOnly ||
xfer == EXS_JPEG2000Multicomponent ) {
return false;
}
// search for mandatory tags
DcmStack stack;
DcmTagKey searchKey;
unsigned int group = 0x0028; // samples per pixel
unsigned int elem = 0x0002; // samples per pixel
searchKey.set(group, elem);
if (dicomFile.search(searchKey, stack, ESM_fromHere, OFTrue) != EC_Normal)
return false;
group = 0x0028; // pixel type
elem = 0x0100; // pixel type
searchKey.set(group, elem);
if (dicomFile.search(searchKey, stack, ESM_fromHere, OFTrue) != EC_Normal)
return false;
return true;
}
DcmElement* DcmtkHelper::searchElement(DcmObject* _object, unsigned int _group, unsigned int _element)
{
try
{
if (_object == NULL)
{
throw DicomLibException("No data!", __FILE__, __LINE__);
}
DcmTagKey searchKey;
searchKey.set(_group, _element);
DcmStack stack;
_object -> search(searchKey, stack, ESM_fromHere, OFTrue);
DcmElement* element = NULL;
element = reinterpret_cast<DcmElement*>(stack.elem(0));
return element;
}
catch (...)
{
throw;
}
}
int AddMetaData::process(GadgetContainerMessage<DcmFileFormat> * m1)
{
DcmFileFormat * dcm = m1->getObjectPtr();
GadgetContainerMessage<std::string> * f;
GadgetContainerMessage<ISMRMRD::MetaContainer>* meta;
if(dcm)
{
f=AsContainerMessage<std::string>(m1->cont());
}
else
{ GERROR("No filename set for DICOM file\n");
return GADGET_FAIL;
}
if(f)
{
meta= AsContainerMessage<ISMRMRD::MetaContainer>(f->cont());
}
else
{
GERROR("No meta data found for DICOM\n");
return GADGET_FAIL;
}
unsigned int BUFSIZE = 1024;
char *buf = new char[BUFSIZE];
const char* checkbuf;
OFCondition status;
DcmTagKey key;
DcmDataset *dataset = dcm->getDataset();
float rescaleIntercept;//= meta->getObjectPtr()->as_double(GADGETRON_IMAGE_SCALE_OFFSET);
float rescaleSlope;//= meta->getObjectPtr()->as_double(GADGETRON_IMAGE_SCALE_RATIO);
static bool studyUIDmade=false;
std::time_t rawtime;
std::time(&rawtime);
std::tm *timeinfo = std::localtime(&rawtime);
if(meta->getObjectPtr()->exists("GADGETRON_IMAGE_SCALE_OFFSET") && meta->getObjectPtr()->exists("GADGETRON_IMAGE_SCALE_RATIO"))
{
rescaleIntercept=meta->getObjectPtr()->as_double(GADGETRON_IMAGE_SCALE_OFFSET);
rescaleIntercept=meta->getObjectPtr()->as_double(GADGETRON_IMAGE_SCALE_OFFSET);
// Window Center
key.set(0x0028, 0x1050);
dataset->remove(key);
// Window Width
key.set(0x0028, 0x1051);
dataset->remove(key);
rescaleIntercept = -1.0*rescaleIntercept*rescaleSlope;
rescaleSlope= 1.0/rescaleSlope;
key.set(0x0028,0x1052);
ACE_OS::snprintf(buf, BUFSIZE, "%f", rescaleIntercept);//
WRITE_DCM_STRING(key, buf);
key.set(0x0028,0x1053);
ACE_OS::snprintf(buf, BUFSIZE, "%f", rescaleSlope);//meta->getObjectPtr()->as_double("Intercept"));
WRITE_DCM_STRING(key, buf);
key.set(0x0028, 0x0030);
ACE_OS::snprintf(buf, BUFSIZE, "%.6f\\%.6f", pixel_spacing_Y, pixel_spacing_X);
WRITE_DCM_STRING(key, buf);
}
key.set(0x0008,0x1030); //Study Description
dataset->findAndGetString(key, checkbuf, false);
if(checkbuf==NULL || !strcmp(checkbuf, "XXXXXXXX"))
{
ACE_OS::snprintf(buf, BUFSIZE, "%s", "Gadgetron^IEV");
WRITE_DCM_STRING(key, buf);
}
key.set(0x0008,0x103E); //Series Description
//if(!dataset->tagExistsWithValue(key))
//{
std::string type;
if(meta)
type=meta->getObjectPtr()->as_str(GADGETRON_DATA_ROLE);
else
type="MRI Images";
ACE_OS::snprintf(buf, BUFSIZE, "%s", type.c_str());
WRITE_DCM_STRING(key, buf);
//}
key.set(0x0020,0x0010);//Study ID
dataset->findAndGetString(key, checkbuf, false);
if(checkbuf==NULL || !strcmp(checkbuf, "XXXXXXXX"))
{
WRITE_DCM_STRING(key, "1");
// be sure to use the same one for all series you generate
}
//Study UID should be created in IEVChannelSumGadget.
key.set(0x0020,0x000D);//Study UID
dataset->findAndGetString(key, checkbuf, false);
if(checkbuf==NULL || !strcmp(checkbuf, "XXXXXXXX"))
{
WRITE_DCM_STRING(key, meta->getObjectPtr()->as_str("StudyInstanceUID"));
// be sure to use the same one for all series you generate
}
std::strftime(buf, 100, "%Y%m%d", timeinfo);
key.set(0x0008,0x0020);//Study Date
dataset->findAndGetString(key, checkbuf, false);
if(checkbuf==NULL || !strcmp(checkbuf, "19000101"))
{
WRITE_DCM_STRING(key, buf);
}
key.set(0x0008,0x0030);//Study Time
dataset->findAndGetString(key, checkbuf, false);
if(checkbuf==NULL || !strcmp(checkbuf, "121212"))
{
WRITE_DCM_STRING(key, buf);
}
key.set(0x0008,0x0021);//Series Date
if(!dataset->tagExistsWithValue(key))
{
WRITE_DCM_STRING(key, buf);
}
key.set(0x0008,0x0012);//Instance Creation Date
if(!dataset->tagExistsWithValue(key))
{
WRITE_DCM_STRING(key, buf);
}
std::strftime(buf, 100, "%H%M%S", timeinfo);
key.set(0x0008,0x0031);//Series Time
if(!dataset->tagExistsWithValue(key))
{
WRITE_DCM_STRING(key, buf);
}
key.set(0x0008,0x0013);//Instance Creation Time
if(!dataset->tagExistsWithValue(key))
{
WRITE_DCM_STRING(key, buf);
}
key.set(0x0018,0x0081);//Echo Time
WRITE_DCM_STRING(key, meta->getObjectPtr()->as_str("TE"));
delete[] buf;
//add try catch
if(-1==this->next()->putq(m1))
{
m1->release();
GERROR("Unable to pass on message\n");
return GADGET_FAIL;
}
return GADGET_OK;
}
