Update deps

Project.toml

@@ -35,15 +35,19 @@ FFTW = "1"
 FileIO = "1"
 FourierTools = "0.4.2"
 HDF5 = "0.17"
+ImageIO = "0.6"
 ImageShow = "0.3"
 IndexFunArrays = "0.2"
+JLD2 = "0.4"
 NDTools = "0.5"
 NNlib = "0.9"
 Optim = "1"
 Parameters = "0.12"
 Plots = "1"
+Pkg = "1"
 RadonKA = "0.6"
 Statistics = "1.10"
+Tar = "1"
 WaveOpticsPropagation = "0.2.0"
 Zygote = "0.6"
 

src/utils.jl

@@ -109,16 +109,16 @@ Plot an intensity histogram of the `printed_object`.
 `target` is the original binary target object.
 `thresholds` should indicate the threshold values for polymerization.
 """
-function plot_intensity_histogram(target, object_printed, thresholds; yscale=:log10)
+function plot_intensity_histogram(target, object_printed, thresholds; yscale=:log10, xlim=(0.0, 1.1))
     # :stephist vs :barhist
     target = Array(target)
     object_printed = Array(object_printed)
     plot_font = "Computer Modern"
     default(fontfamily=plot_font,
 	    linewidth=2, framestyle=:box, label=nothing, grid=false)
     m = maximum(object_printed)
-	plot(object_printed[target .== 0], seriestype=:stephist, bins=(0.0:0.01:m), xlim=(0.0, 1.1), label="dose distribution void", ylabel="voxel count", xlabel="normalized intensity",  ylim=(1, 10000000),  linewidth=1, legend=:topleft, yscale=yscale, size=(400, 300))
-	plot!(object_printed[target .== 1], seriestype=:stephist, bins=(0.0:0.01:m), xlim=(0.0, 1.1), label="dose distribution object", ylabel="voxel count", xlabel="normalized intensity",  ylim=(1, 10000000),  linewidth=1, legend=:topleft, yscale=yscale, size=(400, 300))
+	plot(object_printed[target .== 0], seriestype=:stephist, bins=(0.0:0.01:m), xlim=xlim, label="dose distribution void", ylabel="voxel count", xlabel="normalized intensity",  ylim=(1, 10000000),  linewidth=1, legend=:topleft, yscale=yscale, size=(400, 300))
+	plot!(object_printed[target .== 1], seriestype=:stephist, bins=(0.0:0.01:m),  label="dose distribution object", ylabel="voxel count", xlabel="normalized intensity",  ylim=(1, 10000000),  linewidth=1, legend=:topleft, yscale=yscale, size=(400, 300))
 	plot!([thresholds[1], thresholds[1]], [1, 10000_000], label="lower threshold", linewidth=1)
 	plot!([thresholds[2], thresholds[2]], [1, 10000_000], label="upper threshold", linewidth=1)
 	#plot!([chosen_threshold, chosen_threshold], [1, 30000000], label="chosen threshold", linewidth=3)