vtkScancoCTReader.h

/*=========================================================================
 
  Program: DICOM for VTK
 
  Copyright (c) 2015 David Gobbi
  All rights reserved.
  See Copyright.txt or http://dgobbi.github.io/bsd3.txt 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.
 
=========================================================================*/
#ifndef vtkScancoCTReader_h
#define vtkScancoCTReader_h
 
#include "vtkImageReader2.h"
#include "vtkDICOMModule.h" // For export macro
#include "vtkDICOMConfig.h" // For configuration details
 
#include <iosfwd> // for std::istream
 
//----------------------------------------------------------------------------
class VTKDICOM_EXPORT vtkScancoCTReader : public vtkImageReader2
{
public:
  static vtkScancoCTReader *New();
  vtkTypeMacro(vtkScancoCTReader, vtkImageReader2);
 
  void PrintSelf(ostream& os, vtkIndent indent) VTK_DICOM_OVERRIDE;
 
  const char* GetFileExtensions() VTK_DICOM_OVERRIDE {
    return ".isq .rsq .rad .aim" ; }
 
  const char* GetDescriptiveName() VTK_DICOM_OVERRIDE {
    return "SCANCO MicroCT"; }
 
  int CanReadFile(const char* filename) VTK_DICOM_OVERRIDE;
 
  const char *GetVersion() { return this->Version; }
 
  int GetPatientIndex() { return this->PatientIndex; }
 
  int GetScannerID() { return this->ScannerID; }
 
  const char *GetCreationDate() { return this->CreationDate; }
 
  const char *GetModificationDate() { return this->ModificationDate; }
 
  double GetSliceThickness() { return this->SliceThickness; }
 
  double GetSliceIncrement() { return this->SliceIncrement; }
 
  double GetStartPosition() { return this->StartPosition; }
 
  double GetEndPosition() { return this->EndPosition; }
 
  double GetZPosition() { return this->ZPosition; }
 
  vtkGetVector2Macro(DataRange, double);
 
  double GetMuScaling() { return this->MuScaling; }
 
  int GetNumberOfSamples() { return this->NumberOfSamples; }
 
  int GetNumberOfProjections() { return this->NumberOfProjections; }
 
  double GetScanDistance() { return this->ScanDistance; }
 
  int GetScannerType() { return this->ScannerType; }
 
  double GetSampleTime() { return this->SampleTime; }
 
  int GetMeasurementIndex() { return this->MeasurementIndex; }
 
  int GetSite() { return this->Site; }
 
  double GetReferenceLine() { return this->ReferenceLine; }
 
  int GetReconstructionAlg() { return this->ReconstructionAlg; }
 
  const char *GetPatientName() { return this->PatientName; }
 
  double GetEnergy() { return this->Energy; }
 
  double GetIntensity() { return this->Intensity; }
 
 
  double GetRescaleSlope() { return this->RescaleSlope; }
  double GetRescaleIntercept() { return this->RescaleIntercept; }
 
  int GetRescaleType() { return this->RescaleType; }
 
  const char *GetRescaleUnits() { return this->RescaleUnits; }
 
  double GetMuWater() { return this->MuWater; }
 
  const char *GetCalibrationData() { return this->CalibrationData; }
 
  vtkGetVector3Macro(ScanDimensionsPixels, int);
 
  vtkGetVector3Macro(ScanDimensionsPhysical, double);
 
  void *GetRawHeader() { return this->RawHeader; }
 
protected:
  vtkScancoCTReader();
  ~vtkScancoCTReader() VTK_DICOM_OVERRIDE;
 
  int RequestInformation(
    vtkInformation* request, vtkInformationVector** inputVector,
    vtkInformationVector* outputVector) VTK_DICOM_OVERRIDE;
 
  int RequestData(
    vtkInformation* request, vtkInformationVector** inputVector,
    vtkInformationVector* outputVector) VTK_DICOM_OVERRIDE;
 
  void InitializeHeader();
 
  int ReadISQHeader(std::istream *file, unsigned long bytesRead);
 
  int ReadAIMHeader(std::istream *file, unsigned long bytesRead);
 
 
  static int CheckVersion(const char header[16]);
 
  void DecodeDate(const void *data,
    int& year, int& month, int& day,
    int& hour, int& minute, int& second, int& millis);
 
  static int DecodeInt(const void *data);
 
  static float DecodeFloat(const void *data);
 
  static double StringToDouble(const char *cp, char **cpp=nullptr);
 
  static int StringToInt(const char *cp, char **cpp=nullptr);
 
  static double DecodeDouble(const void *data);
 
 
  static void StripString(char *dest, const char *source, size_t l);
 
  // Header information
  char Version[18];
  char PatientName[42];
  int PatientIndex;
  int ScannerID;
  char CreationDate[32];
  char ModificationDate[32];
  int ScanDimensionsPixels[3];
  double ScanDimensionsPhysical[3];
  double SliceThickness;
  double SliceIncrement;
  double StartPosition;
  double EndPosition;
  double ZPosition;
  double DataRange[2];
  double MuScaling;
  int NumberOfSamples;
  int NumberOfProjections;
  double ScanDistance;
  double SampleTime;
  int ScannerType;
  int MeasurementIndex;
  int Site;
  int ReconstructionAlg;
  double ReferenceLine;
  double Energy;
  double Intensity;
  int RescaleType;
  char RescaleUnits[18];
  char CalibrationData[66];
  double RescaleSlope;
  double RescaleIntercept;
  double MuWater;
  char *RawHeader;
 
  // The compression mode, if any.
  int Compression;
 
private:
#ifdef VTK_DICOM_DELETE
  vtkScancoCTReader(const vtkScancoCTReader&) VTK_DICOM_DELETE;
  void operator=(const vtkScancoCTReader&) VTK_DICOM_DELETE;
#else
  vtkScancoCTReader(const vtkScancoCTReader&) = delete;
  void operator=(const vtkScancoCTReader&) = delete;
#endif
};
 
#endif // vtkScancoCTReader_h