#include "gdcmReader.h"
#include "gdcmMediaStorage.h"
typedef std::set<gdcm::DataElement> DataElementSet;
typedef DataElementSet::const_iterator ConstIterator;
int main(int argc, char *argv [])
{
if( argc < 2 ) return 1;
const char *filename = argv[1];
gdcm::Reader reader;
reader.SetFileName( filename);
if( !reader.Read() )
{
std::cerr << "Could not read: " << filename << std::endl;
return 1;
}
std::stringstream strm;
gdcm::File &file = reader.GetFile();
gdcm::DataSet &ds = file.GetDataSet();
gdcm::FileMetaInformation &fmi = file.GetHeader();
gdcm::MediaStorage ms;
ms.SetFromFile(file);
if( ms != gdcm::MediaStorage::MediaStorageDirectoryStorage )
{
std::cout << "This file is not a DICOMDIR" << std::endl;
return 1;
}
if (fmi.FindDataElement( gdcm::Tag (0x0002, 0x0002)))
{ strm.str("");
fmi.GetDataElement( gdcm::Tag (0x0002, 0x0002) ).GetValue().Print(strm);
}
else
{
std::cerr << " Media Storage Sop Class UID not present" << std::cout;
}
if ("1.2.840.10008.1.3.10"!=strm.str())
{
std::cout << "This file is not a DICOMDIR" << std::endl;
return 1;
}
ConstIterator it = ds.GetDES().begin();
for( ; it != ds.GetDES().end(); ++it)
{
if (it->GetTag()==gdcm::Tag (0x0004, 0x1220))
{
const gdcm::DataElement &de = (*it);
gdcm::SmartPointer<gdcm::SequenceOfItems> sqi =de.GetValueAsSQ();
int itemused = 1;
while (itemused<=sqi->GetNumberOfItems())
{
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0004, 0x1430)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0004, 0x1430)).GetValue().Print(strm);
while((strm.str()=="PATIENT")||((strm.str()=="PATIENT ")))
{
std::cout << strm.str() << std::endl;
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0010, 0x0010)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0010, 0x0010)).GetValue().Print(strm);
std::cout << "PATIENT NAME : " << strm.str() << std::endl;
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0010, 0x0020)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0010, 0x0020)).GetValue().Print(strm);
std::cout << "PATIENT ID : " << strm.str() << std::endl;
std::cout << "=========================== " << std::endl;
itemused++;
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0004, 0x1430)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0004, 0x1430)).GetValue().Print(strm);
while((strm.str()=="STUDY")||((strm.str()=="STUDY ")))
{
std::cout << " " << strm.str() << std::endl;
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0020, 0x000d)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0020, 0x000d)).GetValue().Print(strm);
std::cout << " STUDY UID : " << strm.str() << std::endl;
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0008, 0x0020)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0008, 0x0020)).GetValue().Print(strm);
std::cout << " STUDY DATE : " << strm.str() << std::endl;
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0008, 0x1030)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0008, 0x1030)).GetValue().Print(strm);
std::cout << " STUDY DESCRIPTION : " << strm.str() << std::endl;
std::cout << " " << "=========================== " << std::endl;
itemused++;
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0004, 0x1430)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0004, 0x1430)).GetValue().Print(strm);
while((strm.str()=="SERIES")||((strm.str()=="SERIES ")))
{
std::cout << " " << strm.str() << std::endl;
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0020, 0x000e)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0020, 0x000e)).GetValue().Print(strm);
std::cout << " SERIE UID" << strm.str() << std::endl;
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0008, 0x0060)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0008, 0x0060)).GetValue().Print(strm);
std::cout << " SERIE MODALITY" << strm.str() << std::endl;
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0008, 0x103e)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0008, 0x103e)).GetValue().Print(strm);
std::cout << " SERIE DESCRIPTION" << strm.str() << std::endl;
std::cout << " " << "=========================== " << std::endl;
itemused++;
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0004, 0x1430)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0004, 0x1430)).GetValue().Print(strm);
while ((strm.str()=="IMAGE")||((strm.str()=="IMAGE ")))
{
std::cout << " " << strm.str() << std::endl;
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0004, 0x1511)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0004, 0x1511)).GetValue().Print(strm);
std::cout << " IMAGE UID : " << strm.str() << std::endl;
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0004, 0x1500)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0004, 0x1500)).GetValue().Print(strm);
std::cout << " IMAGE PATH : " << strm.str() << std::endl;
if(itemused <sqi->GetNumberOfItems())
{itemused++;
}else{break;}
strm.str("");
if (sqi->GetItem(itemused).FindDataElement(gdcm::Tag (0x0004, 0x1430)))
sqi->GetItem(itemused).GetDataElement(gdcm::Tag (0x0004, 0x1430)).GetValue().Print(strm);
}
}
}
}
itemused++;
}
}
}
return 0;
}