1
$\begingroup$

Within my field a recent study suggested to use the symmetric properties of certain image datasets to improve signal to noise ratio (SNR). I will spare you the details, but in the end one can get a symmetric dataset (as described below) with improved SNR. However the way that one typically computes SNR is to compare the "amount of signal" (sum of the pixels in the central region, more or less), to the overall standard deviation (std).

The properties of the symmetrized data are such that if one computes the std over the entire image (in a standard way), std is indeed lower than in the original image, leading to the improved SNR.

However I am bit skeptical regarding the fact that one should compute the std in the same way as in a "normal image". Indeed, here we have an image with identical pixels according to certain symmetries (around central vertical axis, horizontal vertical axis, both diagonals, and [$\pi/2$] rotations). Many pixels are hence correlated, and I am not sure that simply computing the std across the whole map is a fair assessment of the SNR.

The question is: how do I compute the standard deviation from such a dataset? Should I take into account the fact that many pixels are not independent (correlated)? How?

example symmetric data set (don't mind the central black circle): https://i.stack.imgur.com/TzlAS.png

enter image description here

$\endgroup$
2
  • 1
    $\begingroup$ Welcome to Mathematics SE. Take a tour. You'll find that simple "Here's the statement of my question, solve it for me" posts will be poorly received. What is better is for you to add context (with an edit): What you understand about the problem, what you've tried so far, etc.; something both to show you are part of the learning experience and to help us guide you to the appropriate help. You can consult this link for further guidance. $\endgroup$
    – AnotherUser
    Sep 14, 2022 at 12:45
  • $\begingroup$ do you know what are the symmetries of images? are they the same across different images? $\endgroup$
    – Aksakal
    Sep 21, 2022 at 13:47

1 Answer 1

0
$\begingroup$

The short answer is that it should be Ok to use SNR as you described in certain case as explained below.

Consider this: the simplest way to improve the example image knowing its symmetry would be to calculate the average on 8 repeating slices. One slice goes from the center to the right horizontally and from the center to upper right corner at 45 degrees. So you stack up all 8 slices, and create a new slice whose pixels are average across 8 slices, then create a new image by repeating this new slice 8 times. The resulting image should have lower noise.

The way you defined the noise includes two components: complexity of an image and bona fide noise. If your image was a simple white noise, then the contribution from the complexity would be zero. In your case there's some structure in the background, which contributes to the noise. By using the symmetry, we reduce noise coming from the bona fide noise, and not the complexity of the background.

So, the simple standard deviation of the whole image will reflect the reduction of bona fide noise after image processing. Hence, it seems to me it's an acceptable metric to compare images with the same symmetric background. If the backgrounds are different, then this metric would not account for the complexity of the background. Whether you need to account for this or not depends on your goals.

In your image the symmetry group is probably $D_2$. This is a very interesting question indeed.

background of the signal

If I had to account for the background of the image, then something like Kernel density could be used to extract the background on the slice. Then I'd calculate the noise over the background and use it to get stdev. This way we could compare noise between images of different background complexity.

Consider the variance equation: $$\frac 1 n \sum_i(x_i-\bar x)^2$$ Instead of using a global average $\bar x$, we'll calculate the local average $\mu_i$, perhaps, using Kernel density estimators such as $$\mu_i=\sum_j K(x_i-x_j)x_j$$ Then the local variance would be $$\frac 1 n \sum_i(x_i-\mu_i)^2$$

This equation excludes the complexity of the background that is represented with smoothed pixels $\mu_i$, and makes it possible to fairly compare images of different backgrounds.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.