Stochastic UpResolution

5 replies113 views
Robert Majewski avatar

A number of imaging systems are undersampled for a variety of reasons. For example, a user may want a wide field of view, but the imaging sensor used has a small size. Also, a system that has pixels corresponding to a small solid angle in object space may result in low SNR imagery. In order to compensate, a user may want to bin the pixels or use a lower f number optical system in order to obtain a higher SNR image. These systems will result in a phenomena referred to as aliasing. High spatial frequency components of the imagery will be aliased to lower spatial frequencies due to the way the image is sampled.

The Drizzle algorithm was developed (Williams4 et al. 1996; Fruchter & Hook 1998) in order to compensate for undersampling with the first set of CCD detectors used with the Hubble Space Telescope.

I have attached a PDF thatStochastic UpResolution.pdf discusses another method to remove aliasing with undersampled imagery. This is accomplished by processing a large number of images for which there are random image shifts. The processing method is similar to the Drizzle algorithm1. However, the details are different and its ability to eliminate aliasing is proved in the appendix of this paper.

Well Written Helpful Insightful
Robin Bosshard avatar

Impressive stuff!

You might clarifiy on Figure 3 that Stochastic UpRes has been applied, or even have a comparison like figures 1 & 2. It would also be helpful, for a more practical context, how many images where used to produce figure 3, since the number in your proof (figure 5) is rather high, while fig. 2 comes from a suprisingly low number of images.

If you put that in a PI Script, I’ll ditch Drizzle for it!

Clear Skies!

Helpful Insightful Respectful Engaging Supportive
Tony Gondola avatar

The first test result in the paper looks impressive. I’m not sure it’s increasing resolution but certainly improves contrast. Hopefully the script makers out there will give it a closer look.

Well Written
Robert Majewski avatar

Robin

In figure 1 only 1 image was taken. In figure 2 128 images were taken. In figure 3, M1 the crab 1344 15 second exposures were taken. The image dither used is obtained by the fact that the mount tracking is not quite perfect. I get it free of charge.

https://www.astrobin.com/8ys0uu/

My software is in a Macintosh complied c application. I do not know how to write a PI script. Maybe other astrobin members would know how to do that.

The image shift measurement algorithm is described in this PDF

Image Autoguiding Theory3.pdf

Robin Bosshard avatar

Robert Majewski · Jan 31, 2026, 10:53 PM

In figure 3, M1 the crab 1344 15 second exposures were taken

Thanks Robert, that’s reasonably practical number, compared to th 10k for fig.5 (which I would not want to put through WBPP for a comparison😉).

Robert Majewski avatar

Robin Bosshard · Jan 31, 2026 at 11:01 PM

Robert Majewski · Jan 31, 2026, 10:53 PM

In figure 3, M1 the crab 1344 15 second exposures were taken

Thanks Robert, that’s reasonably practical number, compared to th 10k for fig.5 (which I would not want to put through WBPP for a comparison😉).

The 10 K value is just something that was easy to do in Mathematica. I was for a one dimensional function. However the basic mathematics is the same.

Bob