- Hnd_header hnd;
- FILE *fp;
-
- fp = fopen (m_FileName.c_str(), "rb");
- if (fp == NULL)
- itkGenericExceptionMacro(<< "Could not open file (for reading): " << m_FileName);
-
- fread ((void *) hnd.sFileType, sizeof(char), 32, fp);
- fread ((void *) &hnd.FileLength, sizeof(uint32_t), 1, fp);
- fread ((void *) hnd.sChecksumSpec, sizeof(char), 4, fp);
- fread ((void *) &hnd.nCheckSum, sizeof(uint32_t), 1, fp);
- fread ((void *) hnd.sCreationDate, sizeof(char), 8, fp);
- fread ((void *) hnd.sCreationTime, sizeof(char), 8, fp);
- fread ((void *) hnd.sPatientID, sizeof(char), 16, fp);
- fread ((void *) &hnd.nPatientSer, sizeof(uint32_t), 1, fp);
- fread ((void *) hnd.sSeriesID, sizeof(char), 16, fp);
- fread ((void *) &hnd.nSeriesSer, sizeof(uint32_t), 1, fp);
- fread ((void *) hnd.sSliceID, sizeof(char), 16, fp);
- fread ((void *) &hnd.nSliceSer, sizeof(uint32_t), 1, fp);
- fread ((void *) &hnd.SizeX, sizeof(uint32_t), 1, fp);
- fread ((void *) &hnd.SizeY, sizeof(uint32_t), 1, fp);
- fread ((void *) &hnd.dSliceZPos, sizeof(double), 1, fp);
- fread ((void *) hnd.sModality, sizeof(char), 16, fp);
- fread ((void *) &hnd.nWindow, sizeof(uint32_t), 1, fp);
- fread ((void *) &hnd.nLevel, sizeof(uint32_t), 1, fp);
- fread ((void *) &hnd.nPixelOffset, sizeof(uint32_t), 1, fp);
- fread ((void *) hnd.sImageType, sizeof(char), 4, fp);
- fread ((void *) &hnd.dGantryRtn, sizeof(double), 1, fp);
- fread ((void *) &hnd.dSAD, sizeof(double), 1, fp);
- fread ((void *) &hnd.dSFD, sizeof(double), 1, fp);
- fread ((void *) &hnd.dCollX1, sizeof(double), 1, fp);
- fread ((void *) &hnd.dCollX2, sizeof(double), 1, fp);
- fread ((void *) &hnd.dCollY1, sizeof(double), 1, fp);
- fread ((void *) &hnd.dCollY2, sizeof(double), 1, fp);
- fread ((void *) &hnd.dCollRtn, sizeof(double), 1, fp);
- fread ((void *) &hnd.dFieldX, sizeof(double), 1, fp);
- fread ((void *) &hnd.dFieldY, sizeof(double), 1, fp);
- fread ((void *) &hnd.dBladeX1, sizeof(double), 1, fp);
- fread ((void *) &hnd.dBladeX2, sizeof(double), 1, fp);
- fread ((void *) &hnd.dBladeY1, sizeof(double), 1, fp);
- fread ((void *) &hnd.dBladeY2, sizeof(double), 1, fp);
- fread ((void *) &hnd.dIDUPosLng, sizeof(double), 1, fp);
- fread ((void *) &hnd.dIDUPosLat, sizeof(double), 1, fp);
- fread ((void *) &hnd.dIDUPosVrt, sizeof(double), 1, fp);
- fread ((void *) &hnd.dIDUPosRtn, sizeof(double), 1, fp);
- fread ((void *) &hnd.dPatientSupportAngle, sizeof(double), 1, fp);
- fread ((void *) &hnd.dTableTopEccentricAngle, sizeof(double), 1, fp);
- fread ((void *) &hnd.dCouchVrt, sizeof(double), 1, fp);
- fread ((void *) &hnd.dCouchLng, sizeof(double), 1, fp);
- fread ((void *) &hnd.dCouchLat, sizeof(double), 1, fp);
- fread ((void *) &hnd.dIDUResolutionX, sizeof(double), 1, fp);
- fread ((void *) &hnd.dIDUResolutionY, sizeof(double), 1, fp);
- fread ((void *) &hnd.dImageResolutionX, sizeof(double), 1, fp);
- fread ((void *) &hnd.dImageResolutionY, sizeof(double), 1, fp);
- fread ((void *) &hnd.dEnergy, sizeof(double), 1, fp);
- fread ((void *) &hnd.dDoseRate, sizeof(double), 1, fp);
- fread ((void *) &hnd.dXRayKV, sizeof(double), 1, fp);
- fread ((void *) &hnd.dXRayMA, sizeof(double), 1, fp);
- fread ((void *) &hnd.dMetersetExposure, sizeof(double), 1, fp);
- fread ((void *) &hnd.dAcqAdjustment, sizeof(double), 1, fp);
- fread ((void *) &hnd.dCTProjectionAngle, sizeof(double), 1, fp);
- fread ((void *) &hnd.dCTNormChamber, sizeof(double), 1, fp);
- fread ((void *) &hnd.dGatingTimeTag, sizeof(double), 1, fp);
- fread ((void *) &hnd.dGating4DInfoX, sizeof(double), 1, fp);
- fread ((void *) &hnd.dGating4DInfoY, sizeof(double), 1, fp);
- fread ((void *) &hnd.dGating4DInfoZ, sizeof(double), 1, fp);
- fread ((void *) &hnd.dGating4DInfoTime, sizeof(double), 1, fp);
- fclose (fp);
-
- /* Convert hnd to ITK image information */
- SetNumberOfDimensions(2);
- SetDimensions(0, hnd.SizeX);
- SetDimensions(1, hnd.SizeY);
- SetComponentType(itk::ImageIOBase::UINT);
+ Hnd_header hnd;
+ FILE *fp;
+
+ fp = fopen (m_FileName.c_str(), "rb");
+ if (fp == NULL)
+ itkGenericExceptionMacro(<< "Could not open file (for reading): " << m_FileName);
+
+ fread ((void *) hnd.sFileType, sizeof(char), 32, fp);
+ fread ((void *) &hnd.FileLength, sizeof(uint32_t), 1, fp);
+ fread ((void *) hnd.sChecksumSpec, sizeof(char), 4, fp);
+ fread ((void *) &hnd.nCheckSum, sizeof(uint32_t), 1, fp);
+ fread ((void *) hnd.sCreationDate, sizeof(char), 8, fp);
+ fread ((void *) hnd.sCreationTime, sizeof(char), 8, fp);
+ fread ((void *) hnd.sPatientID, sizeof(char), 16, fp);
+ fread ((void *) &hnd.nPatientSer, sizeof(uint32_t), 1, fp);
+ fread ((void *) hnd.sSeriesID, sizeof(char), 16, fp);
+ fread ((void *) &hnd.nSeriesSer, sizeof(uint32_t), 1, fp);
+ fread ((void *) hnd.sSliceID, sizeof(char), 16, fp);
+ fread ((void *) &hnd.nSliceSer, sizeof(uint32_t), 1, fp);
+ fread ((void *) &hnd.SizeX, sizeof(uint32_t), 1, fp);
+ fread ((void *) &hnd.SizeY, sizeof(uint32_t), 1, fp);
+ fread ((void *) &hnd.dSliceZPos, sizeof(double), 1, fp);
+ fread ((void *) hnd.sModality, sizeof(char), 16, fp);
+ fread ((void *) &hnd.nWindow, sizeof(uint32_t), 1, fp);
+ fread ((void *) &hnd.nLevel, sizeof(uint32_t), 1, fp);
+ fread ((void *) &hnd.nPixelOffset, sizeof(uint32_t), 1, fp);
+ fread ((void *) hnd.sImageType, sizeof(char), 4, fp);
+ fread ((void *) &hnd.dGantryRtn, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dSAD, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dSFD, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dCollX1, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dCollX2, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dCollY1, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dCollY2, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dCollRtn, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dFieldX, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dFieldY, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dBladeX1, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dBladeX2, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dBladeY1, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dBladeY2, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dIDUPosLng, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dIDUPosLat, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dIDUPosVrt, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dIDUPosRtn, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dPatientSupportAngle, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dTableTopEccentricAngle, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dCouchVrt, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dCouchLng, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dCouchLat, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dIDUResolutionX, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dIDUResolutionY, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dImageResolutionX, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dImageResolutionY, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dEnergy, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dDoseRate, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dXRayKV, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dXRayMA, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dMetersetExposure, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dAcqAdjustment, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dCTProjectionAngle, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dCTNormChamber, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dGatingTimeTag, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dGating4DInfoX, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dGating4DInfoY, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dGating4DInfoZ, sizeof(double), 1, fp);
+ fread ((void *) &hnd.dGating4DInfoTime, sizeof(double), 1, fp);
+ fclose (fp);
+
+ /* Convert hnd to ITK image information */
+ SetNumberOfDimensions(2);
+ SetDimensions(0, hnd.SizeX);
+ SetDimensions(1, hnd.SizeY);
+ SetSpacing(0, hnd.dIDUResolutionX);
+ SetSpacing(1, hnd.dIDUResolutionY);
+ SetOrigin(0, -0.5*(hnd.SizeX-1)*hnd.dIDUResolutionX); //SR: assumed centered
+ SetOrigin(1, -0.5*(hnd.SizeY-1)*hnd.dIDUResolutionY); //SR: assumed centered
+ SetComponentType(itk::ImageIOBase::UINT);
+
+ /* Store important meta information in the meta data dictionary */
+ itk::EncapsulateMetaData<double>(this->GetMetaDataDictionary(), "dCTProjectionAngle", hnd.dCTProjectionAngle);
+ //itk::ExposeMetaData<double>( this->GetMetaDataDictionary(), &(hnd.dCTProjectionAngle), "dCTProjectionAngle");