1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136
|
/*=========================================================================
Program: GDCM (Grassroots DICOM). A DICOM library
Copyright (c) 2006-2011 Mathieu Malaterre
All rights reserved.
See Copyright.txt or http://gdcm.sourceforge.net/Copyright.html for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notice for more information.
=========================================================================*/
#include "gdcmImageRegionReader.h"
#include "gdcmImageHelper.h"
#include "gdcmFilename.h"
#include "gdcmBoxRegion.h"
#include "gdcmTesting.h"
static int TestImageRegionRead(const char* filename, bool verbose = false)
{
if( verbose )
std::cerr << "Reading: " << filename << std::endl;
gdcm::ImageRegionReader reader;
gdcm::Filename fn( filename );
// DMCPACS_ExplicitImplicit_BogusIOP.dcm is very difficult to handle since
// we need to read 3 attribute to detect the "well known" bug. However the third
// attribute is (b500,b700) which make ReadUpToTag(7fe0,0010) fails...
if( strcmp(fn.GetName(), "DMCPACS_ExplicitImplicit_BogusIOP.dcm" ) == 0
|| strcmp(fn.GetName(), "SC16BitsAllocated_8BitsStoredJ2K.dcm" ) == 0 // mismatch pixel format in JPEG 2000 vs DICOM
|| strcmp(fn.GetName(), "SC16BitsAllocated_8BitsStoredJPEG.dcm" ) == 0 // mismatch pixel format in JPEG vs DICOM
|| strcmp(fn.GetName(), "PHILIPS_Gyroscan-12-Jpeg_Extended_Process_2_4.dcm" ) == 0 // bogus JPEG cannot be streamed
|| strcmp(fn.GetName(), "US-YBR_FULL_422-EVRLE.dcm" ) == 0 // FIXME TODO
// FIXME: we should be able to handle those at some point:
|| strcmp(fn.GetName(), "JPEGNote_empty.dcm" ) == 0
|| strcmp(fn.GetName(), "JPEGNote_missing.dcm" ) == 0
|| strcmp(fn.GetName(), "JPEGNote_bogus.dcm" ) == 0
)
{
std::cerr << "Skipping impossible file: " << filename << std::endl;
return 0;
}
reader.SetFileName( filename );
bool canReadInformation = reader.ReadInformation();
if (!canReadInformation)
{
//std::cerr << "Cannot ReadInformation: " << filename << std::endl;
return 0; //unable to read tags as expected.
}
int res = 0;
//std::cout << reader.GetFile().GetDataSet() << std::endl;
std::vector<unsigned int> dims =
gdcm::ImageHelper::GetDimensionsValue(reader.GetFile());
gdcm::BoxRegion box;
box.SetDomain(0, dims[0] - 1, 0, dims[1] - 1, 0, dims[2] - 1);
reader.SetRegion( box );
size_t len = reader.ComputeBufferLength();
if( !len )
{
std::cerr << "No length for: " << filename << std::endl;
return 1;
}
std::vector<char> vbuffer;
vbuffer.resize( len );
char* buffer = vbuffer.data();
bool b = reader.ReadIntoBuffer(buffer, len);
if( !b )
{
std::cerr << "Could not ReadIntoBuffer: " << filename << std::endl;
return 1;
}
#if 0
std::ofstream of( "/tmp/dd.raw", std::ios::binary );
of.write( buffer, len );
of.close();
#endif
const char *ref = gdcm::Testing::GetMD5FromFile(filename);
char digest[33];
gdcm::Testing::ComputeMD5(buffer, len, digest);
// FIXME: PC=1
if( strcmp(fn.GetName(), "ACUSON-24-YBR_FULL-RLE-b.dcm" ) == 0) ref = "2d7a28cae6c3b3183284d1b4ae08307f";
if( strcmp(fn.GetName(), "ACUSON-24-YBR_FULL-RLE.dcm" ) == 0) ref = "429f31f0b70bd515b3feeda5dea5eac0";
if( verbose )
{
if(ref) std::cout << "ref=" << ref << std::endl;
std::cout << "md5=" << digest << std::endl;
}
if( !ref )
{
// new regression image needs a md5 sum
std::cout << "Missing md5 " << digest << " for: " << filename << std::endl;
res = 1;
}
else if( strcmp(digest, ref) != 0 )
{
std::cerr << "Problem reading image from: " << filename << std::endl;
std::cerr << "Found " << digest << " instead of " << ref << std::endl;
res = 1;
}
return res;
}
int TestImageRegionReader1(int argc, char *argv[])
{
if( argc == 2 )
{
const char *filename = argv[1];
return TestImageRegionRead(filename, true);
}
// else
// First of get rid of warning/debug message
gdcm::Trace::DebugOff();
gdcm::Trace::WarningOff();
int r = 0, i = 0;
const char *filename;
const char * const *filenames = gdcm::Testing::GetFileNames();
while( (filename = filenames[i]) )
{
r += TestImageRegionRead(filename);
++i;
}
return r;
}
|