@@ -53,7 +53,7 @@ First, we import all necessary packages.
5353
5454``` {code-cell} ipython3
5555# Uncomment the next line to install dependencies if needed.
56- # !pip install tqdm numpy matplotlib astropy photutils>=2.0 astroquery>=0.4.10 sep>=1.4 firefly_client>=3.2
56+ # !pip install tqdm numpy matplotlib astropy ' photutils>=2.0' ' astroquery>=0.4.10' ' sep>=1.4' ' firefly_client>=3.2'
5757```
5858
5959``` {code-cell} ipython3
@@ -159,7 +159,6 @@ Irsa.list_collections(servicetype='SIA')
159159
160160Here we use the collection * euclid_ero* , containing the Euclid ERO images. We first create a ` SkyCoord ` object and then query the SIA.
161161
162-
163162``` {code-cell} ipython3
164163image_tab = Irsa.query_sia(pos=(coord, search_radius), collection='euclid_ero')
165164print("Number of images available: {}".format(len(image_tab)))
@@ -319,7 +318,6 @@ mask = np.isnan(img)
319318
320319Next, we compute the background statistics what will be used by ` sep ` to extract the sources above a certain threshold.
321320
322-
323321``` {code-cell} ipython3
324322mean, median, std = sigma_clipped_stats(img, sigma=3.0)
325323print(np.array((mean, median, std)))
@@ -359,7 +357,6 @@ psfphot = PSFPhotometry(psf_model,
359357 finder = DAOStarFinder(fwhm=0.1, threshold=3.*std, exclude_border=True), # not needed because we are using fixed initial positions.
360358 aperture_radius = 4,
361359 progress_bar = True)
362-
363360```
364361
365362After initiating the ` PSFPhotometry ` object, we can run the PSF photometry measurement. This can take a while (typically between 1 and 2 minutes).
@@ -448,7 +445,6 @@ mask2 = np.isnan(img2)
448445
449446... and we also get some statistics on the sky background.
450447
451-
452448``` {code-cell} ipython3
453449mean2, median2, std2 = sigma_clipped_stats(img2, sigma=3.0)
454450print(np.array((mean2, median2, std2)))
@@ -521,7 +517,6 @@ Gaia.ROW_LIMIT = -1
521517
522518Next, we request the Gaia catalog around the position of the globular cluster. We use the same size as the cutout size.
523519
524-
525520``` {code-cell} ipython3
526521gaia_objects = Gaia.query_object_async(coordinate=coord, radius = cutout_size/2)
527522print("Number of Gaia stars found: {}".format(len(gaia_objects)))
0 commit comments