import sunpy.map
import astropy.units as u
from sunpy.data.sample import AIA_171_IMAGE
amap = sunpy.map.Map(AIA_171_IMAGE)
amap.plot(clip_interval=(1.25, 99.5)*u.percent)