optimising the dull job of image stacking

The cosmos at our fingertips.
User avatar
stephen63
Science Officer
Posts: 304
Joined: Sat Jan 12, 2013 2:53 am
Location: Pa
Contact:

Re: optimising the dull job of image stacking

Post by stephen63 » Wed Oct 23, 2013 5:28 pm

OK, Chris
I guess read noise IS the limiting factor. So, after the images are stacked, as long as the SNR is at least 3 then it can be separated from the noise.

User avatar
Chris Peterson
Abominable Snowman
Posts: 18573
Joined: Wed Jan 31, 2007 11:13 pm
Location: Guffey, Colorado, USA
Contact:

Re: optimising the dull job of image stacking

Post by Chris Peterson » Wed Oct 23, 2013 5:36 pm

stephen63 wrote:OK, Chris
I guess read noise IS the limiting factor. So, after the images are stacked, as long as the SNR is at least 3 then it can be separated from the noise.
S/N > 3 is commonly used as a target in scientific imaging. For aesthetic imaging, people usually want much larger values.
Chris

*****************************************
Chris L Peterson
Cloudbait Observatory
https://www.cloudbait.com

User avatar
Nitpicker
Inverse Square
Posts: 2692
Joined: Fri Sep 20, 2013 2:39 am
Location: S27 E153

Re: optimising the dull job of image stacking

Post by Nitpicker » Wed Oct 23, 2013 10:42 pm

Chris Peterson wrote:S/N > 3 is commonly used as a target in scientific imaging. For aesthetic imaging, people usually want much larger values.
No wonder I still think my humble photos still have some slight scientific merit. I always try to remind myself that the universe is a noisy place.

User avatar
Beyond
500 Gigaderps
Posts: 6889
Joined: Tue Aug 04, 2009 11:09 am
Location: BEYONDER LAND

Re: optimising the dull job of image stacking

Post by Beyond » Wed Oct 23, 2013 11:24 pm

Nitpicker wrote:
Chris Peterson wrote:S/N > 3 is commonly used as a target in scientific imaging. For aesthetic imaging, people usually want much larger values.
No wonder I still think my humble photos still have some slight scientific merit. I always try to remind myself that the universe is a noisy place.
Yeah :!: It's a good thing 'space' is a hard vacuum, or we wouldn't be able to sleep at night. :lol2:
To find the Truth, you must go Beyond.

User avatar
geckzilla
Ocular Digitator
Posts: 9180
Joined: Wed Sep 12, 2007 12:42 pm
Location: Modesto, CA
Contact:

Re: optimising the dull job of image stacking

Post by geckzilla » Thu Oct 24, 2013 4:31 pm

Chris Peterson wrote:
geckzilla wrote:Coincidentally, I am messing around with four pretty bad exposures from the HLA with low S/N ratio and a crapload of cosmic rays. I'm getting some interesting results by taking four different possible combinations of pairs of them and using minimum which I presume is some kind of floor function and then creating a median stack out of those four combinations. I thought it might be a good way to recover some of the signal loss from using only a minimum stack.
What software are you using? Because median processing increases noise, the method has largely been replaced with more complex stacking algorithms, which do a better job of eliminating single frame artifacts while preserving the overall S/N similar to summing. Sigma mask, sigma clip, and similar algorithms are generally better choices.
I have been thinking about this thread as well as this post since yesterday. I downloaded CCDStack since it offers these more complex algorithms and tried them out. After a crash course with their tutorial, I recreated my homebaked Photoshop method along with two sum stacks using STD sigma data rejection and Poisson sigma data rejection. All results looked "good" as far as this particular data goes (it's the faint, outer shells of the Calabash Nebula. It's hopelessly noisy and chock full of CRs.) Interestingly to me, the Poisson sigma data rejection looked very similar to my home baked idea of pairs of minimum stacks combined into a median stack.

Visually, it's hard for me to determine the S/N so I used the software calculation. Poisson sigma sum came in last at 0.31, my silly median minimum second at 0.34, and STD sigma sum on top at 0.39 for the same rectangle around the object for each.

The bad part about the STD sigma sum stack is that it also left quite a few pieces of cosmic rays for whatever reason, so even though the data overall is less noisy, for aesthetic purposes it might be harder to work with. Then again, maybe a day worth of learning and tinkering with it wasn't enough for me to know how to fine tune it.

PS - I hope you don't mind me interjecting this in your thread, Nitpicker.
Just call me "geck" because "zilla" is like a last name.

User avatar
Chris Peterson
Abominable Snowman
Posts: 18573
Joined: Wed Jan 31, 2007 11:13 pm
Location: Guffey, Colorado, USA
Contact:

Re: optimising the dull job of image stacking

Post by Chris Peterson » Thu Oct 24, 2013 4:46 pm

geckzilla wrote:The bad part about the STD sigma sum stack is that it also left quite a few pieces of cosmic rays for whatever reason, so even though the data overall is less noisy, for aesthetic purposes it might be harder to work with. Then again, maybe a day worth of learning and tinkering with it wasn't enough for me to know how to fine tune it.
When it comes to rejecting transient artifacts, it helps to have more frames. You don't usually have that luxury with HST data.
Chris

*****************************************
Chris L Peterson
Cloudbait Observatory
https://www.cloudbait.com

User avatar
geckzilla
Ocular Digitator
Posts: 9180
Joined: Wed Sep 12, 2007 12:42 pm
Location: Modesto, CA
Contact:

Re: optimising the dull job of image stacking

Post by geckzilla » Thu Oct 24, 2013 5:01 pm

Chris Peterson wrote:
geckzilla wrote:The bad part about the STD sigma sum stack is that it also left quite a few pieces of cosmic rays for whatever reason, so even though the data overall is less noisy, for aesthetic purposes it might be harder to work with. Then again, maybe a day worth of learning and tinkering with it wasn't enough for me to know how to fine tune it.
When it comes to rejecting transient artifacts, it helps to have more frames. You don't usually have that luxury with HST data.
This one actually had 7 exposures which is a little more than usual. One of them was so short it seemed to degrade the whole stack even with sum so I left it out. It needs a couple more hour long exposures and some good luck with the CRs smooth it out. No wonder they had to point the thing so long for the ultra deep fields, which it seems like they are going to do again. Going to watch the Hubble Hangout in three hours to find out more. :)
Just call me "geck" because "zilla" is like a last name.

Post Reply