1 /*=========================================================================
4 Module: $RCSfile: TestWriteSimple.cxx,v $
6 Date: $Date: 2005/01/08 23:18:31 $
7 Version: $Revision: 1.11 $
9 Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
10 l'Image). All rights reserved. See Doc/License.txt or
11 http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details.
13 This software is distributed WITHOUT ANY WARRANTY; without even
14 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
15 PURPOSE. See the above copyright notices for more information.
17 =========================================================================*/
20 * Write a dicom file from nothing
21 * The written image is 256x256, 8 bits, unsigned char
22 * The image content is a horizontal grayscale from
25 #include "gdcmHeader.h"
33 int sizeX; // Size X of the image
34 int sizeY; // Size Y of the image
35 int sizeZ; // Size Z of the image
36 int components; // Number of components for a pixel
37 int componentSize; // Component size (in bits : 8, 16)
38 int componentUse ; // Component size (in bits)
39 int sign; // Sign of components
40 char writeMode; // Write mode
42 // - 'e' : Explicit VR
43 // - 'i' : Implicit VR
47 {256, 256, 1, 1, 8, 8, 0, 'a'},
48 {256, 256, 1, 1, 8, 8, 0, 'e'},
49 {256, 256, 1, 1, 8, 8, 0, 'i'},
51 {512, 256, 1, 1, 8, 8, 0, 'a'},
52 {512, 256, 1, 1, 8, 8, 0, 'e'},
53 {512, 256, 1, 1, 8, 8, 0, 'i'},
55 {256, 512, 1, 1, 8, 8, 0, 'a'},
56 {256, 512, 1, 1, 8, 8, 0, 'e'},
57 {256, 512, 1, 1, 8, 8, 0, 'i'},
59 {256, 512, 1, 1, 16, 16, 0, 'a'},
60 {256, 512, 1, 1, 16, 16, 0, 'e'},
61 {256, 512, 1, 1, 16, 16, 0, 'i'},
62 {256, 512, 1, 1, 16, 16, 0, 'a'},
63 {256, 512, 1, 1, 16, 16, 0, 'e'},
64 {256, 512, 1, 1, 16, 16, 0, 'i'},
65 {0, 0, 1, 1, 8, 8, 0, 'i'} // to find the end
68 int WriteSimple(Image &img)
70 std::string fileName = "TestWriteSimple.dcm";
72 // Step 1 : Create the header of the image
74 gdcm::Header *header = new gdcm::Header();
75 std::ostringstream str;
80 header->ReplaceOrCreate(str.str(),0x0028,0x0011); // Columns
84 header->ReplaceOrCreate(str.str(),0x0028,0x0010); // Rows
90 header->ReplaceOrCreate(str.str(),0x0028,0x0008); // Number of Frames
95 str << img.componentSize;
96 header->ReplaceOrCreate(str.str(),0x0028,0x0100); // Bits Allocated
99 str << img.componentUse;
100 header->ReplaceOrCreate(str.str(),0x0028,0x0101); // Bits Stored
103 str << img.componentSize - 1;
104 header->ReplaceOrCreate(str.str(),0x0028,0x0102); // High Bit
106 // Set the pixel representation
109 header->ReplaceOrCreate(str.str(),0x0028,0x0103); // Pixel Representation
111 // Set the samples per pixel
113 str << img.components;
114 header->ReplaceOrCreate(str.str(),0x0028,0x0002); // Samples per Pixel
116 if( !header->IsReadable() )
118 std::cout << "Failed\n"
119 << " Prepared image isn't readable\n";
125 // Step 2 : Create the output image
127 if( img.componentSize%8 > 0 )
129 img.componentSize += 8-img.componentSize%8;
131 size_t size = img.sizeX * img.sizeY * img.sizeZ
132 * img.components * img.componentSize / 8;
133 unsigned char *imageData = new unsigned char[size];
135 // FIXME : find a best heuristic to create the image
136 unsigned char *tmp = imageData;
137 for(int k=0;k<img.sizeZ;k++)
139 for(int j=0;j<img.sizeY;j++)
141 for(int i=0;i<img.sizeX;i++)
143 for(int c=0;c<img.components;c++)
146 if( img.componentSize>8 )
150 tmp += img.componentSize * img.components/8;
156 // Step 3 : Create the file of the image
158 gdcm::File *file = new gdcm::File(header);
159 file->SetImageData(imageData,size);
161 // Step 4 : Set the writting mode and write the image
164 file->SetWriteModeToRaw();
165 switch (img.writeMode)
167 case 'a' : // Write an ACR file
168 file->SetWriteTypeToAcr();
171 case 'e' : // Write a DICOM Explicit VR file
172 file->SetWriteTypeToDcmExplVR();
175 case 'i' : // Write a DICOM Implicit VR file
176 file->SetWriteTypeToDcmImplVR();
180 std::cout << "Failed\n"
181 << " Write mode '"<<img.writeMode<<"' is undefined\n";
189 if( !file->Write(fileName) )
191 std::cout << "Failed\n"
192 << "File in unwrittable\n";
200 // Step 5 : Read the written image
202 gdcm::File* reread = new gdcm::File( fileName );
203 if( !reread->GetHeader()->IsReadable() )
205 std::cerr << "Failed" << std::endl
206 << "Test::TestReadWriteReadCompare: Could not reread image "
207 << "written:" << fileName << std::endl;
214 // Step 6 : Compare to the written image
216 size_t dataSizeWritten = reread->GetImageDataSize();
217 uint8_t* imageDataWritten = reread->GetImageData();
219 // Test the image size
220 if (header->GetXSize() != reread->GetHeader()->GetXSize() ||
221 header->GetYSize() != reread->GetHeader()->GetYSize() ||
222 header->GetZSize() != reread->GetHeader()->GetZSize())
224 std::cout << "Failed" << std::endl
225 << " X Size differs: "
226 << "X: " << header->GetXSize() << " # "
227 << reread->GetHeader()->GetXSize() << " | "
228 << "Y: " << header->GetYSize() << " # "
229 << reread->GetHeader()->GetYSize() << " | "
230 << "Z: " << header->GetZSize() << " # "
231 << reread->GetHeader()->GetZSize() << std::endl;
240 // Test the data size
241 if (size != dataSizeWritten)
243 std::cout << "Failed" << std::endl
244 << " Pixel areas lengths differ: "
245 << size << " # " << dataSizeWritten << std::endl;
254 // Test the data's content
255 if (int res = memcmp(imageData, imageDataWritten, size) !=0)
258 std::cout << "Failed" << std::endl
259 << " Pixel differ (as expanded in memory)." << std::endl;
268 std::cout << "OK" << std::endl;
278 int TestWriteSimple(int argc, char* argv[])
282 std::cerr << "usage: \n"
283 << argv[0] << " (without parameters) " << std::endl
290 while( Images[i].sizeX>0 && Images[i].sizeY>0 )
292 ret += WriteSimple(Images[i]);