Spring_2024 Updates

Lab_2/.ipynb_checkpoints/Lab 2-checkpoint.ipynb

@@ -0,0 +1,649 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "# <p style=\"text-align: center;\">PHYS 134L Spring 2022 Lab 2</p>"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "<div class=\"alert alert-block alert-danger\"><b>Due date:</b> Sunday, April 17th, 2022 by 11:59pm, submitted through Gradescope.</div>"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "## Names: "
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Enter your name and your partner's name here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Read through this entire lab before you start. This lab will introduce you to the idea of astronomical coordinate systems, sidereal time, focal length and image scale. To complete this lab you should have already read textbook Chapter 1 and Section 2.2. You may also find it useful to read the Wikipedia entries for \"Right ascension\", \"Declination,\" and \"Sidereal time.\""
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "\n",
+    "## <p style=\"text-align: center;\">Part 1: Celestial Coordinates</p>"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Objects in the sky have ''equatorial'' coordinates just as objects on the Earth. But instead of ''longitude'' for the\n",
+    "east-west coordinate, astronomers say ''right ascension'' (or ''RA''), and instead of ''latitude'' for the north-south coordinate, they say ''declination'' (or ''$\\delta$''). Right ascension is typically measured in hours, not degrees. Dividing the full circle of the celestial equator into 24 hours, each hour corresponds to 15 degrees.\n"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "We'll begin by revisiting the same datafile that we used in Lab 1: ```object.fits```.\n",
+    "**In ds9, open \\texttt{object.fits}.  Click on ''file/display\\_header'' to display the header.** The entries ''RA'' and ''Dec'' give the coordinates in the sky where the telescope was pointing while taking the image. Also notice the entry\n",
+    "commented ``start time of the observation.'' This is the time (in Universal Time, or UT) when the shutter opened for the CCD image. UT corresponds to the local time on the prime meridian (longitude = 0 degrees), which passes through the observatory at Greenwich, England. This time zone is 8 hours ahead of Pacific Std Time, and 7 hours ahead of Pacific Daylight Time.\n"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**On a certain day, the Sun crossed the meridian at 12:00 noon UT in Greenwich, England. Santa Barbara lies at 119.7 degrees W longitude. On that same day, what is the UT time when the Sun crosses the meridian in Santa Barbara?**"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**Using ds9 to get information from the headers, list the RA, $\\delta$, date, and UT for the file ```object.fits```**"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**On a sheet of paper, sketch a map of the part of the sky where the image was taken large enough to include a couple of bright, named stars (please label them).  Check out the [Stellarium website](https://stellarium-web.org/) if you need some help. Note that RA is defined so that as the Earth turns, the RA of objects on the meridian increases with time. Draw your map with N up and E to the left (as it would appear if you were facing the southern horizon). Attach the drawing to the end of this lab report before you submit it to gradescope**"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**Does RA increase to the right, or to the left?**"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Positions on the Earth stay fixed in the planet's rotating reference frame.  For example, as time passes, Los Angeles always remains about 100 miles, or 160 km, east-southeast of us here in Santa Barbara.  Positions in the sky are different.  The stars remain basically fixed in their current positions in the sky (they do move across the celestial sphere, but very slowly).  As the Earth rotates about its axis, the entire sky appears to rotate from the perspective of someone fixed on the Earth's surface.\n"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "\n",
+    "## <p style=\"text-align: center;\">Part 2: Sidereal Time</p>"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "The word ''sidereal'' means ''with respect to the stars.'' The current Local Sidereal Time (LST) is the value of the RA in the equatorial\n",
+    "coordinate system that is crossing your meridian at the moment. Since the coordinates of stars are essentially constant over very long times,\n",
+    "at a given LST you will always find the stars in the same apparent positions in the sky. Sidereal time is not the same as solar time (which we normally use) -- at a given solar time (such as noon), we find the {\\it Sun} in the same position, not the stars. Because the Earth orbits the Sun once per year, the sidereal day is about 4 minutes shorter than the solar day. Thus, measuring by solar time, a given star rises and sets about 4 minutes earlier every day."
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**Explain, with words, a sketch (appended to the end of the lab) and/or a calculation, why the solar day is about 4 minutes longer than the sidereal day.**"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answere here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**You are in Goleta (long = 120 degrees W, lat = +34:30) and the LST is 07:30. At the same instant, what is the LST in New York City (long = 74 degrees W, lat = +40:00)? (hint: remember 15 degrees per hour)**"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**What bright star (bright as in ''with a name, not a catalog number'') is most nearly overhead in Goleta at midnght tonight? (Check a star map on the web or the celestial globe.)**"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**What bright star is the most nearly overhead in Goleta over the whole year?**"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**You are in Urumqi, China (long = 87 degrees E, lat = +44:00) and the LST is 07:30. What bright star is most nearly tonight?**"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**What bright star is the most nearly overhead in Urumqi over the whole year?**"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Now we will explore the use of ```astropy``` to calculate local sidereal time and the apparent positions of objects in the sky.  Here we'll continue to use the file ```object.fits``` from the first lab. Recall from last time that (after importing the standard libraries, such as \\texttt{numpy, matplotlib, astropy, etc.}) you can read the image using\n",
+    "```\n",
+    "image = fits.open(imfile)[0].data\n",
+    "```\n",
+    "where ```imfile``` is the name of the file we want to open. As we discussed briefly in the introduction to fits files last lab, indexing the file by ```[0]``` opens the *first* file extension (this file only has one). "
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "To read the header in python, you simply use \n",
+    "```\n",
+    "header = fits.open(imfile)[0].header\n",
+    "```\n",
+    "\n",
+    "The header is read in as a python dictionary, which you now have experience with from the Python Bootcamp.\n",
+    "\n",
+    "**Open the file here and print out the contents of the header**"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "#Your Code here"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "One of the first listings should read\n",
+    "```\n",
+    "\\MJD-OBS =    56365.29978605324 / [UTC days] Start date/time (Modified Julian Dat\n",
+    "```\n",
+    "\n",
+    "There are three components to this.  The first, ```MJD-OBS```, is called the keyword.  It is how we access the information.  The second component, ```56365.29978605324``` is the value.  The third component is a comment.  You may create a new variable and set it equal to a header value using, e.g.,\n",
+    "```\n",
+    "mjd = header['MJD-OBS']}\n",
+    "```"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Now we'll try a bit of astropy's machinery to deal with coordinates and observations.  We'll use the astropy [coordinates](https://docs.astropy.org/en/stable/coordinates/index.html)coordinates, [time](https://docs.astropy.org/en/stable/time/index.html) and [units](https://docs.astropy.org/en/stable/units/index.html) packages. Check out their documentation pages for some examples on how to use them. **To use them, we first need to import the relevant bits, so put these lines in a notebook cell:**\n",
+    "```\n",
+    "from astropy import coordinates, time\n",
+    "from astropy import units as u\n",
+    "```"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Next, we want to create a coordinate object with the sky position of the observation, a location object with the location of our observation, and a time object.  We do that as follows (with ```imfile``` being the path to ````object.fits``` as you did last time):\n",
+    "```\n",
+    "header = fits.open(imfile)[0].header \n",
+    "radec = coordinates.SkyCoord(header['RA'], header['DEC'], unit=(u.hourangle, u.deg)) \n",
+    "location = coordinates.EarthLocation(lon=header['LONGITUD'], lat=header['LATITUDE']) \n",
+    "t = time.Time(header['MJD-OBS'], format='mjd')\n",
+    "```"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**Copy and paste the above code into the cell below to create a time object ```t```.**"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "#Your code here"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "An observer needs to know both their location and the time in order to compute the position of the stars.  These two objects together define the reference frame:\n",
+    "```\n",
+    "frame_altaz = coordinates.AltAz(obstime=t, location=location)\n",
+    "```\n",
+    "\n",
+    "MJD, the modified Julian date, is a unit of time, the (fractional) number of days since a certain date (midnight on November 17, 1858).  For example, Tuesday, Oct 5, 2021 at noon UTC is MJD 59492.50000.  Tuesday at 9:30 Pacific, (16:30 UTC; the start of class) is MJD 59492.68750, etc.  We are using an altitude-azimuth reference frame, which gives the cardinal direction (east-west-north-south) toward which I need to point my telescope, and the elevation or altitude, how high an angle above the horizon I should point.  With a location on the Earth, a time of observation, and a location on the sky, we can compute other quantities.  "
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**Now compute the local sidereal time at the time of observation.**  This calls a *method* of our time *object*.  We also need to specify our longitude on the Earth.\n",
+    "```\n",
+    "lst = t.sidereal\\_time('apparent', header['LONGITUD'])\n",
+    "```\n",
+    "**Compare this computed value to the local sidereal time that you can find in the fits header.**\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "#Your code here"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "\n",
+    "\n",
+    "**Now compute the altitude, or the elevation, of our telescope pointing.** We'll transform to the altitude-azimuth frame that we defined a few lines prior:\n",
+    "```\n",
+    "coord_altaz = radec.transform_to(frame_altaz)\n",
+    "altitude = coord_altaz.alt\n",
+    "```\n",
+    "\n",
+    "In that last line, ```coord_altaz.alt`` fetches the value ```alt`` held within the object ```coord_altaz```. \n",
+    "**Compare this altitude with the value you obtain from the header.**"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "#Your code here"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**Now is your turn to dig into the code a bit. Don't forget to read up on the documentation for how to use the astropy packages.**"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**What was the altitude one hour before this image was taken (hint: what was the mjd one hour before)?**"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "#Your code here"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*You answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**What will be the altitude two hours after the time this image was taken?**"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "#Your code here"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**Approximately when (local time)did the object pass through zenith (maximum elevation), and what was the local sidereal time when this happened?** (hint: recall the relationship between RA and local sidereal time)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "#Your code here"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "(Skip these next two questions and jump to the next section if you're short on time)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**Approximately how long after this image was taken did the object reach an elevation of 40 $^\\circ$?**"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {
+    "vscode": {
+     "languageId": "plaintext"
+    }
+   },
+   "outputs": [],
+   "source": [
+    "#Your code here"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**If you have time and are feeling ambitious, make a plot of the elevation vs. time from when the object rises to when it sets. (optional, not part of grade)**"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {
+    "vscode": {
+     "languageId": "plaintext"
+    }
+   },
+   "outputs": [],
+   "source": [
+    "#Your code here"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "## <p style=\"text-align: center;\">Part 3: Focal Length and Image Scale</p>"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "The size of an astronomical image on the CCD detector depends on the effective focal length (usually abbreviated ''focal length'') of\n",
+    "the telescope. Here is a link to a quick primary on focal length: [Focal length and f/# explained](https://www.paragon-press.com/lens/lenchart.html). This part of the lab will use a fits file called ```cluster.fits``` that should have been downloaded to your JupyterHub account when you clicked the link for this notebook, but it can also be found on the Lab 2 tab on the Gauchospace site. "
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**To get quantitative, an angle of 1 arcsec on the sky maps into a distance $d = fl / 206264.80$ measured in the focal plane of a\n",
+    "telescope having focal length $fl$.\n",
+    "Explain the origin of the magic number 206264.80. If you get stuck here, read the the primer above and think about the relationship between radians and arcseconds.**"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**The header of the image ```cluster.fits``` has a keyword ```FOCALLEN``` that gives an estimate of the telescope’s\n",
+    "focal length. The mean angular diameter of the Moon is about $\\frac{1}{2}$ degree, or 1800 arcsec. What is the diameter of an image of the Moon\n",
+    "made by this telescope, in mm?**"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {
+    "vscode": {
+     "languageId": "plaintext"
+    }
+   },
+   "outputs": [],
+   "source": [
+    "#Your code here"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**Other keywords give the widths and the number of detector pixels, in the x- and y- directions. What is the (x,y) size of the detector as\n",
+    "projected on the sky, in arcsec? Will the entire image of the full Moon fit on it?**"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "#Your code here"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "How do we know if the header value for the image scale of the telescope is correct? We will measure it. **Use ds9 to open the image ```cluster.fits```. Click on ''scale'' and ''zscale'' so you can see faint stars. Display the header, and write down the RA, $\\delta$ coordinates to which the telescope is pointing.**\n"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "*Your answer here*"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "After doing a Y-invert (see below), your image should look like this, except for the little circles and numbers, and the color scheme I used to make the numbers show up better.\n"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "![test](M67_circled.jpeg)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.7.3"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 2
+}