Observatorio de la Ballona

Many years ago I came across a reading in Ezekiel that really struck me. The imagery was powerful and disturbing. I recently came across it again in a daily Mass reading in Magnificat. It is Ezekiel 37:1-14.

The hand of the LORD came upon me, and he led me out in the spirit of the LORD and set me in the center of the plain, which was now filled with bones.
He made me walk among them in every direction so that I saw how many they were on the surface of the plain. How dry they were!
He asked me: Son of man, can these bones come to life? “Lord GOD,” I answered, “you alone know that.”
Then he said to me: Prophesy over these bones, and say to them: Dry bones, hear the word of the LORD!
Thus says the Lord GOD to these bones: See! I will bring spirit into you, that you may come to life.
I will put sinews upon you, make flesh grow over you, cover you with skin, and put spirit in you so that you may come to life and know that I am the LORD.
I prophesied as I had been told, and even as I was prophesying I heard a noise; it was a rattling as the bones came together, bone joining bone.
I saw the sinews and the flesh come upon them, and the skin cover them, but there was no spirit in them.
Then he said to me: Prophesy to the spirit, prophesy, son of man, and say to the spirit: Thus says the Lord GOD: From the four winds come, O spirit, and breathe into these slain that they may come to life.
I prophesied as he told me, and the spirit came into them; they came alive and stood upright, a vast army.
Then he said to me: Son of man, these bones are the whole house of Israel. They have been saying, “Our bones are dried up, our hope is lost, and we are cut off.”
Therefore, prophesy and say to them: Thus says the Lord GOD: O my people, I will open your graves and have you rise from them, and bring you back to the land of Israel.
Then you shall know that I am the LORD, when I open your graves and have you rise from them, O my people!
I will put my spirit in you that you may live, and I will settle you upon your land; thus you shall know that I am the LORD. I have promised, and I will do it, says the LORD.

This time, I was more struck by the image of dry bones coming to life. A real feeling of resurrection and rebirth. This is certainly a powerful reading.

A couple of weeks back, during the last quarter Moon, I was able to get some decent data of M101, the Pinwheel Galaxy. M101 is a spiral galaxy in Ursa Major and is 27 million light years away.

The image is from 160 minutes of LRGB data. That’s 90 minutes of luminance, and 23.3 minutes each of red, green and blue, binned 2×2. LRGB imaging takes advantage of the fact that we perceive most of the detail in an image from the black and white, or luminance part of the image and less detail from the color. I obtained 90 minutes of high-resolution black and white data and combined it with lower resolution color data to produce the image. The color was obtained by binning the pixels, or adding four pixels to create one. This allows for more data collected in a shorter period of time but at half the resolution.

Here is the result. Click on the picture to go to the gallery where you can see a full-sized version of the image.

160 Minutes LRGB (90:20:20:20 subs 300:200:200:200×18:7:7:7)

It has been three months since clear skies, the new Moon, and being in Lake Riverside have converged on the same day and allowed for imaging. These factors came together on May 15th. I collimated the C-11 back at the full Moon and was ready to go. (Dew foiled my attempt to get even any Ha data that full Moon night but that is another story.)

I used my general approach for finding targets. I use The Sky’s database manager data wizard to query objects above 40 degrees or so. Moving from object to object I record my interest, the transit time, and the availability of a guide star. Up until this weekend I have tried to capture all of my data on one side of the meridian, preferably the east. This simplifies the taking of flats (with the ST-10 the camera rotates 180 degrees to maintain framing after a meridian flip) and the east has a better view from the western pier in the observatory.

I chose to image NGC 5033, a spiral galaxy in the constellation Canes Venatici. Is looks very interesting and has a very bright star that can be used for a guide star. Astronomical twilight would occur at about 9:15 PM PDT and the galaxy would transit at 10:32 PM. I hoped to get data both before and after transit, that being the best time to image as the object is highest in the sky.

I had my usual troubles with autoguiding and initial start-up but got imaging going by 9:40 or so. In the first of many little glitches, I had set CCDCommander to send me a text message when the first set of luminance images were done, but the message didn’t get sent because Zone Alarm had interrupted the request to send e-mail. Security software can be annoying. I was in the house watching Midway with the family, so I didn’t go back to get the next series started and lost 20 minutes of imaging time before the meridian crossing.

I planned on getting 50 minutes of luminance data and 20 minutes each of red, green, and blue. After having trouble on my last imaging outing with dark pixels on the color images, I wanted to make sure I had at least six sub-exposures and adequately dithered the frames—six frames to enable a sigma reject algorithm and dithering to move the dark pixels. So my luminance exposures were 300 seconds and my color were 200 seconds. I had other little issues with losing sync on the meridian flip and FocusMax acquire star not returning exactly to the original framing but overall things went fine the rest of the night. These problems would be solved if I were able to spend several days in a row in the observatory—a good goal for the future. I was able to use the new window shade to take flats after putting the white lights on the dimmer for the rope lights.

I ended up with 110 minutes of LRGB data, 50 minutes of luminance and 20 each of red, green, and blue. The RGB images were binned 2×2. I performed my data reduction, alignment, and combination in CCDStack. To simplify alignment, I process the luminance first, then save a binned copy of the final image. I then use that binned image as the master to align each of the sets of color images. This means that frames are only re-sampled once for alignment and that all combined frames are fully aligned.

Once I was processing the data it was clear that I really had inadequate imaging time for this object. At magnitude 10, it is quite dim. When the image was processed to bring detail to visibility, it was quite noisy and grainy. So I did the best that I could with it.

I did most of the processing in PixInsight. I have a standard approach when starting on an image that makes the later combine steps easier. The basic idea is to get a standard crop of the images that will work across all frames, eliminating any artifacts around the edges from the combine / integration process. The RGB frames are scaled up two times to match the luminance. Then, using the dynamic crop tool, I select a crop that will include solid data from all frames. I save an instance of the crop settings on the desktop and then save the process icons. This allows me to start and re-start work on each frame from a consistent point.

On the luminance, I used both the histogram tool and curves to bring out the details. A very nice feature of the histogram tool in PixInsight is the ability to see the number of pixels you have clipped so there is no guesswork. Generally I take several quick passes through the data before settling on a final processing approach. On this image, I took a copy of the luminance and applied a four-level HDR wavelet transform, recombining it with the original using PixelMath. This brought out the dust lanes in the center of the galaxy. I did need to create a star mask to preventing ringing when doing the HDR wavelet transform.

I tried to keep things simple on the color data. I used a channel combine to create an single RGB image, then used histogram to eyeball the black points across the colors. There is are some new tools in PixInsight that really helped the color. The first is Background Neutralization. In my experimentation with this tool, it will allow you to drop out the noisy background in a color image. I used the tool in truncate mode to eliminate most of the background data in the RGB frames. This is consistent with what RickJ suggested in the BAUT Astrophotography Forum. There is also a color calibration tool that appears to correct your color balance. I’ll have to work on understanding that.

My processing approach on images is to create solid luminance and RGB images. I then use channel extraction to pull the individual red, green, and blue components out of the RGB image. The final combine step is to use LRGB
Combination to create the final LRGB image. I touched up color and contrast and then moved to Photoshop.

In Photoshop, I used clone stamp and healing brush to fix a couple egregious color issues. I also used Noel Carboni’s Astronomy Tools actions to reduce noise in the image and tweak the star size. The final step was to add a layer processed with a high pass filter, masked to reveal only a few key locations, and merged with “soft light” blending mode to highlight some details.

Here is the final image. Still a bit noisy and needing more data, but I think it looks OK. Click on the image to go to the gallery and a full sized image. Comments and suggestions are always welcome.

We hosted a Cub Scout troop tonight at our house and the observatory (in Los Angeles). Some friends of ours son is a member of the troop and they asked us to help with a night of astronomy focused on getting the troop members their Astronomy Belt Loop.

It was great fun and I think a great success. We had 6 Scouts, parents, and a sister as our guests. After they completed their opening ceremonies, my wife went through the specific learning about astronomy that is required for the belt loop (what is a planet, etc.). Once it was dark — well, time for a digression. It is dark a bit late this time of year because of the early start of daylight savings. It is starting too early. End of digression.

I had opened the observatory at the start of the meeting. Good hoots and hollers from the scouts at the opening. Once it was dark and the lessons were complete, we started observing.

Mars was first and was not too exciting. In the 8″ SCT and a 9mm Nagler (225X magnification), it was just a small disk with almost no features. We were pushing the seeing in any case.

The next target was M42, the Nebula in Orion. The Trapezium was a clear four stars and there was a nice haze of nebulosity around it. This got a few aahs from the Scouts and other guests.

We then moved on to the Moon. This turned out to be the hit of the night. In the 35mm plossl, it was a full disk view with the terminator right in the middle of the Moon. Great detail in the maria and the peaks. Lots of oohs and aahs. I then went back to the 9mm Nagler for a closer look and this brought some serious study by one of the Scouts. Aimed at the north pole of the Moon, we had a view that showed the shape of the Moon with lots of crater detail. A big hit.

We had reached 8:40 pm or so, so the troop left except for our friends and their son. We chatted a bit and were looking at the constellations (we could pick out Orion, Gemini, Leo, Canus Major [Sirius anyway], and the Big Dipper). Then I saw a bright object at about 25 degrees and realized we had a jackpot view — Saturn!

Saturn was low but we could see the rings well and two moons clearly, Rhea and Titan. Saturn is always a big hit. So we finished the night on a high note.

It was really nice to be able to share both the knowledge of Astronomy and the experience of looking through the telescope with a bright and eager group of kids. Well behaved too. They were all very thankful. And I am very thankful that I had the ability to share the view of the sky with them.

Back over the holidays, we had a clear night but the Moon was out — almost full. After taking some Ha images of M1, I decided to try the Moon itself. It turns out that with the C-11, the new 0.75 focal reducer from Optec, and an Ha filter, a 0.11 second exposure is not blooming or over saturated on the ST-10.

So I took three series of images, one of 20 and two of 40 exposures. February’s Sky and Telescope has a very good article Lunar Photography with AviStack. Based on the advice there, I used AviStack to align and stack the images. I then took them into Photoshop where I used the Photomerge feature to create a mosaic of the three images. Here is the raw image.

Raw Moon Panorama

The next group of steps were performed in PixInsight. I will go through how I used mask and several of the image enhancement functions to produce my final image.

The first step in processing was to change the areas that weren’t part of the Mosaic from white to black. To do that, I cloned my original image by dragging the name tab onto the desktop. I used the histogram tool on the cloned image, moving the black point almost all the way to the top. I then inverted the image, producing a white background with the missing corners black. Using Pixel Match, I created a new image that was the minimum of the original image and the black corners.

The next thing I needed was a mask that would allow me to apply deconvolution, A Trous Wavelets, and HDR Wavelets to the image without creating artifacts around the edge of the Moon. I cloned my new image and reduced it to almost pure black and white using the Histogram tool, moving the white and black points to just above the background sky glow. This gave me a sharp-edged mask, but what I needed was a smoother edged mask. Without a gaussian blur function, I used A Trous Wavelets to blur the image. Here are the settings I used to blur the mask (click on the image for a larger view):

Wavelets Setting to Blur Mask

This produced the following mask image. The basic idea is that you set wavelets to many levels, drop the lower levels and slightly increase the bias near the top level.

Moon Mask

With this mask set on my main image, I went through key image enhancement steps with minor curves adjustments in between each. My first step was deconvolution. I had found that larger standard deviation settings brought out noise in the image, with single pixels leaping in brightness. After trial and error, I settled on a standard deviation of 1.5 and 80 iterations of Regularized Richardson-Lucy deconvolution. Here are the settings I used.

This had the effect of subtly sharpening the details of the Lunar surface. While this was good, I wanted to expose more of the details in this image. So at this point I turned to A Trous Wavelets. I didn’t want to over-process the image (although some might argue that doing wavelets on top of deconvolution is already that) but I wanted more detail. I applied the following, somewhat subtle, A Trous Wavelet transformation.

From there I wanted the global contrast change that you get from HDR Wavelet transformation. Note that for all of the prior sharpenings, I have been using my mask. For this transformation, I needed a sharper edge to the mask and even that left behind some brightness I needed to fix. Here is the HDR Wavelet setting I used. The number of layers is one greater than the default of 6. This matched the scale where I wanted the contrast adjustment.

This process did leave a minor artifact. The mask left a somewhat bright area around the lower edge of the Moon. Rather than tweak the mask again, after applying a very light noise reduction using ACDNR, I took the image into Photoshop. In Photoshop, I burned the lower edge of the Moon very slightly. I also added black shapes around the missing mosaic area to hid the edge artifacts from the A Trous wavelet processing.

This is the final image:

Back in October, I took some data of the Bubble Nebula (NGC 7635) and open cluster October Salt and Pepper (NGC 7634). I took 180 minutes of data:

• 90 minutes Hydrogen Alpha (Ha) in 9, 10-minute sub-exposures unbinned
• 30 minutes Luminance in 6, 5-minute sub-exposures unbinned
• 60 minutes, 20 each of Red, Green, and Blue in 2, 10-minute sub-exposures binned 2×2.

I’ll start with the finish. Here is the final image. (Click on the image for a larger version in the gallery.)

While I do think this is a nice image, it does have the pale salmon reds typical of Ha LRGB processing in basic combine method I used. My basic approach was to make an Ha-Luminance image which was then combined with the RGB to produce the final image. This started off all right as the HaL image looks fairly good. I created the image by adding 3 times the Ha to the Luminance using Pixel Math (most processing was done in PixInsight). I used both HDR Wavelets and A Trous Wavelets to sharpen up the image. Here is the HaL image used to create the final.

I then created an RGB image, adding 3 times the Ha to the red before combining with blue and green. I processed the image, stretching with histogram and curves to end up with a final RGB image. I brought both the Ha-Luminance and the RGB into Photoshop, pasting the Ha-Luminance on top of the RGB and changing the blending mode of the pasted layer to luminance. There was a lot of noise in the color image, with hot pixel color spots from each different color. Selecting the bottom, color layer, I used the spot healing brush tool to clean up the bad spots. This worked very effectively underneath the luminance layer, since only the spurious color was removed, not the details in the image. I now regret that I didn’t save a separate layer version of the image at this stage. Given the results, I figured I would go back to the beginning (and I still may).

I spent a fair amount of time tweaking the color balance, and passed the image through PixInsight’s GREYCstoration noise reduction algorithm.

Since the original processing, I have tried PixInsight’s “A New Approach to Combination of Broadband and Narrowband Data,” but my color data binned 2×2 is too noisy for it to work. I may yet try Rob Gendler’s approach. At this point, I think I am tired of working on this image.

One lesson learned from this is that when one is doing narrowband imaging, it is better to get full resolution color. This provides better data to use when creating the final image, and allows the luminance combine that creates the peachy / salmon red tones. I got better results on the Heart Nebula, but I’m not done processing that one so it will be another post.

I did try a highly stretched Ha-Luminance image to see what detail I could get out of the image. This is stretched very hard and it has brought out the noise, but there is a lot of detail in the image. This is created from the maximum of 5 * Ha or the Luminance image in Pixel Math, then stretched with curves and some blackpoint adjustment in the histogram tool. I did average in an HDR Wavelet and A Trous Wavelet to enhance the detail in the nebula. Click for a bigger image as the thumbnail is pretty small.

I received an e-mail from Carla Befera, who is doing national PR for 400 Years of the Telescope, a show coming up on PBS. She asked me to put in a plug for the show so my readers would know about it. It looks really good. It may have aired in your area as the general release was April 10, but it is on this Thursday, April 16 at 7pm on KCET in Los Angeles.

Here is the blurb from the press release:

This visually stunning program chronicles a sweeping journey, from 1609, when Galileo revealed mankind’s place in the galaxy, to today’s thrilling quests to discover new worlds in the universe. Narrated by NOVA’s Neil deGrasse Tyson, the compelling program takes viewers on an adventure through the heavens and around the globe, visiting the world’s leading astronomers, cosmologists and observatories.

I did let Carla know that I’m not sure how many readers I have. I am tickled pink that I even showed up on the radar! And I am certainly willing to put in a plug for what looks to be an excellent program.

UPDATE

It appears that the show will be on KCET on Thursday April 30th at 9pm. The earlier broadcast date is for “KCET World,” whatever that is. The KCET website is full of bugs too. Links don’t work, information is incomplete. Poorly executed. But this still looks like a good show.

This is a new version of an image I took back in January 2007. I spent two nights capturing data, but I have never been satisfied with the original processing. I have finally have gotten around to trying again.

On preliminary comment in self-defense here. This is a difficult object, even under dark skies. It is diffuse, with the broader areas quite dim. I was imaging in West Los Angeles, with very light-polluted skies. Much of the galaxy is only slightly brighter than the background sky glow in the city. I may have taken eight hours of data, but seeing was not that great and the signal was barely above the background.

Here is the original image. It has many flaws. First, it is mirror-flipped so the orientation is incorrect. It should have been cropped before it was processes to eliminate the lines all around the edges. It is stretched too hard and looks over processed.

Here is the revised image.

I was much less aggressive on stretching the image. Image reduction (darks and flats) was done in Maxim DL. I aligned and stacked the images in CCDStack. I found that the bicubic b-spline resampling produced the least noise in the align process. I used mean combination to combine the subframes. Technique-wise, I created a master luminance frame and then used that to align the R/G/B frames.

Most of the processing was done in PixInsight. I did most of the enhancement work was done on the luminance image. I cropped the image, rescaled it to use the full dynamic range. This was required since the signal was very low and I had combined with a mean combine rather than a sum. A sum would have produced higher totals, but the net result is the same as a mean. I stretched the image with historgram, but did not clip it but for a small amount at the lower brightness levels.

I did try to do most of the processing on a linear version of the image, per the advice of the experts from PixInsight. While I was able to more easily get a star mask built, HDR Wavelet and other functions just didn’t work well.

With the image stretched appropriately, I copied it and did an HDR wavelet transform with a star mask applied. This was averaged with the original image. That was essentially the final luminance image.

For the RGB frames, I did a similar stretch process, attempting to create similar-looking histograms for each one. I then combined all frames into a color image using LRGB combination. I brought that image into Photoshop to adjust the color balance and then used Russ Croman’s Gradient Xterminator to flatten the background. Finally, back into PixInsight for some brightness touch up and noise reduction.

The exposure is LRGB 8 hours, 10 minutes total exposure, L = 115 min (23×5 min) + RGB 375 minutes, RGB = 125 min (25×5 min), all unbinned. The equipment: C-11 @f5.95, CGE mount, ST-10, CFW-8A, AO-7, & Hutech LPS Filter.

AIC 2008 wrapped up just over a week ago. It was a good conference, fairly different in content and organization than prior years. The core of the conference did not focus all that much on astrophotography, but rather on general space science, publishing, and 3-D rendering. I suppose it was good to get general theory rather than cookbook demonstrations as it has been in the past. From an astrophotography standpoint, the sessions on Friday were quite good. This balanced the more general content of Saturday and Sunday with smaller group sessions on tools and techniques.

As always, it was good to see friends from prior conferences and talk about imaging, observatories, and equipment. The vendor section was quite good, with more vendors than last year. I am now again interested in pursuing Hyperstar imaging given a very convincing discussion with the folks at Starizona. My current dream scope is a Planewake CDK 12.5. But the observatory comes first.

Here are some photos of the conference. This first is a shot of the main conference area.

AIC 2008 Meeting Room

Bud Guinn and Venkata model the 3-D glasses we were given for the 3-D rendering presentation.

Astro-Physics, Planewave Instruments, and ASA telescopes provided an impressive display.

The Chronos Mount is an intriguing product. It uses harmonic drives and is an equatorial mount that does not require a meridian flip.

I’m already looking forward to next year’s meeting. I would like to see some presentations on the math behind image processing.

08:00 PST — Start of the final day at AIC. Last night included a demo of Light Buckets, an on-line telescope service. Charges add up quickly as it is from $80 to $175 per hour. Next year’s conference is set at the end of October. It will be larger, they have reserved the upstairs room. Bigger and better.

Ray Gralak Advanced Image Combine Techniques CCD usage, starting in 90s. CCD stacking led to great dynamic range, DDP made it visible, but artifacts in the image came out. These are cosmic ray hits, hot/cold pixels, bad columns, satellite trails, plane trails, and asteroids. Four raw image types: Bias, Dark, Flat, and Light frames. Dark and Bias frames may change over time. Data collection tips: Take all frames at the same temperature, take the same number of darks as lights, use light frame duration for darks, dither your lights and flats, flat frames should be at or near focus and at the same orientation of the light frames, a flat for each filter. Workflow: Create master Bias, master dark, subtract bias from each flat, normalize and combine flat frames, subtract master dark from light, apply flat, align and combine.

Combine methods: Average, median, min / max clip, sigma, SDM. Average, good SNR, but artifacts remain. Median, best noise rejection, but lower SNR. Min/Max, rejects most artifacts, but leaves some, must have >6 images. Sigma Clip, strong noise rejection, requires >10 images to work best, does make errors to reject good or not reject bad. SDM, calculates mean and median, then looks at standard deviation of all pixels, if STD is above a certain number, it uses median, otherwise uses mean. Multiply a sigma factor to select “certain number,” if sigma factor is zero, all pixels are rejected and the mean is used. Lost of detail on settings (I took a picture).

09:00 PST — Wolfgang Promper Making the Most of Your Site This is a topic for me, about imaging from an urban location. I get it! We should image where we can do it more often. Normalize background: Equal pixel counts on background in each channel. Interesting, he cuts out of Maxim and pastes into Photoshop. Gradient Xterminator as primary gradient removal, then manual tweaking of the image. High pass and threshold to find the stars (SSRO has instructions). Add stars as new color channel with cut and paste to use the channel as a selection, expand and feather selection by 2 pixels, 2 pixel minimum filter to reduce size of stars. Use extensive selective color to balance nebular color. DDP 100 lower than auto. “After a lot of processing and versions you often find that the first version was the best.”

09:32 PST — Break

10:00 PST — Door Prizes! Art won a telecompressor!

10:30 PST — Conference Over