Function main

Examples/DicomSeriesFromArray/DicomSeriesFromArray.cxx:93–238

}

int
main(int argc, char * argv[])
{
  if (argc < 2)
  {
    std::cout << "Usage: DicomSeriesFromArray <output_directory>" << std::endl;
    return 1;
  }

  // Set pixel type sitkInt16 or sitkFloat64
  sitk::PixelIDValueEnum pixelType = sitk::sitkInt16;

  // Create image from random data array based on the pixel type
  std::vector<unsigned int> imageSize = { 512, 512, 256 };
  size_t                    totalPixels = imageSize[0] * imageSize[1] * imageSize[2];

  std::random_device rd;
  std::mt19937       gen(rd());
  sitk::Image        newImg;

  if (pixelType == sitk::sitkInt16)
  {
    std::uniform_int_distribution<short> dis(-1000, 1000);
    std::vector<short>                   randomData(totalPixels);
    for (size_t i = 0; i < totalPixels; ++i)
    {
      randomData[i] = dis(gen);
    }
    newImg = sitk::ImportAsInt16(randomData.data(), imageSize);
  }
  else if (pixelType == sitk::sitkFloat64)
  {
    std::uniform_real_distribution<> dis(-1000.0, 1000.0);
    std::vector<double>              randomData(totalPixels);
    for (size_t i = 0; i < totalPixels; ++i)
    {
      randomData[i] = dis(gen);
    }
    newImg = sitk::ImportAsFloat64(randomData.data(), imageSize);
  }
  else
  {
    std::cerr << "Unsupported pixel type" << std::endl;
    return 1;
  }

  newImg.SetSpacing({ 2.5, 3.5, 4.5 });

  // Write the 3D image as a series
  // IMPORTANT: There are many DICOM tags that need to be updated when you modify
  //            an original image. This is a delicate operation and requires
  //            knowledge of the DICOM standard. This example only modifies some.
  //            For a more complete list of tags that need to be modified see:
  //                  http://gdcm.sourceforge.net/wiki/index.php/Writing_DICOM
  //            If it is critical for your work to generate valid DICOM files,
  //            It is recommended to use David Clunie's Dicom3tools to validate
  //            the files:
  //                  http://www.dclunie.com/dicom3tools.html