Author Topic: Jupiter Wavelet Processing using MLT  (Read 931 times)

Offline sandconp

  • Newcomer
  • Posts: 24
Jupiter Wavelet Processing using MLT
« on: 2019 June 20 18:25:12 »
I did some imaging of Jupiter and processed the image using Autostakker and went right into PixInsight hoping I could do the equivalent of Wavelet Processing in Registax using MLT. I first performed Unsharp Mask and then Curves but the output is kind of grainy. I did not use MLT because I wasn't sure how to do it.


Any suggestions if I am doing the steps in the proper order?

Here is a link to my TIF file if anybody would care to take a stab at it.

https://www.dropbox.com/s/ouffz6q555liwab/Jup_234740_g6_ap51_conv.tif?dl=0

Offline Juan Conejero

  • PTeam Member
  • PixInsight Jedi Grand Master
  • ********
  • Posts: 7111
    • http://pixinsight.com/
Re: Jupiter Wavelet Processing using MLT
« Reply #1 on: 2019 June 21 03:34:03 »
Unsharp mask should never be used to process these images. It is more of a 'finishing' tool to improve the presentation of already processed, nonlinear images (e.g., before printing), than an actual image processing and analysis tool.

The first problem with this image, which prevents you from achieving any reasonably good result, is channel misalignment. The red and blue channels are displaced by several pixels with respect to the green channel, especially on the vertical axis. Fortunately, the image includes a satellite, which is a star-like object that can be used to fit PSF functions. So the first step is calculating the required offsets to align the red and blue channels with respect to green:


Click on the image for a full-size version.

From the cx and cy PSF parameters computed for each channel we can know the required offsets in pixels, which we can apply to the RGB image to achieve an accurate channel alignment:


Click on the image for a full-size version.

Now as you always should expect in PixInsight, there are several ways to skin this cat. The most obvious one, i.e. deconvolution, is usually the least efficient option in these cases, mainly because we normally have no idea of what the actual PSF of the image looks like. Since we have a satellite on this particular image, it is a very special case where trying out deconvolution makes sense, so here we go with a synthetic PSF generated with the DynamicPSF tool:


Click on the image for a full-size version.

The wavelet-regularized Van Cittert algorithm is normally the most efficient option for deconvolution of planetary and lunar images, where one can enjoy high signal levels. The strong noise in this image forces us to apply aggressive wavelet regularization parameters. Despite the efficiency of these algorithms, noise is clearly the main limiting factor in this case. Don't worry about the isolated bright and dark noise pixels that remain after deconvolution; they are outliers (also referred to as salt and pepper noise in the image processing jargon) that can be removed very easily with a simple median filter (MorphologicalTransformation).

A nice result, but we want more, right? An efficient and flexible way to process these images is by means of multiscale analysis tools, and we have quite a few of these in PixInsight. Here is a quick go with the starlet transform (a.k.a. à trous wavelet transform):


Click on the image for a full-size version.

The main problem with linear transforms is ringing, and planetary images are no exception. A nonlinear algorithm like the multiscale median transform or the median-wavelet transform is essentially ringing-free:


Click on the image for a full-size version.

Note how the median-wavelet transform preserves the brightness distribution over the entire planetary disk, while the starlet transform does not as a result of bright ringing generated at the disk edge. Along with this large-scale advantage, median-based transforms can enhance more small-scale structures, also as a nice side-effect of the absence of ringing. A drawback of these algorithms is noise intensification, which can lead to generation of artifacts if not properly controlled with noise reduction parameters.

Finally, a different route to process these images is applying restoration filters. In particular, the constrained least squares filtering algorithm is an excellent choice for restoration of lunar and planetary images, especially when one can work with high signal-to-noise ratios. The main problem here is that the PSF of the image is usually unknown, but this is not the case with this image, so here we go:


Click on the image for a full-size version.

I hope this can be of some help for you. It had been a long time since I played with planetary images (years, probably, I don't know), so this has been kind of a refreshing change! :)


Juan Conejero
PixInsight Development Team
http://pixinsight.com/

Offline acmalko

  • Member
  • *
  • Posts: 54
Re: Jupiter Wavelet Processing using MLT
« Reply #2 on: 2019 June 21 03:50:20 »
Hi,

I personnally use Pixinsight just after AS!3
I didn't knew how to align channel but thanks Juan, here I used Registax but you can tick the right case in AS!3 for a good result too

(some grammatical errors could appear, english is not my native language)

-channel alignment Registax
-Restoration filter
-Deconvolution
-TGVDenoise
-Color calibration
-Saturation etc...




Offline Juan Conejero

  • PTeam Member
  • PixInsight Jedi Grand Master
  • ********
  • Posts: 7111
    • http://pixinsight.com/
Re: Jupiter Wavelet Processing using MLT
« Reply #3 on: 2019 June 21 04:34:47 »
Very nice, more subtle and pleasant result than mine. Definitely, I am not a planetary processor by any stretch of the imagination :)

You have applied restoration and deconvolution in unconventional ways, guessing a PSF by trial-error work. Despite that, the result is nice and plausible visually, and of course not more 'by eye' than an arbitrary multiscale transformation. Nice work.

Here is the result of my median-wavelet transform with the correct color balance:



Compared to yours, my version is definitely in the firing line. That's what happens when somebody trained to process deep-sky data dares playing with planets!  :cheesy:
Juan Conejero
PixInsight Development Team
http://pixinsight.com/

Offline sandconp

  • Newcomer
  • Posts: 24
Re: Jupiter Wavelet Processing using MLT
« Reply #4 on: 2019 June 21 05:23:26 »
Hi,

I personnally use Pixinsight just after AS!3
I didn't knew how to align channel but thanks Juan, here I used Registax but you can tick the right case in AS!3 for a good result too

(some grammatical errors could appear, english is not my native language)

-channel alignment Registax
-Restoration filter
-Deconvolution
-TGVDenoise
-Color calibration
-Saturation etc...


So did you use Registax in addition to PixInsight?

Offline sandconp

  • Newcomer
  • Posts: 24
Re: Jupiter Wavelet Processing using MLT
« Reply #5 on: 2019 June 21 05:27:26 »
Juan,

Thank you so much for your information.  I will try to repeat what you did this weekend.

Offline acmalko

  • Member
  • *
  • Posts: 54
Re: Jupiter Wavelet Processing using MLT
« Reply #6 on: 2019 June 21 08:23:17 »
Hi,

I personnally use Pixinsight just after AS!3
I didn't knew how to align channel but thanks Juan, here I used Registax but you can tick the right case in AS!3 for a good result too

(some grammatical errors could appear, english is not my native language)

-channel alignment Registax
-Restoration filter
-Deconvolution
-TGVDenoise
-Color calibration
-Saturation etc...


So did you use Registax in addition to PixInsight?

No, I only use Registax as a fast way to align rgb channels. I very dislike this software.

Don't forget the "RGB align" option in Autostakkert

Offline acmalko

  • Member
  • *
  • Posts: 54
Re: Jupiter Wavelet Processing using MLT
« Reply #7 on: 2019 June 21 08:29:22 »
Very nice, more subtle and pleasant result than mine. Definitely, I am not a planetary processor by any stretch of the imagination :)

You have applied restoration and deconvolution in unconventional ways, guessing a PSF by trial-error work. Despite that, the result is nice and plausible visually, and of course not more 'by eye' than an arbitrary multiscale transformation. Nice work.

Here is the result of my median-wavelet transform with the correct color balance:



Compared to yours, my version is definitely in the firing line. That's what happens when somebody trained to process deep-sky data dares playing with planets!  :cheesy:

Thanks :)
Of course, I tried many parameters to achieve a good result using RestorationFilter wich is the main process.
Deconvolution or MMT are not so good for me

Offline sandconp

  • Newcomer
  • Posts: 24
Re: Jupiter Wavelet Processing using MLT
« Reply #8 on: 2019 June 21 21:40:31 »
I was able to duplicate what you did my using AS!2 and PixInsight without using Registax !

Thanks !

I am attaching my new image.