]> Creatis software - gdcm.git/blob - WriteDicomSimple.cxx
97bfe6dc14c8a6ee3802ba56b8134ae28c190d41
[gdcm.git] / WriteDicomSimple.cxx
1 /*=========================================================================
2                                                                                 
3   Program:   gdcm
4   Module:    $RCSfile: WriteDicomSimple.cxx,v $
5   Language:  C++
6   Date:      $Date: 2005/01/21 11:40:53 $
7   Version:   $Revision: 1.8 $
8                                                                                 
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.
12                                                                                 
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.
16                                                                                 
17 =========================================================================*/
18
19 /**
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 
23  * 
24  */
25 #include "gdcmFile.h"
26 #include "gdcmFileHelper.h"
27
28 #include <iostream>
29 #include <sstream>
30
31 // Image size
32 #define SIZE_X          256
33 #define SIZE_Y          256
34 // Number of components in the image (3 for RGB)
35 #define COMPONENT       1
36 // Size of each component (in byte)
37 #define COMPONENT_SIZE  1
38 // Window / Level
39 #define COLOR_WINDOW    256
40 #define COLOR_LEVEL     128
41
42 int main(int argc, char* argv[])
43 {
44    if (argc < 3) 
45    {
46       std::cerr << "usage: \n" 
47                 << argv[0] << " Output Mode " << std::endl 
48                 << "Output : output file name\n"
49                 << "Mode : \n"
50                 << "   a : ACR, produces a file named Output.ACR\n"
51                 << "   e : DICOM Explicit VR, produces a file named Output.E.DCM\n"
52                 << "   i : DICOM Implicit VR, produces a file named Output.I.DCM\n"
53                 << "   r : RAW, produces a file named Output.RAW\n"
54                 << std::endl;
55       return 1;
56    }
57
58
59 // Step 1 : Create the header of the image
60    gdcm::File *header = new gdcm::File();
61    std::ostringstream str;
62
63    // Set the image size
64    str.str("");
65    str << SIZE_X;
66    header->ReplaceOrCreate(str.str(),0x0028,0x0011); // Columns
67
68    str.str("");
69    str << SIZE_Y;
70    header->ReplaceOrCreate(str.str(),0x0028,0x0010); // Rows
71
72    // Set the pixel type
73    str.str("");
74    str << COMPONENT_SIZE * 8;
75    header->ReplaceOrCreate(str.str(),0x0028,0x0100); // Bits Allocated
76    header->ReplaceOrCreate(str.str(),0x0028,0x0101); // Bits Stored
77
78    str.str("");
79    str << COMPONENT_SIZE * 8 - 1;
80    header->ReplaceOrCreate(str.str(),0x0028,0x0102); // High Bit
81
82    // Set the pixel representation
83    str.str("");
84    str << "0"; // Unsigned
85    header->ReplaceOrCreate(str.str(),0x0028,0x0103); // Pixel Representation
86
87    // Set the samples per pixel
88    str.str("");
89    str << COMPONENT;
90    header->ReplaceOrCreate(str.str(),0x0028,0x0002); // Samples per Pixel
91
92    // Set the Window / Level
93    str.str("");
94    str << COLOR_WINDOW;
95    header->ReplaceOrCreate(str.str(),0x0028,0x1051); // Window Width
96    str.str("");
97    str << COLOR_LEVEL;
98    header->ReplaceOrCreate(str.str(),0x0028,0x1050); // Window Center
99
100    if( !header->IsReadable() )
101    {
102       std::cerr << "-------------------------------\n"
103                 << "Error while creating the file\n"
104                 << "This file is considered to be not readable\n";
105
106       return 1;
107    }
108
109 // Step 2 : Create the output image
110    size_t size = SIZE_X * SIZE_Y * COMPONENT * COMPONENT_SIZE;
111    unsigned char *imageData = new unsigned char[size];
112
113    unsigned char *tmp = imageData;
114    for(int j=0;j<SIZE_Y;j++)
115    {
116       for(int i=0;i<SIZE_X;i++)
117       {
118          for(int c=0;c<COMPONENT;c++)
119          {
120             *tmp = j;
121             tmp += COMPONENT_SIZE; 
122          }
123       }
124    }
125
126 // Step 3 : Create the file of the image
127    gdcm::FileHelper *file = new gdcm::FileHelper(header);
128    file->SetImageData(imageData,size);
129
130 // Step 4 : Set the writting mode and write the image
131    std::string fileName = argv[1]; 
132    std::string mode = argv[2];
133
134    file->SetWriteModeToRaw();
135    switch (mode[0])
136    {
137       case 'a' : // Write an ACR file
138          fileName += ".ACR";
139          file->SetWriteTypeToAcr();
140          std::cout << "Write ACR" << std::endl
141                    << "File :" << fileName << std::endl;
142          break;
143
144       case 'e' : // Write a DICOM Explicit VR file
145          fileName += ".E.DCM";
146          file->SetWriteTypeToDcmExplVR();
147          std::cout << "Write DCM Explicit VR" << std::endl
148                    << "File :" << fileName << std::endl;
149          break;
150
151       case 'i' : // Write a DICOM Implicit VR file
152          fileName += ".I.DCM";
153          file->SetWriteTypeToDcmImplVR();
154          std::cout << "Write DCM Implicit VR" << std::endl
155                    << "File :" << fileName << std::endl;
156          break;
157
158       case 'r' : // Write a RAW file
159          fileName += ".RAW";
160          file->WriteRawData(fileName);
161          std::cout << "Write Raw" << std::endl
162                    << "File :" << fileName << std::endl;
163
164          delete[] imageData;
165          delete file;
166          delete header;
167          return 0;
168
169       default :
170          std::cout << "-------------------------------\n"
171                    << "Write mode undefined...\n"
172                    << "No file written\n";
173
174          delete[] imageData;
175          delete file;
176          delete header;
177          return 1;
178    }
179
180    if( !file->Write(fileName) )
181    {
182       std::cout << "-------------------------------\n"
183                    << "Error when writting the file " << fileName << "\n"
184                 << "No file written\n";
185    }
186
187    delete[] imageData;
188    delete file;
189    delete header;
190
191    return 0;
192 }
193