plot residuals

[aguinot]

Implementation of the plots on the focal plane. We use the defined parameter n_ccd to split the stars per CCDs. This is hard-coded but the parameter has to be defined so it should be OK.

Features:

• Config file
  • parse the display of the focal plane (handle up-side-down CCDs)
  • parse all the plot we want to make (possibility to make operation on the columns)
  • plot parameters
    • nx
    • ny
• Handle "any" focal plane display (only square-ish at the moment)
• compute any residuals requested
• possibility to change the estimator to get the value in each bins (default: mean but could be len, median, ...)
  • Possibility to do it per plot
• Handle already existing run and skip computation (if requested) or only compute missing ones
  • Save information about the plots
  • Re-use already computed maps and only do the plotting
• plotting fonction
• Personnalisation of the plots
  • change cmap
  • force vmin, vmax

[aguinot]

At the moment the computation is done with numpy. I did not manage to get a good solution using dask. I will need to investigate further how dask works and how to correctly setup a dask Client and Cluster.

[aguinot]

Find a better name since we can do more than just PSF residuals. Maybe: focal plane plots?

[aguinot]

First implementation:

• 92M stars
• 40 CCDs
• (20, 40) binning
• 8 plots

Run time: 3min45 (on laptop)
🎉