This fixes compatibility for imshow plots with array data that is a unit-aware ndarray subclass, for example data from yt.units or astropy.units.

Take for example the following script:

import numpy as np
import matplotlib.pyplot as plt
from yt.units import km

arr = np.random.random((400, 400))*km

plt.imshow(a)
plt.show()

This produces the following traceback right now: https://gist.github.com/ngoldbaum/27effb41173859132fed08a0ad5485ee

This is not isolated to yt's units - if one tries the same thing with astropy's units the same issue arises. Pint does not have this issue because Pint's Quantity class is not an ndarray subclass, so when it is converted to a masked array in matplotlib.image, the units are stripped.

This function seems to expect result.data to be an instance of the base ndarray class. This patch just enforces that expectation explicitly.