Extracting unique colors of an overlay AxesImage in matplotlib

Question

I plot two arrays with the following code:

import numpy as np
import matplotlib.pyplot as plt

M1 = ([1,      2,      3,      np.nan], 
      [4,      5,