hi sir,
I am new to cryoSPARC. I have started to learn recently. Kindly help me understand the difference between symmetric, white, and colored noise models.? When should we impose this restriction- during ab initio or nonuniform refinement?
Please explain to me.
Hi @BhawnaMishra, is there a problem in your data analysis that you think relates to the noise model? I’ve rarely needed to touch these kinds of parameters and rarely see it discussed. E.g.,
Hi @BhawnaMishra, and welcome to cryoEM!
To begin with, I agree with @drichman — we have found that the default noise models are almost always the right choice, and I wouldn’t generally recommend changing them from their default values.
However, I would still be happy to to explain the difference in these noise models, and what we use noise models for generally.
What do we mean when we say noise has a color?
We can define “noise” as a series of values with no information content. Put another way, noise is a series of random values. However, defining exactly what we mean by “random” is not as obvious as it first seems!
We can describe different types of noise (randomness) by their “color”.
An image of white noise from Wikipedia
White noise is equally noisy at all frequencies (more formally: white noise has equal power over all frequency bands). This type of noise sounds like an old untuned TV or radio, and is not so common in natural systems. It is called white noise by analogy to white light, which has an equal amount of all the visible frequencies.
An image of pink noise from Wikipedia
Colored noise has a different “amount” of randomness at different frequencies. For instance, pink noise has more low-frequency noise than high frequency noise. It sounds like rushing water, and is more commonly observed in natural systems.
If you compare the image of pink noise to the image of white noise, you see larger regions of the image that have the same average brightness. This is because pink noise has more variation at low frequencies than at high frequencies!
In cryoEM, we generally use a colored noise system, since we expect that our images have more noise at the higher frequencies. We also often observe a very noisy peak in the medium frequencies, the source of which we do not yet fully understand.
What is symmetric noise?
Our images are 2D, so we need to generate 2D noise as well. However, we have no reason to expect that our images will be noisier in, say, the vertical direction than the horizontal direction. In symmetric noise we therefore enforce rotational symmetry on the noise, such that noise of a given frequency in direction A will have the same intensity (that is, the same level of “noisiness”) as the same frequency in direction B.
In 2D plots of noise, this appears as rings of the same intensity, such as the right-hand side of the plot below image. However, with a symmetric noise model, we can just plot a single line from the center of this plot outward, since all such lines will have the same profile. This rotational average is what’s plotted on the left-hand side of the plot below.
Why do we need a noise model in the first place?
You may be wondering what the noise model is used for at all — why do we need to know how noisy our images are, if there’s nothing we can do to remove that noise? The answer is that we can use our noise model to tell us how trustworthy our images are at a given frequency.
For instance, say our noise model tells us that the images are particularly noisy at 10 Å for some reason. During alignment, we can downweight the contribution of the particle images at 10 Å, since we know they are not particularly reliable at that frequency. This may improve the overall quality of the alignment, since we’re only using frequencies we’re more confident in.
I hope that answers any questions you might have, but please feel free to let me know if there’s something I didn’t quite explain clearly!
If we have more high frequency noise on average, is the noise model effectively acting as a fine-tuned low-pass filter for purposes of alignment?
Great question @ccgauvin94!
I want to first be very clear that we do not actually use the noise model to filter the particle images at any point — this is all just to build a mental model of what it’s doing.
With that out of the way, your model is more-or-less correct! It’s not a low-pass filter per se, since there is the mid-resolution noise I mentioned earlier (i.e., the noise is not monotonically increasing), and the response is not necessarily 1.0 anywhere over the frequency domain, but that doesn’t mean it’s a bad mental model.
It’s probably best to think about this as a prior on the error, e.g., “I expect the error due to noise at frequency K to be higher than at frequency L”. We can factor this prior into the score of a given alignment, so that an alignment with more error at frequency K is treated as a better alignment than one with the same error at frequency L. So if in this example K was always greater than L, you could certainly think of it as a tuned low-pass filter during alignment.
The most formally-correct way to think about this is that the noise model is a filter applied to the residual, i.e., the difference between the reference projection and the particle image, in Fourier space. We modulate the error values in each Fourier shell by the amount of noise in that shell before picking the best pose.
Dear sir,
In the Noise Model, generally, the default parameter is symmetric. But when we click there it shows it can be changed to white, or colored are well depending on the data. I have tried to find the explanation for that but I was unsuccessful. But I am keen to understand them in detail. Kindly also quote a link where I can read about them. Also I would like you to put some shine on the fact that is it required to change those parameters or we should keep it default only?
Thankyou so much for your explanation sir
Sir I request you to please tell how to interpret the graph? What is x and y axis? What does this sphere indicate? Also how to downweight the contribution of noise at what angstron? What job tells this? Which parameter?
Hi @BhawnaMishra!
I assume we’re still talking about the noise model graph above? If so, you can interpret the circular graphs as, essentially, the same as the 1-D version to their left – we’re modelling how noisy the data is at each given frequency. A cooler color (more blue) means less noisy, while a warmer color (more red) means noisier. I wouldn’t worry too much about directly interpreting these graphs, to be honest with you – as long as it looks more or less smooth, it’s probably fine!
The X- and Y-axes are in inverse Å. Essentially, the color of a pixel further from the center represents the noise a higher frequency/resolution than a pixel closer to the center. They may look spherical because of the coloring, but there is no real 3-D aspect.
The noise is handled automatically in all CryoSPARC jobs, so there is no specific job you need to run to handle it!
If I misinterpreted any of your questions, please let me know! I’m happy to elaborate on anything I misunderstood
sir, I wanted to say that since I want to understand cryosparc in-depth, I am having difficulty understanding what this sigma and error indicate, and why the sigma is orange outside and the error is blue outside. I am sorry to ask this every time but I need to know.
Hi @BhawnaMishra! Don’t be sorry for asking questions, that’s what we’re here for!
I want to begin by saying that this noise model plot is not expected to be useful for data processing – it’s an internal diagnostic for our use. I personally never look at these noise model plots.
With that out of the way, the difference in color of these two plots is not meaningful. We do not compute or consider Fourier components outside the current resolution. In error, they are set to 0; in sigma, they are set to 1, but these pixels are not used for anything in either case, so the color has no meaning. Sorry if that’s confusing!