default plotting function now pops up figures

Author: davidssmith

src/plotting.jl

@@ -99,29 +99,34 @@ function makeplots(mat::Dict; outdir::String="results")
   title("arterial input function, \$C_p\$")
   savefig("$outdir/aif.pdf")
 
+  figure()
   clf()
   imshow(mat["SER"], interpolation="nearest", cmap="cubehelix", vmin=0, vmax=10)
   colorbar()
   title("signal enhancement ratio")
   savefig("$outdir/ser.pdf")
 
+  figure()
   clf()
   imshow(mask, interpolation="nearest", cmap="gray")
   title("mask")
   savefig("$outdir/mask.pdf")
 
+  figure()
   clf()
   imshow(R1map, interpolation="nearest", cmap="cubehelix", vmin=0, vmax=5)
   title("\$R_1\$ relaxation rate (s\$^{-1}\$)")
   colorbar(ticks=[0:5])
   savefig("$outdir/R1.pdf")
 
+  figure()
   clf()
   imshow(S0map, interpolation="nearest", cmap="gray")
   title("\$S_0\$")
   colorbar()
   savefig("$outdir/S0.pdf")
 
+  figure()
   clf()
   Ct = squeeze(maximum(Ct,1),1)
   x = oplot2(clamp(Ct, 0.0, 5.0), back, mask)
@@ -130,27 +135,31 @@ function makeplots(mat::Dict; outdir::String="results")
   colorbar()
   savefig("$outdir/Ct.pdf")
 
+  figure()
   clf()
   x = oplot2(clamp(Kt, 0.0, 1.0), back, mask)
   imshow(x, interpolation="nearest", cmap="jet", vmin=0, vmax=1)
   title("\$K^\\mathrm{trans}\$ (min\$^{-1}\$)")
   colorbar(ticks=[0:2:10]/10.0)
   savefig("$outdir/Kt.pdf")
 
+  figure()
   clf()
   x = oplot2(clamp(ve, 0.0, 1.0), back, mask)
   imshow(x, interpolation="nearest", cmap="jet", vmin=0, vmax=1)
   title("\$v_e\$")
   colorbar(ticks=[0,0.2,0.4,0.6,0.8,1])
   savefig("$outdir/ve.pdf")
 
+  figure()
   clf()
   x = oplot2(clamp(vp, 0.0, 1.0), back, mask)
   imshow(x, interpolation="nearest", cmap="jet", vmin=0, vmax=1)
   title("\$v_p\$")
   colorbar(ticks=[0,0.2,0.4,0.6,0.8,1])
   savefig("$outdir/vp.pdf")
 
+  figure()
   clf()
   kep = clamp(Kt./ve, 0.0, 10.0)
   x = oplot2(kep, back, mask)
@@ -159,12 +168,14 @@ function makeplots(mat::Dict; outdir::String="results")
   colorbar(ticks=[0,2,4,6,8,10])
   savefig("$outdir/kep.pdf")
 
+  figure()
   clf()
   imshow(modelmap, interpolation="nearest", cmap="cubehelix")
   title("model used")
   colorbar(ticks=[0,1,2,3])
   savefig("$outdir/modelmap.pdf")
 
+  figure()
   clf()
   x = oplot2(clamp(100*resid, 0.0, 1.0), back, mask)
   imshow(x, interpolation="nearest", cmap="jet", vmin=0, vmax=1.0)