Go Back   IceInSpace > Beginners Start Here > Beginners Talk
Register FAQ Calendar Today's Posts Search

Reply
 
Thread Tools Rate Thread
  #1  
Old 26-06-2006, 04:43 PM
davidpretorius's Avatar
davidpretorius
lots of eyes on you!

davidpretorius is offline
 
Join Date: Jun 2005
Location: Launceston Tasmania
Posts: 7,381
Jupiter image files to test your astra image skills on

Hi Guys

http://www.precons.com/iceinspace/ga...r_060426_1147/

grab some of the .tiff files and then split in astraimage and then have a play with deconvolution and unsharp mask and fft etc.

there are various stack sizes and wavelet settings.

Ask questions and post your results to see how different deconvolution works?
Reply With Quote
  #2  
Old 26-06-2006, 05:03 PM
janoskiss's Avatar
janoskiss (Steve H)
Registered User

janoskiss is offline
 
Join Date: Jan 2005
Location: Sale, VIC
Posts: 6,033
there are quite a lot of files there. pick one or two for us Davo.
Reply With Quote
  #3  
Old 26-06-2006, 05:15 PM
Robert_T's Avatar
Robert_T
aiming for 2nd Halley's

Robert_T is offline
 
Join Date: Oct 2005
Location: Australia
Posts: 2,959
Dave, is FFT (what is it?) something that's available on the mid-range Astra Image 2.0 or only on the full blown version?

cheers,
Reply With Quote
  #4  
Old 26-06-2006, 05:21 PM
janoskiss's Avatar
janoskiss (Steve H)
Registered User

janoskiss is offline
 
Join Date: Jan 2005
Location: Sale, VIC
Posts: 6,033
FFT = Fast Fourier Transform. It provides an efficient way or performing a convolution operations (which in turn are used in deconvolution algorithms).
Reply With Quote
  #5  
Old 26-06-2006, 07:28 PM
davidpretorius's Avatar
davidpretorius
lots of eyes on you!

davidpretorius is offline
 
Join Date: Jun 2005
Location: Launceston Tasmania
Posts: 7,381
in the demo version, to the right of VC deconvolutions
Reply With Quote
  #6  
Old 27-06-2006, 01:23 AM
asimov's Avatar
asimov (John)
Planet photographer

asimov is offline
 
Join Date: May 2005
Location: Bundaberg
Posts: 8,819
Cheers Davo, just what I need...more decon practice!
Reply With Quote
  #7  
Old 27-06-2006, 01:49 AM
asimov's Avatar
asimov (John)
Planet photographer

asimov is offline
 
Join Date: May 2005
Location: Bundaberg
Posts: 8,819
VC deconvolution 5 iterations/0.6 curves FFT edit on R & G. Not the result I was looking for...Back in for another go!
Attached Thumbnails
Click for full-size image (daves two joops.jpg)
110.0 KB43 views
Reply With Quote
  #8  
Old 01-07-2006, 12:45 AM
janoskiss's Avatar
janoskiss (Steve H)
Registered User

janoskiss is offline
 
Join Date: Jan 2005
Location: Sale, VIC
Posts: 6,033
I did not use Astra Image but here is my effort anyway. I am at home sick today so this was a nice way to pass the time, keep my mind off my pounding headache.

For the most part I used a piece of free software called ImageJ, and the FFTJ and DeconvolutionJ plugins. I use it in Linux but it runs on just about anything (it's java based).

For best results deconvolution should be applied before any other image processing operations. So I chose a this starting image (see first attached pic): 1__ST149_WV10_10_58_91_118_186.tif, because it roughly looks like it has the degree of blurring expected for the telescope used (which will blur out of details smaller than about 4-7 pixels depending on colour), and does not show any visible signs of processing, i.e., sharpening (unlike e.g. 1__ST149_WV10_10_58_138_293_455.tif , which is a bit sharper but shows a bit of a white halo around the moon shadow).

For the deconvolution operations I estimated the point spread function (PSF) based on the following parameters:

* diameter of Jupiter on image = 412 pixels (measured with kruler off the screen, +/- 4 pixels)
* angular size of Jupiter = 40.1 arc sec (assume picture taken on 25/06/06)
=> we have pixel size = 0.09733 arc sec
* Telescope aperture = 250mm (10" Dob)
* Central obstruction = 65mm (info about 10" Dob in previous posts on IIS)
* wavelength = 475nm (blue), 510nm (green), 650nm (red)

I wrote a short C++ program to calculate the PSF with Bessel function J1 from the GNU Scientific Library. I generated 3 PSFs, one for each colour. (See attached image - PSFs are shown at four times actual size and gamma adjusted to show ring structure.)

Decomposed the image into RGB channels and deconvoluted each using the corresponding PSF. The FFTJ plugin is a one-step direct method of deconvolution. There is a "regularisation parameter" that needs to be chosen appropiately (for best compromise between noise and sharpness). I chose these values by inspection: red: 0.002, green: 0.00005, blue: 0.001, which gave good detail plus a bit of speckle noise.

Then:
- RGB merge
- Despeckle: ImageJ has a very good parameter-free despeckle filter and cleans up the deconvolution noise very nicely.
- Level tweaks: black for each channel separately set to darkest pixel on moon shadow, white to global max (all channels together)
- Gamma: 0.67 0.70 0.60 (RGB) - by eye
- correct for slight purple-green fringing by realigning colours: shift red channel 2/3 pixel up and 1/3 pixel right (blow up image + move + shrink image) - see attached pic.

Result shown in last attached image. I'm quite pleased with the result for my first attempt at this sort of thing.

PS. I kept getting this strange error when trying to open the original tiffs with any image editors I tried: unknown field with tag 317 (0x13d) encountered. Fixed with a pass-through raw2tiff.
Attached Thumbnails
Click for full-size image (jup-01.jpg)
76.1 KB26 views
Click for full-size image (psf.png)
61.4 KB22 views
Click for full-size image (chroma.png)
23.1 KB30 views
Click for full-size image (jup-final.jpg)
60.9 KB35 views

Last edited by janoskiss; 01-07-2006 at 08:39 AM.
Reply With Quote
  #9  
Old 01-07-2006, 01:07 AM
Harpspitfire
Registered User

Harpspitfire is offline
 
Join Date: Jun 2006
Location: pittsburgh pa
Posts: 268
bit off topic here, but is there any download for a deconvultion file anywhere without having to have it in a whole processing program? i tried iris and couldnt even figure out the program, something else (i forget) but it had the L-R decon, van client?, and another kind- ended up with jup images that werent any better at all- thx, john
Reply With Quote
  #10  
Old 01-07-2006, 01:29 AM
C.
Registered User

C. is offline
 
Join Date: Jan 2006
Location: South Australia
Posts: 34
What a good idea for those of us who haven't started imaging yet. I just used my regular image program and had a play. Thank you for the opportunity David. I hope I can load it ok.

Janoskiss I hope you are feeling better now.
Attached Thumbnails
Click for full-size image (jup-01.jpg)
76.1 KB16 views
Click for full-size image (jup-02.jpg)
62.2 KB22 views
Reply With Quote
  #11  
Old 01-07-2006, 10:09 AM
janoskiss's Avatar
janoskiss (Steve H)
Registered User

janoskiss is offline
 
Join Date: Jan 2005
Location: Sale, VIC
Posts: 6,033
John, not sure what you mean. If you are after dedicated deconv software, they tend to be highly specialised and more complicated than integrated deconv routines in imaging programs.

I uploaded the PSFs I used for deconvolution to http://www.iceinspace.com.au/uploads/PSF.zip in case anyone else wants to have a go using these.

Once you have the PSFs, deconvolution is very simple to do with ImageJ. You need to get the program from the Downloads page, and the "FFTJ and DeconvolutionJ" plugin from the Plugins page. To deconvolute a greyscale image:

* open greyscale image to be deconvoluted (use original 16-bit tiff)
* open corresponding PSF file (also 16-bit tiff)
* increase canvas sizes to nearest power of two (if they are not already, in this case 512x512; use Image -> Adjust -> Canvas size)
* Plugins -> FFTJ -> DeconvolutionJ:
. - select blurred image to be deconv'ed and PSF image
. - enter a value for regularisation parameter or leave at default (can experiment here for best results)
. - you might like to do a Process -> Noise -> Despeckle on the deconv'ed image

One thing to note about the PSFs here is that they assume that the diffraction limit of the telescope is the main source of blur in the image, which is likely incorrect. Other sources of blurring, e.g. due to poor seeing, can be thought of as convolution operations too, and can be corrected for using deconvolution. Then one would need to estimate the PSF from the blurred image itself.
Reply With Quote
  #12  
Old 01-07-2006, 12:32 PM
asimov's Avatar
asimov (John)
Planet photographer

asimov is offline
 
Join Date: May 2005
Location: Bundaberg
Posts: 8,819
You certainly know your stuff Steve. I tend not to look into the technical side of it too much; I know next to nothing of the terms & so forth, & basically have no concept of what I'm doing. I do like to muck about in Astraimage to see if a better final image can be derived from it. Great job having never done it before ! Your effort is a lot smoother looking than mine.
Reply With Quote
  #13  
Old 01-07-2006, 12:32 PM
davidpretorius's Avatar
davidpretorius
lots of eyes on you!

davidpretorius is offline
 
Join Date: Jun 2005
Location: Launceston Tasmania
Posts: 7,381
thats awesome steve, excellent stuff!

more stuff to read up on!
Reply With Quote
  #14  
Old 01-07-2006, 01:15 PM
Harpspitfire
Registered User

Harpspitfire is offline
 
Join Date: Jun 2006
Location: pittsburgh pa
Posts: 268
thx for the info steve- after looking at everything thats way beyond my processing skills!- i was just hoping to get a better final image- guess ill just have to play with photoshop for now- john
Reply With Quote
  #15  
Old 01-07-2006, 01:33 PM
asimov's Avatar
asimov (John)
Planet photographer

asimov is offline
 
Join Date: May 2005
Location: Bundaberg
Posts: 8,819
Well if your anything like me John (& it sounds like you are) get yourself Astraimage & just have a play around. http://precons.com/iceinspace/gallery/ai25max.zip
Reply With Quote
  #16  
Old 01-07-2006, 01:34 PM
asimov's Avatar
asimov (John)
Planet photographer

asimov is offline
 
Join Date: May 2005
Location: Bundaberg
Posts: 8,819
BTW nice job C
Reply With Quote
Reply

Bookmarks


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

Forum Jump


All times are GMT +10. The time is now 07:36 AM.

Powered by vBulletin Version 3.8.7 | Copyright ©2000 - 2024, Jelsoft Enterprises Ltd.
Advertisement
Testar
Advertisement
Bintel
Advertisement