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

PanelOffset for FFS

Browse Source
This commit is contained in:
Proton local user
2023-05-27 00:46:19 +02:00
parent 98ada50379
commit 8a34450119
3 changed files with 20 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
reg23Topograms/itkDTRrecon/DRTMetaInformation.cpp
View File

@ -94,6 +94,8 @@ DRTImageMetaInformation(){
this->m_OriginLPS.Fill(0.);
this->m_PanelOffsetLPS.Fill(0.);
}
void

6
reg23Topograms/itkDTRrecon/DRTMetaInformation.h
View File

@ -197,6 +197,9 @@ public:
itkSetMacro(ProjectionOriginLPSZero,PointType);
itkGetMacro(ProjectionOriginLPSZero,PointType);
itkSetMacro(PanelOffsetLPS,PointType);
itkGetMacro(PanelOffsetLPS,PointType);
itkSetMacro(ProjectionAngleLPS,double);
itkGetMacro(ProjectionAngleLPS,double);
@ -225,7 +228,8 @@ protected:
//m_Origin,
m_OriginLPS,
m_ProjectionOriginLPS,
m_ProjectionOriginLPSZero;
m_ProjectionOriginLPSZero,
m_PanelOffsetLPS;
double
m_ProjectionAngleLPS,

15
reg23Topograms/itkDTRrecon/itkImageProcessor.cpp
View File

@ -1252,7 +1252,8 @@ void itkImageProcessor::InitializeProjector()
interpolator1->SetThreshold(
m_DRTGeometryMetaInfo->GetIntensityThreshold() );
interpolator1->SetPanelOffset(
m_DRTGeometryMetaInfo->GetPanel1Offset());
m_DRTImage1MetaInfo->GetPanelOffsetLPS()[0]);
// m_DRTGeometryMetaInfo->GetPanel1Offset());
interpolator1->SetTransform(transform1);
interpolator1->Initialize();
@ -1289,7 +1290,8 @@ void itkImageProcessor::InitializeProjector()
interpolator2->SetThreshold(
m_DRTGeometryMetaInfo->GetIntensityThreshold() );
interpolator2->SetPanelOffset(
m_DRTGeometryMetaInfo->GetPanel2Offset());
m_DRTImage2MetaInfo->GetPanelOffsetLPS()[0]);
// m_DRTGeometryMetaInfo->GetPanel2Offset());
interpolator2->SetTransform(transform2);
interpolator2->Initialize();
@ -1571,6 +1573,11 @@ void itkImageProcessor::UpdateProjectionGeometryMeta(){
PanelOffsetIEC1[1] = 0.0;
PanelOffsetIEC1[2] = 0.0;
m_DRTImage1MetaInfo->SetPanelOffsetLPS(
m_CTMetaInfo->ConvertIECPointToLPSPoint(
PanelOffsetIEC1));
// This HAS TO be calculated from the nominal isocenter.
m_DRTImage1MetaInfo->SetOriginLPS(
CalcDRTImageOrigin(
@ -1671,6 +1678,10 @@ void itkImageProcessor::UpdateProjectionGeometryMeta(){
PanelOffsetIEC2[1] = 0.0;
PanelOffsetIEC2[2] = 0.0;
m_DRTImage2MetaInfo->SetPanelOffsetLPS(
m_CTMetaInfo->ConvertIECPointToLPSPoint(
PanelOffsetIEC2));
m_DRTImage2MetaInfo->SetOriginLPS(
CalcDRTImageOrigin(