summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorErich Eckner <git@eckner.net>2019-10-24 15:50:53 +0200
committerErich Eckner <git@eckner.net>2019-10-24 15:50:53 +0200
commit3bc35c993d34c03c947db13d43af9f3eed79e02c (patch)
tree1b646c04621436a05d1e80d0df156d846f7b0101
parent513a1805cb8dbc883ee3bc245d1107162b68ca63 (diff)
downloadepost-3bc35c993d34c03c947db13d43af9f3eed79e02c.tar.xz
werteunit.inc: Werte sind nun zentriert bei Radontransformation - wir alternieren die Vorzeichen
-rw-r--r--werteunit.inc13
1 files changed, 7 insertions, 6 deletions
diff --git a/werteunit.inc b/werteunit.inc
index 21a3bd9..fe71749 100644
--- a/werteunit.inc
+++ b/werteunit.inc
@@ -190,16 +190,17 @@ begin
iYG:=floor(iY);
iYF:=iY-iYG;
// alternierendes Vorzeichen außen (oY) <=> Ergebnisse zentrieren
- werte[oX+oY*params.xSteps]:=(1-2*byte(odd(oY))) * (
+ // alternierendes Vorzeichen innen (iXG,iYG) <=> Eingangsdaten zentriert
+ werte[oX+oY*params.xSteps]:=(1-2*byte(odd(oY+iXG+iYG))) * (
qu^.reBei2DDoResSmi(iXG,iYG)*(1-iXF)*(1-iYF)
- + qu^.reBei2DDoResSmi(iXG+1,iYG)*iXF*(1-iYF)
- + qu^.reBei2DDoResSmi(iXG,iYG+1)*(1-iXF)*iYF
+ - qu^.reBei2DDoResSmi(iXG+1,iYG)*iXF*(1-iYF)
+ - qu^.reBei2DDoResSmi(iXG,iYG+1)*(1-iXF)*iYF
+ qu^.reBei2DDoResSmi(iXG+1,iYG+1)*iXF*iYF);
if (oY>0) and (2*oY<params.tSiz) then
- werte[oX+(params.tSiz-oY)*params.xSteps]:=(1-2*byte(odd(oY))) * (
+ werte[oX+(params.tSiz-oY)*params.xSteps]:=(1-2*byte(odd(oY+iXG+iYG))) * (
qu^.imBei2DDoResSmi(iXG,iYG)*(1-iXF)*(1-iYF)
- + qu^.imBei2DDoResSmi(iXG+1,iYG)*iXF*(1-iYF)
- + qu^.imBei2DDoResSmi(iXG,iYG+1)*(1-iXF)*iYF
+ - qu^.imBei2DDoResSmi(iXG+1,iYG)*iXF*(1-iYF)
+ - qu^.imBei2DDoResSmi(iXG,iYG+1)*(1-iXF)*iYF
+ qu^.imBei2DDoResSmi(iXG+1,iYG+1)*iXF*iYF);
end;
end;