SigmaPoisson     1592 magick/gem.c         poisson=exp(-SigmaPoisson*QuantumScale*pixel);
SigmaPoisson     1598 magick/gem.c         noise=(double) (QuantumRange*i/SigmaPoisson);