An interesting problem recently developed during my latest post-processing of the Medusa nebula. My standard process for processing images is to use the WBPP script in PixInsight, as it has really been improved over the years and automates the post-processing tasks in a nice, convenient way. I’ve been using it for awhile now, with 5 recent processing tasks completed with no major issues. But with the Medusa nebula things turned out a little strange.
For the Medusa I managed to get 50 Ha, 32 Oii and 17 Sii subs over a period of two evenings. I took flats and darkflats (25 each) for each of the two sessions/evenings as per my usual work flow. Using the WBPP script, with a master dark of 300sec for the NB subs and the flats and darkflats taken during the evenings I processed the set of light data. Examining the final stacked images I noticed that the Ha image seemed to have very little dynamic range.
Suspecting the Cosmetic Correction that I had applied I re-ran the script with no CC. The results were better, but still not up to standard. I ran the readout probe over the image taking some readings in the background and then in the brightest part of the nebula. The background was about 13 ADU, with the brightest regions around 66 ADU. I thought that maybe 300 seconds was insufficient time to get a reasonable level of data, and that maybe the case, but the 5 minute subs did show some fairly good range in the individual subs (as the following table shows).
Similar probing on the Oiii and Sii stacks showed 224-392 for the Oiii stack and 118-220 for the Sii. If 300 seconds was good enough for Oiii and Sii, it should have been just fine for the Ha data. Generally, the Ha signal is the stronger, and it should have been for this object.
I checked the parameters on WBPP and they seemed fine. Nothing stood out; typical settings, similar to what I usually set for my runs.
PI’s Statistics process yielded the following data:
Single sub Exp:300 Gain:100 O:50 Temp: -10
Filter Mean Std Dev Min Max
Ha 508.0 27.5 338 62448
Oiii 519.2 30.8 352 56090
Sii 505.8 28.1 337 64355
Processed in WBPP
Filter Mean Std Dev Min Max Stack count
Ha 7.8 11.7 2.2 5212.7 50
Oiii 242.1 143.2 131.4 64138.7 32
Sii 128.6 153.4 52.0 65535.0 17
The sample subs seemed fine with almost identical statistics. However, the calibrated masters showed significant differences. The Oiii and Sii data was somewhat expected, but the Ha data seemed way off the mark.
I then decided to take the data and run the calibration manually, without using WBPP script. The results were very interesting to say the least. The Oiii and Sii data were very close to each other even though the results using WBPP showed a difference. But the real interesting part was the fact that the Ha data was also very close to both the Oiii and Sii data. The stacked master means were very close to the value of the single sample subs for all three filters.
Integration (manual)
Filter Mean Std Dev Min Max Stack count
Ha 507.5 11.4 496.8 5186.2 50
Oiii 521.4 14.0 433.0 6057.3 32
Sii 506.2 10.1 343.7 4370.7 17
Now I know WBPP may do some additional processing that I may not have done when I calibrated and integrated manually, but the Ha data was at least an order of magnitude better.
I am at a loss as to what is going on here. I was able to complete the post-processing of the Medusa data using the manually calibrated subs and the final result was very respectable. Even the WBPP data was usable, but required much more processing finesse.
Any idea on what I am seeing here?
For the Medusa I managed to get 50 Ha, 32 Oii and 17 Sii subs over a period of two evenings. I took flats and darkflats (25 each) for each of the two sessions/evenings as per my usual work flow. Using the WBPP script, with a master dark of 300sec for the NB subs and the flats and darkflats taken during the evenings I processed the set of light data. Examining the final stacked images I noticed that the Ha image seemed to have very little dynamic range.
Suspecting the Cosmetic Correction that I had applied I re-ran the script with no CC. The results were better, but still not up to standard. I ran the readout probe over the image taking some readings in the background and then in the brightest part of the nebula. The background was about 13 ADU, with the brightest regions around 66 ADU. I thought that maybe 300 seconds was insufficient time to get a reasonable level of data, and that maybe the case, but the 5 minute subs did show some fairly good range in the individual subs (as the following table shows).
Similar probing on the Oiii and Sii stacks showed 224-392 for the Oiii stack and 118-220 for the Sii. If 300 seconds was good enough for Oiii and Sii, it should have been just fine for the Ha data. Generally, the Ha signal is the stronger, and it should have been for this object.
I checked the parameters on WBPP and they seemed fine. Nothing stood out; typical settings, similar to what I usually set for my runs.
PI’s Statistics process yielded the following data:
Single sub Exp:300 Gain:100 O:50 Temp: -10
Filter Mean Std Dev Min Max
Ha 508.0 27.5 338 62448
Oiii 519.2 30.8 352 56090
Sii 505.8 28.1 337 64355
Processed in WBPP
Filter Mean Std Dev Min Max Stack count
Ha 7.8 11.7 2.2 5212.7 50
Oiii 242.1 143.2 131.4 64138.7 32
Sii 128.6 153.4 52.0 65535.0 17
The sample subs seemed fine with almost identical statistics. However, the calibrated masters showed significant differences. The Oiii and Sii data was somewhat expected, but the Ha data seemed way off the mark.
I then decided to take the data and run the calibration manually, without using WBPP script. The results were very interesting to say the least. The Oiii and Sii data were very close to each other even though the results using WBPP showed a difference. But the real interesting part was the fact that the Ha data was also very close to both the Oiii and Sii data. The stacked master means were very close to the value of the single sample subs for all three filters.
Integration (manual)
Filter Mean Std Dev Min Max Stack count
Ha 507.5 11.4 496.8 5186.2 50
Oiii 521.4 14.0 433.0 6057.3 32
Sii 506.2 10.1 343.7 4370.7 17
Now I know WBPP may do some additional processing that I may not have done when I calibrated and integrated manually, but the Ha data was at least an order of magnitude better.
I am at a loss as to what is going on here. I was able to complete the post-processing of the Medusa data using the manually calibrated subs and the final result was very respectable. Even the WBPP data was usable, but required much more processing finesse.
Any idea on what I am seeing here?