@@ -233,20 +233,6 @@ print(hdu_mer_irsa.info())
233233# hdu_mer_irsa.writeto(download_path+'./MER_image_VIS.fits', overwrite=True)
234234```
235235
236- ``` {code-cell} ipython3
237- ## Great! Now we have the data, lets extract just the primary image
238-
239- im_mer_irsa=hdu_mer_irsa[0].data
240- ```
241-
242- ``` {code-cell} ipython3
243- ## Plot a quick simple plot to show the full MER image, large FOV!
244-
245- plt.imshow(im_mer_irsa, cmap='gray', origin='lower',
246- norm=ImageNormalize(im_mer_irsa, interval=PercentileInterval(99.9), stretch=AsinhStretch()))
247- colorbar = plt.colorbar()
248- ```
249-
250236## 4. Overplot the catalog on the MER mosaic image
251237
252238``` {code-cell} ipython3
@@ -272,17 +258,6 @@ df_g_irsa['y_pix']=xy_irsa[1]
272258df_g_irsa
273259```
274260
275- ``` {code-cell} ipython3
276- ## Plot MER catalog sources on the Euclid VIS image
277-
278- plt.imshow(im_mer_irsa, cmap='gray', origin='lower', norm=ImageNormalize(im_mer_irsa, interval=PercentileInterval(99.9), stretch=LogStretch()))
279- colorbar = plt.colorbar()
280- plt.scatter(df_g_irsa["x_pix"], df_g_irsa["y_pix"], s=36, facecolors='none', edgecolors='red')
281-
282- plt.title('Galaxies between z = 1.4 and 1.6')
283- plt.show()
284- ```
285-
286261## Pull the spectra on the top brightest source based on object ID
287262
288263``` {code-cell} ipython3
@@ -373,7 +348,7 @@ new_hdu.header.update(cutout_data.wcs.to_header())
373348## Plot a quick simple plot to show the cutout on the galaxy
374349
375350plt.imshow(new_hdu.data, cmap='gray', origin='lower',
376- norm=ImageNormalize(im_mer_irsa , interval=PercentileInterval(99.9), stretch=AsinhStretch()))
351+ norm=ImageNormalize(new_hdu.data , interval=PercentileInterval(99.9), stretch=AsinhStretch()))
377352colorbar = plt.colorbar()
378353```
379354
0 commit comments