Simple script using Comrade.jl to generate an image model for ngEHT challenge 1 data.
Usage: julia -i imaging.jl
Sampling is expected to take on the order of several hours.
HMC samples are saved into the hchain
variable. Images are centered after sampling. Plots and displays the MAP, mean, std dev, and fractional error (std dev/mean) images.