Save intermediate results of mutli-resolution recon

[dyang37]

Per Rakesh's request, we should add an argument in recon() function to save the intermediate results of multi-resolution recon. Specifically, the interpolated results from lower resolution should be saved. In the case where the experiment gets interrupted, the interploated lower resolution image can be used as an init image to re-launch the experiment.

Proposed interface:

def recon(sino, angles, dist_source_detector, magnification,
          ...
          lib_path=__lib_path, save_init_images=0,
          ...
          ):

• lib_path (str, optional): [Default=~/.cache/mbircone] Path to directory containing library of forward projection matrices.

• save_init_images (int, optional): [Default=0] If (save_init_images>0) and (max_resolutions>0), the init_image for all resolution spaces except for the lowest one will be saved. Each init_image is generated from the interpolation of the recon result from the lower resolution space. Each init_image will be saved as a .npy file to /init_images/.npy, where < sysmatrix_name > is the name of the system matrix corresponding to the geometry parameters of the current resolution space.

[dyang37]

Implementation available at branch feature/save_intermediate.

I can make a PR once we get group agreement on the interface design.

[cabouman]

Diyu, I don't think it is a good idea to save ppt's or binaries in general on the GitHub. I'm looking at your email now.

[dyang37]

Diyu, I don't think it is a good idea to save ppt's or binaries in general on the GitHub. I'm looking at your email now.

Got it. Thanks! I deleted the ppt slides from the reply above.