70 likes | 248 Views
Usage of WATERMODULE A routine for estimating myocardial blood flow based on Iida model developed for MunichHeart. Chunlei Han Turku PET center Turku, Finland November, 2002 chunlei.han@tyks.fi. INPUTPUT. Time: one-dimension array with n-element
E N D
Usage of WATERMODULE A routine for estimating myocardial blood flow based on Iida model developed for MunichHeart Chunlei Han Turku PET center Turku, Finland November, 2002 chunlei.han@tyks.fi
INPUTPUT Time: one-dimension array with n-element LV: one-dimension array with n-element allTACs: two-dimension array (m*n) StudyType: string, “rest” or “stress”, to indicate study type method: calculation method, ‘regress’,’svdc’, ‘constrained_min’ ‘powell’, ‘pfdmin’ and ‘curvefit’. Default is ‘constrained_min’ cProtocol: string, “FreeAlfa” or “FixAlfa”, to indicate calculation with or without fixing alfa, default is “FreeAlfa” alfa_Rest_whole: scalar, rest study alfa value for average ROIs. alfa_Rest_regional: one-dimension array with n-element
Input argument structure Time 0.2 0.3 … .. 5.00 LV 1.2 10.3 … 334.4 ... ... 30.0 allROIs ROI1 ROI2 …… ROIm 0.34 …… ……. 20.5 … 490.5 ... ... 25.5 …… …… Alfa_Rest_regional 0.2 0.3 … 0.5 ... … 0.4 Alfa_Rest_whole=0.33
OUTPUT Output is a structure with fields of whole_Va: Va for average ROIs whole_f: f for average ROIs whole_alfa: alfa for average ROIs cInputfunction: calculated aterial input function regional.Va: Va for all ROIs regional.f: f for all ROIs regional.alfa: alfa for all ROIs regional.fittedcurve: fitted curves for all ROIs regional.fitquality: number 0-7 to indicate of fit goodness, only meaningful in “constrained_min”
Call Sequence General WATERMODULE(time, LV, allTACs, StudyType, method, cProtocol, alfa_Rest_whole, alfa_Rest_regional) Rest study WATERMODULE(time, LV, allTACs, StudyType method,) where StudyType=‘rest’ Stress study with free alfa WATERMODULE(time, LV, allTACs, StudyType, method, cProtocol) where StudyType=‘stress’ and cProtocol=‘FreeAlfa’ Stress study with fixed alfa WATERMODULE(time, LV, allTACs, StudyType, method, cProtocol, alfa_Rest_whole, alfa_Rest_regional) where StudyType=‘stress’ and cProtocol=‘FixAlfa’, also for alfa values from rest study
Default values and common block Default values default_p=0.91 default_beta=0.93 default_A=[0.5,1.0,0.5] default_xbnd=[[0,0,0],[1,50,1]] default_tolerance=1.0e-4 default_method=‘constrained_min’ default_Miter_constr=1000 default_Miter_curvefit=20 default_Miter_powell=200 Common block COMMON Xpara, Xin,Yin, FixAlfa_whole,FixAlfa_regional, constant constant={p:0.0,beta:0.0,xbnd_default:fltarr(3,2),A_initial:fltarr(3,1) tolerance:0.0,Miter_constr:0.0,Miter_curvefit:0.0,Miter_powell:0.0}
Function Call Flow watermodule Iida_whole Iida_regional StudyType=‘rest’ or StudyType=‘stress’ and cProtocol=‘FreeAlfa’ StudyType=‘stress’ and cProtocol=‘FixAlfa’ Func_whole_FixAlfa Func_whole_FixAlfa_curvefit Func_whole_FixAlfa_dfpminGRAD Func_whole_FreeAlfa Func_whole_FreeAlfa_curvefit Func_whole_FreeAlfa_dfpminGRAD Func_regional_FreeAlfa Func_regional_FreeAlfa_curvefit Func_regional_FreeAlfa_dfpminGRAD Func_regional_FixAlfa Func_regional_FixAlfa_curvefit Func_regional_FixAlfa_dfpminGRAD