Set Inf or NaN concentration to zero

davidssmith · Aug 23, 2020 · cd81e49 · cd81e49
1 parent 0846fac
commit cd81e49
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/src/science.jl b/src/science.jl
@@ -233,6 +233,7 @@ function fitdata(opts::Dict)
   end
   R1 = r1eff(dcedata, R10, TR, dceflip)
   Ct = tissueconc(R1, R10, relaxivity)
+  @. Ct[!isfinite(Ct)] = 0
   params, resid, modelmap = fitdce(Ct, mask, t, Cp, models=models)
 
   results = Dict()