For the Poisson distribution the mean ${\displaystyle m}$  and variance ${\displaystyle v}$  are not independent: ${\displaystyle m=v}$ . The Anscombe transform^[1]

${\displaystyle A:x\mapsto 2{\sqrt {x+{\tfrac {3}{8}}}}\,}$ 

aims at transforming the data so that the variance is set approximately 1 for large enough mean; for mean zero, the variance is still zero.

It transforms Poissonian data ${\displaystyle x}$  (with mean ${\displaystyle m}$ ) to approximately Gaussian data of mean ${\displaystyle 2{\sqrt {m+{\tfrac {3}{8}}}}-{\tfrac {1}{4\,m^{1/2}}}+O\left({\tfrac {1}{m^{3/2}}}\right)}$  and standard deviation ${\displaystyle 1+O\left({\tfrac {1}{m^{2}}}\right)}$ . This approximation gets more accurate for larger ${\displaystyle m}$ ,^[2] as can be also seen in the figure.

For a transformed variable of the form ${\displaystyle 2{\sqrt {x+c}}}$ , the expression for the variance has an additional term ${\displaystyle {\frac {{\tfrac {3}{8}}-c}{m}}}$ ; it is reduced to zero at ${\displaystyle c={\tfrac {3}{8}}}$ , which is exactly the reason why this value was picked.

When the Anscombe transform is used in denoising (i.e. when the goal is to obtain from ${\displaystyle x}$  an estimate of ${\displaystyle m}$ ), its inverse transform is also needed in order to return the variance-stabilized and denoised data ${\displaystyle y}$  to the original range. Applying the algebraic inverse

${\displaystyle A^{-1}:y\mapsto \left({\frac {y}{2}}\right)^{2}-{\frac {3}{8}}}$ 

usually introduces undesired bias to the estimate of the mean ${\displaystyle m}$ , because the forward square-root transform is not linear. Sometimes using the asymptotically unbiased inverse^[1]

${\displaystyle y\mapsto \left({\frac {y}{2}}\right)^{2}-{\frac {1}{8}}}$ 

mitigates the issue of bias, but this is not the case in photon-limited imaging, for which the exact unbiased inverse given by the implicit mapping^[3]

${\displaystyle \operatorname {E} \left[2{\sqrt {x+{\tfrac {3}{8}}}}\mid m\right]=2\sum _{x=0}^{+\infty }\left({\sqrt {x+{\tfrac {3}{8}}}}\cdot {\frac {m^{x}e^{-m}}{x!}}\right)\mapsto m}$ 

should be used. A closed-form approximation of this exact unbiased inverse is^[4]

${\displaystyle y\mapsto {\frac {1}{4}}y^{2}-{\frac {1}{8}}+{\frac {1}{4}}{\sqrt {\frac {3}{2}}}y^{-1}-{\frac {11}{8}}y^{-2}+{\frac {5}{8}}{\sqrt {\frac {3}{2}}}y^{-3}.}$ 

There are many other possible variance-stabilizing transformations for the Poisson distribution. Bar-Lev and Enis report^[2] a family of such transformations which includes the Anscombe transform. Another member of the family is the Freeman-Tukey transformation^[5]

${\displaystyle A:x\mapsto {\sqrt {x+1}}+{\sqrt {x}}.\,}$ 

A simplified transformation, obtained as the primitive of the reciprocal of the standard deviation of the data, is

${\displaystyle A:x\mapsto 2{\sqrt {x}}\,}$ 

which, while it is not quite so good at stabilizing the variance, has the advantage of being more easily understood. Indeed, from the delta method,

${\displaystyle V[2{\sqrt {x}}]\approx \left({\frac {d(2{\sqrt {m}})}{dm}}\right)^{2}V[x]=\left({\frac {1}{\sqrt {m}}}\right)^{2}m=1}$ .

While the Anscombe transform is appropriate for pure Poisson data, in many applications the data presents also an additive Gaussian component. These cases are treated by a Generalized Anscombe transform^[6] and its asymptotically unbiased or exact unbiased inverses.^[7]

• Variance-stabilizing transformation
• Box–Cox transformation

• Starck, J.-L.; Murtagh, F. (2001), "Astronomical image and signal processing: looking at noise, information and scale", Signal Processing Magazine, IEEE, vol. 18, no. 2, pp. 30–40, Bibcode:2001ISPM...18...30S, doi:10.1109/79.916319, S2CID 13210703