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

Related to alredy closed Issue #97 - added update on ISO. Also taking out Panel Offset completely from Projection leads to some results. Makes sense. It is now commented on the code but will be kept in case of future analysis

Browse Source
This commit is contained in:
Proton local user
2023-06-07 11:37:17 +02:00
parent 005f42aa55
commit 4e2a7d87b1
1 changed files with 21 additions and 21 deletions
Download Patch File Download Diff File Expand all files Collapse all files

42
reg23Topograms/itkDTRrecon/vtkContourTopogramProjectionFilter.cxx
View File

@ -153,7 +153,7 @@ int vtkContourTopogramProjectionFilter::RequestData(
/** correct x coordinate by the ratio of distance from source
* Basic intercept theorem */
if(pp1[2] != 0.){
pp1[0] = (pp1[0] * dSCD / abs(pp1[2])) + dPOffset;
pp1[0] = (pp1[0] * dSCD / abs(pp1[2])) /*+ dPOffset*/;
} else {
vtkDebugMacro("point at SCD distance from panel. Skipping lateral correction for this point.");
}
@ -168,26 +168,26 @@ int vtkContourTopogramProjectionFilter::RequestData(
pp2[1] -= dProjectionLPSOffset [1];
pp2[2] -= dProjectionLPSOffset [2];
// Workaround to account for PA-LAT HFS-FFS
switch (this->m_ProjectionOrientation) {
case itk::PA:
pp2[0] += dPOffset;
break;
case itk::LAT:
switch (this->m_PatientOrientation) {
case itk::FFS:
pp2[1] += dPOffset;
break;
case itk::HFS:
pp2[1] -= dPOffset;
break;
default:
break;
}
break;
default:
break;
}
// // Workaround to account for PA-LAT HFS-FFS
// switch (this->m_ProjectionOrientation) {
// case itk::PA:
// pp2[0] += dPOffset;
// break;
// case itk::LAT:
// switch (this->m_PatientOrientation) {
// case itk::FFS:
// pp2[1] += dPOffset;
// break;
// case itk::HFS:
// pp2[1] -= dPOffset;
// break;
// default:
// break;
// }
// break;
// default:
// break;
// }
// std::cout << ii << " " << pp2[0] << " " << pp2[1] << " " << pp2[2] << std::endl;
PrjPoints->InsertPoint(ii,pp2);