fattori_g/reg23Topograms
0
0
Fork 0
Code Issues 1 Pull Requests Actions Packages Projects Releases Wiki Activity

Merge branch 'scoutRTUIDevelTemp' into 'ScoutRTUIDevel'

Browse Source 
Scout rtui devel temp into Devel

See merge request cpt_bioeng/drt!24
This commit is contained in:
via_r
2023-05-03 14:19:18 +00:00
parent 97674c6b4b 8cdb53731f
commit 24e7435376
2 changed files with 93 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

81
reg23Topograms/itkDTRrecon/itkImageProcessor.cpp
View File

@ -471,7 +471,7 @@ int itkImageProcessor::load3DSerieFromFiles( std::vector<std::string> v_fnames){
imageSeriesReader3D->SetImageIO(dicomIO);
imageSeriesReader3D->SetFileNames(v_sortedFnames);
//imageSeriesReader3D->SetNumberOfWorkUnits(8);
imageSeriesReader3D->SetNumberOfWorkUnits(8);
imageSeriesReader3D->ForceOrthogonalDirectionOff(); // properly read CTs with gantry tilt
//imageSeriesReader3D->SetSpacingWarningRelThreshold();
@ -1362,6 +1362,85 @@ void itkImageProcessor::CalculateExternalUserTransform(TransformType::Pointer tr
m_TransformMetaInfo->SetUserRotations(LPStoIEC_Directions * rotationUser);
}
itkImageProcessor::TransformType::Pointer
itkImageProcessor::CalculateInternalTransformV2(
ImageType3D::PointType m_ImportOffsetLPS,
ImageType3D::PointType m_RotationOffsetIEC,
ImageType3D::PointType m_TranslationUserIEC,
ImageType3D::PointType m_RotationUserIEC,
ImageType3D::PointType m_ProjectionCenterLPS,
ImageType3D::PointType m_RTIsocenterLPS,
InternalImageType::DirectionType m_IECtoLPSDirections
){
const double dtr = (atan(1.0) * 4.0) / 180.0;
ImageType3D::PointType m_ROffsetLPS =
m_IECtoLPSDirections * m_RotationOffsetIEC;
ImageType3D::PointType m_TUserLPS =
m_IECtoLPSDirections * m_TranslationUserIEC;
ImageType3D::PointType m_RUserLPS =
m_IECtoLPSDirections * m_RotationUserIEC;
TransformType::OutputVectorType translation;
translation[0] = m_TUserLPS[0];
translation[1] = m_TUserLPS[1];
translation[2] = m_TUserLPS[2];
ImageType3D::PointType m_TransformRotations;
m_TransformRotations[0] = m_ROffsetLPS[0] + m_RUserLPS[0] ;
m_TransformRotations[1] = m_ROffsetLPS[1] + m_RUserLPS[1] ;
m_TransformRotations[2] = m_ROffsetLPS[2] + m_RUserLPS[2] ;
ImageType3D::PointType m_RTIsocenterLPSNEG;
m_RTIsocenterLPSNEG[0] = -m_RTIsocenterLPS[0] + m_ProjectionCenterLPS[0];
m_RTIsocenterLPSNEG[1] = -m_RTIsocenterLPS[1] + m_ProjectionCenterLPS[0];
m_RTIsocenterLPSNEG[2] = -m_RTIsocenterLPS[2] + m_ProjectionCenterLPS[0];
// Map user to the projection center
TransformType::Pointer m_IsoTransformZeroAtCTOrigin =
MapTransformToNewOrigin (
m_RTIsocenterLPSNEG,
translation,
m_TransformRotations
);
TransformType::Pointer m_OutputTransform =
TransformType::New();
m_OutputTransform ->SetComputeZYX(true);
m_OutputTransform ->SetIdentity();
TransformType::OutputVectorType translation_Out;
translation_Out[0] =
m_IsoTransformZeroAtCTOrigin->GetTranslation()[0] +
m_ImportOffsetLPS[0] +
m_ProjectionCenterLPS[0];
translation_Out[1] =
m_IsoTransformZeroAtCTOrigin->GetTranslation()[1] +
m_ImportOffsetLPS[1] +
m_ProjectionCenterLPS[1];
translation_Out[2] =
m_IsoTransformZeroAtCTOrigin->GetTranslation()[2] +
m_ImportOffsetLPS[2] +
m_ProjectionCenterLPS[2];
m_OutputTransform->SetTranslation(
translation_Out);
m_OutputTransform->SetRotation(
m_IsoTransformZeroAtCTOrigin->GetAngleX(),
m_IsoTransformZeroAtCTOrigin->GetAngleY(),
m_IsoTransformZeroAtCTOrigin->GetAngleZ());
m_OutputTransform->SetCenter(
m_ProjectionCenterLPS);
return m_OutputTransform;
}
itkImageProcessor::TransformType::Pointer
itkImageProcessor::CalculateInternalTransformV2(

15
reg23Topograms/itkDTRrecon/itkImageProcessor.h
View File

@ -237,7 +237,6 @@ private:
using ResampleFilterType = itk::ResampleImageFilter<InternalImageType, InternalImageType>;
/** Optimizer which tries to find the minimn (Powell Optimizer)*/
// using OptimizerType = itk::PowellOptimizer;
using AmoebaOptimizerType = itk::AmoebaOptimizer;
/** Optimizer which samples the whol space */
using ExhaustiveOptimizerType = itk::ExhaustiveOptimizer;
@ -266,7 +265,9 @@ private:
using ITKtoVTKFilterType = itk::ImageToVTKImageFilter<OutputImageType>;
void
CalculateExternalUserTransform(TransformType::Pointer transform, DRTImageMetaInformation::Pointer imageMetaInfo);
CalculateExternalUserTransform(
TransformType::Pointer transform,
DRTImageMetaInformation::Pointer imageMetaInfo);
// TransformType::Pointer
// CalculateInternalTransform(
@ -293,6 +294,16 @@ private:
InternalImageType::DirectionType m_IECtoLPSDirections
);
TransformType::Pointer
CalculateInternalTransformV2(
ImageType3D::PointType m_ImportOffsetLPS,
ImageType3D::PointType m_RotationOffsetIEC,
ImageType3D::PointType m_TranslationUserIEC,
ImageType3D::PointType m_RotationUserIEC,
ImageType3D::PointType m_ProjectionCenterLPS,
ImageType3D::PointType m_RTIsocenterLPS,
InternalImageType::DirectionType m_IECtoLPSDirections
);
TransformType::Pointer
transform1,