British Astronomical Association
Supporting amateur astronomers since 1890

Secondary menu

Main menu

BAA Community Member Pages Paul Leyland
Image search

An experiment in super-resolution, part 2

Image

Expand image

About this observation
Observer
Paul Leyland
Time of observation
11/11/2018 - 01:55
Object
ES 1128
Observing location
Tacande Observatory
Equipment
0.4m Dilworth-Relay
SBIG-8 with AO-7 autoguider
No filter
Exposure
75 seconds in four subs of 5, 10, 20, 40s
Like This Image
Share

https://britastro.org/node/19348 shows my first, quick, and dirty attempt at super-resolution.  Since then I've implemented the CLEAN deconvolution algorithm.  See https://en.wikipedia.org/wiki/CLEAN_(algorithm) for an introduction and links to further reading.  2000 iterations with a gain of 0.1, followed by a 2-pixel radius Gaussian blur of the CLEAN map and the re-addition of the residuals produced the image above.

It's markedly better than the Richardson-Lucy reconstruction produced earlier, even though it started with exactly the same raw image and PSF.  The signal to noise ratio is much better, the limiting magnitude is better and there is clear sky between the components of the double star.

Comments

Xilman's picture

The image is much easier to measure than its predecessor.  The centroids differ in RA by 0.23 pixels and in Declination by 4.93 pixels.  At an image scale of 0.71 arcsec/pixel the separation evaluates to 3.50 arcsec.  The P.A. is now 177°.

Gaia-DR2 gives 3.247 arcsec with a PA of 172°28'.

nickoftime's picture

This is a very good effort and I'm interested in knowing what it is that leads to a difference of ~0.25 arcsec, measure vs Gaia? Is it the uncertainty associated with measuring the distance between star centres or the derivation of the 0.71 arcsec/pixel calibration factor, or both, or something else? 

Xilman's picture

The honest answer is that I don't know.  A hand-waving answer is that it may be related to distortion of the image caused by the deconvolution process.  A closely related hand-waving answer is that perhaps the PSF wasn't sufficiently good.  In any event, deconvolution artefacts are obvious, as shown by the halo in the vicinity of the deconvolved stars.

It might be interesting to try two further approaches, either individually or in concert.

1) Resample the image and the PSF onto a finer grid and deconvolve the result on to that grid.  Note there are a potentially infinite number of ways of resampling.

2) Use other deconvolution algorithms.  Both CLEAN and R-L are cheap and cheerful.  Maximum entropy, for example, may be much better but it is also much more difficult to implement.

The derivation of the 0.71 arcsec/pixel scale is almost certainly not the issue.  Plate-solving of many thousands of images, as well as computations based on simple geometry of focal length and pixel size, indicates that that figure is good to better than 1%.

nickoftime's picture

Thank you for the explanation. Keep chipping away at it!

Xilman's picture

Thanks.  If you, or anyone else, would like to try their hand at super-resolution  the FITS image can be provided. An ex-colleague at Cambridge University has already expressed interest at giving it to his students as an exercise in signal processing.  Until now he's used synthetic data.

Xilman's picture

According to Table 1 of JDSO V8.2 (April 2012) p102, ES1128 has components of magnitude 11.31 and 12.50, separation 2.87 arcsec and PA 178.5°.

Those numbers are much closer to my measurements of the RL-processed image than those from Gaia-DR2.

Not sure what to make of this.

Copyright of all images and other observations submitted to the BAA remains with the owner of the work. Reproduction of the work by third-parties is expressly forbidden without the consent of the copyright holder. For more information, please contact the webmaster.