Use xarray plotting methods with NAIP

[snowman2]

I was going through: https://www.earthdatascience.org/courses/use-data-open-source-python/multispectral-remote-sensing/intro-naip/

And I was thinking the plotting examples may benefit from: http://xarray.pydata.org/en/stable/user-guide/plotting.html

If you have the long_name/units attributes, it pulls the labels from those IIRC.

Just a thought in case you hadn't seen it yet.

[lwasser]

@snowman2 thank you again for all of these issues! we will begin to work on some of these items as it will improve our lessons!! i will spend some more time digging into this in the upcoming weeks!

cc: @nkorinek

[snowman2]

Another place it may be useful: https://www.earthdatascience.org/courses/use-data-open-source-python/hierarchical-data-formats-hdf/open-MODIS-hdf4-files-python/