[1]:
import astropy.units as u
-import matplotlib.pyplot as plt
-import mocpy
-from astropy.coordinates import Angle, SkyCoord
-from astropy.visualization.wcsaxes.frame import EllipticalFrame
-from astropy.wcs import WCS
-from astroquery.mocserver import MOCServer
-from astroquery.vizier import Vizier
-from ipyaladin import Aladin
-from matplotlib import patheffects
-from regions import CircleSkyRegion
+import astropy.units as u
+import matplotlib.pyplot as plt
+import mocpy
+from astropy.coordinates import Angle, SkyCoord
+from astropy.visualization.wcsaxes.frame import EllipticalFrame
+from astropy.wcs import WCS
+from astroquery.mocserver import MOCServer
+from astroquery.vizier import Vizier
+from ipyaladin import Aladin
+from matplotlib import patheffects
+from regions import CircleSkyRegion
-mocpy.__version__
+mocpy.__version__
[2]:
# HST MOC footprint
-# -----------------
-# We want to retrieve all the HST surveys i.e. the HST surveys covering any region of the sky.
-allsky = CircleSkyRegion(SkyCoord(0, 0, unit="deg"), Angle(180, unit="deg"))
-hst_moc = MOCServer.query_region(
- region=allsky,
- # We want a MOCpy object instead of an astropy table
- return_moc=True,
- # The order of the MOC
- max_norder=11,
- # Expression on the ID meta-data
- meta_data="ID=*HST*",
-)
+# HST MOC footprint
+# -----------------
+# We want to retrieve all the HST surveys i.e. the HST surveys covering any region of the sky.
+allsky = CircleSkyRegion(SkyCoord(0, 0, unit="deg"), Angle(180, unit="deg"))
+hst_moc = MOCServer.query_region(
+ region=allsky,
+ # We want a MOCpy object instead of an astropy table
+ return_moc=True,
+ # The order of the MOC
+ max_norder=11,
+ # Expression on the ID meta-data
+ meta_data="ID=*HST*",
+)
-# SDSS9
-# -----
-sdss_moc = MOCServer.find_datasets(
- meta_data="ID=CDS/P/SDSS9/color",
- return_moc=True,
- max_norder=11,
-)
+# SDSS9
+# -----
+sdss_moc = MOCServer.find_datasets(
+ meta_data="ID=CDS/P/SDSS9/color",
+ return_moc=True,
+ max_norder=11,
+)
[3]:
type(sdss_moc)
+type(sdss_moc)
[4]:
%%time
-sdss_and_hst_moc = sdss_moc.intersection(hst_moc)
-sdss_moc.serialize(format="fits")
-sdss_moc.save("sdss_moc.fits", format="fits", overwrite=True)
+%%time
+sdss_and_hst_moc = sdss_moc.intersection(hst_moc)
+sdss_moc.serialize(format="fits")
+sdss_moc.save("sdss_moc.fits", format="fits", overwrite=True)
[5]:
fig = plt.figure(figsize=(20, 10))
+fig = plt.figure(figsize=(20, 10))
-wcs = WCS(
- {
- "naxis": 2,
- "naxis1": 1620,
- "naxis2": 810,
- "crpix1": 810.5,
- "crpix2": 405.5,
- "cdelt1": -0.2,
- "cdelt2": 0.2,
- "ctype1": "RA---AIT",
- "ctype2": "DEC--AIT",
- },
-)
+wcs = WCS(
+ {
+ "naxis": 2,
+ "naxis1": 1620,
+ "naxis2": 810,
+ "crpix1": 810.5,
+ "crpix2": 405.5,
+ "cdelt1": -0.2,
+ "cdelt2": 0.2,
+ "ctype1": "RA---AIT",
+ "ctype2": "DEC--AIT",
+ },
+)
-# Create a matplotlib axe and give it a astropy.wcs.WCS object
-ax = fig.add_subplot(1, 1, 1, projection=wcs, frame_class=EllipticalFrame)
+# Create a matplotlib axe and give it a astropy.wcs.WCS object
+ax = fig.add_subplot(1, 1, 1, projection=wcs, frame_class=EllipticalFrame)
-# Fill the SDSS MOC in red with an opacity of 70%
-sdss_moc.fill(
- ax=ax,
- wcs=wcs,
- edgecolor="r",
- facecolor="r",
- linewidth=0,
- fill=True,
- alpha=0.7,
- label="sdss9 footprint",
-)
-# Draw its perimeter in black
-sdss_moc.border(ax=ax, wcs=wcs, color="black", alpha=0.5)
+# Fill the SDSS MOC in red with an opacity of 70%
+sdss_moc.fill(
+ ax=ax,
+ wcs=wcs,
+ edgecolor="r",
+ facecolor="r",
+ linewidth=0,
+ fill=True,
+ alpha=0.7,
+ label="sdss9 footprint",
+)
+# Draw its perimeter in black
+sdss_moc.border(ax=ax, wcs=wcs, color="black", alpha=0.5)
-# Fill the HST surveys MOC in green with an opacity of 70%
-hst_moc.fill(
- ax=ax,
- wcs=wcs,
- edgecolor="g",
- facecolor="g",
- linewidth=0,
- fill=True,
- alpha=0.7,
- label="hst all surveys footprint",
-)
-# Draw its perimeter in black
-hst_moc.border(ax=ax, wcs=wcs, color="black", alpha=0.5)
+# Fill the HST surveys MOC in green with an opacity of 70%
+hst_moc.fill(
+ ax=ax,
+ wcs=wcs,
+ edgecolor="g",
+ facecolor="g",
+ linewidth=0,
+ fill=True,
+ alpha=0.7,
+ label="hst all surveys footprint",
+)
+# Draw its perimeter in black
+hst_moc.border(ax=ax, wcs=wcs, color="black", alpha=0.5)
-# Fill the intersection MOC in blue
-sdss_and_hst_moc.fill(
- ax=ax,
- wcs=wcs,
- edgecolor="b",
- facecolor="b",
- linewidth=0,
- fill=True,
- alpha=0.7,
- label="intersection",
-)
-# Draw its perimeter in black
-sdss_and_hst_moc.border(ax=ax, wcs=wcs, color="black", alpha=0.5)
+# Fill the intersection MOC in blue
+sdss_and_hst_moc.fill(
+ ax=ax,
+ wcs=wcs,
+ edgecolor="b",
+ facecolor="b",
+ linewidth=0,
+ fill=True,
+ alpha=0.7,
+ label="intersection",
+)
+# Draw its perimeter in black
+sdss_and_hst_moc.border(ax=ax, wcs=wcs, color="black", alpha=0.5)
-# we set the axis off by half a pixel otherwise they are slightly outside the sphere
-# see astropy issue for the half-pixel trick https://github.com/astropy/astropy/issues/10201
-ax.set_xlim(-0.5, 1620 - 0.5)
-ax.set_ylim(-0.5, 810 - 0.5)
-ax.set_aspect(1.0)
+# we set the axis off by half a pixel otherwise they are slightly outside the sphere
+# see astropy issue for the half-pixel trick https://github.com/astropy/astropy/issues/10201
+ax.set_xlim(-0.5, 1620 - 0.5)
+ax.set_ylim(-0.5, 810 - 0.5)
+ax.set_aspect(1.0)
-# Usual matplotlib calls
-plt.title("Using matplotlib to vizualize MOCs")
-plt.xlabel("ra")
-plt.ylabel("dec")
-plt.legend()
-plt.grid(color="black", linestyle="dotted")
-path_effects = [patheffects.withStroke(linewidth=3, foreground="black")]
-ax.coords[0].set_ticklabel(color="white", path_effects=path_effects)
-plt.show()
-plt.close()
+# Usual matplotlib calls
+plt.title("Using matplotlib to vizualize MOCs")
+plt.xlabel("ra")
+plt.ylabel("dec")
+plt.legend()
+plt.grid(color="black", linestyle="dotted")
+path_effects = [patheffects.withStroke(linewidth=3, foreground="black")]
+ax.coords[0].set_ticklabel(color="white", path_effects=path_effects)
+plt.show()
+plt.close()
[ ]:
viz = Vizier(columns=["*", "_RAJ2000", "_DEJ2000"])
-viz.ROW_LIMIT = -1
-# Photometric standard stars (tables II and IV of paper)
-tables = viz.get_catalogs("II/50")
-tables
+viz = Vizier(columns=["*", "_RAJ2000", "_DEJ2000"])
+viz.ROW_LIMIT = -1
+# Photometric standard stars (tables II and IV of paper)
+tables = viz.get_catalogs("II/50")
+tables
[ ]:
table_II50 = tables[0]
-table_II50
+table_II50 = tables[0]
+table_II50
[ ]:
idx_inside = sdss_and_hst_moc.contains_lonlat(
- table_II50["_RAJ2000"].T * u.deg,
- table_II50["_DEJ2000"].T * u.deg,
-)
-sources_inside = table_II50[idx_inside]
-sources_outside = table_II50[~idx_inside]
-sources_inside
+idx_inside = sdss_and_hst_moc.contains_lonlat(
+ table_II50["_RAJ2000"].T * u.deg,
+ table_II50["_DEJ2000"].T * u.deg,
+)
+sources_inside = table_II50[idx_inside]
+sources_outside = table_II50[~idx_inside]
+sources_inside
[ ]:
aladin = Aladin()
-aladin
+aladin = Aladin()
+aladin
[ ]:
aladin.target = "messier 51"
-aladin.fov = 1
-aladin.height = 600
+aladin.target = "messier 51"
+aladin.fov = 1
+aladin.height = 600
[ ]:
aladin.coo_frame = "icrs"
+aladin.coo_frame = "icrs"
[ ]:
aladin.survey = "P/SDSS9/color"
+aladin.survey = "P/SDSS9/color"
[ ]:
aladin.add_moc(
- sdss_and_hst_moc,
- color="blue",
- opacity=0.3,
- fill=True,
- edge=False,
- perimeter=True,
-)
+aladin.add_moc(
+ sdss_and_hst_moc,
+ color="blue",
+ opacity=0.3,
+ fill=True,
+ edge=False,
+ perimeter=True,
+)
[ ]:
# Add astropy source tables to the aladin lite viewer
-aladin.add_table(sources_inside, color="green", shape="rhomb")
-aladin.add_table(sources_outside, color="red", shape="cross")
+# Add astropy source tables to the aladin lite viewer
+aladin.add_table(sources_inside, color="green", shape="rhomb")
+aladin.add_table(sources_outside, color="red", shape="cross")
[ ]:
# change the fov and target
-aladin.target = "13 04 4.193 -03 34 13.54"
-aladin.fov = 11
+# change the fov and target
+aladin.target = "13 04 4.193 -03 34 13.54"
+aladin.fov = 11
© Copyright 2024, Matthieu Baumann <matthieu.baumann@astro.unistra.fr> Thomas Boch <thomas.boch@astro.unistra.fr> Manon Marchand <manon.marchand@astro.unistra.fr> F.-X. Pineau <francois-xavier.pineau@astro.unistra.fr>.
- Last built 26 Jul 2024.
+ Last built 17 Oct 2024.