X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=tools%2FclitkConeBeamProjectImageFilter.txx;h=b4dba6d1f90c7fb2ce96f73f453cfcc00147af9f;hb=8ba069c89cddf44e9bd529b781bb5ddfd8d21adf;hp=173411f5531be7b8bd2057beaaecfb5a5a9d407b;hpb=a26cd8a19e1b9ad8344ab501436045f171a73713;p=clitk.git diff --git a/tools/clitkConeBeamProjectImageFilter.txx b/tools/clitkConeBeamProjectImageFilter.txx index 173411f..b4dba6d 100755 --- a/tools/clitkConeBeamProjectImageFilter.txx +++ b/tools/clitkConeBeamProjectImageFilter.txx @@ -114,7 +114,7 @@ namespace clitk itk::Matrix centeredRotationMatrix = GetCenteredRotationMatrix3D(rotationParameters,transformedCenter); // Compose this rotation with the rigid transform matrix - itk::Matrix finalTransform = centeredRotationMatrix * m_RigidTransformMatrix; + itk::Matrix finalTransform = m_RigidTransformMatrix * centeredRotationMatrix ; // Set the rotation itk::Matrix finalRotation = GetRotationalPartMatrix3D(finalTransform); @@ -169,7 +169,8 @@ namespace clitk // JV -1 seems to correspond better with shearwarp of Simon Rit typename InterpolatorType::InputPointType originOutput; originOutput[0] = m_IsoCenter[0]- (m_SourceToScreen - m_SourceToAxis); - originOutput[1] = m_IsoCenter[1]-static_cast(sizeOuput[1]-1)*spacingOutput[1]/2.0; + DD(m_PanelShift); + originOutput[1] = m_IsoCenter[1]-static_cast(sizeOuput[1]-1)*spacingOutput[1]/2.0 - m_PanelShift; originOutput[2] = m_IsoCenter[2]-static_cast(sizeOuput[2]-1)*spacingOutput[2]/2.0; m_Resampler->SetOutputOrigin( originOutput ); if (m_Verbose)std::cout<<"The origin of the flat panel is at "<< originOutput <<",..."<< std::endl;