Подавление вейвлет коэффициентов байесовским подходом

                double T = 0.1;
                double sum = 0;
                double sko = 0;
                for (int i = 0; i < data.GetLength(0); i++)
                    for (int j = 0; j < data.GetLength(1); j++)
                        sum += data[i, j];
                for (int i = 0; i < data.GetLength(0); i++)
                    for (int j = 0; j < data.GetLength(1); j++)
                        sko = data[i, j] - sum / (data.GetLength(0) * data.GetLength(1));
                sko = sko * sko / data.GetLength(0) * data.GetLength(1);
                double s = 0;
                double ro = 0;
                double M = 0;
                for (int i = 1; i < data.GetLength(0) - 1; i++)
                    for (int j = 1; j < data.GetLength(1) - 1; j++)
                    {
                        s = Math.Sqrt(0.5 * (sko - (
                            (Math.Pow(data[i, j], 2) * 1 / 9 + Math.Pow(data[i + 1, j], 2) * 1 / 9 + Math.Pow(data[i + 1, j + 1], 2) * 1 / 9 + Math.Pow(data[i, j + 1], 2) * 1 / 9 + Math.Pow(data[i - 1, j + 1], 2) * 1 / 9 + Math.Pow(data[i - 1, j], 2) * 1 / 9 + Math.Pow(data[i - 1, j - 1], 2) * 1 / 9 + Math.Pow(data[i, j - 1], 2) * 1 / 9 + Math.Pow(data[i + 1, j - 1], 2) * 1 / 9) -
                            Math.Pow(data[i, j] * 1 / 9 + data[i + 1, j] * 1 / 9 + data[i + 1, j + 1] * 1 / 9 + data[i, j + 1] * 1 / 9 + data[i - 1, j + 1] * 1 / 9 + data[i - 1, j] * 1 / 9 + data[i - 1, j - 1] * 1 / 9 + data[i, j - 1] * 1 / 9 + data[i + 1, j - 1] * 1 / 9, 2))));
                        ro = Math.Exp(-T / s) / 2 * s;
                        M = (data[i + 1, j] + data[i + 1, j + 1] + data[i, j + 1] + data[i - 1, j + 1] + data[i - 1, j] + data[i - 1, j - 1] + data[i, j - 1] + data[i + 1, j - 1]) / 8;
                        data[i, j] = data[i, j] * (1 - 1 / (1 + ro * (
                            ((Math.Exp(T / s) / 2 * s) * Math.Exp(-Math.Abs(data[i, j]) / s)) /
                            ((Math.Exp(T / s) / (2 * s * (Math.Exp(T / s) - 1)) * Math.Exp(-Math.Abs(data[i, j]) / s)))) *
                            (
                            ((Math.Exp(T / s) / 2 * s) * Math.Exp(-Math.Abs(M) / s)) /
                            ((Math.Exp(T / s) / (2 * s * (Math.Exp(T / s) - 1)) * Math.Exp(-Math.Abs(M) / s))))));
                    }