//povprečje
                double povpOrg = 0;
                double povpDec = 0;
                for(int i=0;i<or.Length;i++)
                {
                    povpOrg += or[i];
                    povpDec += de[i];
                }
                povpOrg /= or.Length;
                povpDec /= de.Length;

                //varianca
                double varOrg = 0;
                double varDec = 0;
                double covar = 0;
                for (int i = 0; i < or.Length; i++)
                {
                    varOrg += Math.Pow(or[i] - povpOrg, 2);
                    varDec += Math.Pow(de[i] - povpDec, 2);

                    //covar
                    covar += (or[i] - povpOrg) * (de[i] - povpDec);
                }
                varOrg /= or.Length;
                varDec /= or.Length;

                varOrg = Math.Sqrt(varOrg);
                varDec = Math.Sqrt(varDec);

                covar /= or.Length;

                double ssim = (2 * povpOrg * povpDec + Math.Pow(255 * 0.01,2)) * (2 * covar + Math.Pow(0.03 * 255,2));
                ssim /= (Math.Pow(povpOrg, 2) + Math.Pow(povpDec, 2) + Math.Pow(255 * 0.01, 2)) * (Math.Pow(varDec,2) + Math.Pow(varOrg,2) + Math.Pow(0.03 * 255, 2));