Make use of Matplotlib's projection= framework

This makes the initialization of WCSAxes much easier:

```
import matplotlib.pyplot as plt
from wcsaxes import datasets, WCS

hdu = datasets.fetch_msx_hdu()
wcs = WCS(hdu.header)

fig = plt.figure()
ax = fig.add_subplot(1, 1, 1, projection=wcs)
```

@cdeil - could you review the changes to the docs since you raise this before?

cc @anizami 
Yes it will no longer be necessary once we merge wcsaxes into astropy.