Wednesday, April 2, 2014

Consecutive Evenings

For the first time in this project, I have data from two consecutive evenings!

Plus, it looks like maybe tonight is also going to be clear and calm.  Three nights in a row?  Maybe so!

Here's a quick couple images from last night showing the relative motion of Io and Europa, as well as motion of the Jupiter system itself against the background stars:

Figure 1: The Jupiter system 02 April 2014

It's gonna be sorta neat to watch the photometry of Io and Europa over this evening.  At the start, the photometry will show one object since they're so close.  So the values will be some combination of the two.  As the night goes on, however, Io and Europa separate enough for the software to distinguish two objects.  So the photometry will show a "single object" at a certain brightness, and then that line will split into two separate brightnesses, presumably slightly fainter than the combined brightness, but at pretty much the same distance and position angle.

1800 images were taken last night.  About 1282 have data in them.  962 images have Jupiter in them.  So as long as I'm hitting about 50% with the "Jupiter in the image" stat, I'm ok with my performance as a substitute telescope drive.

Tuesday, April 1, 2014

Clear night and April Events

A very decent session last night tobegin the new month and the new lunar cycle.  The moon was a thin, 30 hour sliver and was nice to see so young.

I'll post more later as I have a number of new things to talk about.

But for now, here's an image from last night (notice I've reversed the greyscale, so black = bright):

Figure 1: The Jupiter system 04:00:52 UTC 01 April 2014
I've numbered five background stars in this image so you can compare it to the AAVSO sky chart of the same area:

Figure 2: AAVSO chart of same region


The crosshair in Figure 2 correspond to RA 06:49:14.5 and Dec +23:12:29.  I mention this because at least according to the ephemeris on the Heaven's Above website, Jupiter should be at RA 06:50:14.5 Dec +23:12:29.  As you can see from the image, Jupiter is not at the location specified by Heaven's Above.

This site seems to do a better job at calculating the correct position.  Here's the same chart showing the calculated position of Jupiter:


Figure 3: A better position of Jupiter marked

So, Heaven's Above -- you're fired!  This makes me question all the other information on this site now.

I'll make another post later today or tomorrow as I have a number of things to talk about.

==== Mid-day Update ====

Wow this is a really great site.  A lot of very useful info for me can be gotten from this site, including positions of the Jovian satellites and how far away they are from Jupiter.  Plus, very cool charts like this one:



April Jupiter moon eclipses that I'll be able to record:

DD  HMM.M
 6  414.0  II EC R
 8  343.8   I EC R
11  329.4 III EC R
15  538.8   I EC R
18  4 7.5 III EC D

There are also some near conjunctions that have become pretty interesting to record.

MMDDYY    Moons   HH:MM
040914    I/E     05:00
042214    I/E/C   03:00
050214    I/G     03:00
050414    I/E     03:00
050514    E/G     03:00
050614    I/G     05:00

Thursday, March 27, 2014

Systematic Errors Yes and No

Does my data contain systematic errors?  You bet, it does.

Do those systematic errors make the data useless?  Not in the least.

The systematic errors are contributing to the overall error, but when the overall error is 2 - 4 percent, I'm not going to worry about them too much.

I spent all of yesterday looking at these errors and, while they're annoying, they appear to be in the raw data which means nothing I can do on the software end is going to change that.  I'm not going to try to explain to my readers what I'm seeing, because unless you've been looking at all of my data as much as I have, you probably won't see it.  However, this is the kind of plot I was looking at all day:

Figure 1: Photometry data from 23 Mar 2014 UTC

Now I guess that might look a bit scattered.  Staring at it a while will allow you to start seeing some patterns.  See them?  But consider the scale on the y-axis, plus the overall scale of the data (zero to 65535).  When I do this, I realize that the data in Figure 1 has an error of 2.18%.  This is well within my tolerances as far as having confidence that I'll be able to see variations in the lightcurves as these moons orbit Jupiter.

So now I can say, for the first time in this project, that I'm confident about the data that I have.  It seems consistent to a very reasonable tolerance.

I now, therefore, move on to more important matters.

Photometry versus Phase Angle

What I'm slowly moving towards is being able to plot phase angle versus photometric value.  When I do this, I should see variations similar to those shown in those USGS plots I did on a previous post.

To calculate the phase angle, I need to know three pieces of information:

  • The distance from the target to Jupiter
  • The direction of motion of the target (toward or away from Jupiter)
  • The side that the target is on (east or west of Jupiter)

I've created a little parameter file that will contain data used for this calculation (as well as the data needed to calculate the airmass of the observations -- needed to correct for atmospheric extinction).  In this case, the data is the distance (in AU) to Jupiter.  I need this in order to properly calculate the radius of the orbits for each moon in pixels.

The phase angle is then defined:

cos (phase angle) = measured distance / radius

I also need to know what "quadrant" the moon is in, and the direction of motion and the side tells me that:

Quad 1: east moving away
Quad 2: east moving toward
Quad 3: west moving away
Quad 4: west moving toward

So with all of this info, I can calculate the phase angle, which of course will go from zero to 360 degrees.  I've decided to define zero degrees as the point in the orbit when the moon is directly behind Jupiter and moving east.  This is an arbitrary definition and I'm not sure if it's the standard definition.  Once I find out what that is, it's likely I'll have to change mine to fit the accepted geometry.  No biggie.  The parameter file I mentioned above will also contain this information.

"Lost" Data

The other problem I'm trying to figure out is that a lot of data is "lost" because Jupiter is out of the field and therefore I can't compute a distance.  I've got an idea about how to deal with this and I'm hoping to try that out sometime very soon.  Callisto will especially suffer from this since it almost always is going to appear pretty far from Jupiter.  But as I said, I've got a possible solution and I'll see if that helps in any way to at least give me a way of looking at all of the data even though it may not be as exact as I would like it.

Data Reduction

I still have a backlog of new data to reduce.  I also need to go back and re-reduce all of my old data so I can start looking at this phase angle stuff.  Once I get the phase angle code running, I look forward to plotting all of this data to start looking for photometric variations.

Tuesday, March 25, 2014

Asteroid Data Hunter, Backlog Jupiter Data, Relative Photometry, Things To Do

Asteroid Data Hunter

I'm still sorta recovering from my marathon session creating my entry for the Asteroid Data Hunter competition.  I should know the results of the reviews later today.  I'm having very strong doubts that my proposal is "good enough", but that decision is out of my hands.  I look forward to seeing all the other entries.

This is a pretty big deal if I win, but for now I'm keeping my excitement to a minimum.

Backlog Jupiter Data

Wow -- I have a backlog of data that I still need to reduce and get prelim analysis.  Never thought I'd be saying such a thing for my own data.  It's just plain old work to reduce data.  Just gotta do it, that's all.

I've modified my data collection technique in that now I'm taking 300 images per "set" rather than 200.  While 846/1800 (48%) images actually have Jupiter in them, because of the geometry of the moons (Callisto being pretty far away), I'd say that something like 85% of the images have SOMETHING in them of value.  So that's something around 1500 new data points.

I'm pretty excited about these last couple of data sets since they contain my first known observations of eclipses.  The one involving Callisto ingress into eclipse on 20 Mar 2014 caught me totally by surprise.  The session from 23 Mar 2014 has an Io egress out of eclipse.   I look forward to plotting each of those light curves to see when and how these events occurred.  I need to go back though all of my data now to see if I've already recorded one of these events without knowing it.

Here are pics of the Io egress event:

Figure 1: Io egress event

The times for these images are:

Upper left: 05:22:04 UTC
Upper right: 05:22:27 UTC
Lower left: 05:24:02 UTC

I'd say my timing error is probably +/- 2 seconds but I'll have to compare the internal clock with the WWV to be certain.

Using the predictions from this site, I note that the time of this event is 05:24:54.  So I see that my timing is off by about three minutes (three minutes early).  Not sure why.  I need to check my internal clock on the computer but I doubt it's drifted that much in just a month.  I also don't know what this 05:24:54 time is -- if it's for the start of the egress, the end of the egress, or the middle.  All I can say at this point is that it appeared fairly suddenly over the course of about one or two minutes.  If my timings are correct, then I'd say that the time from IMCEE is the time of the END of the event (when Io reaches its nominal magnitude and is totally out of the shadow of Jupiter).

I guess there's nothing else I should really do until I get the atmospheric extinction stuff built into the code.  This requires me writing some software that'll do least-squares fitting to a polynomial of an arbitrary order.  I want to write a general fitting piece of software like this so I don't have to come back later an expand on just a simple linear fit.  This requires dusting off some old books and reminding myself how to compute things like determinants and cofactors.  Ah, mathematics.  So really, until I have some working code to do this kind of thin, there's no point in running through the data.  The data will be there and continue to collect it.  Tonight looks like it's gonna be a good night with Io at max elongation and Callisto and Ganymede in nice positions.  Europa is gonna be heading closer to Jupiter, so I'm likely not going to get anything useful outta her.

Relative Photometry

I was staring at a background star the other night in the same FOV as the Jupiter system and remarked to myself that I could be doing relative photometry of each moon with that star.  I can assume for the time being (until I find out what star this is) that the star maintained a constant brightness over the period of the session.  This means I can compute the ratio (of the difference) between the background star and the moons to see how they vary in brightness over time.  This kind of measurement mostly eliminates the need for extinction correction.  But it also pretty much eliminates all the noise from the system, leaving just the raw photometry.

Things To Do

Gadzooks!  So what first?  Do the extinction software work, or focus on relative photometry?  I guess I'm leaning more towards the latter because it probably involves a bit less work.  Linear algebra was never my strong point, although as I recall it was easy to encode.  I'll work on both today, for sure.  Just not sure which one first.  Probably the relative photometry part.

I've also started looking at geometry.  What I need to start doing asap is creating phase plots.  In order to do this, I need to know where in its orbit each of these moons are.  The three pieces of information I have about that is location, distance, and direction of motion.  Phase plots are likely the best way I can visualize any kind of variations in brightness.  The jury is still out on whether or not I'm seeing variations in brightness, but that'll come when I look at more data and put several sessions together (which I'm about able to do!).

Much to do and I'll have an update (hopefully) sometime later today.  For now I'm just trying to distract my brain from just focusing on the Asteroid Data Hunter results.

Thursday, March 20, 2014

Callisto Eclipse, Asteroid Data Hunter

A surprise event happened during last night's session: Callisto was eclipsed by the shadow of Jupiter -- and I watched it happen.  Very cool.


Going....

Jupiter system 04:18 20 Mar 2014 UTC

Going ....

Jupiter system 04:27 20 Mar 2014 UTC

Gone!

Jupiter system 04:32 20 Mar 2014

More info in this to come in later posts.

At the moment, I'm busy working on my entry for the first phase of the Asteroid Data Hunter challenge.

Monday, March 17, 2014

Data From Yesterday Looks Good

Glancing at some images from last night's session has me anxious to reduce the data and look at the results.  I've reserve any expectations because sofar I've been surprised every time.  The seeing was great and I got probably the best focus sofar. 

I collected 1600 images in total.  841 of those images had Jupiter in them.  Looking at individual image statistics, I calculate that about 1141 images had at least one target.  So that's an efficiency of 71.25%.  Not bad at all.  I look forward to the day when I don't have to measure "observing efficiency".

Here's a sample of the data:

Figure 1: The Jupiter system as seem from Earth at 17 Mar 2014 03:30:24 UTC