updated ARG simulations added

arg/arg_method/radiowater.fit

+FIT1        fit_winp 0.2.0  (c) 2010 by Turku PET Centre
+Date:       2010-09-19 19:19:34
+Data file:  mean.kbq
+Data unit:  kBq/mL
+Time unit:  sec
+Nr of VOIs: 1
+Region        Plane  Start   End     dataNr WSS       parNr Type Parameters
+000001 .      .        0.010 324.493    138 7.35E+003     6 0331 7.180211E-001 1.708856E+001 5.888741E+002 7.143209E-002 8.410205E+001 1.925347E-003

arg/arg_method/readme.txt

+Data for demonstration of ARG method for radiowater studies
+===========================================================
+
+See:
+simul.bat

arg/arg_method/simul.bat

+::
+::  Basic simulations on ARG method
+::
+
+
+@echo off
+setLocal EnableDelayedExpansion
+chcp 1252
+
+@echo.
+@echo Make result folders
+@echo.
+if not exist sim mkdir sim
+del /Q sim\*.* > NUL
+
+@echo.
+@echo Create input blood TAC
+fit2dat -c=0,360,1 radiowater.fit sim\radiowater_6-min.dat
+if %errorlevel% neq 0 goto FAILED
+fit2dat -c=0,90,1 radiowater.fit sim\radiowater_90-sec.dat
+if %errorlevel% neq 0 goto FAILED
+
+@echo Add delay
+tactime sim\radiowater_6-min.dat  +15 sim\radiowater_6-min_delayed.dat
+if %errorlevel% neq 0 goto FAILED
+tactime sim\radiowater_90-sec.dat +15 sim\radiowater_90-sec_delayed.dat
+if %errorlevel% neq 0 goto FAILED
+tacren sim\radiowater_6-min_delayed.dat 1 Blood__
+if %errorlevel% neq 0 goto FAILED
+tacren sim\radiowater_90-sec_delayed.dat 1 Blood__
+if %errorlevel% neq 0 goto FAILED
+
+@echo.
+@echo Simulate high perfusion TACs
+b2t_h2o -nosub -voiname=f=40 sim\radiowater_90-sec_delayed.dat  40 0.8 1 0 1 sim\high-040.dat
+if %errorlevel% neq 0 goto FAILED
+b2t_h2o -nosub -voiname=f=60 sim\radiowater_90-sec_delayed.dat  60 0.8 1 0 1 sim\high-060.dat
+if %errorlevel% neq 0 goto FAILED
+b2t_h2o -nosub -voiname=f=80 sim\radiowater_90-sec_delayed.dat  80 0.8 1 0 1 sim\high-080.dat
+if %errorlevel% neq 0 goto FAILED
+b2t_h2o -nosub -voiname=f=100 sim\radiowater_90-sec_delayed.dat 100 0.8 1 0 1 sim\high-100.dat
+if %errorlevel% neq 0 goto FAILED
+b2t_h2o -nosub -voiname=f=120 sim\radiowater_90-sec_delayed.dat 120 0.8 1 0 1 sim\high-120.dat
+if %errorlevel% neq 0 goto FAILED
+
+@echo.
+@echo Simulate low perfusion TACs
+b2t_h2o -nosub -voiname=f=2 sim\radiowater_6-min_delayed.dat  2 0.8 1 0 1 sim\low-02.dat
+if %errorlevel% neq 0 goto FAILED
+b2t_h2o -nosub -voiname=f=4 sim\radiowater_6-min_delayed.dat  4 0.8 1 0 1 sim\low-04.dat
+if %errorlevel% neq 0 goto FAILED
+b2t_h2o -nosub -voiname=f=6 sim\radiowater_6-min_delayed.dat  6 0.8 1 0 1 sim\low-06.dat
+if %errorlevel% neq 0 goto FAILED
+b2t_h2o -nosub -voiname=f=8 sim\radiowater_6-min_delayed.dat  8 0.8 1 0 1 sim\low-08.dat
+if %errorlevel% neq 0 goto FAILED
+b2t_h2o -nosub -voiname=f=10 sim\radiowater_6-min_delayed.dat 10 0.8 1 0 1 sim\low-10.dat
+if %errorlevel% neq 0 goto FAILED
+
+@echo.
+@echo Plot TACs
+tac2svg -legend=y -xt="Time p.i. (s)" -yt="Concentration (kBq/cc)" -x1=0 -x2=375 -y1=0 sim\radiowater-arg-low.svg -l sim\radiowater_6-min_delayed.dat sim\low-??.dat
+if %errorlevel% neq 0 goto FAILED
+tac2svg -legend=y -xt="Time p.i. (s)" -yt="Concentration (kBq/cc)" -x1=0 -x2=105 -y1=0 sim\radiowater-arg-high.svg -l sim\radiowater_90-sec_delayed.dat sim\high-???.dat
+if %errorlevel% neq 0 goto FAILED
+
+@echo.
+@echo Calculate ARG look-up tables
+arlkup -static=n -nr=500 sim\radiowater_90-sec_delayed.dat 0.8 120 15 90 sim\high.lkup
+if %errorlevel% neq 0 goto FAILED
+arlkup -static=n -nr=500 sim\radiowater_6-min_delayed.dat 0.8 10 15 360 sim\low.lkup
+if %errorlevel% neq 0 goto FAILED
+
+@echo.
+@echo Plot lookup tables
+tac2svg -mt="" -legend=n -xt="AUC 0-360 s" -yt="Blood flow (mL/min/dL)" -x1=0 -y1=0 sim\radiowater-arg-low-lkup.svg -l sim\low.lkup
+if %errorlevel% neq 0 goto FAILED
+tac2svg -mt="" -legend=n -xt="AUC 0-90 s" -yt="Blood flow (mL/min/dL)" -x1=0 -y1=0 sim\radiowater-arg-high-lkup.svg -l sim\high.lkup
+if %errorlevel% neq 0 goto FAILED
+
+
+@echo.
+@echo AUC versus perfusion
+
+@echo Low flow
+
+echo 1 2 > sim\x.dat
+echo 2 4 >> sim\x.dat
+echo 3 6 >> sim\x.dat
+echo 4 8 >> sim\x.dat
+echo 5 10 >> sim\x.dat
+
+interpol -i -x=375 -header=n sim\low-02.dat  > sim\auc-low.dat
+if %errorlevel% neq 0 goto FAILED
+interpol -i -x=375 -header=n sim\low-04.dat >> sim\auc-low.dat
+if %errorlevel% neq 0 goto FAILED
+interpol -i -x=375 -header=n sim\low-06.dat >> sim\auc-low.dat
+if %errorlevel% neq 0 goto FAILED
+interpol -i -x=375 -header=n sim\low-08.dat >> sim\auc-low.dat
+if %errorlevel% neq 0 goto FAILED
+interpol -i -x=375 -header=n sim\low-10.dat >> sim\auc-low.dat
+if %errorlevel% neq 0 goto FAILED
+
+tacsetx sim\auc-low.dat sim\x.dat 1
+if %errorlevel% neq 0 goto FAILED
+
+tac2svg -legend=n -mt="" -yt="AUC 0-360 s" -xt="Blood flow (mL/min/dL)" -x1=0 -y1=0 sim\radiowater-auc-vs-flow-low.svg  sim\auc-low.dat
+if %errorlevel% neq 0 goto FAILED
+
+
+@echo High flow
+
+echo 1 40 > sim\x2.dat
+echo 2 60 >> sim\x2.dat
+echo 3 80 >> sim\x2.dat
+echo 4 100 >> sim\x2.dat
+echo 5 120 >> sim\x2.dat
+
+interpol -i -x=105 -header=n sim\high-040.dat  > sim\auc-high.dat
+if %errorlevel% neq 0 goto FAILED
+interpol -i -x=105 -header=n sim\high-060.dat >> sim\auc-high.dat
+if %errorlevel% neq 0 goto FAILED
+interpol -i -x=105 -header=n sim\high-080.dat >> sim\auc-high.dat
+if %errorlevel% neq 0 goto FAILED
+interpol -i -x=105 -header=n sim\high-100.dat >> sim\auc-high.dat
+if %errorlevel% neq 0 goto FAILED
+interpol -i -x=105 -header=n sim\high-120.dat >> sim\auc-high.dat
+if %errorlevel% neq 0 goto FAILED
+
+tacsetx sim\auc-high.dat sim\x2.dat 1
+if %errorlevel% neq 0 goto FAILED
+
+tac2svg -legend=n -mt="" -yt="AUC 0-90 s" -xt="Blood flow (mL/min/dL)" -x1=0 -y1=0 sim\radiowater-auc-vs-flow-high.svg  sim\auc-high.dat
+if %errorlevel% neq 0 goto FAILED
+
+
+
+
+:PASSED
+@echo.
+@echo ==========================================================================
+@echo All steps fine.
+@echo ==========================================================================
+goto END
+:FAILED
+@echo.
+@echo ==========================================================================
+@echo Failed!
+@echo ==========================================================================
+pause
+:END
+@echo.

arg/arg_vs_va/Va.dat

+1 0.0
+2 0.5
+3 1.0
+4 1.5
+5 2.0
+6 3.0
+7 4.0
+8 5.0
+9 6.0
+10 7.0
+11 8.0

arg/arg_vs_va/radiowater.fit

+FIT1        fit_winp 0.2.0  (c) 2010 by Turku PET Centre
+Date:       2010-09-19 19:19:34
+Data file:  mean.kbq
+Data unit:  kBq/mL
+Time unit:  sec
+Nr of VOIs: 1
+Region        Plane  Start   End     dataNr WSS       parNr Type Parameters
+000001 .      .        0.010 324.493    138 7.35E+003     6 0331 7.180211E-001 1.708856E+001 5.888741E+002 7.143209E-002 8.410205E+001 1.925347E-003

arg/arg_vs_va/readme.txt

+Simulating the bias in ARG method caused by of ignored Va
+=========================================================
+
+See:
+simul.bat

arg/auc_vs_f/radiowater.fit

+FIT1        fit_winp 0.2.0  (c) 2010 by Turku PET Centre
+Date:       2010-09-19 19:19:34
+Data file:  mean.kbq
+Data unit:  kBq/mL
+Time unit:  sec
+Nr of VOIs: 1
+Region        Plane  Start   End     dataNr WSS       parNr Type Parameters
+000001 .      .        0.010 324.493    138 7.35E+003     6 0331 7.180211E-001 1.708856E+001 5.888741E+002 7.143209E-002 8.410205E+001 1.925347E-003

arg/auc_vs_f/readme.txt

+Demonstration of AUC versus perfusion in radiowater studies
+===========================================================
+
+See:
+simul.bat