From 32114111394f7796650041c00bdf8fce552f9b25 Mon Sep 17 00:00:00 2001 From: MelissaGraham Date: Mon, 10 Aug 2026 16:42:09 +0000 Subject: [PATCH] fix omission --- .../301_DP2_observations/301_1_DP2_overview.ipynb | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) diff --git a/DP2/300_Science_demos/301_DP2_observations/301_1_DP2_overview.ipynb b/DP2/300_Science_demos/301_DP2_observations/301_1_DP2_overview.ipynb index d91376ae..0fa4bbad 100644 --- a/DP2/300_Science_demos/301_DP2_observations/301_1_DP2_overview.ipynb +++ b/DP2/300_Science_demos/301_DP2_observations/301_1_DP2_overview.ipynb @@ -463,7 +463,8 @@ "source": [ "for i, name in enumerate(region_names):\n", " ra, dec, rad = regions[name]\n", - " offsets = np.sqrt((visit_table['ra']-ra)**2 + (visit_table['dec']-dec)**2)\n", + " cosdec = np.cos(np.deg2rad(dec))\n", + " offsets = np.sqrt((cosdec*(visit_table['ra']-ra))**2 + (visit_table['dec']-dec)**2)\n", " tx = np.where(offsets < rad)[0]\n", " for filt in filter_names:\n", " field_filt_dist[filt][i] = len(np.where(visit_table['band'][tx] == filt)[0])\n", @@ -519,7 +520,8 @@ "\n", "for i, name in enumerate(region_names):\n", " ra, dec, rad = regions[name]\n", - " offsets = np.sqrt((visit_table['ra']-ra)**2 + (visit_table['dec']-dec)**2)\n", + " cosdec = np.cos(np.deg2rad(dec))\n", + " offsets = np.sqrt((cosdec*(visit_table['ra']-ra))**2 + (visit_table['dec']-dec)**2)\n", " tx = np.where(offsets < rad)[0]\n", " field_epochs['visits'][i] = len(tx)\n", " unique_nights = np.unique(np.floor(visit_table['expMidptMJD'][tx]))\n", @@ -817,7 +819,8 @@ "\n", "for i, name in enumerate(region_names):\n", " ra, dec, rad = regions[name]\n", - " offsets = np.sqrt((joined_table['ra']-ra)**2 + (joined_table['dec']-dec)**2)\n", + " cosdec = np.cos(np.deg2rad(dec))\n", + " offsets = np.sqrt((cosdec*(joined_table['ra']-ra))**2 + (joined_table['dec']-dec)**2)\n", " tx = np.where(offsets < rad)[0]\n", " for filt in filter_names:\n", " fx = np.where(joined_table['band'][tx] == filt)[0]\n",