Module: $RCSfile: gdcmFileHelper.cxx,v $
Language: C++
- Date: $Date: 2005/04/27 09:52:28 $
- Version: $Revision: 1.35 $
+ Date: $Date: 2005/06/24 10:55:59 $
+ Version: $Revision: 1.46 $
Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
l'Image). All rights reserved. See Doc/License.txt or
*/
FileHelper::~FileHelper()
{
- if( PixelReadConverter )
+ if ( PixelReadConverter )
{
delete PixelReadConverter;
}
- if( PixelWriteConverter )
+ if ( PixelWriteConverter )
{
delete PixelWriteConverter;
}
- if( Archive )
+ if ( Archive )
{
delete Archive;
}
- if( SelfHeader )
+ if ( SelfHeader )
{
delete FileInternal;
}
ValEntry *FileHelper::InsertValEntry(std::string const &content,
uint16_t group, uint16_t elem)
{
- return FileInternal->InsertValEntry(content,group,elem);
+ return FileInternal->InsertValEntry(content, group, elem);
}
/**
{
return PixelWriteConverter->GetUserDataSize();
}
-
return PixelReadConverter->GetRGBSize();
}
{
return PixelWriteConverter->GetUserDataSize();
}
-
return PixelReadConverter->GetRawSize();
}
return 0;
}
memcpy( destination,
- (void*)PixelReadConverter->GetRaw(),
+ (void *)PixelReadConverter->GetRaw(),
PixelReadConverter->GetRawSize() );
return PixelReadConverter->GetRawSize();
}
*/
uint8_t* FileHelper::GetLutRGBA()
{
+ if ( PixelReadConverter->GetLutRGBA() ==0 )
+ PixelReadConverter->BuildLUTRGBA();
return PixelReadConverter->GetLutRGBA();
}
+/**
+ * \brief Access to the underlying \ref PixelReadConverter RGBA LUT Item Number
+ */
+int FileHelper::GetLutItemNumber()
+{
+ return PixelReadConverter->GetLutItemNumber();
+}
+
+/**
+ * \brief Access to the underlying \ref PixelReadConverter RGBA LUT Item Size
+ */
+int FileHelper::GetLutItemSize()
+{
+ return PixelReadConverter->GetLutItemSize();
+}
+
/**
* \brief Writes on disk A SINGLE Dicom file
* NO test is performed on processor "Endiannity".
return false;
}
- if( PixelWriteConverter->GetUserData() )
+ if ( PixelWriteConverter->GetUserData() )
{
- fp1.write( (char*)PixelWriteConverter->GetUserData(),
+ fp1.write( (char *)PixelWriteConverter->GetUserData(),
PixelWriteConverter->GetUserDataSize() );
}
else if ( PixelReadConverter->GetRGB() )
{
- fp1.write( (char*)PixelReadConverter->GetRGB(),
+ fp1.write( (char *)PixelReadConverter->GetRGB(),
PixelReadConverter->GetRGBSize());
}
else if ( PixelReadConverter->GetRaw() )
{
- fp1.write( (char*)PixelReadConverter->GetRaw(),
+ fp1.write( (char *)PixelReadConverter->GetRaw(),
PixelReadConverter->GetRawSize());
}
else
/**
* \brief Writes on disk A SINGLE Dicom file,
* using the Implicit Value Representation convention
- * NO test is performed on processor "Endiannity".
+ * NO test is performed on processor "Endianity".
* @param fileName name of the file to be created
* (any already existing file is overwritten)
* @return false if write fails
{
case ImplicitVR:
SetWriteFileTypeToImplicitVR();
- CheckMandatoryElements();
break;
case Unknown: // should never happen; ExplicitVR is the default value
case ExplicitVR:
SetWriteFileTypeToExplicitVR();
- CheckMandatoryElements();
break;
case ACR:
case ACR_LIBIDO:
- // Just to avoid further trouble if user create a file ex-nihilo,
+ // NOTHING is done here just for LibIDO.
+ // Just to avoid further trouble if user creates a file ex-nihilo,
// wants to write it as an ACR-NEMA file,
- // and forget to create any Entry belonging to group 0008
+ // and forgets to create any Entry belonging to group 0008
// (shame on him !)
// We add Recognition Code (RET)
if ( ! FileInternal->GetValEntry(0x0008, 0x0010) )
FileInternal->InsertValEntry("", 0x0008, 0x0010);
SetWriteFileTypeToACR();
- SetWriteFileTypeToImplicitVR();
- CheckMandatoryElements();
+ // SetWriteFileTypeToImplicitVR(); // ACR IS implicit VR !
break;
-
}
+ CheckMandatoryElements();
// --------------------------------------------------------------
// Special Patch to allow gdcm to re-write ACR-LibIDO formated images
/// \todo the best trick would be *change* the recognition code
/// but pb expected if user deals with, e.g. COMPLEX images
- if( WriteType == ACR_LIBIDO )
+ if ( WriteType == ACR_LIBIDO )
{
SetWriteToLibido();
}
}
bool check = CheckWriteIntegrity(); // verifies length
- if(check)
+ if (check)
{
check = FileInternal->Write(fileName,WriteType);
}
*/
bool FileHelper::CheckWriteIntegrity()
{
- if(PixelWriteConverter->GetUserData())
+ if ( PixelWriteConverter->GetUserData() )
{
int numberBitsAllocated = FileInternal->GetBitsAllocated();
if ( numberBitsAllocated == 0 || numberBitsAllocated == 12 )
{
+ gdcmWarningMacro( "numberBitsAllocated changed from "
+ << numberBitsAllocated << " to 16 "
+ << " for consistency purpose" );
numberBitsAllocated = 16;
}
size_t decSize = FileInternal->GetXSize()
- * FileInternal->GetYSize()
- * FileInternal->GetZSize()
- * ( numberBitsAllocated / 8 )
- * FileInternal->GetSamplesPerPixel();
+ * FileInternal->GetYSize()
+ * FileInternal->GetZSize()
+ * FileInternal->GetSamplesPerPixel()
+ * ( numberBitsAllocated / 8 );
size_t rgbSize = decSize;
- if( FileInternal->HasLUT() )
+ if ( FileInternal->HasLUT() )
rgbSize = decSize * 3;
switch(WriteMode)
{
case WMODE_RAW :
- if( decSize!=PixelWriteConverter->GetUserDataSize() )
+ if ( decSize!=PixelWriteConverter->GetUserDataSize() )
{
gdcmWarningMacro( "Data size (Raw) is incorrect. Should be "
<< decSize << " / Found :"
}
break;
case WMODE_RGB :
- if( rgbSize!=PixelWriteConverter->GetUserDataSize() )
+ if ( rgbSize!=PixelWriteConverter->GetUserDataSize() )
{
gdcmWarningMacro( "Data size (RGB) is incorrect. Should be "
<< decSize << " / Found "
}
break;
}
- }
-
+ }
return true;
}
*/
void FileHelper::SetWriteToRaw()
{
- if( FileInternal->GetNumberOfScalarComponents() == 3
- && !FileInternal->HasLUT())
+ if ( FileInternal->GetNumberOfScalarComponents() == 3
+ && !FileInternal->HasLUT() )
{
SetWriteToRGB();
}
else
{
ValEntry *photInt = CopyValEntry(0x0028,0x0004);
- if(FileInternal->HasLUT())
+ if (FileInternal->HasLUT() )
{
photInt->SetValue("PALETTE COLOR ");
}
PixelReadConverter->GetRawSize());
std::string vr = "OB";
- if( FileInternal->GetBitsAllocated()>8 )
+ if ( FileInternal->GetBitsAllocated()>8 )
vr = "OW";
- if( FileInternal->GetBitsAllocated()==24 ) // For RGB ACR files
+ if ( FileInternal->GetBitsAllocated()==24 ) // For RGB ACR files
vr = "OB";
BinEntry *pixel =
CopyBinEntry(GetFile()->GetGrPixel(),GetFile()->GetNumPixel(),vr);
*/
void FileHelper::SetWriteToRGB()
{
- if(FileInternal->GetNumberOfScalarComponents()==3)
+ if ( FileInternal->GetNumberOfScalarComponents()==3 )
{
PixelReadConverter->BuildRGBImage();
ValEntry *photInt = CopyValEntry(0x0028,0x0004);
photInt->SetValue("RGB ");
- if(PixelReadConverter->GetRGB())
+ if ( PixelReadConverter->GetRGB() )
{
PixelWriteConverter->SetReadData(PixelReadConverter->GetRGB(),
PixelReadConverter->GetRGBSize());
}
std::string vr = "OB";
- if( FileInternal->GetBitsAllocated()>8 )
+ if ( FileInternal->GetBitsAllocated()>8 )
vr = "OW";
- if( FileInternal->GetBitsAllocated()==24 ) // For RGB ACR files
+ if ( FileInternal->GetBitsAllocated()==24 ) // For RGB ACR files
vr = "OB";
BinEntry *pixel =
CopyBinEntry(GetFile()->GetGrPixel(),GetFile()->GetNumPixel(),vr);
Archive->Push(0x0028,0x1202);
Archive->Push(0x0028,0x1203);
+ // push out Palette Color Lookup Table UID, if any
+ Archive->Push(0x0028,0x1199);
+
// For old '24 Bits' ACR-NEMA
// Thus, we have a RGB image and the bits allocated = 24 and
// samples per pixels = 1 (in the read file)
- if(FileInternal->GetBitsAllocated()==24)
+ if ( FileInternal->GetBitsAllocated()==24 )
{
ValEntry *bitsAlloc = CopyValEntry(0x0028,0x0100);
bitsAlloc->SetValue("8 ");
Archive->Restore(0x0028,0x1202);
Archive->Restore(0x0028,0x1203);
+ // For the Palette Color Lookup Table UID
+ Archive->Restore(0x0028,0x1203);
+
+
// group 0002 may be pushed out for ACR-NEMA writting purposes
Archive->Restore(0x0002,0x0000);
Archive->Restore(0x0002,0x0001);
ValEntry *oldCol = dynamic_cast<ValEntry *>
(FileInternal->GetDocEntry(0x0028, 0x0011));
- if( oldRow && oldCol )
+ if ( oldRow && oldCol )
{
std::string rows, columns;
{
ValEntry *recCode = dynamic_cast<ValEntry *>
(FileInternal->GetDocEntry(0x0008,0x0010));
- if( recCode )
+ if ( recCode )
{
- if( recCode->GetValue() == "ACRNEMA_LIBIDO_1.1" )
+ if ( recCode->GetValue() == "ACRNEMA_LIBIDO_1.1" )
{
ValEntry *libidoCode = CopyValEntry(0x0008,0x0010);
libidoCode->SetValue("");
DocEntry *oldE = FileInternal->GetDocEntry(group, elem);
ValEntry *newE;
- if( oldE )
+ if ( oldE )
{
newE = new ValEntry(oldE->GetDictEntry());
newE->Copy(oldE);
DocEntry *oldE = FileInternal->GetDocEntry(group, elem);
BinEntry *newE;
- if( oldE )
- if( oldE->GetVR()!=vr )
+ if ( oldE )
+ if ( oldE->GetVR()!=vr )
oldE = NULL;
- if( oldE )
+ if ( oldE )
{
newE = new BinEntry(oldE->GetDictEntry());
newE->Copy(oldE);
void FileHelper::CheckMandatoryElements()
{
// just to remember : 'official' 0002 group
-
+ if ( WriteType != ACR && WriteType != ACR_LIBIDO )
+ {
+ // Group 000002 (Meta Elements) already pushed out
+
//0002 0000 UL 1 Meta Group Length
//0002 0001 OB 1 File Meta Information Version
//0002 0002 UI 1 Media Stored SOP Class UID
// Create them if not found
// Always modify the value
// Push the entries to the archive.
- ValEntry *e_0002_0000 = CopyValEntry(0x0002,0x0000);
+ ValEntry *e_0002_0000 = CopyValEntry(0x0002,0x0000);
e_0002_0000->SetValue("0"); // for the moment
Archive->Push(e_0002_0000);
- BinEntry *e_0002_0001 = CopyBinEntry(0x0002,0x0001, "OB");
+ BinEntry *e_0002_0001 = CopyBinEntry(0x0002,0x0001, "OB");
e_0002_0001->SetBinArea((uint8_t*)Util::GetFileMetaInformationVersion(),
false);
e_0002_0001->SetLength(2);
Archive->Push(e_0002_0001);
// 'Media Stored SOP Class UID'
- ValEntry *e_0002_0002 = CopyValEntry(0x0002,0x0002);
+ ValEntry *e_0002_0002 = CopyValEntry(0x0002,0x0002);
// [Secondary Capture Image Storage]
e_0002_0002->SetValue("1.2.840.10008.5.1.4.1.1.7");
Archive->Push(e_0002_0002);
// 'Media Stored SOP Instance UID'
- ValEntry *e_0002_0003 = CopyValEntry(0x0002,0x0003);
+ ValEntry *e_0002_0003 = CopyValEntry(0x0002,0x0003);
e_0002_0003->SetValue(Util::CreateUniqueUID());
Archive->Push(e_0002_0003);
// 'Implementation Class UID'
- ValEntry *e_0002_0012 = CopyValEntry(0x0002,0x0012);
+ ValEntry *e_0002_0012 = CopyValEntry(0x0002,0x0012);
e_0002_0012->SetValue(Util::CreateUniqueUID());
Archive->Push(e_0002_0012);
// 'Implementation Version Name'
- ValEntry *e_0002_0013 = CopyValEntry(0x0002,0x0013);
- e_0002_0013->SetValue("GDCM 1.0");
+ ValEntry *e_0002_0013 = CopyValEntry(0x0002,0x0013);
+ e_0002_0013->SetValue("GDCM 1.1");
Archive->Push(e_0002_0013);
//'Source Application Entity Title' Not Mandatory
//ValEntry *e_0002_0016 = CopyValEntry(0x0002,0x0016);
// e_0002_0016->SetValue("1.2.840.10008.5.1.4.1.1.7");
// Archive->Push(e_0002_0016);
+ }
// Push out 'LibIDO-special' entries, if any
Archive->Push(0x0028,0x0015);
Archive->Push(e_0028_0100);
}
- // --- Check UID-related Entries ---
+ // Check if user wasn't drunk ;-)
+
+ std::ostringstream s;
+ // check 'Bits Allocated' vs decent values
+ int nbBitsAllocated = FileInternal->GetBitsAllocated();
+ if ( nbBitsAllocated == 0 || nbBitsAllocated > 32)
+ {
+ ValEntry *e_0028_0100 = CopyValEntry(0x0028,0x0100);
+ e_0028_0100->SetValue("16");
+ Archive->Push(e_0028_0100);
+ gdcmWarningMacro("(0028,0100) changed from "
+ << nbBitsAllocated << " to 16 for consistency purpose");
+ nbBitsAllocated = 16;
+ }
+ // check 'Bits Stored' vs 'Bits Allocated'
+ int nbBitsStored = FileInternal->GetBitsStored();
+ if ( nbBitsStored == 0 || nbBitsStored > nbBitsAllocated )
+ {
+ s << nbBitsAllocated;
+ ValEntry *e_0028_0101 = CopyValEntry(0x0028,0x0101);
+ e_0028_0101->SetValue( s.str() );
+ Archive->Push(e_0028_0101);
+ gdcmWarningMacro("(0028,0101) changed from "
+ << nbBitsStored << " to " << nbBitsAllocated
+ << " for consistency purpose" );
+ nbBitsStored = nbBitsAllocated;
+ }
+ // check 'Hight Bit Position' vs 'Bits Allocated' and 'Bits Stored'
+ int highBitPosition = FileInternal->GetHighBitPosition();
+ if ( highBitPosition == 0 ||
+ highBitPosition > nbBitsAllocated-1 ||
+ highBitPosition < nbBitsStored-1 )
+ {
+ ValEntry *e_0028_0102 = CopyValEntry(0x0028,0x0102);
+
+ s << nbBitsStored - 1;
+ e_0028_0102->SetValue( s.str() );
+ Archive->Push(e_0028_0102);
+ gdcmWarningMacro("(0028,0102) changed from "
+ << highBitPosition << " to " << nbBitsAllocated-1
+ << " for consistency purpose");
+ }
+ // --- Check UID-related Entries ---
// If 'SOP Class UID' exists ('true DICOM' image)
// we create the 'Source Image Sequence' SeqEntry
if ( e_0008_0016 != 0 )
{
// Create 'Source Image Sequence' SeqEntry
- SeqEntry *s = new SeqEntry (
+ SeqEntry *sis = new SeqEntry (
Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0008, 0x2112) );
SQItem *sqi = new SQItem(1);
// (we assume 'SOP Instance UID' exists too)
e_0008_1155->SetValue( e_0008_0018->GetValue());
sqi->AddEntry(e_0008_1155);
- s->AddSQItem(sqi,1);
+ sis->AddSQItem(sqi,1);
// temporarily replaces any previous 'Source Image Sequence'
- Archive->Push(s);
+ Archive->Push(sis);
// 'Image Type' (The written image is no longer an 'ORIGINAL' one)
ValEntry *e_0008_0008 = CopyValEntry(0x0008,0x0008);
Archive->Restore(0x0002,0x0001);
Archive->Restore(0x0002,0x0002);
Archive->Restore(0x0002,0x0003);
+ Archive->Restore(0x0002,0x0010);
Archive->Restore(0x0002,0x0012);
Archive->Restore(0x0002,0x0013);
Archive->Restore(0x0002,0x0016);
Archive->Restore(0x0002,0x0100);
Archive->Restore(0x0002,0x0102);
+
+ Archive->Restore(0x0008,0x0012);
+ Archive->Restore(0x0008,0x0013);
+ Archive->Restore(0x0008,0x0016);
+ Archive->Restore(0x0008,0x0018);
+ Archive->Restore(0x0008,0x0060);
+ Archive->Restore(0x0008,0x0070);
+ Archive->Restore(0x0008,0x0080);
+ Archive->Restore(0x0008,0x0090);
+ Archive->Restore(0x0008,0x2112);
+
+ Archive->Restore(0x0010,0x0010);
+ Archive->Restore(0x0010,0x0030);
+ Archive->Restore(0x0010,0x0040);
+
+ Archive->Restore(0x0020,0x000d);
+ Archive->Restore(0x0020,0x000e);
+
}
//-----------------------------------------------------------------------------
// The Raw image migth not be loaded yet:
std::ifstream *fp = FileInternal->OpenFile();
PixelReadConverter->ReadAndDecompressPixelData( fp );
- if(fp)
+ if ( fp )
FileInternal->CloseFile();
raw = PixelReadConverter->GetRaw();
return 0;
}
}
-
return raw;
}