--- a/docs/ipac/irsa/irsa.rst	2024-02-24 00:13:11.717826600 +0800
+++ b/docs/ipac/irsa/irsa.rst	2024-03-10 16:26:14.665398517 +0800
@@ -291,11 +291,10 @@
    from astropy.io import fits
    from astropy.nddata import Cutout2D
    from astropy.wcs import WCS
+   from astropy.table import Table
    import matplotlib.pyplot as plt
    coord = SkyCoord('150.01d 2.2d', frame='icrs')
-   spitzer_images = Irsa.query_sia(pos=(coord, 1 * u.arcmin), collection='spitzer_seip')
-   science_image = spitzer_images[spitzer_images['dataproduct_subtype'] == 'science'][0]
-   with fits.open(science_image['access_url'], use_fsspec=True) as hdul:
+   with fits.open('60095931.60095931-14.IRAC.4.mosaic.fits') as hdul:
         cutout = Cutout2D(hdul[0].section, position=coord, size=2 *
    u.arcmin, wcs=WCS(hdul[0].header))
    plt.imshow(cutout.data, cmap='grey')